An evolution-based framework for describing human gut bacteria.

The human gut microbiome contains many bacterial strains of the same species ('strain-level variants'). Describing strains in a biologically meaningful way rather than purely taxonomically is an important goal but challenging due to the genetic complexity of strain-level variation. Here, we measured patterns of co-evolution across >7,000 strains spanning the bacterial tree-of-life. Using these patterns as a prior for studying hundreds of gut commensal strains that we isolated, sequenced, and metabolically profiled revealed widespread structure beneath the phylogenetic level of species. Defining strains by their co-evolutionary signatures enabled predicting their metabolic phenotypes and engineering consortia from strain genome content alone. Our findings demonstrate a biologically relevant organization to strain-level variation and motivate a new schema for describing bacterial strains based on their evolutionary history.

Statistically learning the functional landscape of microbial communities.

Microbial consortia exhibit complex functional properties in contexts ranging from soils to bioreactors to human hosts. Understanding how community composition determines function is a major goal of microbial ecology. Here we address this challenge using the concept of community-function landscapes-analogues to fitness landscapes-that capture how changes in community composition alter collective function. Using datasets that represent a broad set of community functions, from production/degradation of specific compounds to biomass generation, we show that statistically inferred landscapes quantitatively predict community functions from knowledge of species presence or absence. Crucially, community-function landscapes allow prediction without explicit knowledge of abundance dynamics or interactions between species and can be accurately trained using measurements from a small subset of all possible community compositions. The success of our approach arises from the fact that empirical community-function landscapes appear to be not rugged, meaning that they largely lack high-order epistatic contributions that would be difficult to fit with limited data. Finally, we show that this observation holds across a wide class of ecological models, suggesting community-function landscapes can be efficiently inferred across a broad range of ecological regimes. Our results open the door to the rational design of consortia without detailed knowledge of abundance dynamics or interactions.

Defining hierarchical protein interaction networks from spectral analysis of bacterial proteomes.

Cellular behaviors emerge from layers of molecular interactions: proteins interact to form complexes, pathways, and phenotypes. We show that hierarchical networks of protein interactions can be defined from the statistical pattern of proteome variation measured across thousands of diverse bacteria and that these networks reflect the emergence of complex bacterial phenotypes. Our results are validated through gene-set enrichment analysis and comparison to existing experimentally derived databases. We demonstrate the biological utility of our approach by creating a model of motility in Pseudomonas aeruginosa and using it to identify a protein that affects pilus-mediated motility. Our method, SCALES (Spectral Correlation Analysis of Layered Evolutionary Signals), may be useful for interrogating genotype-phenotype relationships in bacteria.

Evaluating microbiome-directed fibre snacks in gnotobiotic mice and humans.

Changing food preferences brought about by westernization that have deleterious health effects1,2-combined with myriad forces that are contributing to increased food insecurity-are catalysing efforts to identify more nutritious and affordable foods3. Consumption of dietary fibre can help to prevent cardiovascular disease, type 2 diabetes and obesity4-6. A substantial number of reports have explored the effects of dietary fibre on the gut microbial community7-9. However, the microbiome is complex, dynamic and exhibits considerable intra- and interpersonal variation in its composition and functions. The large number of potential interactions between the components of the microbiome makes it challenging to define the mechanisms by which food ingredients affect community properties. Here we address the question of how foods containing different fibre preparations can be designed to alter functions associated with specific components of the microbiome. Because a marked increase in snack consumption is associated with westernization, we formulated snack prototypes using plant fibres from different sustainable sources that targeted distinct features of the gut microbiomes of individuals with obesity when transplanted into gnotobiotic mice. We used these snacks to supplement controlled diets that were consumed by adult individuals with obesity or who were overweight. Fibre-specific changes in their microbiomes were linked to changes in their plasma proteomes indicative of an altered physiological state.

A Microbiota-Directed Food Intervention for Undernourished Children.

BACKGROUND: More than 30 million children worldwide have moderate acute malnutrition. Current treatments have limited effectiveness, and much remains unknown about the pathogenesis of this condition. Children with moderate acute malnutrition have perturbed development of their gut microbiota. METHODS: In this study, we provided a microbiota-directed complementary food prototype (MDCF-2) or a ready-to-use supplementary food (RUSF) to 123 slum-dwelling Bangladeshi children with moderate acute malnutrition between the ages of 12 months and 18 months. The supplementation was given twice daily for 3 months, followed by 1 month of monitoring. We obtained weight-for-length, weight-for-age, and length-for-age z scores and mid-upper-arm circumference values at baseline and every 2 weeks during the intervention period and at 4 months. We compared the rate of change of these related phenotypes between baseline and 3 months and between baseline and 4 months. We also measured levels of 4977 proteins in plasma and 209 bacterial taxa in fecal samples. RESULTS: A total of 118 children (59 in each study group) completed the intervention. The rates of change in the weight-for-length and weight-for-age z scores are consistent with a benefit of MDCF-2 on growth over the course of the study, including the 1-month follow-up. Receipt of MDCF-2 was linked to the magnitude of change in levels of 70 plasma proteins and of 21 associated bacterial taxa that were positively correlated with the weight-for-length z score (P<0.001 for comparisons of both protein and bacterial taxa). These proteins included mediators of bone growth and neurodevelopment. CONCLUSIONS: These findings provide support for MDCF-2 as a dietary supplement for young children with moderate acute malnutrition and provide insight into mechanisms by which this targeted manipulation of microbiota components may be linked to growth. (Supported by the Bill and Melinda Gates Foundation and the National Institutes of Health; ClinicalTrials.gov number, NCT04015999.).

Identifying determinants of bacterial fitness in a model of human gut microbial succession.

Human gut microbiota development has been associated with healthy growth but understanding the determinants of community assembly and composition is a formidable challenge. We cultured bacteria from serially collected fecal samples from a healthy infant; 34 sequenced strains containing 103,102 genes were divided into two consortia representing earlier and later stages in community assembly during the first six postnatal months. The two consortia were introduced alone (singly), or sequentially in different order, or simultaneously into young germ-free mice fed human infant formula. The pattern of fitness of bacterial strains observed across the different colonization conditions indicated that later-phase strains substantially outcompete earlier-phase strains, although four early-phase members persist. Persistence was not determined by order of introduction, suggesting that priority effects are not prominent in this model. To characterize succession in the context of the metabolic potential of consortium members, we performed in silico reconstructions of metabolic pathways involved in carbohydrate utilization and amino acid and B-vitamin biosynthesis, then quantified the fitness (abundance) of strains in serially collected fecal samples and their transcriptional responses to different histories of colonization. Applying feature-reduction methods disclosed a set of metabolic pathways whose presence and/or expression correlates with strain fitness and that enable early-stage colonizers to survive during introduction of later colonizers. The approach described can be used to test the magnitude of the contribution of identified metabolic pathways to fitness in different community contexts, study various ecological processes thought to govern community assembly, and facilitate development of microbiota-directed therapeutics.

A sparse covarying unit that describes healthy and impaired human gut microbiota development.

Characterizing the organization of the human gut microbiota is a formidable challenge given the number of possible interactions between its components. Using a statistical approach initially applied to financial markets, we measured temporally conserved covariance among bacterial taxa in the microbiota of healthy members of a Bangladeshi birth cohort sampled from 1 to 60 months of age. The results revealed an "ecogroup" of 15 covarying bacterial taxa that provide a concise description of microbiota development in healthy children from this and other low-income countries, and a means for monitoring community repair in undernourished children treated with therapeutic foods. Features of ecogroup population dynamics were recapitulated in gnotobiotic piglets as they transitioned from exclusive milk feeding to a fully weaned state consuming a representative Bangladeshi diet.

Origins of Allostery and Evolvability in Proteins: A Case Study.

Proteins display the capacity for adaptation to new functions, a property critical for evolvability. But what structural principles underlie the capacity for adaptation? Here, we show that adaptation to a physiologically distinct class of ligand specificity in a PSD95, DLG1, ZO-1 (PDZ) domain preferentially occurs through class-bridging intermediate mutations located distant from the ligand-binding site. These mutations provide a functional link between ligand classes and demonstrate the principle of "conditional neutrality" in mediating evolutionary adaptation. Structures show that class-bridging mutations work allosterically to open up conformational plasticity at the active site, permitting novel functions while retaining existing function. More generally, the class-bridging phenotype arises from mutations in an evolutionarily conserved network of coevolving amino acids in the PDZ family (the sector) that connects the active site to distant surface sites. These findings introduce the concept that allostery in proteins could have its origins not in protein function but in the capacity to adapt.

Photofragment imaging study of the CH2CCH2OH radical intermediate of the OH+allene reaction.

These velocity map imaging experiments characterize the photolytic generation of one of the two radical intermediates formed when OH reacts via an addition mechanism with allene. The CH2CCH2OH radical intermediate is generated photolytically from the photodissociation of 2-chloro-2-propen-1-ol at 193 nm. Detecting the Cl atoms using [2+1] resonance-enhanced multiphoton ionization evidences an isotropic angular distribution for the Cl+CH2CCH2OH photofragments, a spin-orbit branching ratio for Cl(2P1/2):Cl(2P3/2) of 0.28, and a bimodal recoil kinetic energy distribution. Conservation of momentum and energy allows us to determine from this data the internal energy distribution of the nascent CH2CCH2OH radical cofragment. To assess the possible subsequent decomposition pathways of this highly vibrationally excited radical intermediate, we include electronic structure calculations at the G3//B3LYP level of theory. They predict the isomerization and dissociation transition states en route from the initial CH2CCH2OH radical intermediate to the three most important product channels for the OH+allene reaction expected from this radical intermediate: formaldehyde+C2H3, H+acrolein, and ethene+CHO. We also calculate the intermediates and transition states en route from the other radical adduct, formed by addition of the OH to the center carbon of allene, to the ketene+CH3 product channel. We compare our results to a previous theoretical study of the O+allyl reaction conducted at the CBS-QB3 level of theory, as the two reactions include several common intermediates.