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.

Uncoupling of invasive bacterial mucosal immunogenicity from pathogenicity.

There is the notion that infection with a virulent intestinal pathogen induces generally stronger mucosal adaptive immunity than the exposure to an avirulent strain. Whether the associated mucosal inflammation is important or redundant for effective induction of immunity is, however, still unclear. Here we use a model of auxotrophic Salmonella infection in germ-free mice to show that live bacterial virulence factor-driven immunogenicity can be uncoupled from inflammatory pathogenicity. Although live auxotrophic Salmonella no longer causes inflammation, its mucosal virulence factors remain the main drivers of protective mucosal immunity; virulence factor-deficient, like killed, bacteria show reduced efficacy. Assessing the involvement of innate pathogen sensing mechanisms, we show MYD88/TRIF, Caspase-1/Caspase-11 inflammasome, and NOD1/NOD2 nodosome signaling to be individually redundant. In colonized animals we show that microbiota metabolite cross-feeding may recover intestinal luminal colonization but not pathogenicity. Consequent immunoglobulin A immunity and microbial niche competition synergistically protect against Salmonella wild-type infection.

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.

Microbial abundance, activity and population genomic profiling with mOTUs2.

Metagenomic sequencing has greatly improved our ability to profile the composition of environmental and host-associated microbial communities. However, the dependency of most methods on reference genomes, which are currently unavailable for a substantial fraction of microbial species, introduces estimation biases. We present an updated and functionally extended tool based on universal (i.e., reference-independent), phylogenetic marker gene (MG)-based operational taxonomic units (mOTUs) enabling the profiling of >7700 microbial species. As more than 30% of them could not previously be quantified at this taxonomic resolution, relative abundance estimates based on mOTUs are more accurate compared to other methods. As a new feature, we show that mOTUs, which are based on essential housekeeping genes, are demonstrably well-suited for quantification of basal transcriptional activity of community members. Furthermore, single nucleotide variation profiles estimated using mOTUs reflect those from whole genomes, which allows for comparing microbial strain populations (e.g., across different human body sites).

The Ocean Gene Atlas: exploring the biogeography of plankton genes online.

The Ocean Gene Atlas is a web service to explore the biogeography of genes from marine planktonic organisms. It allows users to query protein or nucleotide sequences against global ocean reference gene catalogs. With just one click, the abundance and location of target sequences are visualized on world maps as well as their taxonomic distribution. Interactive results panels allow for adjusting cutoffs for alignment quality and displaying the abundances of genes in the context of environmental features (temperature, nutrients, etc.) measured at the time of sampling. The ease of use enables non-bioinformaticians to explore quantitative and contextualized information on genes of interest in the global ocean ecosystem. Currently the Ocean Gene Atlas is deployed with (i) the Ocean Microbial Reference Gene Catalog (OM-RGC) comprising 40 million non-redundant mostly prokaryotic gene sequences associated with both Tara Oceans and Global Ocean Sampling (GOS) gene abundances and (ii) the Marine Atlas of Tara Ocean Unigenes (MATOU) composed of >116 million eukaryote unigenes. Additional datasets will be added upon availability of further marine environmental datasets that provide the required complement of sequence assemblies, raw reads and contextual environmental parameters. Ocean Gene Atlas is a freely-available web service at: http://tara-oceans.mio.osupytheas.fr/ocean-gene-atlas/.

Innate immunity restricts Citrobacter rodentium A/E pathogenesis initiation to an early window of opportunity.

Citrobacter rodentium infection is a mouse model for the important human diarrheal infection caused by enteropathogenic E. coli (EPEC). The pathogenesis of both species is very similar and depends on their unique ability to form intimately epithelium-adherent microcolonies, also known as "attachment/effacement" (A/E) lesions. These microcolonies must be dynamic and able to self-renew by continuous re-infection of the rapidly regenerating epithelium. It is unknown whether sustained epithelial A/E lesion pathogenesis is achieved through re-infection by planktonic bacteria from the luminal compartment or local spread of sessile bacteria without a planktonic phase. Focusing on the earliest events as C. rodentium becomes established, we show here that all colonic epithelial A/E microcolonies are clonal bacterial populations, and thus depend on local clonal growth to persist. In wild-type mice, microcolonies are established exclusively within the first 18 hours of infection. These early events shape the ongoing intestinal geography and severity of infection despite the continuous presence of phenotypically virulent luminal bacteria. Mechanistically, induced resistance to A/E lesion de-novo formation is mediated by TLR-MyD88/Trif-dependent signaling and is induced specifically by virulent C. rodentium in a virulence gene-dependent manner. Our data demonstrate that the establishment phase of C. rodentium pathogenesis in vivo is restricted to a very short window of opportunity that determines both disease geography and severity.

D-Alanine-Controlled Transient Intestinal Mono-Colonization with Non-Laboratory-Adapted Commensal E. coli Strain HS.

Soon after birth the mammalian gut microbiota forms a permanent and collectively highly resilient consortium. There is currently no robust method for re-deriving an already microbially colonized individual again-germ-free. We previously developed the in vivo growth-incompetent E. coli K-12 strain HA107 that is auxotrophic for the peptidoglycan components D-alanine (D-Ala) and meso-diaminopimelic acid (Dap) and can be used to transiently associate germ-free animals with live bacteria, without permanent loss of germ-free status. Here we describe the translation of this experimental model from the laboratory-adapted E. coli K-12 prototype to the better gut-adapted commensal strain E. coli HS. In this genetic background it was necessary to complete the D-Ala auxotrophy phenotype by additional knockout of the hypothetical third alanine racemase metC. Cells of the resulting fully auxotrophic strain assembled a peptidoglycan cell wall of normal composition, as long as provided with D-Ala and Dap in the medium, but could not proliferate a single time after D-Ala/Dap removal. Yet, unsupplemented bacteria remained active and were able to complete their cell cycle with fully sustained motility until immediately before autolytic death. Also in vivo, the transiently colonizing bacteria retained their ability to stimulate a live-bacteria-specific intestinal Immunoglobulin (Ig)A response. Full D-Ala auxotrophy enabled rapid recovery to again-germ-free status. E. coli HS has emerged from human studies and genomic analyses as a paradigm of benign intestinal commensal E. coli strains. Its reversibly colonizing derivative may provide a versatile research tool for mucosal bacterial conditioning or compound delivery without permanent colonization.

Microbiota-derived compounds drive steady-state granulopoiesis via MyD88/TICAM signaling.

Neutropenia is probably the strongest known predisposition to infection with otherwise harmless environmental or microbiota-derived species. Because initial swarming of neutrophils at the site of infection occurs within minutes, rather than the hours required to induce "emergency granulopoiesis," the relevance of having high numbers of these cells available at any one time is obvious. We observed that germ-free (GF) animals show delayed clearance of an apathogenic bacterium after systemic challenge. In this article, we show that the size of the bone marrow myeloid cell pool correlates strongly with the complexity of the intestinal microbiota. The effect of colonization can be recapitulated by transferring sterile heat-treated serum from colonized mice into GF wild-type mice. TLR signaling was essential for microbiota-driven myelopoiesis, as microbiota colonization or transferring serum from colonized animals had no effect in GF MyD88(-/-)TICAM1(-/-) mice. Amplification of myelopoiesis occurred in the absence of microbiota-specific IgG production. Thus, very low concentrations of microbial Ags and TLR ligands, well below the threshold required for induction of adaptive immunity, sets the bone marrow myeloid cell pool size. Coevolution of mammals with their microbiota has probably led to a reliance on microbiota-derived signals to provide tonic stimulation to the systemic innate immune system and to maintain vigilance to infection. This suggests that microbiota changes observed in dysbiosis, obesity, or antibiotic therapy may affect the cross talk between hematopoiesis and the microbiota, potentially exacerbating inflammatory or infectious states in the host.