The PATRIC Bioinformatics Resource Center: expanding data and analysis capabilities.

The PathoSystems Resource Integration Center (PATRIC) is the bacterial Bioinformatics Resource Center funded by the National Institute of Allergy and Infectious Diseases (https://www.patricbrc.org). PATRIC supports bioinformatic analyses of all bacteria with a special emphasis on pathogens, offering a rich comparative analysis environment that provides users with access to over 250 000 uniformly annotated and publicly available genomes with curated metadata. PATRIC offers web-based visualization and comparative analysis tools, a private workspace in which users can analyze their own data in the context of the public collections, services that streamline complex bioinformatic workflows and command-line tools for bulk data analysis. Over the past several years, as genomic and other omics-related experiments have become more cost-effective and widespread, we have observed considerable growth in the usage of and demand for easy-to-use, publicly available bioinformatic tools and services. Here we report the recent updates to the PATRIC resource, including new web-based comparative analysis tools, eight new services and the release of a command-line interface to access, query and analyze data.

PATRIC as a unique resource for studying antimicrobial resistance.

The Pathosystems Resource Integration Center (PATRIC, www.patricbrc.org) is designed to provide researchers with the tools and services that they need to perform genomic and other 'omic' data analyses. In response to mounting concern over antimicrobial resistance (AMR), the PATRIC team has been developing new tools that help researchers understand AMR and its genetic determinants. To support comparative analyses, we have added AMR phenotype data to over 15 000 genomes in the PATRIC database, often assembling genomes from reads in public archives and collecting their associated AMR panel data from the literature to augment the collection. We have also been using this collection of AMR metadata to build machine learning-based classifiers that can predict the AMR phenotypes and the genomic regions associated with resistance for genomes being submitted to the annotation service. Likewise, we have undertaken a large AMR protein annotation effort by manually curating data from the literature and public repositories. This collection of 7370 AMR reference proteins, which contains many protein annotations (functional roles) that are unique to PATRIC and RAST, has been manually curated so that it projects stably across genomes. The collection currently projects to 1 610 744 proteins in the PATRIC database. Finally, the PATRIC Web site has been expanded to enable AMR-based custom page views so that researchers can easily explore AMR data and design experiments based on whole genomes or individual genes.

Transcriptomics reveal different metabolic strategies for acid resistance and gamma-aminobutyric acid (GABA) production in select Levilactobacillus brevis strains.

BACKGROUND: Of the many neurotransmitters in humans, gamma-aminobutyric acid (GABA) shows potential for improving several mental health indications such as stress and anxiety. The microbiota-gut-brain axis is an important pathway for GABAergic effects, as microbially-secreted GABA within the gut can affect host mental health outcomes. Understanding the molecular characteristics of GABA production by microbes within the gut can offer insight to novel therapies for mental health. RESULTS: Three strains of Levilactobacillus brevis with syntenous glutamate decarboxylase (GAD) operons were evaluated for overall growth, glutamate utilization, and GABA production in typical synthetic growth media supplemented with monosodium glutamate (MSG). Levilactobacillus brevis Lbr-6108™ (Lbr-6108), formerly known as L. brevis DPC 6108, and Levilactobacillus brevis Lbr-35 ™ (Lbr-35) had similar growth profiles but differed significantly in GABA secretion and acid resistance. Lbr-6108 produced GABA early within the growth phase and produced significantly more GABA than Lbr-35 and the type strain Levilactobacillus brevis ATCC 14869 after the stationary phase. The global gene expression during GABA production at several timepoints was determined by RNA sequencing. The GAD operon, responsible for GABA production and secretion, activated in Lbr-6108 after only 6 h of fermentation and continued throughout the stationary phase. Furthermore, Lbr-6108 activated many different acid resistance mechanisms concurrently, which contribute to acid tolerance and energy production. In contrast, Lbr-35, which has a genetically similar GAD operon, including two copies of the GAD gene, showed no upregulation of the GAD operon, even when cultured with MSG. CONCLUSIONS: This study is the first to evaluate whole transcriptome changes in Levilactobacillus brevis during GABA production in different growth phases. The concurrent expression of multiple acid-resistance mechanisms reveals niche-specific metabolic functionality between common human commensals and highlights the complex regulation of GABA metabolism in this important microbial species. Furthermore, the increased and rapid GABA production of Lbr-6108 highlights the strain's potential as a therapeutic and the overall value of screening microbes for effector molecule output.

Improvements to PATRIC, the all-bacterial Bioinformatics Database and Analysis Resource Center.

The Pathosystems Resource Integration Center (PATRIC) is the bacterial Bioinformatics Resource Center (https://www.patricbrc.org). Recent changes to PATRIC include a redesign of the web interface and some new services that provide users with a platform that takes them from raw reads to an integrated analysis experience. The redesigned interface allows researchers direct access to tools and data, and the emphasis has changed to user-created genome-groups, with detailed summaries and views of the data that researchers have selected. Perhaps the biggest change has been the enhanced capability for researchers to analyze their private data and compare it to the available public data. Researchers can assemble their raw sequence reads and annotate the contigs using RASTtk. PATRIC also provides services for RNA-Seq, variation, model reconstruction and differential expression analysis, all delivered through an updated private workspace. Private data can be compared by 'virtual integration' to any of PATRIC's public data. The number of genomes available for comparison in PATRIC has expanded to over 80 000, with a special emphasis on genomes with antimicrobial resistance data. PATRIC uses this data to improve both subsystem annotation and k-mer classification, and tags new genomes as having signatures that indicate susceptibility or resistance to specific antibiotics.

A strictly monofunctional bacterial hydroxymethylpyrimidine phosphate kinase precludes damaging errors in thiamin biosynthesis.

The canonical kinase (ThiD) that converts the thiamin biosynthesis intermediate hydroxymethylpyrimidine (HMP) monophosphate to the diphosphate can also very efficiently convert free HMP to the monophosphate in prokaryotes, plants, and fungi. This HMP kinase activity enables salvage of HMP, but it is not substrate-specific and so allows toxic HMP analogs and damage products to infiltrate the thiamin biosynthesis pathway. Comparative analysis of bacterial genomes uncovered a gene, thiD2 , that is often fused to the thiamin synthesis gene thiE and could potentially encode a replacement for ThiD. Standalone ThiD2 proteins and ThiD2 fusion domains are small (~130-residues) and do not belong to any previously known protein family. Genetic and biochemical analyses showed that representative standalone and fused ThiD2 proteins catalyze phosphorylation of HMP monophosphate, but not of HMP or its toxic analogs and damage products such as bacimethrin and 5-(hydroxymethyl)-2-methylpyrimidin-4-ol. As strictly monofunctional HMP monophosphate kinases, ThiD2 proteins eliminate a potentially fatal vulnerability of canonical ThiD, at the cost of the ability to reclaim HMP formed by thiamin turnover.

Noncanonical coproporphyrin-dependent bacterial heme biosynthesis pathway that does not use protoporphyrin.

It has been generally accepted that biosynthesis of protoheme (heme) uses a common set of core metabolic intermediates that includes protoporphyrin. Herein, we show that the Actinobacteria and Firmicutes (high-GC and low-GC Gram-positive bacteria) are unable to synthesize protoporphyrin. Instead, they oxidize coproporphyrinogen to coproporphyrin, insert ferrous iron to make Fe-coproporphyrin (coproheme), and then decarboxylate coproheme to generate protoheme. This pathway is specified by three genes named hemY, hemH, and hemQ. The analysis of 982 representative prokaryotic genomes is consistent with this pathway being the most ancient heme synthesis pathway in the Eubacteria. Our results identifying a previously unknown branch of tetrapyrrole synthesis support a significant shift from current models for the evolution of bacterial heme and chlorophyll synthesis. Because some organisms that possess this coproporphyrin-dependent branch are major causes of human disease, HemQ is a novel pharmacological target of significant therapeutic relevance, particularly given high rates of antimicrobial resistance among these pathogens.

High-throughput comparison, functional annotation, and metabolic modeling of plant genomes using the PlantSEED resource.

The increasing number of sequenced plant genomes is placing new demands on the methods applied to analyze, annotate, and model these genomes. Today's annotation pipelines result in inconsistent gene assignments that complicate comparative analyses and prevent efficient construction of metabolic models. To overcome these problems, we have developed the PlantSEED, an integrated, metabolism-centric database to support subsystems-based annotation and metabolic model reconstruction for plant genomes. PlantSEED combines SEED subsystems technology, first developed for microbial genomes, with refined protein families and biochemical data to assign fully consistent functional annotations to orthologous genes, particularly those encoding primary metabolic pathways. Seamless integration with its parent, the prokaryotic SEED database, makes PlantSEED a unique environment for cross-kingdom comparative analysis of plant and bacterial genomes. The consistent annotations imposed by PlantSEED permit rapid reconstruction and modeling of primary metabolism for all plant genomes in the database. This feature opens the unique possibility of model-based assessment of the completeness and accuracy of gene annotation and thus allows computational identification of genes and pathways that are restricted to certain genomes or need better curation. We demonstrate the PlantSEED system by producing consistent annotations for 10 reference genomes. We also produce a functioning metabolic model for each genome, gapfilling to identify missing annotations and proposing gene candidates for missing annotations. Models are built around an extended biomass composition representing the most comprehensive published to date. To our knowledge, our models are the first to be published for seven of the genomes analyzed.

Systematic identification and analysis of frequent gene fusion events in metabolic pathways.

BACKGROUND: Gene fusions are the most powerful type of in silico-derived functional associations. However, many fusion compilations were made when <100 genomes were available, and algorithms for identifying fusions need updating to handle the current avalanche of sequenced genomes. The availability of a large fusion dataset would help probe functional associations and enable systematic analysis of where and why fusion events occur. RESULTS: Here we present a systematic analysis of fusions in prokaryotes. We manually generated two training sets: (i) 121 fusions in the model organism Escherichia coli; (ii) 131 fusions found in B vitamin metabolism. These sets were used to develop a fusion prediction algorithm that captured the training set fusions with only 7 % false negatives and 50 % false positives, a substantial improvement over existing approaches. This algorithm was then applied to identify 3.8 million potential fusions across 11,473 genomes. The results of the analysis are available in a searchable database at http://modelseed.org/projects/fusions/ . A functional analysis identified 3,000 reactions associated with frequent fusion events and revealed areas of metabolism where fusions are particularly prevalent. CONCLUSIONS: Customary definitions of fusions were shown to be ambiguous, and a stricter one was proposed. Exploring the genes participating in fusion events showed that they most commonly encode transporters, regulators, and metabolic enzymes. The major rationales for fusions between metabolic genes appear to be overcoming pathway bottlenecks, avoiding toxicity, controlling competing pathways, and facilitating expression and assembly of protein complexes. Finally, our fusion dataset provides powerful clues to decipher the biological activities of domains of unknown function.

The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST).

In 2004, the SEED (http://pubseed.theseed.org/) was created to provide consistent and accurate genome annotations across thousands of genomes and as a platform for discovering and developing de novo annotations. The SEED is a constantly updated integration of genomic data with a genome database, web front end, API and server scripts. It is used by many scientists for predicting gene functions and discovering new pathways. In addition to being a powerful database for bioinformatics research, the SEED also houses subsystems (collections of functionally related protein families) and their derived FIGfams (protein families), which represent the core of the RAST annotation engine (http://rast.nmpdr.org/). When a new genome is submitted to RAST, genes are called and their annotations are made by comparison to the FIGfam collection. If the genome is made public, it is then housed within the SEED and its proteins populate the FIGfam collection. This annotation cycle has proven to be a robust and scalable solution to the problem of annotating the exponentially increasing number of genomes. To date, >12 000 users worldwide have annotated >60 000 distinct genomes using RAST. Here we describe the interconnectedness of the SEED database and RAST, the RAST annotation pipeline and updates to both resources.

Genomic and experimental evidence for multiple metabolic functions in the RidA/YjgF/YER057c/UK114 (Rid) protein family.

BACKGROUND: It is now recognized that enzymatic or chemical side-reactions can convert normal metabolites to useless or toxic ones and that a suite of enzymes exists to mitigate such metabolite damage. Examples are the reactive imine/enamine intermediates produced by threonine dehydratase, which damage the pyridoxal 5'-phosphate cofactor of various enzymes causing inactivation. This damage is pre-empted by RidA proteins, which hydrolyze the imines before they do harm. RidA proteins belong to the YjgF/YER057c/UK114 family (here renamed the Rid family). Most other members of this diverse and ubiquitous family lack defined functions. RESULTS: Phylogenetic analysis divided the Rid family into a widely distributed, apparently archetypal RidA subfamily and seven other subfamilies (Rid1 to Rid7) that are largely confined to bacteria and often co-occur in the same organism with RidA and each other. The Rid1 to Rid3 subfamilies, but not the Rid4 to Rid7 subfamilies, have a conserved arginine residue that, in RidA proteins, is essential for imine-hydrolyzing activity. Analysis of the chromosomal context of bacterial RidA genes revealed clustering with genes for threonine dehydratase and other pyridoxal 5'-phosphate-dependent enzymes, which fits with the known RidA imine hydrolase activity. Clustering was also evident between Rid family genes and genes specifying FAD-dependent amine oxidases or enzymes of carbamoyl phosphate metabolism. Biochemical assays showed that Salmonella enterica RidA and Rid2, but not Rid7, can hydrolyze imines generated by amino acid oxidase. Genetic tests indicated that carbamoyl phosphate overproduction is toxic to S. enterica cells lacking RidA, and metabolomic profiling of Rid knockout strains showed ten-fold accumulation of the carbamoyl phosphate-related metabolite dihydroorotate. CONCLUSIONS: Like the archetypal RidA subfamily, the Rid2, and probably the Rid1 and Rid3 subfamilies, have imine-hydrolyzing activity and can pre-empt damage from imines formed by amine oxidases as well as by pyridoxal 5'-phosphate enzymes. The RidA subfamily has an additional damage pre-emption role in carbamoyl phosphate metabolism that has yet to be biochemically defined. Finally, the Rid4 to Rid7 subfamilies appear not to hydrolyze imines and thus remain mysterious.

HemQ: An iron-coproporphyrin oxidative decarboxylase for protoheme synthesis in Firmicutes and Actinobacteria.

Genes for chlorite dismutase-like proteins are found widely among heme-synthesizing bacteria and some Archaea. It is now known that among the Firmicutes and Actinobacteria these proteins do not possess chlorite dismutase activity but instead are essential for heme synthesis. These proteins, named HemQ, are iron-coproporphyrin (coproheme) decarboxylases that catalyze the oxidative decarboxylation of coproheme III into protoheme IX. As purified, HemQs do not contain bound heme, but readily bind exogeneously supplied heme with low micromolar affinity. The heme-bound form of HemQ has low peroxidase activity and in the presence of peroxide the bound heme may be destroyed. Thus, it is possible that HemQ may serve a dual role as a decarboxylase in heme biosynthesis and a regulatory protein in heme homeostasis.

Salvage of the thiamin pyrimidine moiety by plant TenA proteins lacking an active-site cysteine.

The TenA protein family occurs in prokaryotes, plants and fungi; it has two subfamilies, one (TenA_C) having an active-site cysteine, the other (TenA_E) not. TenA_C proteins participate in thiamin salvage by hydrolysing the thiamin breakdown product amino-HMP (4-amino-5-aminomethyl-2-methylpyrimidine) to HMP (4-amino-5-hydroxymethyl-2-methylpyrimidine); the function of TenA_E proteins is unknown. Comparative analysis of prokaryote and plant genomes predicted that (i) TenA_E has a salvage role similar to, but not identical with, that of TenA_C and (ii) that TenA_E and TenA_C also have non-salvage roles since they occur in organisms that cannot make thiamin. Recombinant Arabidopsis and maize TenA_E proteins (At3g16990, GRMZM2G080501) hydrolysed amino-HMP to HMP and, far more actively, hydrolysed the N-formyl derivative of amino-HMP to amino-HMP. Ablating the At3g16990 gene in a line with a null mutation in the HMP biosynthesis gene ThiC prevented its rescue by amino-HMP. Ablating At3g16990 in the wild-type increased sensitivity to paraquat-induced oxidative stress; HMP overcame this increased sensitivity. Furthermore, the expression of TenA_E and ThiC genes in Arabidopsis and maize was inversely correlated. These results indicate that TenA_E proteins mediate amidohydrolase and aminohydrolase steps in the salvage of thiamin breakdown products. As such products can be toxic, TenA_E proteins may also pre-empt toxicity.

The alternative coproporphyrinogen III oxidase (CgoN) catalyzes the oxygen-independent conversion of coproporphyrinogen III into coproporphyrin III.

Nature utilizes three distinct pathways to synthesize the essential enzyme cofactor heme. The coproporphyrin III-dependent pathway, predominantly present in Bacillaceae, employs an oxygen-dependent coproporphyrinogen III oxidase (CgoX) that converts coproporphyrinogen III into coproporphyrin III. In this study, we report the bioinformatic-based identification of a gene called ytpQ, encoding a putative oxygen-independent counterpart, which we propose to term CgoN, from Priestia (Bacillus) megaterium. The recombinantly produced, purified, and monomeric YtpQ (CgoN) protein is shown to catalyze the oxygen-independent conversion of coproporphyrinogen III into coproporphyrin III. Minimal non-enzymatic conversion of coproporphyrinogen III was observed under the anaerobic test conditions employed in this study. FAD was identified as a cofactor, and menadione served as an artificial acceptor for the six abstracted electrons, with a KM value of 3.95 µmol/L and a kcat of 0.63 per min for the substrate. The resulting coproporphyrin III, in turn, acts as an effective substrate for the subsequent enzyme of the pathway, the coproporphyrin III ferrochelatase (CpfC). Under aerobic conditions, oxygen directly serves as an electron acceptor, but is replaced by the more efficient action of menadione. An AlphaFold2 model of the enzyme suggests that YtpQ adopts a compact triangular shape consisting of three domains. The N-terminal domain appears to be flexible with respect to the rest of the structure, potentially creating a ligand binding site that opens and closes during the catalytic cycle. A catalytic mechanism similar to the oxygen-independent protoporphyrinogen IX oxidase PgoH1 (HemG), based on the flavin-dependent abstraction of six electrons from coproporphyrinogen III and their potential quinone-dependent transfer to a membrane-localized electron transport chain, is proposed.

Plant B vitamin pathways and their compartmentation: a guide for the perplexed.

The B vitamins and the cofactors derived from them are essential for life. B vitamin synthesis in plants is consequently as crucial to plants themselves as it is to humans and animals, whose B vitamin nutrition depends largely on plants. The synthesis and salvage pathways for the seven plant B vitamins are now broadly known, but certain enzymes and many transporters have yet to be identified, and the subcellular locations of various reactions are unclear. Although very substantial, what is not known about plant B vitamin pathways is regrettably difficult to discern from the literature or from biochemical pathway databases. Nor do databases accurately represent all that is known about B vitamin pathways-above all their compartmentation-because the facts are scattered throughout the literature, and thus hard to piece together. These problems (i) deter discoveries because newcomers to B vitamins cannot see which mysteries still need solving; and (ii) impede metabolic reconstruction and modelling of B vitamin pathways because genes for reactions or transport steps are missing. This review therefore takes a fresh approach to capture current knowledge of B vitamin pathways in plants. The synthesis pathways, key salvage routes, and their subcellular compartmentation are surveyed in depth, and encoded in the SEED database (http://pubseed.theseed.org/seedviewer.cgi?page=PlantGateway) for Arabidopsis and maize. The review itself and the encoded pathways specifically identify enigmatic or missing reactions, enzymes, and transporters. The SEED-encoded B vitamin pathway collection is a publicly available, expertly curated, one-stop resource for metabolic reconstruction and modeling.

Red death in Caenorhabditis elegans caused by Pseudomonas aeruginosa PAO1.

During host injury, Pseudomonas aeruginosa can be cued to express a lethal phenotype within the intestinal tract reservoir-a hostile, nutrient scarce environment depleted of inorganic phosphate. Here we determined if phosphate depletion activates a lethal phenotype in P. aeruginosa during intestinal colonization. To test this, we allowed Caenorhabditis elegans to feed on lawns of P. aeruginosa PAO1 grown on high and low phosphate media. Phosphate depletion caused PAO1 to kill 60% of nematodes whereas no worms died on high phosphate media. Unexpectedly, intense redness was observed in digestive tubes of worms before death. Using a combination of transcriptome analyses, mutants, and reporter constructs, we identified 3 global virulence systems that were involved in the "red death" response of P. aeruginosa during phosphate depletion; they included phosphate signaling (PhoB), the MvfR-PQS pathway of quorum sensing, and the pyoverdin iron acquisition system. Activation of all 3 systems was required to form a red colored PQS+Fe(3+) complex which conferred a lethal phenotype in this model. When pyoverdin production was inhibited in P. aeruginosa by providing excess iron, red death was attenuated in C. elegans and mortality was decreased in mice intestinally inoculated with P. aeruginosa. Introduction of the red colored PQS+Fe(3+) complex into the digestive tube of C. elegans or mouse intestine caused mortality associated with epithelial disruption and apoptosis. In summary, red death in C. elegans reveals a triangulated response between PhoB, MvfR-PQS, and pyoverdin in response to phosphate depletion that activates a lethal phenotype in P. aeruginosa.

Discovery and Characterization of HemQ: an essential heme biosynthetic pathway component.

Here we identify a previously undescribed protein, HemQ, that is required for heme synthesis in Gram-positive bacteria. We have characterized HemQ from Bacillus subtilis and a number of Actinobacteria. HemQ is a multimeric heme-binding protein. Spectroscopic studies indicate that this heme is high spin ferric iron and is ligated by a conserved histidine with the sixth coordination site available for binding a small molecule. The presence of HemQ along with the terminal two pathway enzymes, protoporphyrinogen oxidase (HemY) and ferrochelatase, is required to synthesize heme in vivo and in vitro. Although the exact role played by HemQ remains to be characterized, to be fully functional in vitro it requires the presence of a bound heme. HemQ possesses minimal peroxidase activity, but as a catalase it has a turnover of over 10(4) min(-1). We propose that this activity may be required to eliminate hydrogen peroxide that is generated by each turnover of HemY. Given the essential nature of heme synthesis and the restricted distribution of HemQ, this protein is a potential antimicrobial target for pathogens such as Mycobacterium tuberculosis.

Genomics-driven reconstruction of acinetobacter NAD metabolism: insights for antibacterial target selection.

Enzymes involved in the last steps of NAD biogenesis, nicotinate mononucleotide adenylyltransferase (NadD) and NAD synthetase (NadE), are conserved and essential in most bacterial species and are established targets for antibacterial drug development. Our genomics-based reconstruction of NAD metabolism in the emerging pathogen Acinetobacter baumannii revealed unique features suggesting an alternative targeting strategy. Indeed, genomes of all analyzed Acinetobacter species do not encode NadD, which is functionally replaced by its distant homolog NadM. We combined bioinformatics with genetic and biochemical techniques to elucidate this and other important features of Acinetobacter NAD metabolism using a model (nonpathogenic) strain Acinetobacter baylyi sp. ADP1. Thus, a comparative kinetic characterization of PncA, PncB, and NadV enzymes allowed us to suggest distinct physiological roles for the two alternative, deamidating and nondeamidating, routes of nicotinamide salvage/recycling. The role of the NiaP transporter in both nicotinate and nicotinamide salvage was confirmed. The nondeamidating route was shown to be transcriptionally regulated by an ADP-ribose-responsive repressor NrtR. The NadM enzyme was shown to possess dual substrate specificity toward both nicotinate and nicotinamide mononucleotide substrates, which is consistent with its essential role in all three routes of NAD biogenesis, de novo synthesis as well as the two salvage pathways. The experimentally confirmed unconditional essentiality of nadM provided support for the choice of the respective enzyme as a drug target. In contrast, nadE, encoding a glutamine-dependent NAD synthetase, proved to be dispensable when the nondeamidating salvage pathway functioned as the only route of NAD biogenesis.

Synergistic use of plant-prokaryote comparative genomics for functional annotations.

BACKGROUND: Identifying functions for all gene products in all sequenced organisms is a central challenge of the post-genomic era. However, at least 30-50% of the proteins encoded by any given genome are of unknown or vaguely known function, and a large number are wrongly annotated. Many of these 'unknown' proteins are common to prokaryotes and plants. We set out to predict and experimentally test the functions of such proteins. Our approach to functional prediction integrates comparative genomics based mainly on microbial genomes with functional genomic data from model microorganisms and post-genomic data from plants. This approach bridges the gap between automated homology-based annotations and the classical gene discovery efforts of experimentalists, and is more powerful than purely computational approaches to identifying gene-function associations. RESULTS: Among Arabidopsis genes, we focused on those (2,325 in total) that (i) are unique or belong to families with no more than three members, (ii) occur in prokaryotes, and (iii) have unknown or poorly known functions. Computer-assisted selection of promising targets for deeper analysis was based on homology-independent characteristics associated in the SEED database with the prokaryotic members of each family. In-depth comparative genomic analysis was performed for 360 top candidate families. From this pool, 78 families were connected to general areas of metabolism and, of these families, specific functional predictions were made for 41. Twenty-one predicted functions have been experimentally tested or are currently under investigation by our group in at least one prokaryotic organism (nine of them have been validated, four invalidated, and eight are in progress). Ten additional predictions have been independently validated by other groups. Discovering the function of very widespread but hitherto enigmatic proteins such as the YrdC or YgfZ families illustrates the power of our approach. CONCLUSIONS: Our approach correctly predicted functions for 19 uncharacterized protein families from plants and prokaryotes; none of these functions had previously been correctly predicted by computational methods. The resulting annotations could be propagated with confidence to over six thousand homologous proteins encoded in over 900 bacterial, archaeal, and eukaryotic genomes currently available in public databases.

Discovery of a gene involved in a third bacterial protoporphyrinogen oxidase activity through comparative genomic analysis and functional complementation.

Tetrapyrroles are ubiquitous molecules in nearly all living organisms. Heme, an iron-containing tetrapyrrole, is widely distributed in nature, including most characterized aerobic and facultative bacteria. A large majority of bacteria that contain heme possess the ability to synthesize it. Despite this capability and the fact that the biosynthetic pathway has been well studied, enzymes catalyzing at least three steps have remained "missing" in many bacteria. In the current work, we have employed comparative genomics via the SEED genomic platform, coupled with experimental verification utilizing Acinetobacter baylyi ADP1, to identify one of the missing enzymes, a new protoporphyrinogen oxidase, the penultimate enzyme in heme biosynthesis. COG1981 was identified by genomic analysis as a candidate protein family for the missing enzyme in bacteria that lacked HemG or HemY, two known protoporphyrinogen oxidases. The predicted amino acid sequence of COG1981 is unlike those of the known enzymes HemG and HemY, but in some genomes, the gene encoding it is found neighboring other heme biosynthetic genes. When the COG1981 gene was deleted from the genome of A. baylyi, a bacterium that lacks both hemG and hemY, the organism became auxotrophic for heme. Cultures accumulated porphyrin intermediates, and crude cell extracts lacked protoporphyrinogen oxidase activity. The heme auxotrophy was rescued by the presence of a plasmid-borne protoporphyrinogen oxidase gene from a number of different organisms, such as hemG from Escherichia coli, hemY from Myxococcus xanthus, or the human gene for protoporphyrinogen oxidase.

Prevention of siderophore- mediated gut-derived sepsis due to P. aeruginosa can be achieved without iron provision by maintaining local phosphate abundance: role of pH.

BACKGROUND: During extreme physiological stress, the intestinal tract can be transformed into a harsh environment characterized by regio- spatial alterations in oxygen, pH, and phosphate concentration. When the human intestine is exposed to extreme medical interventions, the normal flora becomes replaced by pathogenic species whose virulence can be triggered by various physico-chemical cues leading to lethal sepsis. We previously demonstrated that phosphate depletion develops in the mouse intestine following surgical injury and triggers intestinal P. aeruginosa to express a lethal phenotype that can be prevented by oral phosphate ([Pi]) supplementation. RESULTS: In this study we examined the role of pH in the protective effect of [Pi] supplementation as it has been shown to be increased in the distal gut following surgical injury. Surgically injured mice drinking 25 mM [Pi] at pH 7.5 and intestinally inoculated with P. aeruginosa had increased mortality compared to mice drinking 25 mM [Pi] at pH 6.0 (p < 0.05). This finding was confirmed in C. elegans. Transcriptional analysis of P. aeruginosa demonstrated enhanced expression of various genes involved in media alkalization at pH 6.0 and a global increase in the expression of all iron-related genes at pH 7.5. Maintaining the pH at 6.0 via phosphate supplementation led to significant attenuation of iron-related genes as demonstrated by microarray and confirmed by QRT-PCR analyses. CONCLUSION: Taken together, these data demonstrate that increase in pH in distal intestine of physiologically stressed host colonized by P. aeruginosa can lead to the expression of siderophore-related virulence in bacteria that can be prevented without providing iron by maintaining local phosphate abundance at pH 6.0. This finding is particularly important as provision of exogenous iron has been shown to have untoward effects when administered to critically ill and septic patients. Given that phosphate, pH, and iron are near universal cues that dictate the virulence status of a broad range of microorganisms relevant to serious gut origin infection and sepsis in critically ill patients, the maintenance of phosphate and pH at appropriate physiologic levels to prevent virulence activation in a site specific manner can be considered as a novel anti-infective therapy in at risk patients.