The rate at which genome sequencing data is accruing demands enhanced methods for functional annotation and metabolism discovery. Solute binding proteins (SBPs) facilitate the transport of the first reactant in a metabolic pathway, thereby constraining the regions of chemical space and the chemistries that must be considered for pathway reconstruction. We describe high-throughput protein production and differential scanning fluorimetry platforms, which enabled the screening of 158 SBPs against a 189 component library specifically tailored for this class of proteins. Like all screening efforts, this approach is limited by the practical constraints imposed by construction of the library, i.e., we can study only those metabolites that are known to exist and which can be made in sufficient quantities for experimentation. To move beyond these inherent limitations, we illustrate the promise of crystallographic- and mass spectrometric-based approaches for the unbiased use of entire metabolomes as screening libraries. Together, our approaches identified 40 new SBP ligands, generated experiment-based annotations for 2084 SBPs in 71 isofunctional clusters, and defined numerous metabolic pathways, including novel catabolic pathways for the utilization of ethanolamine as sole nitrogen source and the use of d-Ala-d-Ala as sole carbon source. These efforts begin to define an integrated strategy for realizing the full value of amassing genome sequence data.

A major challenge in genomics is the knowledge gap between sequence and its encoded function. Gain-of-function methods based on gene overexpression are attractive avenues for phenotype-based functional screens, but are not easily applied in high-throughput across many experimental conditions. Here, we present Dual Barcoded Shotgun Expression Library Sequencing (Dub-seq), a method that uses random DNA barcodes to greatly increase experimental throughput. As a demonstration of this approach, we construct a Dub-seq library with Escherichia coli genomic DNA, performed 155 genome-wide fitness assays in 52 experimental conditions, and identified overexpression phenotypes for 813 genes. We show that Dub-seq data is reproducible, accurately recapitulates known biology, and identifies hundreds of novel gain-of-function phenotypes for E. coli genes, a subset of which we verified with assays of individual strains. Dub-seq provides complementary information to loss-of-function approaches and will facilitate rapid and systematic functional characterization of microbial genomes.

Members of the genus Bifidobacterium are common inhabitants of the human gastrointestinal tract. Previously it was shown that arabino-oligosaccharides (AOS) might act as prebiotics and stimulate the bifidobacterial growth in the gut. However, despite the rapid accumulation of genomic data, the precise mechanisms by which these sugars are utilized and associated transcription control still remain unclear. In the current study, we used a comparative genomic approach to reconstruct arabinose and AOS utilization pathways in over 40 bacterial species belonging to the Bifidobacteriaceae family. The results indicate that the gene repertoire involved in the catabolism of these sugars is highly diverse, and even phylogenetically close species may differ in their utilization capabilities. Using bioinformatics analysis we identified potential DNA-binding motifs and reconstructed putative regulons for the arabinose and AOS utilization genes in the Bifidobacteriaceae genomes. Six LacI-family transcriptional factors (named AbfR, AauR, AauU1, AauU2, BauR1 and BauR2) and a TetR-family regulator (XsaR) presumably act as local repressors for AOS utilization genes encoding various α- or β-L-arabinofuranosidases and predicted AOS transporters. The ROK-family regulator AraU and the LacI-family regulator AraQ control adjacent operons encoding putative arabinose transporters and catabolic enzymes, respectively. However, the AraQ regulator is universally present in all Bifidobacterium species including those lacking the arabinose catabolic genes araBDA, suggesting its control of other genes. Comparative genomic analyses of prospective AraQ-binding sites allowed the reconstruction of AraQ regulons and a proposed binary repression/activation mechanism. The conserved core of reconstructed AraQ regulons in bifidobacteria includes araBDA, as well as genes from the central glycolytic and fermentation pathways (pyk, eno, gap, tkt, tal, galM, ldh). The current study expands the range of genes involved in bifidobacterial arabinose/AOS utilization and demonstrates considerable variations in associated metabolic pathways and regulons. Detailed comparative and phylogenetic analyses allowed us to hypothesize how the identified reconstructed regulons evolved in bifidobacteria. Our findings may help to improve carbohydrate catabolic phenotype prediction and metabolic modeling, while it may also facilitate rational development of novel prebiotics.

Bacteroides thetaiotaomicron, a predominant member of the human gut microbiota, is characterized by its ability to utilize a wide variety of polysaccharides using the extensive saccharolytic machinery that is controlled by an expanded repertoire of transcription factors (TFs). The availability of genomic sequences for multiple Bacteroides species opens an opportunity for their comparative analysis to enable characterization of their metabolic and regulatory networks.

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in bacteria is one of the critical tasks of modern genomics. The Shewanella genus is comprised of metabolically versatile gamma-proteobacteria, whose lifestyles and natural environments are substantially different from Escherichia coli and other model bacterial species. The comparative genomics approaches and computational identification of regulatory sites are useful for the in silico reconstruction of transcriptional regulatory networks in bacteria.

Genome scale annotation of regulatory interactions and reconstruction of regulatory networks are the crucial problems in bacterial genomics. The Lactobacillales order of bacteria collates various microorganisms having a large economic impact, including both human and animal pathogens and strains used in the food industry. Nonetheless, no systematic genome-wide analysis of transcriptional regulation has been previously made for this taxonomic group.

Bacterial response to nitric oxide (NO) is of major importance since NO is an obligatory intermediate of the nitrogen cycle. Transcriptional regulation of the dissimilatory nitric oxides metabolism in bacteria is diverse and involves FNR-like transcription factors HcpR, DNR, and NnrR; two-component systems NarXL and NarQP; NO-responsive activator NorR; and nitrite-sensitive repressor NsrR. Using comparative genomics approaches, we predict DNA-binding motifs for these transcriptional factors and describe corresponding regulons in available bacterial genomes. Within the FNR family of regulators, we observed a correlation of two specificity-determining amino acids and contacting bases in corresponding DNA recognition motif. Highly conserved regulon HcpR for the hybrid cluster protein and some other redox enzymes is present in diverse anaerobic bacteria, including Clostridia, Thermotogales, and delta-proteobacteria. NnrR and DNR control denitrification in alpha- and beta-proteobacteria, respectively. Sigma-54-dependent NorR regulon found in some gamma- and beta-proteobacteria contains various enzymes involved in the NO detoxification. Repressor NsrR, which was previously known to control only nitrite reductase operon in Nitrosomonas spp., appears to be the master regulator of the nitric oxides' metabolism, not only in most gamma- and beta-proteobacteria (including well-studied species such as Escherichia coli), but also in gram-positive Bacillus and Streptomyces species. Positional analysis and comparison of regulatory regions of NO detoxification genes allows us to propose the candidate NsrR-binding motif. The most conserved member of the predicted NsrR regulon is the NO-detoxifying flavohemoglobin Hmp. In enterobacteria, the regulon also includes two nitrite-responsive loci, nipAB (hcp-hcr) and nipC (dnrN), thus confirming the identity of the effector, i.e. nitrite. The proposed NsrR regulons in Neisseria and some other species are extended to include denitrification genes. As the result, we demonstrate considerable interconnection between various nitrogen-oxides-responsive regulatory systems for the denitrification and NO detoxification genes and evolutionary plasticity of this transcriptional network.

A comparative genomic approach was used to reconstruct transcriptional regulation of NAD biosynthesis in bacteria containing orthologs of Bacillus subtilis gene yrxA, a previously identified niacin-responsive repressor of NAD de novo synthesis. Members of YrxA family (re-named here NiaR) are broadly conserved in the Bacillus/Clostridium group and in the deeply branching Fusobacteria and Thermotogales lineages. We analyzed upstream regions of genes associated with NAD biosynthesis to identify candidate NiaR-binding DNA motifs and assess the NiaR regulon content in these species. Representatives of the two distinct types of candidate NiaR-binding sites, characteristic of the Firmicutes and Thermotogales, were verified by an electrophoretic mobility shift assay. In addition to transcriptional control of the nadABC genes, the NiaR regulon in some species extends to niacin salvage (the pncAB genes) and includes uncharacterized membrane proteins possibly involved in niacin transport. The involvement in niacin uptake proposed for one of these proteins (re-named NiaP), encoded by the B. subtilis gene yceI, was experimentally verified. In addition to bacteria, members of the NiaP family are conserved in multicellular eukaryotes, including human, pointing to possible NaiP involvement in niacin utilization in these organisms. Overall, the analysis of the NiaR and NrtR regulons (described in the accompanying paper) revealed mechanisms of transcriptional regulation of NAD metabolism in nearly a hundred diverse bacteria.

We used comparative genomics to investigate the distribution of conserved DNA-binding motifs in the regulatory regions of genes involved in iron and manganese homeostasis in alpha-proteobacteria. Combined with other computational approaches, this allowed us to reconstruct the metal regulatory network in more than three dozen species with available genome sequences. We identified several classes of cis-acting regulatory DNA motifs (Irr-boxes or ICEs, RirA-boxes, Iron-Rhodo-boxes, Fur-alpha-boxes, Mur-box or MRS, MntR-box, and IscR-boxes) in regulatory regions of various genes involved in iron and manganese uptake, Fe-S and heme biosynthesis, iron storage, and usage. Despite the different nature of the iron regulons in selected lineages of alpha-proteobacteria, the overall regulatory network is consistent with, and confirmed by, many experimental observations. This study expands the range of genes involved in iron homeostasis and demonstrates considerable interconnection between iron-responsive regulatory systems. The detailed comparative and phylogenetic analyses of the regulatory systems allowed us to propose a theory about the possible evolution of Fe and Mn regulons in alpha-proteobacteria. The main evolutionary event likely occurred in the common ancestor of the Rhizobiales and Rhodobacterales, where the Fur protein switched to regulating manganese transporters (and hence Fur had become Mur). In these lineages, the role of global iron homeostasis was taken by RirA and Irr, two transcriptional regulators that act by sensing the physiological consequence of the metal availability rather than its concentration per se, and thus provide for more flexible regulation.

The database of Alignable Tight Genomic Clusters (ATGCs) consists of closely related genomes of archaea and bacteria, and is a resource for research into prokaryotic microevolution. Construction of a data set with appropriate characteristics is a major hurdle for this type of studies. With the current rate of genome sequencing, it is difficult to follow the progress of the field and to determine which of the available genome sets meet the requirements of a given research project, in particular, with respect to the minimum and maximum levels of similarity between the included genomes. Additionally, extraction of specific content, such as genomic alignments or families of orthologs, from a selected set of genomes is a complicated and time-consuming process. The database addresses these problems by providing an intuitive and efficient web interface to browse precomputed ATGCs, select appropriate ones and access ATGC-derived data such as multiple alignments of orthologous proteins, matrices of pairwise intergenomic distances based on genome-wide analysis of synonymous and nonsynonymous substitution rates and others. The ATGC database will be regularly updated following new releases of the NCBI RefSeq. The database is hosted by the Genomics Division at Lawrence Berkeley National laboratory and is publicly available at http://atgc.lbl.gov.

The functions of most proteins are yet to be determined. The function of an enzyme is often defined by its interacting partners, including its substrate and product, and its role in larger metabolic networks. Here, we describe a computational method that predicts the functions of orphan enzymes by organizing them into a linear metabolic pathway. Given candidate enzyme and metabolite pathway members, this aim is achieved by finding those pathways that satisfy structural and network restraints implied by varied input information, including that from virtual screening, chemoinformatics, genomic context analysis, and ligand -binding experiments. We demonstrate this integrative pathway mapping method by predicting the L-gulonate catabolic pathway in Haemophilus influenzae Rd KW20. The prediction was subsequently validated experimentally by enzymology, crystallography, and metabolomics. Integrative pathway mapping by satisfaction of structural and network restraints is extensible to molecular networks in general and thus formally bridges the gap between structural biology and systems biology.

Although the impact of medicinal and culinary herbs on health and disease has been studied to varying extents, scarcely little is known about the impact of these herbs on gut microbiota and how such effects might contribute to their health benefits. We applied in vitro anaerobic cultivation of human fecal microbiota followed by 16S rRNA sequencing to study the modulatory effects of 4 culinary spices: Curcuma longa (turmeric), Zingiber officinale (ginger), Piper longum (pipli or long pepper), and Piper nigrum (black pepper). All herbs analyzed possessed substantial power to modulate fecal bacterial communities to include potential prebiotic and beneficial repressive effects. We additionally analyzed the sugar composition of each herb by mass spectrometry and conducted genome reconstruction of 11 relevant sugar utilization pathways, glycosyl hydrolase gene representation, and both butyrate and propionate biosynthesis potential to facilitate our ability to functionally interpret microbiota profiles. Results indicated that sugar composition is not predictive of the taxa responding to each herb; however, glycosyl hydrolase gene representation is strongly modulated by each herb, suggesting that polysaccharide substrates present in herbs provide selective potential on gut communities. Additionally, we conclude that catabolism of herbs by gut communities primarily involves sugar fermentation at the expense of amino acid metabolism. Among the herbs analyzed, only turmeric induced changes in community composition that are predicted to increase butyrate-producing taxa. Our data suggests that substrates present in culinary spices may drive beneficial alterations in gut communities thereby altering their collective metabolism to contribute to the salubrious effects on digestive efficiency and health. These results support the potential value of further investigations in human subjects to delineate whether the metabolism of these herbs contributes to documented and yet to be discovered health benefits.

Many organizations face challenges in managing and analyzing data, especially when relevant datasets arise from multiple sources and methods. Analyzing heterogeneous datasets and additional derived data requires rigorous tracking of their interrelationships and provenance. This task has long been a Grand Challenge of data science and has more recently been formalized in the FAIR principles: that all data objects be Findable, Accessible, Interoperable, and Reusable, both for machines and for people. Adherence to these principles is necessary for proper stewardship of information, for testing regulatory compliance, for measuring the efficiency of processes, and for facilitating reuse of data-analytical frameworks.

Understanding how members of the human gut microbiota prioritize nutrient resources is one component of a larger effort to decipher the mechanisms defining microbial community robustness and resiliency in health and disease. This knowledge is foundational for development of microbiota-directed therapeutics. To model how bacteria prioritize glycans in the gut, germfree mice were colonized with 13 human gut bacterial strains, including seven saccharolytic Bacteroidaceae species. Animals were fed a Western diet supplemented with pea fiber. After community assembly, an inducible CRISPR-based system was used to selectively and temporarily reduce the absolute abundance of Bacteroides thetaiotaomicron or B. cellulosilyticus by 10- to 60-fold. Each knockdown resulted in specific, reproducible increases in the abundances of other Bacteroidaceae and dynamic alterations in their expression of genes involved in glycan utilization. Emergence of these "alternate consumers" was associated with preservation of community saccharolytic activity. Using an inducible system for CRISPR base editing in vitro, we disrupted translation of transporters critical for utilizing dietary polysaccharides in Phocaeicola vulgatus, a B. cellulosilyticus knockdown-responsive taxon. In vitro and in vivo tests of the resulting P. vulgatus mutants allowed us to further characterize mechanisms associated with its increased fitness after knockdown. In principle, the approach described can be applied to study utilization of a range of nutrients and to preclinical efforts designed to develop therapeutic strategies for precision manipulation of microbial communities.

Parkinson's disease (PD) is a neurodegenerative disease characterized by the loss of dopaminergic neurons. Although the etiology of PD remains elusive, it has been hypothesized that initial dysregulation may occur in the gastrointestinal tract and may be accompanied by gut barrier defects. A strong clinical interest in developing therapeutics exists, including for the treatment of gut microbiota and physiology. We previously reported the impact of healthy fecal microbiota anaerobic cultures supplemented with nootropic herbs. Here, we evaluated the effect of nootropic Ayurvedic herbs on fecal microbiota derived from subjects with PD in vitro using 16S rRNA sequencing. The microbiota underwent substantial change in response to each treatment, comparable in magnitude to that observed from healthy subjects. However, the fecal samples derived from each participant displayed unique changes, consistent with a personalized response. We used genome-wide metabolic reconstruction to predict the community's metabolic potential to produce products relevant to PD pathology, including SCFAs, vitamins and amino acid degradation products. These results suggest the potential value of conducting in vitro cultivation and analyses of PD stool samples as a means of prescreening patients to select the medicinal herbs for which that individual is most likely to respond and derive benefit.

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.

The prebiotic potential of nervine herbal medicines has been scarcely studied. We therefore used anaerobic human fecal cultivation to investigate whether medicinal herbs commonly used as treatment in neurological health and disease in Ayurveda and other traditional systems of medicine modulate gut microbiota. Profiling of fecal cultures supplemented with either Kapikacchu, Gotu Kola, Bacopa/Brahmi, Shankhapushpi, Boswellia/Frankincense, Jatamansi, Bhringaraj, Guduchi, Ashwagandha or Shatavari by 16S rRNA sequencing revealed profound changes in diverse taxa. Principal coordinate analysis highlights that each herb drives the formation of unique microbial communities predicted to display unique metabolic potential. The relative abundance of approximately one-third of the 243 enumerated species was altered by all herbs. Additional species were impacted in an herb-specific manner. In this study, we combine genome reconstruction of sugar utilization and short chain fatty acid (SCFA) pathways encoded in the genomes of 216 profiled taxa with monosaccharide composition analysis of each medicinal herb by quantitative mass spectrometry to enhance the interpretation of resulting microbial communities and discern potential drivers of microbiota restructuring. Collectively, our results indicate that gut microbiota engage in both protein and glycan catabolism, providing amino acid and sugar substrates that are consumed by fermentative species. We identified taxa that are efficient amino acid fermenters and those capable of both amino acid and sugar fermentation. Herb-induced microbial communities are predicted to alter the relative abundance of taxa encoding SCFA (butyrate and propionate) pathways. Co-occurrence network analyses identified a large number of taxa pairs in medicinal herb cultures. Some of these pairs displayed related culture growth relationships in replicate cultures highlighting potential functional interactions among medicinal herb-induced taxa.

The central carbon/lactate utilization pathway in the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough, is encoded by the highly conserved operon DVU3025-3033. Our earlier in vitro genome-wide study had suggested a network of four two-component system regulators that target this large operon; however, how these four regulators control this operon was not known. Here, we probe the regulation of the lactate utilization operon with mutant strains and DNA-protein binding assays. We show that the LurR response regulator is required for optimal growth and complete lactate utilization, and that it activates the DVU3025-3033 lactate oxidation operon as well as DVU2451, a lactate permease gene, in the presence of lactate. We show by electrophoretic mobility shift assays that LurR binds to three sites in the upstream region of DVU3025, the first gene of the operon. NrfR, a response regulator that is activated under nitrite stress, and LurR share similar binding site motifs and bind the same sites upstream of DVU3025. The DVU3025 promoter also has a binding site motif (Pho box) that is bound by PhoB, a two-component response regulator activated under phosphate limitation. The lactate utilization operon, the regulator LurR, and LurR binding sites are conserved across the order Desulfovibrionales whereas possible modulation of the lactate utilization genes by additional regulators such as NrfR and PhoB appears to be limited to D. vulgaris.

L-rhamnose (L-Rha) is a deoxy-hexose sugar commonly found in nature. L-Rha catabolic pathways were previously characterized in various bacteria including Escherichia coli. Nevertheless, homology searches failed to recognize all the genes for the complete L-Rha utilization pathways in diverse microbial species involved in biomass decomposition. Moreover, the regulatory mechanisms of L-Rha catabolism have remained unclear in most species. A comparative genomics approach was used to reconstruct the L-Rha catabolic pathways and transcriptional regulons in the phyla Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria, and Thermotogae. The reconstructed pathways include multiple novel enzymes and transporters involved in the utilization of L-Rha and L-Rha-containing polymers. Large-scale regulon inference using bioinformatics revealed remarkable variations in transcriptional regulators for L-Rha utilization genes among bacteria. A novel bifunctional enzyme, L-rhamnulose-phosphate aldolase (RhaE) fused to L-lactaldehyde dehydrogenase (RhaW), which is not homologous to previously characterized L-Rha catabolic enzymes, was identified in diverse bacteria including Chloroflexi, Bacilli, and Alphaproteobacteria. By using in vitro biochemical assays we validated both enzymatic activities of the purified recombinant RhaEW proteins from Chloroflexus aurantiacus and Bacillus subtilis. Another novel enzyme of the L-Rha catabolism, L-lactaldehyde reductase (RhaZ), was identified in Gammaproteobacteria and experimentally validated by in vitro enzymatic assays using the recombinant protein from Salmonella typhimurium. C. aurantiacus induced transcription of the predicted L-Rha utilization genes when L-Rha was present in the growth medium and consumed L-Rha from the medium. This study provided comprehensive insights to L-Rha catabolism and its regulation in diverse Bacteria.

The MocR-subfamily transcription factors (MocR-TFs) characterized by the GntR-family DNA-binding domain and aminotransferase-like sensory domain are broadly distributed among certain lineages of Bacteria. Characterized MocR-TFs bind pyridoxal 5'-phosphate (PLP) and control transcription of genes involved in PLP, gamma aminobutyric acid (GABA) and taurine metabolism via binding specific DNA operator sites. To identify putative target genes and DNA binding motifs of MocR-TFs, we performed comparative genomics analysis of over 250 bacterial genomes. The reconstructed regulons for 825 MocR-TFs comprise structural genes from over 200 protein families involved in diverse biological processes. Using the genome context and metabolic subsystem analysis we tentatively assigned functional roles for 38 out of 86 orthologous groups of studied regulators. Most of these MocR-TF regulons are involved in PLP metabolism, as well as utilization of GABA, taurine and ectoine. The remaining studied MocR-TF regulators presumably control genes encoding enzymes involved in reduction/oxidation processes, various transporters and PLP-dependent enzymes, for example aminotransferases. Predicted DNA binding motifs of MocR-TFs are generally similar in each orthologous group and are characterized by two to four repeated sequences. Identified motifs were classified according to their structures. Motifs with direct and/or inverted repeat symmetry constitute the majority of inferred DNA motifs, suggesting preferable TF dimerization in head-to-tail or head-to-head configuration. The obtained genomic collection of in silico reconstructed MocR-TF motifs and regulons in Bacteria provides a basis for future experimental characterization of molecular mechanisms for various regulators in this family.