Flexible genomic island conservation across freshwater and marine Methylophilaceae.

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Neoadjuvant systemic oncolytic vesicular stomatitis virus is safe and may enhance long-term survivorship in dogs with naturally occurring osteosarcoma.

Osteosarcoma is a devastating bone cancer that disproportionally afflicts children, adolescents, and young adults. Standard therapy includes surgical tumor resection combined with multiagent chemotherapy, but many patients still suffer from metastatic disease progression. Neoadjuvant systemic oncolytic virus (OV) therapy has the potential to improve clinical outcomes by targeting primary and metastatic tumor sites and inducing durable antitumor immune responses. Here we describe the first evaluation of neoadjuvant systemic therapy with a clinical-stage recombinant oncolytic vesicular stomatitis virus (VSV), VSV-IFNß-NIS, in naturally occurring cancer, specifically appendicular osteosarcoma in companion dogs. Canine osteosarcoma has a similar natural disease history as its human counterpart. VSV-IFNß-NIS was administered prior to standard of care surgical resection, permitting microscopic and genomic analysis of tumors. Treatment was well-tolerated and a "tail" of long-term survivors (∼35%) was apparent in the VSV-treated group, a greater proportion than observed in two contemporary control cohorts. An increase in tumor inflammation was observed in VSV-treated tumors and RNA-seq analysis showed that all the long-term responders had increased expression of a T cell anchored immune gene cluster. We conclude that neoadjuvant VSV-IFNß-NIS is safe and may increase long-term survivorship in dogs with naturally occurring osteosarcoma, particularly those that exhibit pre-existing antitumor immunity.

Draft genome sequence of the BAL58 Betaproteobacteria representative strain LSUCC0117.

Here, we present the draft genome sequence of strain LSUCC0117, a representative of the abundant aquatic BAL58 Betaproteobacteria group which we isolated from a coastal site in the northern Gulf of Mexico. The genome is estimated at over 99% complete, with a genome size of 2,687,225 bp.

Genome sequences of four agarolytic bacteria from the Bacteroidia and Gammaproteobacteria.

Here we present the genomes of four marine agarolytic bacteria belonging to the Bacteroidota and Proteobacteria. Two genomes are closed and two are in draft form, but all are at least 99% complete and offer new opportunities to study agar-degradation in marine bacteria.

The AEGEAN-169 clade of bacterioplankton is synonymous with SAR11 subclade V (HIMB59) and metabolically distinct.

Bacterioplankton of the SAR11 clade are the most abundant marine microorganisms and consist of numerous subclades spanning order-level divergence (Pelagibacterales). The assignment of the earliest diverging subclade V (a.k.a. HIMB59) to the Pelagibacterales is highly controversial, with multiple recent phylogenetic studies placing them completely separate from SAR11. Other than through phylogenomics, subclade V has not received detailed examination due to limited genomes from this group. Here, we assessed the ecogenomic characteristics of subclade V to better understand the role of this group in comparison to the Pelagibacterales. We used a new isolate genome, recently released single-amplified genomes and metagenome-assembled genomes, and previously established SAR11 genomes to perform a comprehensive comparative genomics analysis. We paired this analysis with the recruitment of metagenomes spanning the open ocean, coastal, and brackish systems. Phylogenomics, average amino acid identity, and 16S rRNA gene phylogeny indicate that SAR11 subclade V is synonymous with the ubiquitous AEGEAN-169 clade and support the contention that this group represents a taxonomic family. AEGEAN-169 shared many bulk genome qualities with SAR11, such as streamlining and low GC content, but genomes were generally larger. AEGEAN-169 had overlapping distributions with SAR11 but was metabolically distinct from SAR11 in its potential to transport and utilize a broader range of sugars as well as in the transport of trace metals and thiamin. Thus, regardless of the ultimate phylogenetic placement of AEGEAN-169, these organisms have distinct metabolic capacities that likely allow them to differentiate their niche from canonical SAR11 taxa. IMPORTANCE One goal of marine microbiologists is to uncover the roles various microorganisms are playing in biogeochemical cycles. Success in this endeavor relies on differentiating groups of microbes and circumscribing their relationships. An early-diverging group (subclade V) of the most abundant bacterioplankton, SAR11, has recently been proposed as a separate lineage that does not share a most recent common ancestor. But beyond phylogenetics, little has been done to evaluate how these organisms compare with SAR11. Our work leverages dozens of new genomes to demonstrate the similarities and differences between subclade V and SAR11. In our analysis, we also establish that subclade V is synonymous with a group of bacteria established from 16S rRNA gene sequences, AEGEAN-169. Subclade V/AEGEAN-169 has clear metabolic distinctions from SAR11 and their shared traits point to remarkable convergent evolution if they do not share a most recent common ancestor.

Neoadjuvant systemic oncolytic vesicular stomatitis virus is safe and may enhance long-term survivorship in dogs with naturally occurring osteosarcoma.

Osteosarcoma is a devastating bone cancer that disproportionally afflicts children, adolescents, and young adults. Standard therapy includes surgical tumor resection combined with multiagent chemotherapy, but many patients still suffer from metastatic disease progression. Neoadjuvant systemic oncolytic virus (OV) therapy has the potential to improve clinical outcomes by targeting primary and metastatic tumor sites and inducing durable antitumor immune responses. Here we described the first evaluation of neoadjuvant systemic therapy with a clinical-stage recombinant oncolytic Vesicular stomatitis virus (VSV), VSV-IFNß-NIS, in naturally occurring cancer, specifically appendicular osteosarcoma in companion dogs. Canine osteosarcoma has a similar natural disease history as its human counterpart. VSV-IFNß-NIS was administered prior to standard of care surgical resection, permitting microscopic and genomic analysis of tumors. Treatment was well-tolerated and a 'tail' of long-term survivors (~35%) was apparent in the VSV-treated group, a greater proportion than observed in two contemporary control cohorts. An increase in tumor inflammation was observed in VSV-treated tumors and RNAseq analysis showed that all the long-term responders had increased expression of a T-cell anchored immune gene cluster. We conclude that neoadjuvant VSV-IFNß-NIS is safe and may increase long-term survivorship in dogs with naturally occurring osteosarcoma, particularly those that exhibit pre-existing antitumor immunity.

Limits of entrainment of circadian neuronal networks.

Circadian rhythmicity lies at the center of various important physiological and behavioral processes in mammals, such as sleep, metabolism, homeostasis, mood changes, and more. Misalignment of intrinsic neuronal oscillations with the external day-night cycle can disrupt such processes and lead to numerous disorders. In this work, we computationally determine the limits of circadian synchronization to external light signals of different frequency, duty cycle, and simulated amplitude. Instead of modeling circadian dynamics with generic oscillator models (e.g., Kuramoto-type), we use a detailed computational neuroscience model, which integrates biomolecular dynamics, neuronal electrophysiology, and network effects. This allows us to investigate the effect of small drug molecules, such as Longdaysin, and connect our results with experimental findings. To combat the high dimensionality of such a detailed model, we employ a matrix-free approach, while our entire algorithmic pipeline enables numerical continuation and construction of bifurcation diagrams using only direct simulation. We, thus, computationally explore the effect of heterogeneity in the circadian neuronal network, as well as the effect of the corrective therapeutic intervention of Longdaysin. Last, we employ unsupervised learning to construct a data-driven embedding space for representing neuronal heterogeneity.

Ecophysiology and genomics of the brackish water adapted SAR11 subclade IIIa.

The Order Pelagibacterales (SAR11) is the most abundant group of heterotrophic bacterioplankton in global oceans and comprises multiple subclades with unique spatiotemporal distributions. Subclade IIIa is the primary SAR11 group in brackish waters and shares a common ancestor with the dominant freshwater IIIb (LD12) subclade. Despite its dominance in brackish environments, subclade IIIa lacks systematic genomic or ecological studies. Here, we combine closed genomes from new IIIa isolates, new IIIa MAGS from San Francisco Bay (SFB), and 460 highly complete publicly available SAR11 genomes for the most comprehensive pangenomic study of subclade IIIa to date. Subclade IIIa represents a taxonomic family containing three genera (denoted as subgroups IIIa.1, IIIa.2, and IIIa.3) that had distinct ecological distributions related to salinity. The expansion of taxon selection within subclade IIIa also established previously noted metabolic differentiation in subclade IIIa compared to other SAR11 subclades such as glycine/serine prototrophy, mosaic glyoxylate shunt presence, and polyhydroxyalkanoate synthesis potential. Our analysis further shows metabolic flexibility among subgroups within IIIa. Additionally, we find that subclade IIIa.3 bridges the marine and freshwater clades based on its potential for compatible solute transport, iron utilization, and bicarbonate management potential. Pure culture experimentation validated differential salinity ranges in IIIa.1 and IIIa.3 and provided detailed IIIa cell size and volume data. This study is an important step forward for understanding the genomic, ecological, and physiological differentiation of subclade IIIa and the overall evolutionary history of SAR11.

Evaluation of coated platelets, a subset of highly procoagulant platelets, in healthy dogs and dogs with neoplasia.

OBJECTIVE: To determine if dogs with neoplasia produce more coated platelets, a subpopulation of activated platelets generated by dual stimulation with thrombin and convulxin, a glycoprotein VI agonist, than healthy control dogs. ANIMALS: Client-owned dogs diagnosed with lymphoma (n = 19) or solid tumors (14) and healthy control dogs (14). PROCEDURES: Platelets were stimulated ex vivo with thrombin and convulxin. Flow cytometry was used to quantify the percentage of coated platelets based on high levels of surface fibrinogen. To compare the percentage of coated platelets between the three groups, an ANOVA was performed followed by pairwise 95% confidence intervals (CI) adjusted for multiple comparisons using Tukey's method. RESULTS: We observed a greater mean percentage of coated platelets in dogs with solid tumors, compared with healthy control dogs, by 10.9 percentage points (95% CI: -1.0, 22.8), and a mean percentage of coated platelets in dogs with lymphoma that was less than healthy control dogs by 0.3 percentage points (95% CI: -11.4, 10.8). CLINICAL RELEVANCE: This study provides the first data-based evidence that dogs with solid tumors may have a greater mean coated platelet percentage when compared with healthy control dogs, although there is overlap between groups. Further studies are needed investigating coated platelets in specific subsets of neoplasia and investigating additional mechanisms of hypercoagulability in dogs with neoplasia.

Draft Genome Sequence of the Marine Flavobacteriaceae sp. Strain LSUCC0859.

A new marine Flavobacteriaceae sp. strain, LSUCC0859, was isolated off the coast of Louisiana with artificial seawater via high-throughput dilution-to-extinction (DTE) cultivation. The 2,168,862-bp genome sequence provides opportunities to investigate the biology of a poorly understood lineage within the Bacteroidetes.

Metagenomics of Antarctic Marine Sediment Reveals Potential for Diverse Chemolithoautotrophy.

The microbial biogeochemical processes occurring in marine sediment in Antarctica remain underexplored due to limited access. Further, these polar habitats are unique, as they are being exposed to significant changes in their climate. To explore how microbes drive biogeochemistry in these sediments, we performed a shotgun metagenomic survey of marine surficial sediment (0 to 3 cm of the seafloor) collected from 13 locations in western Antarctica and assembled 16 high-quality metagenome assembled genomes for focused interrogation of the lifestyles of some abundant lineages. We observe an abundance of genes from pathways for the utilization of reduced carbon, sulfur, and nitrogen sources. Although organotrophy is pervasive, nitrification and sulfide oxidation are the dominant lithotrophic pathways and likely fuel carbon fixation via the reverse tricarboxylic acid and Calvin cycles. Oxygen-dependent terminal oxidases are common, and genes for reduction of oxidized nitrogen are sporadically present in our samples. Our results suggest that the underlying benthic communities are well primed for the utilization of settling organic matter, which is consistent with findings from highly productive surface water. Despite the genetic potential for nitrate reduction, the net catabolic pathway in our samples remains aerobic respiration, likely coupled to the oxidation of sulfur and nitrogen imported from the highly productive Antarctic water column above. IMPORTANCE The impacts of climate change in polar regions, like Antarctica, have the potential to alter numerous ecosystems and biogeochemical cycles. Increasing temperature and freshwater runoff from melting ice can have profound impacts on the cycling of organic and inorganic nutrients between the pelagic and benthic ecosystems. Within the benthos, sediment microbial communities play a critical role in carbon mineralization and the cycles of essential nutrients like nitrogen and sulfur. Metagenomic data collected from sediment samples from the continental shelf of western Antarctica help to examine this unique system and document the metagenomic potential for lithotrophic metabolisms and the cycles of both nitrogen and sulfur, which support not only benthic microbes but also life in the pelagic zone.

Model Systems to Study the Chronic, Polymicrobial Infections in Cystic Fibrosis: Current Approaches and Exploring Future Directions.

A recent workshop titled "Developing Models to Study Polymicrobial Infections," sponsored by the Dartmouth Cystic Fibrosis Center (DartCF), explored the development of new models to study the polymicrobial infections associated with the airways of persons with cystic fibrosis (CF). The workshop gathered 35+ investigators over two virtual sessions. Here, we present the findings of this workshop, summarize some of the challenges involved with developing such models, and suggest three frameworks to tackle this complex problem. The frameworks proposed here, we believe, could be generally useful in developing new model systems for other infectious diseases. Developing and validating new approaches to study the complex polymicrobial communities in the CF airway could open windows to new therapeutics to treat these recalcitrant infections, as well as uncovering organizing principles applicable to chronic polymicrobial infections more generally.

Modelling the functional roles of synaptic and extra-synaptic γ-aminobutyric acid receptor dynamics in circadian timekeeping.

In the suprachiasmatic nucleus (SCN), γ-aminobutyric acid (GABA) is a primary neurotransmitter. GABA can signal through two types of GABAA receptor subunits, often referred to as synaptic GABAA (gamma subunit) and extra-synaptic GABAA (delta subunit). To test the functional roles of these distinct GABAA in regulating circadian rhythms, we developed a multicellular SCN model where we could separately compare the effects of manipulating GABA neurotransmitter or receptor dynamics. Our model predicted that blocking GABA signalling modestly increased synchrony among circadian cells, consistent with published SCN pharmacology. Conversely, the model predicted that lowering GABAA receptor density reduced firing rate, circadian cell fraction, amplitude and synchrony among individual neurons. When we tested these predictions, we found that the knockdown of delta GABAA reduced the amplitude and synchrony of clock gene expression among cells in SCN explants. The model further predicted that increasing gamma GABAA densities could enhance synchrony, as opposed to increasing delta GABAA densities. Overall, our model reveals how blocking GABAA receptors can modestly increase synchrony, while increasing the relative density of gamma over delta subunits can dramatically increase synchrony. We hypothesize that increased gamma GABAA density in the winter could underlie the tighter phase relationships among SCN cells.

Interrogation of the perturbed gut microbiota in gouty arthritis patients through in silico metabolic modeling.

Recent studies have shown perturbed gut microbiota associated with gouty arthritis, a metabolic disease characterized by an imbalance between uric acid production and excretion. To mechanistically investigate altered microbiota metabolism associated with gout disease, 16S rRNA gene amplicon sequence data from stool samples of gout patients and healthy controls were computationally analyzed through bacterial community metabolic models. Patient-specific community models constructed with the metagenomics modeling pipeline, mgPipe, were used to perform k-means clustering of samples according to their metabolic capabilities. The clustering analysis generated statistically significant partitioning of samples into a Bacteroides-dominated, high gout cluster and a Faecalibacterium-elevated, low gout cluster. The high gout cluster was predicted to allow elevated synthesis of the amino acids D-alanine and L-alanine and byproducts of branched-chain amino acid catabolism, while the low gout cluster allowed higher production of butyrate, the sulfur-containing amino acids L-cysteine and L-methionine, and the L-cysteine catabolic product H2S. By expanding the capabilities of mgPipe to provide taxa-level resolution of metabolite exchange rates, acetate, D-lactate and succinate exchanged from Bacteroides to Faecalibacterium were predicted to enhance butyrate production in the low gout cluster. Model predictions suggested that sulfur-containing amino acid metabolism generally and H2S more specifically could be novel gout disease markers.

Ecophysiology of the Cosmopolitan OM252 Bacterioplankton (Gammaproteobacteria).

Among the thousands of species that comprise marine bacterioplankton communities, most remain functionally obscure. One key cosmopolitan group in this understudied majority is the OM252 clade of Gammaproteobacteria. Although frequently found in sequence data and even previously cultured, the diversity, metabolic potential, physiology, and distribution of this clade has not been thoroughly investigated. Here, we examined these features of OM252 bacterioplankton using a newly isolated strain and genomes from publicly available databases. We demonstrated that this group constitutes a globally distributed novel genus ("Candidatus Halomarinus"), sister to Litoricola, comprising two subclades and multiple distinct species. OM252 organisms have small genomes (median, 2.21 Mbp) and are predicted obligate aerobes capable of alternating between chemoorganoheterotrophic and chemolithotrophic growth using reduced sulfur compounds as electron donors. Subclade I genomes encode genes for the Calvin-Benson-Bassham cycle for carbon fixation. One representative strain of subclade I, LSUCC0096, had extensive halotolerance and a mesophilic temperature range for growth, with a maximum rate of 0.36 doublings/h at 35°C. Cells were curved rod/spirillum-shaped, ∼1.5 by 0.2 µm. Growth yield on thiosulfate as the sole electron donor under autotrophic conditions was roughly one-third that of heterotrophic growth, even though calculations indicated similar Gibbs energies for both catabolisms. These phenotypic data show that some "Ca. Halomarinus" organisms can switch between serving as carbon sources or sinks and indicate the likely anabolic cost of lithoautotrophic growth. Our results thus provide new hypotheses about the roles of these organisms in global biogeochemical cycling of carbon and sulfur. IMPORTANCE Marine microbial communities are teeming with understudied taxa due to the sheer numbers of species in any given sample of seawater. One group, the OM252 clade of Gammaproteobacteria, has been identified in gene surveys from myriad locations, and one isolated organism has even been genome sequenced (HIMB30). However, further study of these organisms has not occurred. Using another isolated representative (strain LSUCC0096) and publicly available genome sequences from metagenomic and single-cell genomic data sets, we examined the diversity within the OM252 clade and the distribution of these taxa in the world's oceans, reconstructed the predicted metabolism of the group, and quantified growth dynamics in LSUCC0096. Our results generate new knowledge about the previously enigmatic OM252 clade and point toward the importance of facultative chemolithoautotrophy for supporting some clades of ostensibly "heterotrophic" taxa.

Mechanisms driving genome reduction of a novel Roseobacter lineage.

Members of the marine Roseobacter group are key players in the global carbon and sulfur cycles. While over 300 species have been described, only 2% possess reduced genomes (mostly 3-3.5 Mbp) compared to an average roseobacter (>4 Mbp). These taxonomic minorities are phylogenetically diverse but form a Pelagic Roseobacter Cluster (PRC) at the genome content level. Here, we cultivated eight isolates constituting a novel Roseobacter lineage which we named 'CHUG'. Metagenomic and metatranscriptomic read recruitment analyses showed that CHUG members are globally distributed and active in marine pelagic environments. CHUG members possess some of the smallest genomes (~2.6 Mb) among all known roseobacters, but they do not exhibit canonical features of typical bacterioplankton lineages theorized to have undergone genome streamlining processes, like higher coding density, fewer paralogues and rarer pseudogenes. While CHUG members form a genome content cluster with traditional PRC members, they show important differences. Unlike other PRC members, neither the relative abundances of CHUG members nor their relative gene expression levels are correlated with chlorophyll a concentration across the global samples. CHUG members cannot utilize most phytoplankton-derived metabolites or synthesize vitamin B12, a key metabolite mediating the roseobacter-phytoplankton interactions. This combination of features is evidence for the hypothesis that CHUG members may have evolved a free-living lifestyle decoupled from phytoplankton. This ecological transition was accompanied by the loss of signature genes involved in roseobacter-phytoplankton symbiosis, suggesting that relaxation of purifying selection owing to lifestyle shift is likely an important driver of genome reduction in CHUG.

Metabolic Modeling to Interrogate Microbial Disease: A Tale for Experimentalists.

The explosion of microbiome analyses has helped identify individual microorganisms and microbial communities driving human health and disease, but how these communities function is still an open question. For example, the role for the incredibly complex metabolic interactions among microbial species cannot easily be resolved by current experimental approaches such as 16S rRNA gene sequencing, metagenomics and/or metabolomics. Resolving such metabolic interactions is particularly challenging in the context of polymicrobial communities where metabolite exchange has been reported to impact key bacterial traits such as virulence and antibiotic treatment efficacy. As novel approaches are needed to pinpoint microbial determinants responsible for impacting community function in the context of human health and to facilitate the development of novel anti-infective and antimicrobial drugs, here we review, from the viewpoint of experimentalists, the latest advances in metabolic modeling, a computational method capable of predicting metabolic capabilities and interactions from individual microorganisms to complex ecological systems. We use selected examples from the literature to illustrate how metabolic modeling has been utilized, in combination with experiments, to better understand microbial community function. Finally, we propose how such combined, cross-disciplinary efforts can be utilized to drive laboratory work and drug discovery moving forward.

One versus Many: Polymicrobial Communities and the Cystic Fibrosis Airway.

Culture-independent studies have revealed that chronic lung infections in persons with cystic fibrosis (pwCF) are rarely limited to one microbial species. Interactions among bacterial members of these polymicrobial communities in the airways of pwCF have been reported to modulate clinically relevant phenotypes. Furthermore, it is clear that a single polymicrobial community in the context of CF airway infections cannot explain the diversity of clinical outcomes. While large 16S rRNA gene-based studies have allowed us to gain insight into the microbial composition and predicted functional capacities of communities found in the CF lung, here we argue that in silico approaches can help build clinically relevant in vitro models of polymicrobial communities that can in turn be used to experimentally test and validate computationally generated hypotheses. Furthermore, we posit that combining computational and experimental approaches will enhance our understanding of mechanisms that drive microbial community function and identify new therapeutics to target polymicrobial infections.

Computational modeling of the gut microbiota reveals putative metabolic mechanisms of recurrent Clostridioides difficile infection.

Approximately 30% of patients who have Clostridioides difficile infection (CDI) will suffer at least one incident of reinfection. While the underlying causes of CDI recurrence are poorly understood, interactions between C. difficile and commensal gut bacteria are thought to play an important role. In this study, an in silico pipeline was used to process 16S rRNA gene amplicon sequence data of 225 stool samples from 93 CDI patients into sample-specific models of bacterial community metabolism. Clustered metabolite production rates generated from post-diagnosis samples generated a high Enterobacteriaceae abundance cluster containing disproportionately large numbers of recurrent samples and patients. This cluster was predicted to have significantly reduced capabilities for secondary bile acid synthesis but elevated capabilities for aromatic amino acid catabolism. When applied to 16S sequence data of 40 samples from fecal microbiota transplantation (FMT) patients suffering from recurrent CDI and their stool donors, the community modeling method generated a high Enterobacteriaceae abundance cluster with a disproportionate large number of pre-FMT samples. This cluster also was predicted to exhibit reduced secondary bile acid synthesis and elevated aromatic amino acid catabolism. Collectively, these in silico predictions suggest that Enterobacteriaceae may create a gut environment favorable for C. difficile spore germination and/or toxin synthesis.

Pseudomonas aeruginosa reverse diauxie is a multidimensional, optimized, resource utilization strategy.

Pseudomonas aeruginosa is a globally-distributed bacterium often found in medical infections. The opportunistic pathogen uses a different, carbon catabolite repression (CCR) strategy than many, model microorganisms. It does not utilize a classic diauxie phenotype, nor does it follow common systems biology assumptions including preferential consumption of glucose with an 'overflow' metabolism. Despite these contradictions, P. aeruginosa is competitive in many, disparate environments underscoring knowledge gaps in microbial ecology and systems biology. Physiological, omics, and in silico analyses were used to quantify the P. aeruginosa CCR strategy known as 'reverse diauxie'. An ecological basis of reverse diauxie was identified using a genome-scale, metabolic model interrogated with in vitro omics data. Reverse diauxie preference for lower energy, nonfermentable carbon sources, such as acetate or succinate over glucose, was predicted using a multidimensional strategy which minimized resource investment into central metabolism while completely oxidizing substrates. Application of a common, in silico optimization criterion, which maximizes growth rate, did not predict the reverse diauxie phenotypes. This study quantifies P. aeruginosa metabolic strategies foundational to its wide distribution and virulence including its potentially, mutualistic interactions with microorganisms found commonly in the environment and in medical infections.