Gut bacteria that prevent growth impairments transmitted by microbiota from malnourished children

Undernourished children exhibit impaired development of their gut microbiota. Transplanting microbiota from 6- and 18-month-old healthy or undernourished Malawian donors into young germ-free mice that were fed a Malawian diet revealed that immature microbiota from undernourished infants and children transmit impaired growth phenotypes. The representation of several age-discriminatory taxa in recipient animals correlated with lean body mass gain; liver, muscle, and brain metabolism; and bone morphology. Mice were cohoused shortly after receiving microbiota from healthy or severely stunted and underweight infants; age- and growth-discriminatory taxa from the microbiota of the former were able to invade that of the latter, which prevented growth impairments in recipient animals. Adding two invasive species, Ruminococcus gnavus and Clostridium symbiosum, to the microbiota from undernourished donors also ameliorated growth and metabolic abnormalities in recipient animals. These results provide evidence that microbiota immaturity is causally related to undernutrition and reveal potential therapeutic targets and agents.

RegPrecise 3.0--a resource for genome-scale exploration of transcriptional regulation in bacteria

Background

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in prokaryotes is one of the critical tasks of modern genomics. Bacteria from different taxonomic groups, whose lifestyles and natural environments are substantially different, possess highly diverged transcriptional regulatory networks. The comparative genomics approaches are useful for in silico reconstruction of bacterial regulons and networks operated by both transcription factors (TFs) and RNA regulatory elements (riboswitches).

Description

RegPrecise (http://regprecise.lbl.gov) is a web resource for collection, visualization and analysis of transcriptional regulons reconstructed by comparative genomics. We significantly expanded a reference collection of manually curated regulons we introduced earlier. RegPrecise 3.0 provides access to inferred regulatory interactions organized by phylogenetic, structural and functional properties. Taxonomy-specific collections include 781 TF regulogs inferred in more than 160 genomes representing 14 taxonomic groups of Bacteria. TF-specific collections include regulogs for a selected subset of 40 TFs reconstructed across more than 30 taxonomic lineages. Novel collections of regulons operated by RNA regulatory elements (riboswitches) include near 400 regulogs inferred in 24 bacterial lineages. RegPrecise 3.0 provides four classifications of the reference regulons implemented as controlled vocabularies: 55 TF protein families; 43 RNA motif families; ~150 biological processes or metabolic pathways; and ~200 effectors or environmental signals. Genome-wide visualization of regulatory networks and metabolic pathways covered by the reference regulons are available for all studied genomes. A separate section of RegPrecise 3.0 contains draft regulatory networks in 640 genomes obtained by an conservative propagation of the reference regulons to closely related genomes.

Conclusions

RegPrecise 3.0 gives access to the transcriptional regulons reconstructed in bacterial genomes. Analytical capabilities include exploration of: regulon content, structure and function; TF binding site motifs; conservation and variations in genome-wide regulatory networks across all taxonomic groups of Bacteria. RegPrecise 3.0 was selected as a core resource on transcriptional regulation of the Department of Energy Systems Biology Knowledgebase, an emerging software and data environment designed to enable researchers to collaboratively generate, test and share new hypotheses about gene and protein functions, perform large-scale analyses, and model interactions in microbes, plants, and their communities.

Effects of microbiota-directed foods in gnotobiotic animals and undernourished children

To examine the contributions of impaired gut microbial community development to childhood undernutrition, we combined metabolomic and proteomic analyses of plasma samples with metagenomic analyses of fecal samples to characterize the biological state of Bangladeshi children with severe acute malnutrition (SAM) as they transitioned, after standard treatment, to moderate acute malnutrition (MAM) with persistent microbiota immaturity. Host and microbial effects of microbiota-directed complementary food (MDCF) prototypes targeting weaning-phase bacterial taxa underrepresented in SAM and MAM microbiota were characterized in gnotobiotic mice and gnotobiotic piglets colonized with age- and growth-discriminatory bacteria. A randomized, double-blind controlled feeding study identified a lead MDCF that changes the abundances of targeted bacteria and increases plasma biomarkers and mediators of growth, bone formation, neurodevelopment, and immune function in children with MAM.

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.

Micronutrient Requirements and Sharing Capabilities of the Human Gut Microbiome

The human gut microbiome harbors a diverse array of metabolic pathways contributing to its development and homeostasis via a complex web of diet-dependent metabolic interactions within the microbial community and host. Genomics-based reconstruction and predictive modeling of these interactions would provide a framework for diagnostics and treatment of dysbiosis-related syndromes via rational selection of therapeutic prebiotics and dietary nutrients. Of particular interest are micronutrients, such as B-group vitamins, precursors of indispensable metabolic cofactors, that are produced de novo by some gut bacteria (prototrophs) but must be provided exogenously in the diet for many other bacterial species (auxotrophs) as well as for the mammalian host. Cross-feeding of B vitamins between prototrophic and auxotrophic species is expected to strongly contribute to the homeostasis of microbial communities in the distal gut given the efficient absorption of dietary vitamins in the upper gastrointestinal tract. To confidently estimate the balance of microbiome micronutrient biosynthetic capabilities and requirements using available genomic data, we have performed a subsystems-based reconstruction of biogenesis, salvage and uptake for eight B vitamins (B1, B2, B3, B5, B6, B7, B9, and B12) and queuosine (essential factor in tRNA modification) over a reference set of 2,228 bacterial genomes representing 690 cultured species of the human gastrointestinal microbiota. This allowed us to classify the studied organisms with respect to their pathway variants and infer their prototrophic vs. auxotrophic phenotypes. In addition to canonical vitamin pathways, several conserved partial pathways were identified pointing to alternative routes of syntrophic metabolism and expanding a microbial vitamin "menu" by several pathway intermediates (vitamers) such as thiazole, quinolinate, dethiobiotin, pantoate. A cross-species comparison was applied to assess the extent of conservation of vitamin phenotypes at distinct taxonomic levels (from strains to families). The obtained reference collection combined with 16S rRNA gene-based phylogenetic profiles was used to deduce phenotype profiles of the human gut microbiota across in two large cohorts. This analysis provided the first estimate of B-vitamin requirements, production and sharing capabilities in the human gut microbiome establishing predictive phenotype profiling as a new approach to classification of microbiome samples. Future expansion of our reference genomic collection of metabolic phenotypes will allow further improvement in coverage and accuracy of predictive phenotype profiling of the human microbiome.

A novel antibiotic class targeting the lipopolysaccharide transporter

Carbapenem-resistant Acinetobacter baumannii (CRAB) has emerged as a major global pathogen with limited treatment options1. No new antibiotic chemical class with activity against A. baumannii has reached patients in over 50 years1. Here we report the identification and optimization of tethered macrocyclic peptide (MCP) antibiotics with potent antibacterial activity against CRAB. The mechanism of action of this molecule class involves blocking the transport of bacterial lipopolysaccharide from the inner membrane to its destination on the outer membrane, through inhibition of the LptB2FGC complex. A clinical candidate derived from the MCP class, zosurabalpin (RG6006), effectively treats highly drug-resistant contemporary isolates of CRAB both in vitro and in mouse models of infection, overcoming existing antibiotic resistance mechanisms. This chemical class represents a promising treatment paradigm for patients with invasive infections due to CRAB, for whom current treatment options are inadequate, and additionally identifies LptB2FGC as a tractable target for antimicrobial drug development.

Duodenal Microbiota in Stunted Undernourished Children with Enteropathy

BACKGROUND

Environmental enteric dysfunction (EED) is an enigmatic disorder of the small intestine that is postulated to play a role in childhood undernutrition, a pressing global health problem. Defining the incidence of this disorder, its pathophysiological features, and its contribution to impaired linear and ponderal growth has been hampered by the difficulty in directly sampling the small intestinal mucosa and microbial community (microbiota).

METHODS

In this study, among 110 young children (mean age, 18 months) with linear growth stunting who were living in an urban slum in Dhaka, Bangladesh, and had not benefited from a nutritional intervention, we performed endoscopy in 80 children who had biopsy-confirmed EED and available plasma and duodenal samples. We quantified the levels of 4077 plasma proteins and 2619 proteins in duodenal biopsy samples obtained from these children. The levels of bacterial strains in microbiota recovered from duodenal aspirate from each child were determined with the use of culture-independent methods. In addition, we obtained 21 plasma samples and 27 fecal samples from age-matched healthy children living in the same area. Young germ-free mice that had been fed a Bangladeshi diet were colonized with bacterial strains cultured from the duodenal aspirates.

RESULTS

Of the bacterial strains that were obtained from the children, the absolute levels of a shared group of 14 taxa (which are not typically classified as enteropathogens) were negatively correlated with linear growth (length-for-age z score, r = -0.49; P = 0.003) and positively correlated with duodenal proteins involved in immunoinflammatory responses. The representation of these 14 duodenal taxa in fecal microbiota was significantly different from that in samples obtained from healthy children (P<0.001 by permutational multivariate analysis of variance). Enteropathy of the small intestine developed in gnotobiotic mice that had been colonized with cultured duodenal strains obtained from children with EED.

CONCLUSIONS

These results provide support for a causal relationship between growth stunting and components of the small intestinal microbiota and enteropathy and offer a rationale for developing therapies that target these microbial contributions to EED. (Funded by the Bill and Melinda Gates Foundation and others; ClinicalTrials.gov number, NCT02812615.).

B-Vitamin Sharing Promotes Stability of Gut Microbial Communities

Cross-feeding on intermediary and end-point metabolites plays an important role in the dynamic interactions of host-associated microbial communities. While gut microbiota possess inherent resilience to perturbation, variations in the intake of certain nutrients may lead to changes in the community composition with potential consequences on host physiology. Syntrophic relationships and mutualism at the level of major carbon and energy sources have been documented, however, relatively little is known about metabolic interactions involving micronutrients, such as B-vitamins, biosynthetic precursors of essential cofactors in the mammalian host and numerous members of the gut microbiota alike. In silico genomic reconstruction and prediction of community-wide metabolic phenotypes for eight major B-vitamins (B1, B2, B3, B5, B6, B7, B9, and B12), suggests that a significant fraction of microbial gut communities (>20% by abundance) are represented by auxotrophic species whose viability is strictly dependent on acquiring one or more B-vitamins from diet and/or prototrophic microbes via committed salvage pathways. Here, we report the stability of gut microbiota using humanized gnotobiotic mice and in vitro anaerobic fecal culture in the context of extreme variations of dietary B-vitamin supply as revealed by phylotype-to-phenotype prediction from 16S rRNA profiling and metabolomic measurements. The observed nearly unaltered relative abundance of auxotrophic species in gut communities in the face of diet or media lacking B-vitamins or containing them in great excess (∼30-fold above normal) points to a strong contribution of metabolic cooperation (B-vitamin exchange and sharing) to the stability of gut bacterial populations.

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.

Genomic reconstruction of the transcriptional regulatory network in Bacillus subtilis

The adaptation of microorganisms to their environment is controlled by complex transcriptional regulatory networks (TRNs), which are still only partially understood even for model species. Genome scale annotation of regulatory features of genes and TRN reconstruction are challenging tasks of microbial genomics. We used the knowledge-driven comparative-genomics approach implemented in the RegPredict Web server to infer TRN in the model Gram-positive bacterium Bacillus subtilis and 10 related Bacillales species. For transcription factor (TF) regulons, we combined the available information from the DBTBS database and the literature with bioinformatics tools, allowing inference of TF binding sites (TFBSs), comparative analysis of the genomic context of predicted TFBSs, functional assignment of target genes, and effector prediction. For RNA regulons, we used known RNA regulatory motifs collected in the Rfam database to scan genomes and analyze the genomic context of new RNA sites. The inferred TRN in B. subtilis comprises regulons for 129 TFs and 24 regulatory RNA families. First, we analyzed 66 TF regulons with previously known TFBSs in B. subtilis and projected them to other Bacillales genomes, resulting in refinement of TFBS motifs and identification of novel regulon members. Second, we inferred motifs and described regulons for 28 experimentally studied TFs with previously unknown TFBSs. Third, we discovered novel motifs and reconstructed regulons for 36 previously uncharacterized TFs. The inferred collection of regulons is available in the RegPrecise database (http://regprecise.lbl.gov/) and can be used in genetic experiments, metabolic modeling, and evolutionary analysis.

Transcriptional Regulation of Carbohydrate Utilization Pathways in the Bifidobacterium Genus

Bifidobacteria, which represent common commensals of mammalian gut, are believed to have positive effects on human health. The influence of certain non-digestible carbohydrates (and their use as so-called prebiotics) on growth and metabolic activity of bifidobacteria is of increasing interest; however, mechanisms of transcriptional control of carbohydrate metabolism are poorly understood in these species. We used a comparative genomics approach to reconstruct carbohydrate utilization pathways and transcriptional regulons in 10 Bifidobacterium genomes. Analysis of regulatory gene regions revealed candidate DNA motifs and reconstructed regulons for 268 transcription factors from the LacI, ROK, DeoR, AraC, GntR, and TetR families that form 64 orthologous groups of regulators. Most of the reconstructed regulons are local and control specific catabolic pathways for host- and diet-derived glycans and monosaccharides. Mosaic distributions of many of these local regulators across Bifidobacterium species correlate with distribution of corresponding catabolic pathways. In contrast, the maltose, galactose, sucrose, and fructose regulons, as well as a novel global LacI-family regulator that is predicted to control the central carbohydrate metabolism and arabinose catabolism genes, are universally present in all 10 studied bifidobacteria. A novel group of TetR-family regulators presumably controls the glucoside and galactoside utilization pathways. Paralogs of the ribose repressor RbsR control the pyrimidine nucleoside utilization genes. Multiple paralogs of the maltose regulator MalR co-regulate large sets of genes involved in maltodextrin utilization. The inferred metabolic regulons provide new insights on diverse carbohydrate utilization networks in bifidobacteria that can be employed in metabolic modeling, phenotype prediction and the rational development of novel prebiotics.

Tetracenomycin X inhibits translation by binding within the ribosomal exit tunnel

The increase in multi-drug resistant pathogenic bacteria is making our current arsenal of clinically used antibiotics obsolete, highlighting the urgent need for new lead compounds with distinct target binding sites to avoid cross-resistance. Here we report that the aromatic polyketide antibiotic tetracenomycin (TcmX) is a potent inhibitor of protein synthesis, and does not induce DNA damage as previously thought. Despite the structural similarity to the well-known translation inhibitor tetracycline, we show that TcmX does not interact with the small ribosomal subunit, but rather binds to the large subunit, within the polypeptide exit tunnel. This previously unappreciated binding site is located adjacent to the macrolide-binding site, where TcmX stacks on the noncanonical basepair formed by U1782 and U2586 of the 23S ribosomal RNA. Although the binding site is distinct from the macrolide antibiotics, our results indicate that like macrolides, TcmX allows translation of short oligopeptides before further translation is blocked.

Control of proteobacterial central carbon metabolism by the HexR transcriptional regulator: a case study in Shewanella oneidensis

Bacteria exploit multiple mechanisms for controlling central carbon metabolism (CCM). Thus, a bioinformatic analysis together with some experimental data implicated the HexR transcriptional factor as a global CCM regulator in some lineages of Gammaproteobacteria operating as a functional replacement of the Cra regulator characteristic of Enterobacteriales. In this study, we combined a large scale comparative genomic reconstruction of HexR-controlled regulons in 87 species of Proteobacteria with the detailed experimental analysis of the HexR regulatory network in the Shewanella oneidensis model system. Although nearly all of the HexR-controlled genes are associated with CCM, remarkable variations were revealed in the scale (from 1 to 2 target operons in Enterobacteriales up to 20 operons in Aeromonadales) and gene content of HexR regulons between 11 compared lineages. A predicted 17-bp pseudo-palindrome with a consensus tTGTAATwwwATTACa was confirmed as a HexR-binding motif for 15 target operons (comprising 30 genes) by in vitro binding assays. The negative effect of the key CCM intermediate, 2-keto-3-deoxy-6-phosphogluconate, on the DNA-regulator complex formation was verified. A dual mode of HexR action on various target promoters, repression of genes involved in catabolic pathways and activation of gluconeogenic genes, was for the first time predicted by the bioinformatic analysis and experimentally verified by changed gene expression pattern in S. oneidensis ΔhexR mutant. Phenotypic profiling revealed the inability of this mutant to grow on lactate or pyruvate as a single carbon source. A comparative metabolic flux analysis of wild-type and mutant strains of S. oneidensis using [(13)C]lactate labeling and GC-MS analysis confirmed the hypothesized HexR role as a master regulator of gluconeogenic flux from pyruvate via the transcriptional activation of phosphoenolpyruvate synthase (PpsA).

Combined Prebiotic and Microbial Intervention Improves Oral Cholera Vaccination Responses in a Mouse Model of Childhood Undernutrition

Undernourished children in low-income countries often exhibit poor responses to oral vaccination. Perturbed microbiota development is linked to undernutrition, but whether and how microbiota changes affect vaccine responsiveness remains unclear. Here, we show that gnotobiotic mice colonized with microbiota from undernourished Bangladeshi children and fed a Bangladeshi diet exhibited microbiota-dependent differences in mucosal IgA responses to oral vaccination with cholera toxin (CT). Supplementation with a nutraceutical consisting of spirulina, amaranth, flaxseed, and micronutrients augmented CT-IgA production. Mice initially colonized with a microbiota associated with poor CT responses exhibited improved immunogenicity upon invasion of bacterial taxa from cagemates colonized with a more "responsive" microbiota. Additionally, a consortium of five cultured bacterial invaders conferred augmented CT-IgA responses in mice fed the supplemented diet and colonized with the "hypo-responsive" community. These results provide preclinical proof-of-concept that diet and microbiota influence mucosal immune responses to CT vaccination and identify a candidate synbiotic formulation.

Comparative genomics of metabolic capacities of regulons controlled by cis-regulatory RNA motifs in bacteria

Background

In silico comparative genomics approaches have been efficiently used for functional prediction and reconstruction of metabolic and regulatory networks. Riboswitches are metabolite-sensing structures often found in bacterial mRNA leaders controlling gene expression on transcriptional or translational levels.

An increasing number of riboswitches and other cis-regulatory RNAs have been recently classified into numerous RNA families in the Rfam database. High conservation of these RNA motifs provides a unique advantage for their genomic identification and comparative analysis.

Results

A comparative genomics approach implemented in the RegPredict tool was used for reconstruction and functional annotation of regulons controlled by RNAs from 43 Rfam families in diverse taxonomic groups of Bacteria. The inferred regulons include ~5200 cis-regulatory RNAs and more than 12000 target genes in 255 microbial genomes. All predicted RNA-regulated genes were classified into specific and overall functional categories. Analysis of taxonomic distribution of these categories allowed us to establish major functional preferences for each analyzed cis-regulatory RNA motif family. Overall, most RNA motif regulons showed predictable functional content in accordance with their experimentally established effector ligands. Our results suggest that some RNA motifs (including thiamin pyrophosphate and cobalamin riboswitches that control the cofactor metabolism) are widespread and likely originated from the last common ancestor of all bacteria. However, many more analyzed RNA motifs are restricted to a narrow taxonomic group of bacteria and likely represent more recent evolutionary innovations.

Conclusions

The reconstructed regulatory networks for major known RNA motifs substantially expand the existing knowledge of transcriptional regulation in bacteria. The inferred regulons can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. The obtained genome-wide collection of reference RNA motif regulons is available in the RegPrecise database (http://regprecise.lbl.gov/).

Bioremediation of a Common Product of Food Processing by a Human Gut Bacterium

Dramatic increases in processed food consumption represent a global health threat. Maillard reaction products (MRPs), which are common in processed foods, form upon heat-induced reaction of amino acids with reducing sugars and include advanced glycation end products with deleterious health effects. To examine how processed foods affect the microbiota, we fed gnotobiotic mice, colonized with 54 phylogenetically diverse human gut bacterial strains, defined sugar-rich diets containing whey as the protein source or a matched amino acid mixture. Whey or ϵ-fructoselysine, an MRP in whey and many processed foods, selectively increases Collinsella intestinalis absolute abundance and induces Collinsella expression of genomic loci directing import and metabolism of ϵ-fructoselysine to innocuous products. This locus is repressed by glucose in C. aerofaciens, whose abundance decreases with whey, but is not repressed in C. intestinalis. Identifying gut organisms responding to and degrading potentially harmful processed food components has implications for food science, microbiome science, and public health.

N-acetylgalactosamine utilization pathway and regulon in proteobacteria: genomic reconstruction and experimental characterization in Shewanella

We used a comparative genomics approach to reconstruct the N-acetyl-d-galactosamine (GalNAc) and galactosamine (GalN) utilization pathways and transcriptional regulons in Proteobacteria. The reconstructed GalNAc/GalN utilization pathways include multiple novel genes with specific functional roles. Most of the pathway variations were attributed to the amino sugar transport, phosphorylation, and deacetylation steps, whereas the downstream catabolic enzymes in the pathway were largely conserved. The predicted GalNAc kinase AgaK, the novel variant of GalNAc-6-phosphate deacetylase AgaA(II) and the GalN-6-phosphate deaminase AgaS from Shewanella sp. ANA-3 were validated in vitro using individual enzymatic assays and reconstitution of the three-step pathway. By using genetic techniques, we confirmed that AgaS but not AgaI functions as the main GalN-6-P deaminase in the GalNAc/GalN utilization pathway in Escherichia coli. Regulons controlled by AgaR repressors were reconstructed by bioinformatics in most proteobacterial genomes encoding GalNAc pathways. Candidate AgaR-binding motifs share a common sequence with consensus CTTTC that was found in multiple copies and arrangements in regulatory regions of aga genes. This study provides comprehensive insights into the common and distinctive features of the GalNAc/GalN catabolism and its regulation in diverse Proteobacteria.

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.

Comparative genomics of DtxR family regulons for metal homeostasis in Archaea

The DtxR family consists of metal-dependent transcription factors (DtxR-TFs) that regulate the expression of genes involved in metal homeostasis in the cell. The majority of characterized DtxR-TFs belong to Bacteria. In the current work, we applied a comparative genomics approach to predict DNA-binding sites and reconstruct regulons for DtxR-TFs in Archaea. As a result, we inferred 575 candidate binding sites for 139 DtxR-TFs in 77 genomes from 15 taxonomic orders. Novel DNA motifs of archaeal DtxR-TFs that have a common palindromic structure were classified into 10 distinct groups. By combining functional regulon reconstructions with phylogenetic analysis, we selected 28 DtxR-TF clades and assigned them metal specificities and regulator names. The reconstructed FetR (ferrous iron), MntR (manganese), and ZntR (zinc) regulons largely contain known or putative metal uptake transporters from the FeoAB, NRAMP, ZIP, and TroA families. A novel family of putative iron transporters (named Irt), including multiple FetR-regulated paralogs, was identified in iron-oxidizing Archaea from the Sulfolobales order. The reconstructed DtxR-TF regulons were reconciled with available transcriptomics data in Archaeoglobus, Halobacterium, and Thermococcus spp.

Bioactive glycans in a microbiome-directed food for children with malnutrition

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12-18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Genomic Reconstruction of Carbohydrate Utilization Capacities in Microbial-Mat Derived Consortia

Two nearly identical unicyanobacterial consortia (UCC) were previously isolated from benthic microbial mats that occur in a heliothermal saline lake in northern Washington State. Carbohydrates are a primary source of carbon and energy for most heterotrophic bacteria. Since CO₂ is the only carbon source provided, the cyanobacterium must provide a source of carbon to the heterotrophs. Available genomic sequences for all members of the UCC provide opportunity to investigate the metabolic routes of carbon transfer between autotroph and heterotrophs. Here, we applied a subsystem-based comparative genomics approach to reconstruct carbohydrate utilization pathways and identify glycohydrolytic enzymes, carbohydrate transporters and pathway-specific transcriptional regulators in 17 heterotrophic members of the UCC. The reconstructed metabolic pathways include 800 genes, near a one-fourth of which encode enzymes, transporters and regulators with newly assigned metabolic functions resulting in discovery of novel functional variants of carbohydrate utilization pathways. The in silico analysis revealed the utilization capabilities for 40 carbohydrates and their derivatives. Two Halomonas species demonstrated the largest number of sugar catabolic pathways. Trehalose, sucrose, maltose, glucose, and beta-glucosides are the most commonly utilized saccharides in this community. Reconstructed regulons for global regulators HexR and CceR include central carbohydrate metabolism genes in the members of Gammaproteobacteria and Alphaproteobacteria, respectively. Genomics analyses were supplemented by experimental characterization of metabolic phenotypes in four isolates derived from the consortia. Measurements of isolate growth on the defined medium supplied with individual carbohydrates confirmed most of the predicted catabolic phenotypes. Not all consortia members use carbohydrates and only a few use complex polysaccharides suggesting a hierarchical carbon flow from cyanobacteria to each heterotroph. In summary, the genomics-based identification of carbohydrate utilization capabilities provides a basis for future experimental studies of carbon flow in UCC.

Comparative genomics and evolution of transcriptional regulons in Proteobacteria

Comparative genomics approaches are broadly used for analysis of transcriptional regulation in bacterial genomes. In this work, we identified binding sites and reconstructed regulons for 33 orthologous groups of transcription factors (TFs) in 196 reference genomes from 21 taxonomic groups of Proteobacteria. Overall, we predict over 10 600 TF binding sites and identified more than 15 600 target genes for 1896 TFs constituting the studied orthologous groups of regulators. These include a set of orthologues for 21 metabolism-associated TFs from Escherichia coli and/or Shewanella that are conserved in five or more taxonomic groups and several additional TFs that represent non-orthologous substitutions of the metabolic regulators in some lineages of Proteobacteria. By comparing gene contents of the reconstructed regulons, we identified the core, taxonomy-specific and genome-specific TF regulon members and classified them by their metabolic functions. Detailed analysis of ArgR, TyrR, TrpR, HutC, HypR and other amino-acid-specific regulons demonstrated remarkable differences in regulatory strategies used by various lineages of Proteobacteria. The obtained genomic collection of in silico reconstructed TF regulons contains a large number of new regulatory interactions that await future experimental validation. The collection provides a framework for future evolutionary studies of transcriptional regulatory networks in Bacteria. It can be also used for functional annotation of putative metabolic transporters and enzymes that are abundant in the reconstructed regulons.

Comparative genomics of transcriptional regulation of methionine metabolism in Proteobacteria

Methionine metabolism and uptake genes in Proteobacteria are controlled by a variety of RNA and DNA regulatory systems. We have applied comparative genomics to reconstruct regulons for three known transcription factors, MetJ, MetR, and SahR, and three known riboswitch motifs, SAH, SAM-SAH, and SAM_alpha, in ∼ 200 genomes from 22 taxonomic groups of Proteobacteria. We also identified two novel regulons: a SahR-like transcription factor SamR controlling various methionine biosynthesis genes in the Xanthomonadales group, and a potential RNA regulatory element with terminator-antiterminator mechanism controlling the metX or metZ genes in beta-proteobacteria. For each analyzed regulator we identified the core, taxon-specific and genome-specific regulon members. By analyzing the distribution of these regulators in bacterial genomes and by comparing their regulon contents we elucidated possible evolutionary scenarios for the regulation of the methionine metabolism genes in Proteobacteria.

Global Expansion of Linezolid-Resistant Coagulase-Negative Staphylococci

Coagulase-negative staphylococci (CoNS) for a long time were considered avirulent constituents of the human and warm-blooded animal microbiota. However, at present, S. epidermidis, S. haemolyticus, and S. hominis are recognized as opportunistic pathogens. Although linezolid is not registered for the treatment of CoNS infections, it is widely used off-label, promoting emergence of resistance. Bioinformatic analysis based on maximum-likelihood phylogeny and Bayesian clustering of the CoNS genomes obtained in the current study and downloaded from public databases revealed the existence of international linezolid-resistant lineages, each of which probably had a common predecessor. Linezolid-resistant S. epidermidis sequence-type (ST) 2 from Russia, France, and Germany formed a compact group of closely related genomes with a median pairwise single nucleotide polymorphism (SNP) difference of fewer than 53 SNPs, and a common ancestor of this lineage appeared in 1998 (1986-2006) before introduction of linezolid in practice. Another compact group of linezolid-resistant S. epidermidis was represented by ST22 isolates from France and Russia with a median pairwise SNP difference of 40; a common ancestor of this lineage appeared in 2011 (2008-2013). Linezolid-resistant S. hominis ST2 from Russia, Germany, and Brazil also formed a group with a high-level genome identity with median 25.5 core-SNP differences; the appearance of the common progenitor dates to 2003 (1996-2012). Linezolid-resistant S. hominis isolates from Russia demonstrated associated resistance to teicoplanin. Analysis of a midpoint-rooted phylogenetic tree of the group confirmed the genetic proximity of Russian and German isolates; Brazilian isolates were phylogenetically distant. repUS5-like plasmids harboring cfr were detected in S. hominis and S. haemolyticus.