Genetic determinants of in vivo fitness and diet responsiveness in multiple human gut Bacteroides

Diet shapes host and gut microbe fitness The human gut microbiota is hugely diverse, with many strain variants having a multiplicity of effects on host metabolism and immunity. To define some of these functions, Wu et al. made libraries of mutants of Bacteroides species known for their capacity to process otherwise intractable dietary fiber. Germ-free mice colonized with defined gut microbiota communities containing the mutants were fed specific diets containing different ratios of fat and fiber. Genes, strains, and species were identified that were associated with specific metabolic pathways. The community responses to dietary shifts were manipulated in an attempt to characterize species for their probiotic or therapeutic potential. Science, this issue 10.1126/science.aac5992> To design probiotics, gut microbe fitness determinants and niches were characterized and responses to dietary changes monitored. INTRODUCTION Relatively little is known about the genetic factors that allow members of the human gut microbiota to occupy their niches. Identification of these factors is important for understanding mechanisms that determine microbiota assembly and perturbation through diet, disease, and clinical treatments. Discovery of these factors should enable new approaches for intervening therapeutically in the functional properties of the human gut microbiota. We present a generalizable approach by which to identify fitness determinants for multiple bacterial strains simultaneously in a model human gut microbiota, obtain gene-level characterization of responses to diet change, and design prebiotics for precision microbiota manipulation. RATIONALE We developed a method—multi-taxon INsertion Sequencing (INSeq)—for monitoring the behavior of tens of thousands of transposon (Tn) mutants of multiple bacterial species and strains simultaneously in the guts of gnotobiotic mice. We focused on four prominent human gut Bacteroides: one strain of B. cellulosilyticus, one strain of B. ovatus, and two strains of B. thetaiotaomicron. INSeq libraries, each composed of 87,000 to 167,000 isogenic Tn mutant strains, were produced (single site of Tn insertion per mutant strain; a total of 11 to 26 Tn insertions represented in the library per gene; and 82 to 92% genes covered per genome). The four mutant libraries were introduced into germ-free mice together with 11 wild-type species commonly present in the human gut microbiota. Animals were given a diet rich in fat and simple sugars but devoid of complex polysaccharides [diet 1 (D1)] or one rich in plant polysaccharides and low in fat (D2), either monotonously or in the sequence D1-D2-D1 or D2-D1-D2. Wecalculated a “fitness index” for each gene on the basis of the relative abundance of its INSeq reads in the fecal or cecal microbiota compared with the input library. In vivo INSeq data were correlated with INSeq data generated from organisms cultured under defined in vitro conditions; microbial RNA-seq profiling of the community’s metatranscriptome; and reconstructions of metabolic pathways, regulons, and polysaccharide utilization loci. On the basis of the results, we designed a prebiotic intervention. RESULTS Multi-taxon INSeq (i) provided a digital readout of the remarkably consistent pattern of community assembly; (ii) identified shared as well as species-, strain-, and diet-specific fitness determinants associated with a variety of metabolic or nutrient processing pathways, including those involving amino acids, carbohydrates, and vitamins/cofactors; (iii) enabled quantitative gene-level measurement of the resilience of the responses to diet perturbations; (iv) revealed that arabinoxylan, the most common hemicellulose in cereals, could be used to deliberately manipulate the representation of Bacteroides cellulosilyticus; and (v) defined the niche adjustments of this and the other Bacteroides to arabinoxylan supplementation of the high-fat diet. CONCLUSION In principle, the approach described can be used to obtain a more comprehensive understanding of how host genotype, diet, physiologic, metabolic, and immune factors, as well as pathologic states, affect niches in gut and nongut habitats, as well as to facilitate development of therapeutic interventions for modifying community structure/function. Identification of a prebiotic that increases the abundance of B. cellulosilyticus. (Left) The four mutant libraries were pooled together with 11 other phylogenetically diverse wild-type strains, and this consortium, representing an artificial human gut microbiota, was introduced into germ-free mice. Community assembly, the effects of diet, and recovery from diet oscillations were characterized at a community, strain, and gene level in these gnotobiotic animals by use of multi-taxon INSeq. (Middle) Multi-taxon INSeq revealed an arabinoxylan utilization locus in B. cellulosilyticus that is critical for the organism’s fitness in the high-fat/simple-sugar diet (D1) context but not in the D2 context. A homologous arabinoxylan utilization locus in B. ovatus was not a fitness determinant with either diet. (Right) Consistent with this finding, supplementation of drinking water with arabinoxylan in mice consuming D1 selectively increased the abundance of B. cellulosilyticus but not B. ovatus. Libraries of tens of thousands of transposon mutants generated from each of four human gut Bacteroides strains, two representing the same species, were introduced simultaneously into gnotobiotic mice together with 11 other wild-type strains to generate a 15-member artificial human gut microbiota. Mice received one of two distinct diets monotonously, or both in different ordered sequences. Quantifying the abundance of mutants in different diet contexts allowed gene-level characterization of fitness determinants, niche, stability, and resilience and yielded a prebiotic (arabinoxylan) that allowed targeted manipulation of the community. The approach described is generalizable and should be useful for defining mechanisms critical for sustaining and/or approaches for deliberately reconfiguring the highly adaptive and durable relationship between the human gut microbiota and host in ways that promote wellness.

[1]  J. Gordon,et al.  Coordinate Regulation of Glycan Degradation and Polysaccharide Capsule Biosynthesis by a Prominent Human Gut Symbiont , 2009, The Journal of Biological Chemistry.

[2]  M. Gelfand,et al.  Comparative Genomics of the Vitamin B12 Metabolism and Regulation in Prokaryotes* , 2003, Journal of Biological Chemistry.

[3]  J. W. Campbell,et al.  Experimental Determination and System Level Analysis of Essential Genes in Escherichia coli MG1655 , 2003, Journal of bacteriology.

[4]  B. Henrissat,et al.  Glycan complexity dictates microbial resource allocation in the large intestine , 2015, Nature Communications.

[5]  B. Haas,et al.  A Catalog of Reference Genomes from the Human Microbiome , 2010, Science.

[6]  R. Geyer,et al.  Structural characterization of N-glycans from the freshwater snail Biomphalaria glabrata cross-reacting with Schistosoma mansoni glycoconjugates. , 2007, Glycobiology.

[7]  Bernard Henrissat,et al.  Effects of Diet on Resource Utilization by a Model Human Gut Microbiota Containing Bacteroides cellulosilyticus WH2, a Symbiont with an Extensive Glycobiome , 2013, PLoS biology.

[8]  R. Mackie,et al.  Xylan degradation, a metabolic property shared by rumen and human colonic Bacteroidetes , 2011, Molecular microbiology.

[9]  A A Mironov,et al.  [Software for analyzing bacterial genomes]. , 2000, Molekuliarnaia biologiia.

[10]  William J. Riehl,et al.  RegPrecise 3.0 – A resource for genome-scale exploration of transcriptional regulation in bacteria , 2013, BMC Genomics.

[11]  A. A. Mironov,et al.  Software for analysis of bacterial genomes , 2000, Molecular Biology.

[12]  J. Mullikin,et al.  SSAHA: a fast search method for large DNA databases. , 2001, Genome research.

[13]  J. Gordon,et al.  Message from a human gut symbiont: sensitivity is a prerequisite for sharing. , 2004, Trends in microbiology.

[14]  J. Clemente,et al.  The Long-Term Stability of the Human Gut Microbiota , 2013 .

[15]  Aaron A. Klammer,et al.  Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data , 2013, Nature Methods.

[16]  Abigail A. Salyers,et al.  Characterization of Four Outer Membrane Proteins Involved in Binding Starch to the Cell Surface ofBacteroides thetaiotaomicron , 2000, Journal of bacteriology.

[17]  W. Huber,et al.  which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets , 2011 .

[18]  J. Faith,et al.  Predicting a Human Gut Microbiota’s Response to Diet in Gnotobiotic Mice , 2011, Science.

[19]  N. Perna,et al.  progressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement , 2010, PloS one.

[20]  D. Haller,et al.  Transcriptome analysis of Enterococcus faecalis toward its adaption to surviving in the mouse intestinal tract , 2014, Archives of Microbiology.

[21]  Sara M. Johnson,et al.  Sap Transporter Mediated Import and Subsequent Degradation of Antimicrobial Peptides in Haemophilus , 2011, PLoS pathogens.

[22]  T. Cullen,et al.  Antimicrobial peptide resistance mediates resilience of prominent gut commensals during inflammation , 2015, Science.

[23]  Dmitry A Rodionov,et al.  Comparative genomic reconstruction of transcriptional regulatory networks in bacteria. , 2007, Chemical reviews.

[24]  Steven J. M. Jones,et al.  Circos: an information aesthetic for comparative genomics. , 2009, Genome research.

[25]  P. Bork,et al.  A human gut microbial gene catalogue established by metagenomic sequencing , 2010, Nature.

[26]  Rob Knight,et al.  Identifying genetic determinants needed to establish a human gut symbiont in its habitat. , 2009, Cell host & microbe.

[27]  Bernard Henrissat,et al.  Recognition and Degradation of Plant Cell Wall Polysaccharides by Two Human Gut Symbionts , 2011, PLoS biology.

[28]  Adam Godzik,et al.  Polysaccharides utilization in human gut bacterium Bacteroides thetaiotaomicron: comparative genomics reconstruction of metabolic and regulatory networks , 2013, BMC Genomics.

[29]  J. Gordon,et al.  Gnotobiotic mouse model of phage–bacterial host dynamics in the human gut , 2013, Proceedings of the National Academy of Sciences.

[30]  J. Gordon,et al.  Metabolic niche of a prominent sulfate-reducing human gut bacterium , 2013, Proceedings of the National Academy of Sciences.

[31]  A. Salyers,et al.  Characterization of four outer membrane proteins that play a role in utilization of starch by Bacteroides thetaiotaomicron , 1997, Journal of bacteriology.

[32]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[33]  C. S. Holling Resilience and Stability of Ecological Systems , 1973 .

[34]  P. Degnan,et al.  Human gut microbes use multiple transporters to distinguish vitamin B₁₂ analogs and compete in the gut. , 2014, Cell host & microbe.

[35]  Felipe Garbelini Marques,et al.  Dual diaminopimelate biosynthesis pathways in Bacteroides fragilis and Clostridium thermocellum. , 2011, Biochimica et biophysica acta.

[36]  L. Ursell,et al.  Genetically dictated change in host mucus carbohydrate landscape exerts a diet-dependent effect on the gut microbiota , 2013, Proceedings of the National Academy of Sciences.

[37]  J. Gordon,et al.  Identifying microbial fitness determinants by insertion sequencing using genome-wide transposon mutant libraries , 2011, Nature Protocols.

[38]  A. Salyers,et al.  A Bacteroides thetaiotaomicron outer membrane protein that is essential for utilization of maltooligosaccharides and starch , 1996, Journal of bacteriology.

[39]  J. Gordon,et al.  Mucosal glycan foraging enhances fitness and transmission of a saccharolytic human gut bacterial symbiont. , 2008, Cell host & microbe.

[40]  Lynn K. Carmichael,et al.  A Genomic View of the Human-Bacteroides thetaiotaomicron Symbiosis , 2003, Science.

[41]  Milton H. Saier,et al.  The Transporter Classification Database , 2013, Nucleic Acids Res..

[42]  Benjamin P. Westover,et al.  Glycan Foraging in Vivo by an Intestine-Adapted Bacterial Symbiont , 2005, Science.

[43]  T. Snijders Multivariate Statistics and Matrices in Statistics , 1995 .

[44]  M. Gelfand,et al.  Comparative Genomics of Thiamin Biosynthesis in Procaryotes , 2002, The Journal of Biological Chemistry.

[45]  Naryttza N. Diaz,et al.  The Subsystems Approach to Genome Annotation and its Use in the Project to Annotate 1000 Genomes , 2005, Nucleic acids research.