Metabolic network analysis of pre-ASD newborns and 5-year-old children with autism spectrum disorder.

Classical metabolomic and new metabolic network methods were used to study the developmental features of autism spectrum disorder (ASD) in newborns (n = 205) and 5-year-old children (n = 53). Eighty percent of the metabolic impact in ASD was caused by 14 shared biochemical pathways that led to decreased anti-inflammatory and antioxidant defenses, and to increased physiologic stress molecules like lactate, glycerol, cholesterol, and ceramides. CIRCOS plots and a new metabolic network parameter, V ° net, revealed differences in both the kind and degree of network connectivity. Of 50 biochemical pathways and 450 polar and lipid metabolites examined, the developmental regulation of the purine network was most changed. Purine network hub analysis revealed a 17-fold reversal in typically developing children. This purine network reversal did not occur in ASD. These results revealed previously unknown metabolic phenotypes, identified new developmental states of the metabolic correlation network, and underscored the role of mitochondrial functional changes, purine metabolism, and purinergic signaling in autism spectrum disorder.

A validated pangenome-scale metabolic model for the Klebsiella pneumoniae species complex.

The Klebsiella pneumoniae species complex (KpSC) is a major source of nosocomial infections globally with high rates of resistance to antimicrobials. Consequently, there is growing interest in understanding virulence factors and their association with cellular metabolic processes for developing novel anti-KpSC therapeutics. Phenotypic assays have revealed metabolic diversity within the KpSC, but metabolism research has been neglected due to experiments being difficult and cost-intensive. Genome-scale metabolic models (GSMMs) represent a rapid and scalable in silico approach for exploring metabolic diversity, which compile genomic and biochemical data to reconstruct the metabolic network of an organism. Here we use a diverse collection of 507 KpSC isolates, including representatives of globally distributed clinically relevant lineages, to construct the most comprehensive KpSC pan-metabolic model to date, KpSC pan v2. Candidate metabolic reactions were identified using gene orthology to known metabolic genes, prior to manual curation via extensive literature and database searches. The final model comprised a total of 3550 reactions, 2403 genes and can simulate growth on 360 unique substrates. We used KpSC pan v2 as a reference to derive strain-specific GSMMs for all 507 KpSC isolates, and compared these to GSMMs generated using a prior KpSC pan-reference (KpSC pan v1) and two single-strain references. We show that KpSC pan v2 includes a greater proportion of accessory reactions (8.8 %) than KpSC pan v1 (2.5 %). GSMMs derived from KpSC pan v2 also generate more accurate growth predictions, with high median accuracies of 95.4 % (aerobic, n=37 isolates) and 78.8 % (anaerobic, n=36 isolates) for 124 matched carbon substrates. KpSC pan v2 is freely available at https://github.com/kelwyres/KpSC-pan-metabolic-model, representing a valuable resource for the scientific community, both as a source of curated metabolic information and as a reference to derive accurate strain-specific GSMMs. The latter can be used to investigate the relationship between KpSC metabolism and traits of interest, such as reservoirs, epidemiology, drug resistance or virulence, and ultimately to inform novel KpSC control strategies.

Predicting the primary infection source of Escherichia coli bacteremia using virulence-associated genes.

PURPOSE: To investigate the role of E. coli virulence-associated genes (VAGs) in predicting urinary tract infection (UTI) as the source of bacteremia in two distinct hospital populations, one with a large general catchment area and one dominated by referrals. METHODS: E. coli bacteremias identified at Department of Clinical Microbiology (DCM), Hvidovre Hospital and DCM, Rigshospitalet in the Capital Region of Denmark from October to December 2018. Using whole genome sequencing (WGS), we identified 358 VAGs from 224 E. coli bacteremia. For predictive analysis, VAGs were paired with clinical source of UTI from local bacteremia databases. RESULTS: VAGs strongly predicting of UTI as primary infection source of bacteremia were primarily found within the pap gene family. papX (PPV 96%, sensitivity 54%) and papGII (PPV 93%, sensitivity 56%) were found highly predictive, but showed low sensitivities. The strength of VAG predictions of UTI as source varied significantly between the two hospital populations. VAGs had weaker predictions in the tertiary referral center (Rigshospitalet), a disparity likely stemming from differences in patient population and department specialization. CONCLUSION: WGS data was used to predict the primary source of E. coli bacteremia and is an attempt on a new and different type of infection source identification. Genomic data showed potential to be utilized to predict the primary source of infection; however, discrepancy between the best performing profile of VAGs between acute care hospitals and tertiary hospitals makes it difficult to implement in clinical practice.

Metabolic features of treatment-refractory major depressive disorder with suicidal ideation.

Peripheral blood metabolomics was used to gain chemical insight into the biology of treatment-refractory Major Depressive Disorder with suicidal ideation, and to identify individualized differences for personalized care. The study cohort consisted of 99 patients with treatment-refractory major depressive disorder and suicidal ideation (trMDD-SI n = 52 females and 47 males) and 94 age- and sex-matched healthy controls (n = 48 females and 46 males). The median age was 29 years (IQR 22-42). Targeted, broad-spectrum metabolomics measured 448 metabolites. Fibroblast growth factor 21 (FGF21) and growth differentiation factor 15 (GDF15) were measured as biomarkers of mitochondrial dysfunction. The diagnostic accuracy of plasma metabolomics was over 90% (95%CI: 0.80-1.0) by area under the receiver operator characteristic (AUROC) curve analysis. Over 55% of the metabolic impact in males and 75% in females came from abnormalities in lipids. Modified purines and pyrimidines from tRNA, rRNA, and mRNA turnover were increased in the trMDD-SI group. FGF21 was increased in both males and females. Increased lactate, glutamate, and saccharopine, and decreased cystine provided evidence of reductive stress. Seventy-five percent of the metabolomic abnormalities found were individualized. Personalized deficiencies in CoQ10, flavin adenine dinucleotide (FAD), citrulline, lutein, carnitine, or folate were found. Pathways regulated by mitochondrial function dominated the metabolic signature. Peripheral blood metabolomics identified mitochondrial dysfunction and reductive stress as common denominators in suicidal ideation associated with treatment-refractory major depressive disorder. Individualized metabolic differences were found that may help with personalized management.

Global pathogenomic analysis identifies known and candidate genetic antimicrobial resistance determinants in twelve species.

Surveillance programs for managing antimicrobial resistance (AMR) have yielded thousands of genomes suited for data-driven mechanism discovery. We present a workflow integrating pangenomics, gene annotation, and machine learning to identify AMR genes at scale. When applied to 12 species, 27,155 genomes, and 69 drugs, we 1) find AMR gene transfer mostly confined within related species, with 925 genes in multiple species but just eight in multiple phylogenetic classes, 2) demonstrate that discovery-oriented support vector machines outperform contemporary methods at recovering known AMR genes, recovering 263 genes compared to 145 by Pyseer, and 3) identify 142 AMR gene candidates. Validation of two candidates in E. coli BW25113 reveals cases of conditional resistance: ΔcycA confers ciprofloxacin resistance in minimal media with D-serine, and frdD V111D confers ampicillin resistance in the presence of ampC by modifying the overlapping promoter. We expect this approach to be adaptable to other species and phenotypes.

Bactabolize is a tool for high-throughput generation of bacterial strain-specific metabolic models.

Metabolic capacity can vary substantially within a bacterial species, leading to ecological niche separation, as well as differences in virulence and antimicrobial susceptibility. Genome-scale metabolic models are useful tools for studying the metabolic potential of individuals, and with the rapid expansion of genomic sequencing there is a wealth of data that can be leveraged for comparative analysis. However, there exist few tools to construct strain-specific metabolic models at scale. Here, we describe Bactabolize, a reference-based tool which rapidly produces strain-specific metabolic models and growth phenotype predictions. We describe a pan reference model for the priority antimicrobial-resistant pathogen, Klebsiella pneumoniae, and a quality control framework for using draft genome assemblies as input for Bactabolize. The Bactabolize-derived model for K. pneumoniae reference strain KPPR1 performed comparatively or better than currently available automated approaches CarveMe and gapseq across 507 substrate and 2317 knockout mutant growth predictions. Novel draft genomes passing our systematically defined quality control criteria resulted in models with a high degree of completeness (≥99% genes and reactions captured compared to models derived from matched complete genomes) and high accuracy (mean 0.97, n=10). We anticipate the tools and framework described herein will facilitate large-scale metabolic modelling analyses that broaden our understanding of diversity within bacterial species and inform novel control strategies for priority pathogens.

Rev1 deficiency induces a metabolic shift in MEFs that can be manipulated by the NAD+ precursor nicotinamide riboside.

Replication stress, caused by Rev1 deficiency, is associated with mitochondrial dysfunction, and metabolic stress. However, the overall metabolic alterations and possible interventions to rescue the deficits due to Rev1 loss remain unclear. Here, we report that loss of Rev1 leads to intense changes in metabolites and that this can be manipulated by NAD + supplementation. Autophagy decreases in Rev1-/- mouse embryonic fibroblasts (MEFs) and can be restored by supplementing the NAD+ precursor nicotinamide riboside (NR). The abnormal mitochondrial morphology in Rev1-/- MEFs can be partially reversed by NR supplementation, which also protects the mitochondrial cristae from rotenone-induced degeneration. In nematodes rev-1 deficiency causes sensitivity to oxidative stress but this cannot be rescued by NR supplementation. In conclusion, Rev1 deficiency leads to metabolic dysregulation of especially lipid and nucleotide metabolism, impaired autophagy, and mitochondrial anomalies, and all of these phenotypes can be improved by NR replenishment in MEFs.

CyuR is a Dual Regulator for L-Cysteine Dependent Antimicrobial Resistance in Escherichia coli.

Hydrogen sulfide (H 2 S), mainly produced from L-cysteine (Cys), renders bacteria highly resistant to oxidative stress. This mitigation of oxidative stress was suggested to be an important survival mechanism to achieve antimicrobial resistance (AMR) in many pathogenic bacteria. CyuR (known as DecR or YbaO) is a recently characterized Cys-dependent transcription regulator, responsible for the activation of the cyuAP operon and generation of hydrogen sulfide from Cys. Despite its potential importance, the regulatory network of CyuR remains poorly understood. In this study, we investigated the roles of the CyuR regulon in a Cys-dependent AMR mechanism in E. coli strains. We found: 1) Cys metabolism has a significant role in AMR and its effect is conserved in many E. coli strains, including clinical isolates; 2) CyuR negatively controls the expression of mdlAB encoding a transporter that exports antibiotics such as cefazolin and vancomycin; 3) CyuR binds to a DNA sequence motif 'GAAwAAATTGTxGxxATTTsyCC' in the absence of Cys, confirmed by an in vitro binding assay; and 4) CyuR may regulate 25 additional genes as suggested by in silico motif scanning and transcriptome sequencing. Collectively, our findings expanded the understanding of the biological roles of CyuR relevant to antibiotic resistance associated with Cys.

Whole-genome sequences from wild-type and laboratory-evolved strains define the alleleome and establish its hallmarks.

The genomic diversity across strains of a species forms the genetic basis for differences in their behavior. A large-scale assessment of sequence variation has been made possible by the growing availability of strain-specific whole-genome sequences (WGS) and with the advent of large-scale databases of laboratory-acquired mutations. We define the Escherichia coli "alleleome" through a genome-scale assessment of amino acid (AA) sequence diversity in open reading frames across 2,661 WGS from wild-type strains. We observe a highly conserved alleleome enriched in mutations unlikely to affect protein function. In contrast, 33,000 mutations acquired in laboratory evolution experiments result in more severe AA substitutions that are rarely achieved by natural selection. Large-scale assessment of the alleleome establishes a method for the quantification of bacterial allelic diversity, reveals opportunities for synthetic biology to explore novel sequence space, and offers insights into the constraints governing evolution.

Comprehensive whole genome sequencing with hybrid assembly of multi-drug resistant Candida albicans isolate causing cerebral abscess.

Comprehensive whole genome sequencing (WGS) with hybrid assembly of a multi-drug resistant (MDR) Candida albicans (CA) isolate causing cerebral abscess was performed using Illumina paired end and Oxford Nanopore long read technologies. The innovative technologies utilized here enabled us to resolve fragmented assemblies, and implement comprehensive and detailed genomic analyses involved in antifungal resistance of Candida spp. Functionally important genes (MDR1, CDR2 and SQN2) involved in antifungal resistance were identified and a phylogenetic analysis of the clinical isolate was performed. Additionally, our clinical isolate was found to share 4 single nucleotide polymorphisms with two other sequenced strains of MDR C. auris (381 and 386) including translation elongation factor EF1α and EF3, ATPase activity associated proteins, and the lysine tRNA ligase.

Historical biomonitoring of pollution trends in the North Pacific using archived samples from the Continuous Plankton Recorder Survey.

First started in 1931, the Continuous Plankton Recorder (CPR) Survey is the longest-running and most geographically extensive marine plankton sampling program in the world. This pilot study investigates the feasibility of biomonitoring the spatiotemporal trends of marine pollution using archived CPR samples from the North Pacific. We selected specimens collected from three different locations (British Columbia Shelf, Northern Gulf of Alaska, and Aleutian Shelf) in the North Pacific between 2002 and 2020. Comprehensive profiling of the plankton chemical exposome was conducted using liquid and gas chromatography coupled with tandem mass spectrometry (LC-MS/MS and GC-MS/MS). Our results show that phthalates, plasticizers, persistent organic pollutants (POPs), pesticides, pharmaceuticals, and personal care products were present in the plankton exposome, and that many of these pollutants have decreased in amount over the last two decades, which was most pronounced for tri-n-butyl phosphate. In addition, the plankton exposome differed significantly by regional human activities, with the most polluted samples coming from the nearshore area. Exposome-wide association analysis revealed that bioaccumulation of environmental pollutants was highly correlated with the biomass of different plankton taxa. Overall, this study demonstrates that exposomic analysis of archived samples from the CPR Survey is effective for long-term biomonitoring of the spatial and temporal trends of environmental pollutants in the marine environment.

Metabolomic and exposomic biomarkers of risk of future neurodevelopmental delay in human milk.

BACKGROUND: The chemical composition of human milk has long-lasting effects on brain development. We examined the prognostic value of the human milk metabolome and exposome in children with the risk of neurodevelopmental delay (NDD). METHODS: This retrospective cohort study included 82 mother-infant pairs (40 male and 42 female infants). A total of 59 milk samples were from mothers with typically developing children and 23 samples were from mothers of children at risk. Milk samples were collected before 9 months of age (4.6 ± 2.5 months, mean ± SD). Neurocognitive development was assessed by maternal report at 14.2 ± 3.1 months using the Ages and Stages Questionnaires-2. RESULTS: Metabolome and exposome profiling identified 453 metabolites and 61 environmental chemicals in milk. Machine learning tools identified changes in deoxysphingolipids, phospholipids, glycosphingolipids, plasmalogens, and acylcarnitines in the milk of mothers with children at risk for future delay. A predictive classifier had a diagnostic accuracy of 0.81 (95% CI: 0.66-0.96) for females and 0.79 (95% CI: 0.62-0.94) for males. CONCLUSIONS: Once validated in larger studies, the chemical analysis of human milk might be added as an option in well-baby checks to help identify children at risk of NDD before the first symptoms appear. IMPACT: Maternal milk for infants sampled before 9 months of age contained sex-specific differences in deoxysphingolipids, sphingomyelins, plasmalogens, phospholipids, and acylcarnitines that predicted the risk of neurodevelopmental delay at 14.2 months of age. Once validated, this early biosignature in human milk might be incorporated into well-baby checks and help to identify infants at risk so early interventions might be instituted before the first symptoms appear.

High-quality genome-scale metabolic network reconstruction of probiotic bacterium Escherichia coli Nissle 1917.

BACKGROUND: Escherichia coli Nissle 1917 (EcN) is a probiotic bacterium used to treat various gastrointestinal diseases. EcN is increasingly being used as a chassis for the engineering of advanced microbiome therapeutics. To aid in future engineering efforts, our aim was to construct an updated metabolic model of EcN with extended secondary metabolite representation. RESULTS: An updated high-quality genome-scale metabolic model of EcN, iHM1533, was developed based on comparison with 55 E. coli/Shigella reference GEMs and manual curation, including expanded secondary metabolite pathways (enterobactin, salmochelins, aerobactin, yersiniabactin, and colibactin). The model was validated and improved using phenotype microarray data, resulting in an 82.3% accuracy in predicting growth phenotypes on various nutrition sources. Flux variability analysis with previously published 13C fluxomics data validated prediction of the internal central carbon fluxes. A standardised test suite called Memote assessed the quality of iHM1533 to have an overall score of 89%. The model was applied by using constraint-based flux analysis to predict targets for optimisation of secondary metabolite production. Modelling predicted design targets from across amino acid metabolism, carbon metabolism, and other subsystems that are common or unique for influencing the production of various secondary metabolites. CONCLUSION: iHM1533 represents a well-annotated metabolic model of EcN with extended secondary metabolite representation. Phenotype characterisation and the iHM1533 model provide a better understanding of the metabolic capabilities of EcN and will help future metabolic engineering efforts.

Experimental Evolution Reveals Unifying Systems-Level Adaptations but Diversity in Driving Genotypes.

Genotype-fitness maps of evolution have been well characterized for biological components, such as RNA and proteins, but remain less clear for systems-level properties, such as those of metabolic and transcriptional regulatory networks. Here, we take multi-omics measurements of 6 different E. coli strains throughout adaptive laboratory evolution (ALE) to maximal growth fitness. The results show the following: (i) convergence in most overall phenotypic measures across all strains, with the notable exception of divergence in NADPH production mechanisms; (ii) conserved transcriptomic adaptations, describing increased expression of growth promoting genes but decreased expression of stress response and structural components; (iii) four groups of regulatory trade-offs underlying the adjustment of transcriptome composition; and (iv) correlates that link causal mutations to systems-level adaptations, including mutation-pathway flux correlates and mutation-transcriptome composition correlates. We thus show that fitness landscapes for ALE can be described with two layers of causation: one based on system-level properties (continuous variables) and the other based on mutations (discrete variables). IMPORTANCE Understanding the mechanisms of microbial adaptation will help combat the evolution of drug-resistant microbes and enable predictive genome design. Although experimental evolution allows us to identify the causal mutations underlying microbial adaptation, it remains unclear how causal mutations enable increased fitness and is often explained in terms of individual components (i.e., enzyme rate) as opposed to biological systems (i.e., pathways). Here, we find that causal mutations in E. coli are linked to systems-level changes in NADPH balance and expression of stress response genes. These systems-level adaptation patterns are conserved across diverse E. coli strains and thus identify cofactor balance and proteome reallocation as dominant constraints governing microbial adaptation.

Genome-scale metabolic network reconstructions of diverse Escherichia strains reveal strain-specific adaptations.

Bottom-up approaches to systems biology rely on constructing a mechanistic basis for the biochemical and genetic processes that underlie cellular functions. Genome-scale network reconstructions of metabolism are built from all known metabolic reactions and metabolic genes in a target organism. A network reconstruction can be converted into a mathematical format and thus lend itself to mathematical analysis. Genome-scale models (GEMs) of metabolism enable a systems approach to characterize the pan and core metabolic capabilities of the Escherichia genus. In this work, GEMs were constructed for 222 representative strains of Escherichia across HC1100 levels spanning the known Escherichia phylogeny. The models were used to study Escherichia metabolic diversity and speciation on a large scale. The results show that unique strain-specific metabolic capabilities correspond to different species and nutrient niches. This work is a first step towards a curated reconstruction of pan-Escherichia metabolism. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Estimated Roles of the Carrier and the Bacterial Strain When Methicillin-Resistant Staphylococcus aureus Decolonization Fails: a Case-Control Study.

Methicillin-resistant Staphylococcus aureus (MRSA) is a common bacterial pathogen that frequently colonizes healthy individuals, with potential to cause invasive infection. In Denmark, to keep the prevalence low, MRSA carriers are recommended to undergo decolonization treatments, but achieving decolonization is challenging. Knowledge about the factors contributing to decolonization is scarce. We aimed to identify bacterial genome and clinical factors influencing MRSA decolonization. We identified all new MRSA patients above 2 years of age within the Hvidovre catchment area, Copenhagen, Denmark, in 2017 and 2018. Carriers were defined as chronic carriers (cases) if they were MRSA positive after two or more treatments and as nonchronic carriers (controls) if they were MRSA free after the first or second treatment. Using whole-genome sequencing (WGS), we constructed a pangenome of bacterial strains. With the incorporation of bacterial genome and clinical patient data, machine learning and multivariate analyses were performed to determine the factors associated with decolonization. A total of 477 MRSA carriers were included. An age of ≥13 years was significantly associated with nonchronic carriage. We identified 278 bacterial genetic features that were statistically significantly associated with chronic carriage (P < 0.05 by Fisher's exact test). Chronic MRSA carriage was predicted with 68% accuracy using a combination of bacterial genome data and patient clinical data. Decolonization success is multifactorial. Apart from the 68% predicted accuracy found in this study, we estimate that the remaining 32% is a result of host factors and microbiome composition. IMPORTANCE Carriage of methicillin-resistant Staphylococcus aureus (MRSA) and other multiresistant bacteria is a prerequisite for infection and transmission. Successful decolonization treatment removes these risks. We aimed to identify bacterial genome and host clinical factors that influence MRSA decolonization to estimate the roles of the carrier and the bacterial strain, respectively, when decolonization fails. The long-term goal, beyond this study, is to optimize decolonization success, minimize MRSA transmission, and, ultimately, improve the quality of life of MRSA carriers.

Pangenome analysis of Enterobacteria reveals richness of secondary metabolite gene clusters and their associated gene sets.

In silico genome mining provides easy access to secondary metabolite biosynthetic gene clusters (BGCs) encoding the biosynthesis of many bioactive compounds, which are the basis for many important drugs used in human medicine. However, the association between BGCs and other functions encoded in the genomes of producers have remained elusive. Here, we present a systems biology workflow that integrates genome mining with a detailed pangenome analysis for detecting genes associated with a particular BGC. We analyzed 3,889 enterobacterial genomes and found 13,266 BGCs, represented by 252 distinct BGC families and 347 additional singletons. A pangenome analysis revealed 88 genes putatively associated with a specific BGC coding for the colon cancer-related colibactin that code for diverse metabolic and regulatory functions. The presented workflow opens up the possibility to discover novel secondary metabolites, better understand their physiological roles, and provides a guide to identify and analyze BGC associated gene sets.

Cerebrospinal fluid and plasma metabolomics of acute endurance exercise.

Metabolomics has emerged as a powerful new tool in precision medicine. No studies have yet been published on the metabolomic changes in cerebrospinal fluid (CSF) produced by acute endurance exercise. CSF and plasma were collected from 19 young active adults (13 males and 6 females) before and 60 min after a 90-min monitored outdoor run. The median age, BMI, and VO2 max of subjects was 25 years (IQR 22-31), 23.2 kg/m2 (IQR 21.7-24.5), and 47 ml/kg/min (IQR 38-51), respectively. Targeted, broad-spectrum metabolomics was performed by liquid chromatography, tandem mass spectrometry (LC-MS/MS). In the CSF, purines and pyrimidines accounted for 32% of the metabolic impact after acute endurance exercise. Branch chain amino acids, amino acid neurotransmitters, fatty acid oxidation, phospholipids, and Krebs cycle metabolites traceable to mitochondrial function accounted for another 52% of the changes. A narrow but important channel of metabolic communication was identified between the brain and body by correlation network analysis. By comparing these results to previous work in experimental animal models, we found that over 80% of the changes in the CSF correlated with a cascade of mitochondrial and metabolic changes produced by ATP signaling. ATP is released as a co-neurotransmitter and neuromodulator at every synapse studied to date. By regulating brain mitochondrial function, ATP release was identified as an early step in the kinetic cascade of layered benefits produced by endurance exercise.

Transmission of Klebsiella strains and plasmids within and between grey-headed flying fox colonies.

The grey-headed flying fox (Pteropus poliocephalus) is an endemic Australian fruit bat, known to carry zoonotic pathogens. We recently showed they harbour bacterial pathogen Klebsiella pneumoniae and closely related species in the K. pneumoniae species complex (KpSC); however, the dynamics of KpSC transmission and gene flow within flying fox colonies are poorly understood. High-resolution genome comparisons of 39 KpSC isolates from grey-headed flying foxes identified five putative strain transmission clusters (four intra- and one inter-colony). The instance of inter-colony strain transmission of K. africana was found between two flying fox populations within flying distance, indicating either direct or indirect transmission through a common food/water source. All 11 plasmids identified within the KpSC isolates showed 73% coverage (mean) and ≥95% identity to human-associated KpSC plasmids, indicating gene flow between human clinical and grey-headed flying fox isolates. Along with strain transmission, inter-species horizontal plasmid transmission between K. pneumoniae and Klebsiella africana was also identified within a flying fox colony. Finally, genome-scale metabolic models were generated to predict and compare substrate usage to previously published KpSC models, from human and environmental sources. These models indicated no distinction on the basis of metabolic capabilities. Instead, metabolic capabilities were consistent with population structure and ST/lineage.

A systems approach discovers the role and characteristics of seven LysR type transcription factors in Escherichia coli.

Although Escherichia coli K-12 strains represent perhaps the best known model bacteria, we do not know the identity or functions of all of their transcription factors (TFs). It is now possible to systematically discover the physiological function of TFs in E. coli BW25113 using a set of synergistic methods; including ChIP-exo, growth phenotyping, conserved gene clustering, and transcriptome analysis. Among 47 LysR-type TFs (LTFs) found on the E. coli K-12 genome, many regulate nitrogen source utilization or amino acid metabolism. However, 19 LTFs remain unknown. In this study, we elucidated the regulation of seven of these 19 LTFs: YbdO, YbeF, YcaN, YbhD, YgfI, YiaU, YneJ. We show that: (1) YbdO (tentatively re-named CitR) regulation has an effect on bacterial growth at low pH with citrate supplementation. CitR is a repressor of the ybdNM operon and is implicated in the regulation of citrate lyase genes (citCDEFG); (2) YgfI (tentatively re-named DhfA) activates the dhaKLM operon that encodes the phosphotransferase system, DhfA is involved in formate, glycerol and dihydroxyacetone utilization; (3) YiaU (tentatively re-named LpsR) regulates the yiaT gene encoding an outer membrane protein, and waaPSBOJYZU operon is also important in determining cell density at the stationary phase and resistance to oxacillin microaerobically; (4) YneJ, re-named here as PtrR, directly regulates the expression of the succinate-semialdehyde dehydrogenase, Sad (also known as YneI), and is a predicted regulator of fnrS (a small RNA molecule). PtrR is important for bacterial growth in the presence of L-glutamate and putrescine as nitrogen/energy sources; and (5) YbhD and YcaN regulate adjacent y-genes on the genome. We have thus established the functions for four LTFs and identified the target genes for three LTFs.