Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome.

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms. VIDEO ABSTRACT.

mBARq: a versatile and user-friendly framework for the analysis of DNA barcodes from transposon insertion libraries, knockout mutants, and isogenic strain populations.

MOTIVATION: DNA barcoding has become a powerful tool for assessing the fitness of strains in a variety of studies, including random transposon mutagenesis screens, attenuation of site-directed mutants, and population dynamics of isogenic strain pools. However, the statistical analysis, visualization, and contextualization of the data resulting from such experiments can be complex and require bioinformatic skills. RESULTS: Here, we developed mBARq, a user-friendly tool designed to simplify these steps for diverse experimental setups. The tool is seamlessly integrated with an intuitive web app for interactive data exploration via the STRING and KEGG databases to accelerate scientific discovery. AVAILABILITY AND IMPLEMENTATION: The tool is implemented in Python. The source code is freely available (https://github.com/MicrobiologyETHZ/mbarq) and the web app can be accessed at: https://microbiomics.io/tools/mbarq-app.

A multiproducer microbiome generates chemical diversity in the marine sponge Mycale hentscheli.

Bacterial specialized metabolites are increasingly recognized as important factors in animal-microbiome interactions: for example, by providing the host with chemical defenses. Even in chemically rich animals, such compounds have been found to originate from individual members of more diverse microbiomes. Here, we identified a remarkable case of a moderately complex microbiome in the sponge host Mycale hentscheli in which multiple symbionts jointly generate chemical diversity. In addition to bacterial pathways for three distinct polyketide families comprising microtubule-inhibiting peloruside drug candidates, mycalamide-type contact poisons, and the eukaryotic translation-inhibiting pateamines, we identified extensive biosynthetic potential distributed among a broad phylogenetic range of bacteria. Biochemical data on one of the orphan pathways suggest a previously unknown member of the rare polytheonamide-type cytotoxin family as its product. Other than supporting a scenario of cooperative symbiosis based on bacterial metabolites, the data provide a rationale for the chemical variability of M. hentscheli and could pave the way toward biotechnological peloruside production. Most bacterial lineages in the compositionally unusual sponge microbiome were not known to synthesize bioactive metabolites, supporting the concept that microbial dark matter harbors diverse producer taxa with as yet unrecognized drug discovery potential.

Assessing microbiome population dynamics using wild-type isogenic standardized hybrid (WISH)-tags.

Microbiomes feature recurrent compositional structures under given environmental conditions. However, these patterns may conceal diverse underlying population dynamics that require intrastrain resolution. Here we developed a genomic tagging system, termed wild-type isogenic standardized hybrid (WISH)-tags, that can be combined with quantitative polymerase chain reaction and next-generation sequencing for microbial strain enumeration. We experimentally validated the performance of 62 tags and showed that they can be differentiated with high precision. WISH-tags were introduced into model and non-model bacterial members of the mouse and plant microbiota. Intrastrain priority effects were tested using one species of isogenic barcoded bacteria in the murine gut and the Arabidopsis phyllosphere, both with and without microbiota context. We observed colonization resistance against late-arriving strains of Salmonella Typhimurium in the mouse gut, whereas the phyllosphere accommodated Sphingomonas latecomers in a manner proportional to their presence at the late inoculation timepoint. This demonstrates that WISH-tags are a resource for deciphering population dynamics underlying microbiome assembly across biological systems.

Indole prevents Escherichia coli cell division by modulating membrane potential.

Indole is a bacterial signalling molecule that blocks E. coli cell division at concentrations of 3-5 mM. We have shown that indole is a proton ionophore and that this activity is key to the inhibition of division. By reducing the electrochemical potential across the cytoplasmic membrane of E. coli, indole deactivates MinCD oscillation and prevents formation of the FtsZ ring that is a prerequisite for division. This is the first example of a natural ionophore regulating a key biological process. Our findings have implications for our understanding of membrane biology, bacterial cell cycle control and potentially for the design of antibiotics that target the cell membrane.

Metabolic interaction models recapitulate leaf microbiota ecology.

Resource allocation affects the structure of microbiomes, including those associated with living hosts. Understanding the degree to which this dependency determines interspecies interactions may advance efforts to control host-microbiome relationships. We combined synthetic community experiments with computational models to predict interaction outcomes between plant-associated bacteria. We mapped the metabolic capabilities of 224 leaf isolates from Arabidopsis thaliana by assessing the growth of each strain on 45 environmentally relevant carbon sources in vitro. We used these data to build curated genome-scale metabolic models for all strains, which we combined to simulate >17,500 interactions. The models recapitulated outcomes observed in planta with >89% accuracy, highlighting the role of carbon utilization and the contributions of niche partitioning and cross-feeding in the assembly of leaf microbiomes.

Indole inhibition of ColE1 replication contributes to stable plasmid maintenance.

In the absence of active partitioning, strict control of plasmid copy number is required to minimise the possibility of plasmid loss at bacterial cell division. An important cause of multicopy plasmid instability is the formation of plasmid dimers by recombination and their subsequent proliferation by over-replication in a process known as the dimer catastrophe. This leads to the formation of dimer-only cells in which plasmid copy number is substantially lower than in cells containing only monomers, and which have a greatly increased probability of plasmid loss at division. The accumulation of dimers triggers the synthesis of the regulatory small RNA, Rcd, which stimulates tryptophanase and increases the production of indole. This, in turn, inhibits Escherichia coli cell division. The Rcd checkpoint hypothesis proposes that delaying cell division allows time for the relatively slow conversion of plasmid dimers to monomers by Xer-cer site-specific recombination. In the present work we have re-evaluated this hypothesis and concluded that a cell division block is insufficient to prevent the dimer catastrophe. Plasmid replication must also be inhibited. In vivo experiments have shown that indole, when added exogenously to a broth culture of E. coli does indeed stop plasmid replication as well as cell division. We have also shown that indole inhibits the activity of DNA gyrase in vitro and propose that this is the mechanism by which plasmid replication is blocked. The simultaneous effects of upon growth, cell division and DNA replication in E. coli suggest that indole acts as a true cell cycle regulator.

Genome sequences of Rhizopogon roseolus, Mariannaea elegans, Myrothecium verrucaria, and Sphaerostilbella broomeana and the identification of biosynthetic gene clusters for fungal peptide natural products.

In recent years, a variety of fungal cyclic peptides with interesting bioactivities have been discovered. For many of these peptides, the biosynthetic pathways are unknown and their elucidation often holds surprises. The cyclic and backbone N-methylated omphalotins from Omphalotus olearius were recently shown to constitute a novel class (borosins) of ribosomally synthesized and posttranslationally modified peptides, members of which are produced by many fungi, including species of the genus Rhizopogon. Other recently discovered fungal peptide macrocycles include the mariannamides from Mariannaea elegans and the backbone N-methylated verrucamides and broomeanamides from Myrothecium verrucaria and Sphaerostilbella broomeana, respectively. Here, we present draft genome sequences of four fungal species Rhizopogon roseolus, Mariannaea elegans, Myrothecium verrucaria, and Sphaerostilbella broomeana. We screened these genomes for precursor proteins or gene clusters involved in the mariannamide, verrucamide, and broomeanamide biosynthesis including a general screen for borosin-producing precursor proteins. While our genomic screen for potential ribosomally synthesized and posttranslationally modified peptide precursor proteins of mariannamides, verrucamides, broomeanamides, and borosins remained unsuccessful, antiSMASH predicted nonribosomal peptide synthase gene clusters that may be responsible for the biosynthesis of mariannamides, verrucamides, and broomeanamides. In M. verrucaria, our antiSMASH search led to a putative NRPS gene cluster with a predicted peptide product of 20 amino acids, including multiple nonproteinogenic isovalines. This cluster likely encodes a member of the peptaibols, an antimicrobial class of peptides previously isolated primarily from the Genus Trichoderma. The nonribosomal peptide synthase gene clusters discovered in our screenings are promising candidates for future research.

The plant NADPH oxidase RBOHD is required for microbiota homeostasis in leaves.

The plant microbiota consists of a multitude of microorganisms that can affect plant health and fitness. However, it is currently unclear how the plant shapes its leaf microbiota and what role the plant immune system plays in this process. Here, we evaluated Arabidopsis thaliana mutants with defects in different parts of the immune system for an altered bacterial community assembly using a gnotobiotic system. While higher-order mutants in receptors that recognize microbial features and in defence hormone signalling showed substantial microbial community alterations, the absence of the plant NADPH oxidase RBOHD caused the most pronounced change in the composition of the leaf microbiota. The rbohD knockout resulted in an enrichment of specific bacteria. Among these, we identified Xanthomonas strains as opportunistic pathogens that colonized wild-type plants asymptomatically but caused disease in rbohD knockout plants. Strain dropout experiments revealed that the lack of RBOHD unlocks the pathogenicity of individual microbiota members driving dysbiosis in rbohD knockout plants. For full protection, healthy plants require both a functional immune system and a microbial community. Our results show that the NADPH oxidase RBOHD is essential for microbiota homeostasis and emphasizes the importance of the plant immune system in controlling the leaf microbiota.

High throughput sequencing provides exact genomic locations of inducible prophages and accurate phage-to-host ratios in gut microbial strains.

BACKGROUND: Temperate phages influence the density, diversity and function of bacterial populations. Historically, they have been described as carriers of toxins. More recently, they have also been recognised as direct modulators of the gut microbiome, and indirectly of host health and disease. Despite recent advances in studying prophages using non-targeted sequencing approaches, methodological challenges in identifying inducible prophages in bacterial genomes and quantifying their activity have limited our understanding of prophage-host interactions. RESULTS: We present methods for using high-throughput sequencing data to locate inducible prophages, including those previously undiscovered, to quantify prophage activity and to investigate their replication. We first used the well-established Salmonella enterica serovar Typhimurium/p22 system to validate our methods for (i) quantifying phage-to-host ratios and (ii) accurately locating inducible prophages in the reference genome based on phage-to-host ratio differences and read alignment alterations between induced and non-induced prophages. Investigating prophages in bacterial strains from a murine gut model microbiota known as Oligo-MM12 or sDMDMm2, we located five novel inducible prophages in three strains, quantified their activity and showed signatures of lateral transduction potential for two of them. Furthermore, we show that the methods were also applicable to metagenomes of induced faecal samples from Oligo-MM12 mice, including for strains with a relative abundance below 1%, illustrating its potential for the discovery of inducible prophages also in more complex metagenomes. Finally, we show that predictions of prophage locations in reference genomes of the strains we studied were variable and inconsistent for four bioinformatic tools we tested, which highlights the importance of their experimental validation. CONCLUSIONS: This study demonstrates that the integration of experimental induction and bioinformatic analysis presented here is a powerful approach to accurately locate inducible prophages using high-throughput sequencing data and to quantify their activity. The ability to generate such quantitative information will be critical in helping us to gain better insights into the factors that determine phage activity and how prophage-bacteria interactions influence our microbiome and impact human health. Video abstract.

A general non-self response as part of plant immunity.

Plants, like other multicellular lifeforms, are colonized by microorganisms. How plants respond to their microbiota is currently not well understood. We used a phylogenetically diverse set of 39 endogenous bacterial strains from Arabidopsis thaliana leaves to assess host transcriptional and metabolic adaptations to bacterial encounters. We identified a molecular response, which we termed the general non-self response (GNSR) that involves the expression of a core set of 24 genes. The GNSR genes are not only consistently induced by the presence of most strains, they also comprise the most differentially regulated genes across treatments and are predictive of a hierarchical transcriptional reprogramming beyond the GNSR. Using a complementary untargeted metabolomics approach we link the GNSR to the tryptophan-derived secondary metabolism, highlighting the importance of small molecules in plant-microbe interactions. We demonstrate that several of the GNSR genes are required for resistance against the bacterial pathogen Pseudomonas syringae. Our results suggest that the GNSR constitutes a defence adaptation strategy that is consistently elicited by diverse strains from various phyla, contributes to host protection and involves secondary metabolism.

Identification, heterologous production and bioactivity of lentinulin A and dendrothelin A, two natural variants of backbone N-methylated peptide macrocycle omphalotin A.

Backbone N-methylation and macrocyclization improve the pharmacological properties of peptides by enhancing their proteolytic stability, membrane permeability and target selectivity. Borosins are backbone N-methylated peptide macrocycles derived from a precursor protein which contains a peptide α-N-methyltransferase domain autocatalytically modifying the core peptide located at its C-terminus. Founding members of borosins are the omphalotins from the mushroom Omphalotus olearius (omphalotins A-I) with nine out of 12 L-amino acids being backbone N-methylated. The omphalotin biosynthetic gene cluster codes for the precursor protein OphMA, the protease prolyloligopeptidase OphP and other proteins that are likely to be involved in other post-translational modifications of the peptide. Mining of available fungal genome sequences revealed the existence of highly homologous gene clusters in the basidiomycetes Lentinula edodes and Dendrothele bispora. The respective borosins, referred to as lentinulins and dendrothelins are naturally produced by L. edodes and D. bispora as shown by analysis of respective mycelial extracts. We produced all three homologous peptide natural products by coexpression of OphMA hybrid proteins and OphP in the yeast Pichia pastoris. The recombinant peptides differ in their nematotoxic activity against the plant pathogen Meloidogyne incognita. Our findings pave the way for the production of borosin peptide natural products and their potential application as novel biopharmaceuticals and biopesticides.

Complex general stress response regulation in Sphingomonas melonis Fr1 revealed by transcriptional analyses.

The general stress response (GSR) represents an important trait to survive in the environment by leading to multiple stress resistance. In alphaproteobacteria, the GSR is under the transcriptional control of the alternative sigma factor EcfG. Here we performed transcriptome analyses to investigate the genes controlled by EcfG of Sphingomonas melonis Fr1 and the plasticity of this regulation under stress conditions. We found that EcfG regulates genes for proteins that are typically associated with stress responses. Moreover, EcfG controls regulatory proteins, which likely fine-tune the GSR. Among these, we identified a novel negative GSR feedback regulator, termed NepR2, on the basis of gene reporter assays, phenotypic analyses, and biochemical assays. Transcriptional profiling of signaling components upstream of EcfG under complex stress conditions showed an overall congruence with EcfG-regulated genes. Interestingly however, we found that the GSR is transcriptionally linked to the regulation of motility and biofilm formation via the single domain response regulator SdrG and GSR-activating histidine kinases. Altogether, our findings indicate that the GSR in S. melonis Fr1 underlies a complex regulation to optimize resource allocation and resilience in stressful and changing environments.

Synthetic microbiota reveal priority effects and keystone strains in the Arabidopsis phyllosphere.

Multicellular organisms, including plants, are colonized by microorganisms, some of which are beneficial to growth and health. The assembly rules for establishing plant microbiota are not well understood, and neither is the extent to which their members interact. We conducted drop-out and late introduction experiments by inoculating Arabidopsis thaliana with synthetic communities from a resource of 62 native bacterial strains to test how arrival order shapes community structure. As a read-out we tracked the relative abundance of all strains in the phyllosphere of individual plants. Our results showed that community assembly is historically contingent and subject to priority effects. Missing strains could, to various degrees, invade an already established microbiota, which was itself resistant and remained largely unaffected by latecomers. Additionally, our results indicate that individual strains of Proteobacteria (Sphingomonas, Rhizobium) and Actinobacteria (Microbacterium, Rhodococcus) have the greatest potential to affect community structure as keystone species.