One-shot 13 C15 N-metabolic flux analysis for simultaneous quantification of carbon and nitrogen flux

Optimizing in vitro spherulation cues in the fungal pathogen Coccidioides

Coccidioides spp. are part of a group of thermally dimorphic fungal pathogens, which grow as filamentous cells (hyphae) in the soil and transform to a different morphology upon inhalation into the host. The Coccidioides host form, the spherule, is unique and highly undercharacterized due to both technical and biocontainment challenges. Each spherule arises from an environmental spore (arthroconidium), matures, and develops hundreds of internal endospores, which are released from the spherule upon rupture. Each endospore can then go on to form another spherule in a cycle called spherulation. One of the foremost technical challenges has been reliably growing spherules in culture without the formation of contaminating hyphae and consistently inducing endospore release from spherules. Here, we present optimization of in vitro spherule growth and endospore release, by closely controlling starting cell density in the culture, using freshly harvested arthroconidia, and decreasing the concentration of multiple salts in spherulation media. We developed a minimal medium to test spherule growth on various carbon and nitrogen sources. We defined a critical role for the dispersant Tamol in both early spherule formation and prevention of the accumulation of a visible film around spherules. Finally, we examined how the conditions under which arthroconidia are generated influence their transcriptome and subsequent development into spherules, demonstrating that this is an important variable to control when designing spherulation experiments. Together, our data reveal multiple strategies to optimize in vitro spherulation growth, enabling characterization of this virulence-relevant morphology.IMPORTANCECoccidioides spp. are thermally dimorphic fungal pathogens found in the Southwest United States, Mexico, Central America, and South America. Coccidioides can infect both immunocompetent and immunocompromised people and can cause a devastating disseminated infection, including meningitis, with 30% mortality despite all currently available treatments. In this work, we tackle one of the current largest technical barriers to studying the fungus Coccidioides: reliably growing its host form in vitro. Our work is impactful because we have created a set of foundational tools for the burgeoning field of Coccidioides pathogenesis research. We have carefully optimized conditions that allow the development of Coccidioides in vitro into its pathogenic form. This work will open up many lines of investigation into the molecules that underlie Coccidioides pathogenesis.

Catalase activity deficiency sensitizes multidrug-resistant Mycobacterium tuberculosis to the ATP synthase inhibitor bedaquiline

Multidrug-resistant tuberculosis (MDR-TB), defined as resistance to the first-line drugs isoniazid and rifampin, is a growing source of global mortality and threatens global control of tuberculosis disease. The diarylquinoline bedaquiline has recently emerged as a highly efficacious drug against MDR-TB and kills Mycobacterium tuberculosis by inhibiting mycobacterial ATP synthase. However, the mechanisms underlying bedaquiline's efficacy against MDR-TB remain unknown. Here we investigate bedaquiline hyper-susceptibility in drug-resistant Mycobacterium tuberculosis using systems biology approaches. We discovered that MDR clinical isolates are commonly sensitized to bedaquiline. This hypersensitization is caused by several physiological changes induced by deficient catalase activity. These include enhanced accumulation of reactive oxygen species, increased susceptibility to DNA damage, induction of sensitizing transcriptional programs, and metabolic repression of several biosynthetic pathways. In this work we demonstrate how resistance-associated changes in bacterial physiology can mechanistically induce collateral antimicrobial drug sensitivity and reveal druggable vulnerabilities in antimicrobial resistant pathogens.

Enhanced catabolism of glycine betaine and derivatives provides improved osmotic stress protection in Methylorubrum extorquens PA1

Integration of metabolites into the overall metabolic network of a cell requires careful coordination dependent upon the ultimate usage of the metabolite. Different stoichiometric needs, and thus pathway fluxes, must exist for compounds destined for diverse uses, such as carbon sources, nitrogen sources, or stress-protective agents. Herein, we expand upon our previous work that highlighted the nature of glycine betaine (GB) metabolism in Methylobacteria to examine the utilization of GB-derivative compounds dimethylglycine (DMG) and sarcosine into Methylorubrum extorquens in different metabolic capacities, including as sole nitrogen and/or carbon sources. We isolated gain-of-function mutations that allowed M. extorquens PA1 to utilize dimethylglycine as a carbon source and dimethylglycine and sarcosine as nitrogen source. Characterization of mutants demonstrated selection for variants of the AraC-like regulator Mext_3735 that confer constitutive expression of the GB metabolic gene cluster, allowing direct utilization of the downstream GB derivatives. Finally, among the distinct isolates examined, we found that catabolism of the osmoprotectant used for selection (GB or dimethylglycine) enhanced osmotic stress resistance provided in the presence of that particular osmolyte. Thus, access to the carbon and nitrogen and osmoprotective effects of GB and DMG are made readily accessible through adaptive mutations. In M. extorquens PA1, the limitations to exploiting this group of compounds appear to exist predominantly at the levels of gene regulation and functional activity, rather than being constrained by transport or toxicity.IMPORTANCEOsmotic stress is a common challenge for bacteria colonizing the phyllosphere, where glycine betaine (GB) can be found as a prevalent osmoprotectant. Though Methylorubrum extorquens PA1 cannot use GB or its demethylation products, dimethylglycine (DMG) and sarcosine, as a sole carbon source, utilization is highly selectable via single nucleotide changes for both GB and DMG growth. The innate inability to use these compounds is due to limited flux through steps in the pathway and regulatory constraints. Herein, the characterization of the transcriptional regulator, Mext_3735 (GbdR), expands our understanding of the various roles in which GB derivatives can be used in M. extorquens PA1. Interestingly, increased catabolism of GB and derivatives does not interfere with, but rather improves, the ability of cells to thrive under increased salt stress conditions, suggesting that metabolic flux improves stress tolerance rather than providing a distinct tension between uses.

Optimizing in vitro spherulation cues in the fungal pathogen Coccidioides

Coccidioides spp . are part of a group of thermally dimorphic fungal pathogens, which grow as filamentous cells (hyphae) in the soil and transform to a different morphology upon inhalation into the host. The Coccidioides host form, the spherule, is unique and highly under characterized due to both technical and biocontainment challenges. Each spherule arises from an environmental spore (arthroconidium), matures, and develops hundreds of internal endospores, which are released from the spherule upon rupture. Each endospore can then go on to form another spherule in a cycle called spherulation. One of the foremost technical challenges has been reliably growing spherules in culture without the formation of contaminating hyphae, and consistently inducing endospore release from spherules. Here, we present optimization of in vitro spherule growth and endospore release, by closely controlling starting cell density in the culture, using freshly-harvested arthroconidia, and decreasing the concentration of multiple salts in spherulation media. We developed a minimal media to test spherule growth on various carbon and nitrogen sources. We defined a critical role for the dispersant Tamol in both early spherule formation and prevention of the accumulation of a visible film around spherules. Finally, we examined how the conditions under which arthroconidia are generated influence their transcriptome and subsequent development into spherules, demonstrating that this is an important variable to control when designing spherulation experiments. Together, our data reveal multiple strategies to optimize in vitro spherulation growth, enabling characterization of this virulence-relevant morphology.

Isotopically nonstationary metabolic flux analysis of plants: recent progress and future opportunities

Metabolic flux analysis (MFA) is a valuable tool for quantifying cellular phenotypes and to guide plant metabolic engineering. By introducing stable isotopic tracers and employing mathematical models, MFA can quantify the rates of metabolic reactions through biochemical pathways. Recent applications of isotopically nonstationary MFA (INST-MFA) to plants have elucidated nonintuitive metabolism in leaves under optimal and stress conditions, described coupled fluxes for fast-growing algae, and produced a synergistic multi-organ flux map that is a first in MFA for any biological system. These insights could not be elucidated through other approaches and show the potential of INST-MFA to correct an oversimplified understanding of plant metabolism.

Rethinking 13C-metabolic flux analysis - The Bayesian way of flux inference

Metabolic reaction rates (fluxes) play a crucial role in comprehending cellular phenotypes and are essential in areas such as metabolic engineering, biotechnology, and biomedical research. The state-of-the-art technique for estimating fluxes is metabolic flux analysis using isotopic labelling (13C-MFA), which uses a dataset-model combination to determine the fluxes. Bayesian statistical methods are gaining popularity in the field of life sciences, but the use of 13C-MFA is still dominated by conventional best-fit approaches. The slow take-up of Bayesian approaches is, at least partly, due to the unfamiliarity of Bayesian methods to metabolic engineering researchers. To address this unfamiliarity, we here outline similarities and differences between the two approaches and highlight particular advantages of the Bayesian way of flux analysis. With a real-life example, re-analysing a moderately informative labelling dataset of E. coli, we identify situations in which Bayesian methods are advantageous and more informative, pointing to potential pitfalls of current 13C-MFA evaluation approaches. We propose the use of Bayesian model averaging (BMA) for flux inference as a means of overcoming the problem of model uncertainty through its tendency to assign low probabilities to both, models that are unsupported by data, and models that are overly complex. In this capacity, BMA resembles a tempered Ockham's razor. With the tempered razor as a guide, BMA-based 13C-MFA alleviates the problem of model selection uncertainty and is thereby capable of becoming a game changer for metabolic engineering by uncovering new insights and inspiring novel approaches.

SIMPEL: using stable isotopes to elucidate dynamics of context specific metabolism

The capacity to leverage high resolution mass spectrometry (HRMS) with transient isotope labeling experiments is an untapped opportunity to derive insights on context-specific metabolism, that is difficult to assess quantitatively. Tools are needed to comprehensively mine isotopologue information in an automated, high-throughput way without errors. We describe a tool, Stable Isotope-assisted Metabolomics for Pathway Elucidation (SIMPEL), to simplify analysis and interpretation of isotope-enriched HRMS datasets. The efficacy of SIMPEL is demonstrated through examples of central carbon and lipid metabolism. In the first description, a dual-isotope labeling experiment is paired with SIMPEL and isotopically nonstationary metabolic flux analysis (INST-MFA) to resolve fluxes in central metabolism that would be otherwise challenging to quantify. In the second example, SIMPEL was paired with HRMS-based lipidomics data to describe lipid metabolism based on a single labeling experiment. Available as an R package, SIMPEL extends metabolomics analyses to include isotopologue signatures necessary to quantify metabolic flux.

Bayesian 13C-Metabolic Flux Analysis of Parallel Tracer Experiments in Granulocytes: A Directional Shift within the Non-Oxidative Pentose Phosphate Pathway Supports Phagocytosis

The pentose phosphate pathway (PPP) plays a key role in the cellular regulation of immune function; however, little is known about the interplay of metabolic adjustments in granulocytes, especially regarding the non-oxidative PPP. For the determination of metabolic mechanisms within glucose metabolism, we propose a novel set of measures for 13C-metabolic flux analysis based on ex vivo parallel tracer experiments ([1,2-13C]glucose, [U-13C]glucose, [4,5,6-13C]glucose) and gas chromatography-mass spectrometry labeling measurements of intracellular metabolites, such as sugar phosphates and their fragments. A detailed constraint analysis showed that the permission range for net and irreversible fluxes was limited to a three-dimensional space. The overall workflow, including its Bayesian flux estimation, resulted in precise flux distributions and pairwise confidence intervals, some of which could be represented as a line due to the strength of their correlation. The principal component analysis that was enabled by these behaviors comprised three components that explained 99.6% of the data variance. It showed that phagocytic stimulation reversed the direction of non-oxidative PPP net fluxes from ribose-5-phosphate biosynthesis toward glycolytic pathways. This process was closely associated with the up-regulation of the oxidative PPP to promote the oxidative burst.

KatG catalase deficiency confers bedaquiline hyper-susceptibility to isoniazid resistant Mycobacterium tuberculosis

Multidrug-resistant tuberculosis (MDR-TB) is a growing source of global mortality and threatens global control of tuberculosis (TB) disease. The diarylquinoline bedaquiline (BDQ) recently emerged as a highly efficacious drug against MDR-TB, defined as resistance to the first-line drugs isoniazid (INH) and rifampin. INH resistance is primarily caused by loss-of-function mutations in the catalase KatG, but mechanisms underlying BDQ's efficacy against MDR-TB remain unknown. Here we employ a systems biology approach to investigate BDQ hyper-susceptibility in INH-resistant Mycobacterium tuberculosis . We found hyper-susceptibility to BDQ in INH-resistant cells is due to several physiological changes induced by KatG deficiency, including increased susceptibility to reactive oxygen species and DNA damage, remodeling of transcriptional programs, and metabolic repression of folate biosynthesis. We demonstrate BDQ hyper-susceptibility is common in INH-resistant clinical isolates. Collectively, these results highlight how altered bacterial physiology can impact drug efficacy in drug-resistant bacteria.

Elucidating dynamic anaerobe metabolism with HRMAS 13C NMR and genome-scale modeling

Anaerobic microbial metabolism drives critical functions within global ecosystems, host-microbiota interactions, and industrial applications, yet remains ill-defined. Here we advance a versatile approach to elaborate cellular metabolism in obligate anaerobes using the pathogen Clostridioides difficile, an amino acid and carbohydrate-fermenting Clostridia. High-resolution magic angle spinning nuclear magnetic resonance (NMR) spectroscopy of C. difficile, grown with fermentable 13C substrates, informed dynamic flux balance analysis (dFBA) of the pathogen's genome-scale metabolism. Analyses identified dynamic recruitment of oxidative and supporting reductive pathways, with integration of high-flux amino acid and glycolytic metabolism at alanine's biosynthesis to support efficient energy generation, nitrogen handling and biomass generation. Model predictions informed an approach leveraging the sensitivity of 13C NMR spectroscopy to simultaneously track cellular carbon and nitrogen flow from [U-13C]glucose and [15N]leucine, confirming the formation of [13C,15N]alanine. Findings identify metabolic strategies used by C. difficile to support its rapid colonization and expansion in gut ecosystems.