Bacillus subtilis engineered for topical delivery of an antifungal agent

Fungal skin infections are a common condition affecting 20-25 percent of the world population. While these conditions are treatable with regular application of an antifungal medication, we sought to develop a more convenient, longer-lasting topical antifungal platform that could increase patient adherence to treatment regimens by using Bacillus subtilis, a naturally antifungal bacteria found on the skin, for drug production and delivery. In this study, we engineered B. subtilis for increased production of the antifungal lipopeptide iturin A by overexpression of the pleiotropic regulator DegQ. The engineered strain had an over 200% increase in iturin A production as detected by HPLC, accompanied by slower growth but the same terminal cell density as determined by absorbance measurements of liquid culture. In an in vitro antifungal assay, we found that despite its higher iturin A production, the engineered strain was less effective at reducing the growth of a plug of the pathogenic fungus Trichophyton mentagrophytes on an agar plate compared to the parent strain. The reduced efficacy of the engineered strain may be explained by its reduced growth rate, which highlights the need to address trade-offs between titers (e.g. measured drug production) and other figures of merit (e.g. growth rate) during metabolic engineering.

Complex Dependence of Escherichia coli-based Cell-Free Expression on Sonication Energy During Lysis

Cell lysis─by sonication or bead beating, for example─is a key step in preparing extracts for cell-free expression systems. To create high protein-production capacity extracts, standard practice is to lyse cells sufficiently to thoroughly disrupt the membrane and thus extract expression machinery but without degrading that machinery. Here, we investigate the impact of different sonication energy inputs on the protein-production capacity of Escherichia coli extracts. While the existence of operator-specific optimal sonication energy inputs is widely known, our findings show that the sonication energy input that yields maximal protein output from a given expression template may depend on plasmid concentration, transcriptional and translational features (e.g., promoter), and other expression vector components (e.g., origin of replication). These results indicate that sonication protocols cannot be standardized to a single optimum, suggest strategies for improving protein yields, and more broadly highlight the need for better metrics and protocols for characterizing cell extracts.

Plasmid Crosstalk in Cell-Free Expression Systems

Although cell-free protein expression has been widely used for the synthesis of single proteins, cell-free synthetic biology has rapidly expanded to new, more complex applications. One such application is the prototyping or implementation of complex genetic networks involving the expression of multiple proteins at precise ratios, often from different plasmids. However, expression of multiple proteins from multiple plasmids may inadvertently result in unexpected, off-target changes to the levels of the proteins being expressed, a phenomenon termed plasmid crosstalk. Here, we show that the effects of plasmid crosstalk─even at the qualitative level of increases vs decreases in protein expression─depend on the concentration of plasmids in the reaction and the type of transcriptional machinery involved in the expression. This crosstalk can have a significant impact on genetic circuitry function and even interpretation of simple experimental results and thus should be taken into consideration during the development of cell-free applications.

A Cell-Free Biosensor for Assessment of Hyperhomocysteinemia

Hyperhomocysteinemia─a condition characterized by elevated levels of homocysteine in the blood─is associated with multiple health conditions including folate deficiency and birth defects, but there are no convenient, low-cost methods to measure homocysteine in plasma. A cell-free biosensor that harnesses the native homocysteine sensing machinery of Escherichia coli bacteria could satisfy the need for a detection platform with these characteristics. Here, we describe our efforts to engineer a cell-free biosensor for point-of-care, low-cost assessment of homocysteine status. This biosensor can detect physiologically relevant concentrations of homocysteine in plasma with a colorimetric output visible to the naked eye in under 1.5 h, making it a fast, convenient tool for point-of-use diagnosis and monitoring of hyperhomocysteinemia and related health conditions.

Introduction to Engineering Biology: A Conceptual Framework for Teaching Synthetic Biology

As the impacts of engineering biology grow, it is important to introduce the field early and in an accessible way. However, teaching engineering biology poses challenges, such as limited representation of the field in widely used scientific textbooks or curricula, and the interdisciplinary nature of the subject. We have created an adaptable curriculum module that can be used by anyone to teach the basic principles and applications of engineering biology. The module consists of a versatile, concept-based slide deck designed by experts across engineering biology to cover key topic areas. Starting with the design, build, test, and learn cycle, the slide deck covers the framework, core tools, and applications of the field at an undergraduate level. The module is available for free on a public website and can be used in a stand-alone fashion or incorporated into existing curricular materials. Our aim is that this modular, accessible slide deck will improve the ease of teaching current engineering biology topics and increase public engagement with the field.

Short Activators and Repressors of RNA Toehold Switches

RNA toehold switches are a widely used class of molecule to detect specific RNA "trigger" sequences, but their design, intended function, and characterization to date leave it unclear whether they can function properly with triggers shorter than 36 nucleotides. Here, we explore the feasibility of using standard toehold switches with 23-nucleotide truncated triggers. We assess the crosstalk of different triggers with significant homology and identify a highly sensitive trigger region where just one mutation from the consensus trigger sequence can reduce switch activation by 98.6%. However, we also find that triggers with as many as seven mutations outside of this region can still lead to 5-fold induction of the switch. We also present a new approach using 18- to 22-nucleotide triggers as translational repressors for toehold switches and assess the off-target regulation for this strategy as well. The development and characterization of these strategies could help enable applications like microRNA sensors, where well-characterized crosstalk between sensors and detection of short target sequences are critical.

Towards inferring absolute concentrations from relative abundance in time-course GC-MS metabolomics data

Metabolomics, the large-scale study of metabolites, has significant appeal as a source of information for metabolic modeling and other scientific applications. One common approach for measuring metabolomics data is gas chromatography-mass spectrometry (GC-MS). However, GC-MS metabolomics data are typically reported as relative abundances, precluding their use with approaches and tools where absolute concentrations are necessary. While chemical standards can be used to help provide quantification, their use is time-consuming, expensive, or even impossible due to their limited availability. The ability to infer absolute concentrations from GC-MS metabolomics data without chemical standards would have significant value. We hypothesized that when analyzing time-course metabolomics datasets, the mass balances of metabolism and other biological information could provide sufficient information towards inference of absolute concentrations. To demonstrate this, we developed and characterized MetaboPAC, a computational framework that uses two approaches-one based on kinetic equations and another using biological heuristics-to predict the most likely response factors that allow translation between relative abundances and absolute concentrations. When used to analyze noiseless synthetic data generated from multiple types of kinetic rate laws, MetaboPAC performs significantly better than negative control approaches when 20% of kinetic terms are known a priori. Under conditions of lower sampling frequency and high noise, MetaboPAC is still able to provide significant inference of concentrations in 3 of 4 models studied. This provides a starting point for leveraging biological knowledge to extract concentration information from time-course intracellular GC-MS metabolomics datasets, particularly for systems that are well-studied and have partially known kinetic structures.

Rapid and Finely-Tuned Expression for Deployable Sensing Applications

Organisms from across the tree of life have evolved highly efficient mechanisms for sensing molecules of interest using biomolecular machinery that can in turn be quite valuable for the development of biosensors. However, purification of such machinery for use in in vitro biosensors is costly, while the use of whole cells as in vivo biosensors often leads to long sensor response times and unacceptable sensitivity to the chemical makeup of the sample. Cell-free expression systems overcome these weaknesses by removing the requirements associated with maintaining living sensor cells, allowing for increased function in toxic environments and rapid sensor readout at a production cost that is often more reasonable than purification. Here, we focus on the challenge of implementing cell-free protein expression systems that meet the stringent criteria required for them to serve as the basis for field-deployable biosensors. Fine-tuning expression to meet these requirements can be achieved through careful selection of the sensing and output elements, as well as through optimization of reaction conditions via tuning of DNA/RNA concentrations, lysate preparation methods, and buffer conditions. Through careful sensor engineering, cell-free systems can continue to be successfully used for the production of tightly regulated, rapidly expressing genetic circuits for biosensors.

MaHPIC malaria systems biology data from Plasmodium cynomolgi sporozoite longitudinal infections in macaques

Plasmodium cynomolgi causes zoonotic malarial infections in Southeast Asia and this parasite species is important as a model for Plasmodium vivax and Plasmodium ovale. Each of these species produces hypnozoites in the liver, which can cause relapsing infections in the blood. Here we present methods and data generated from iterative longitudinal systems biology infection experiments designed and performed by the Malaria Host-Pathogen Interaction Center (MaHPIC) to delve deeper into the biology, pathogenesis, and immune responses of P. cynomolgi in the Macaca mulatta host. Infections were initiated by sporozoite inoculation. Blood and bone marrow samples were collected at defined timepoints for biological and computational experiments and integrative analyses revolving around primary illness, relapse illness, and subsequent disease and immune response patterns. Parasitological, clinical, haematological, immune response, and -omic datasets (transcriptomics, proteomics, metabolomics, and lipidomics) including metadata and computational results have been deposited in public repositories. The scope and depth of these datasets are unprecedented in studies of malaria, and they are projected to be a F.A.I.R., reliable data resource for decades.

Harnessing Escherichia coli's Native Machinery for Detection of Vitamin C (Ascorbate) Deficiency

Vitamin C (l-ascorbate) deficiency is a global public health issue most prevalent in resource-limited regions, creating a need for an inexpensive detection platform. Here, we describe efforts to engineer whole-cell and cell-free ascorbate biosensors. Both sensors used the protein UlaR, which binds to a metabolite of ascorbate and regulates transcription. The whole-cell sensor could detect lower, physiologically relevant concentrations of ascorbate, which we attributed to intact functionality of a phosphotransferase system (PTS) that transports ascorbate across the cell membrane and phosphorylates it to form UlaR's ligand. We used multiple strategies to enhance cell-free PTS functionality (which has received little previous attention), improving the cell-free sensor's performance, but the whole-cell sensor remained more sensitive. These efforts demonstrated an advantage of whole-cell sensors for detection of molecules─like ascorbate─transformed by a PTS, but also proof of principle for cell-free sensors requiring membrane-bound components like the PTS. In addition, the cell-free sensor was functional in plasma, setting the stage for future implementation of ascorbate sensors for clinically relevant biofluids in field-deployable formats.

Nucleic acid partitioning in PEG-Ficoll protocells

The phase separation of aqueous polymer solutions is a widely used method for producing self-assembled, membraneless droplet protocells. Non-ionic synthetic polymers forming an aqueous two-phase system (ATPS) have been shown to reliably form protocells that, when equipped with biological materials, are useful for applications such as analyte detection. Previous characterization of an ATPS-templated protocell did not investigate the effects of its biological components on phase stability. Here we report the phase diagram of a PEG 35k-Ficoll 400k-water ATPS at baseline and in the presence of necessary protocell components. Because the stability of an ATPS can be sensitive to small changes in composition, which in turn impacts solute partitioning, we present partitioning data of a variety of nucleic acids in response to protocell additives. The results show that the additives-particularly a mixture of salts and small organic molecules-have profound positive effects on ATPS stability and nucleic acid partitioning, both of which significantly contribute to protocell function. Our data uncovers several new areas of optimization for future protocell engineering.

Systems biology-based analysis of cell-free systems

Cell-free expression systems are becoming increasingly widely used due to their diverse applications in biotechnology. Despite this rapid expansion in adoption, many aspects of cell-free systems remain surprisingly poorly understood. Systems biology approaches make it possible to characterize cell-free systems deeply and broadly to better understand their underlying complexity. Here, we review recent systems biology studies that have provided insight into cell-free systems. We focus on characterization of the cell-free proteome, including its dependence on preparation protocol and host strain, as well as the cell-free metabolome and the relationship of endogenous metabolism to system performance. We conclude by highlighting promising future research directions.

Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework

Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, Linear Kinetics-Dynamic Flux Balance Analysis (LK-DFBA), that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this framework's success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that when genetic perturbations were implemented in the systems, the optimal constraint approach typically remained the same as for the wild-type regardless of the model topology or parameterization, indicating that just a single wild-type dataset could allow identification of the ideal constraint to enable model predictivity for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems.

Metabolomics Analysis of Cell-Free Expression Systems Using Gas Chromatography-Mass Spectrometry

Metabolomics is the systems-scale study of the biochemical intermediates of metabolism. This approach has great potential to uncover how metabolic intermediates are used and generated in cell-free expression systems, something that is to date not well understood. Here, we present a detailed metabolomics protocol for characterization of the small molecules in cell-free systems. We specifically focus on the analysis of Escherichia coli lysate-based cell-free systems using gas chromatography coupled to mass spectrometry. Measuring and monitoring the metabolic changes in cell-free systems can provide insight into the ways that metabolites affect the productivity of cell-free reactions, ultimately allowing for more informed engineering and optimization efforts for cell-free systems.

Point-of-Care Analyte Quantification and Digital Readout via Lysate-Based Cell-Free Biosensors Interfaced with Personal Glucose Monitors

Field-deployable diagnostics based on cell-free systems have advanced greatly, but on-site quantification of target analytes remains a challenge. Here we demonstrate that Escherichia coli lysate-based cell-free biosensors coupled to a personal glucose monitor (PGM) can enable on-site analyte quantification, with the potential for straightforward reconfigurability to diverse types of analytes. We show that analyte-responsive regulators of transcription and translation can modulate the production of the reporter enzyme β-galactosidase, which in turn converts lactose into glucose for PGM quantification. Because glycolysis is active in the lysate and would readily deplete converted glucose, we decoupled enzyme production and glucose conversion to increase the end point signal output. However, this lysate metabolism did allow for one-pot removal of glucose present in complex samples (like human serum) without confounding target quantification. Taken together, our results show that integrating lysate-based cell-free biosensors with PGMs enables accessible target detection and quantification at the point of need.

Dramatic transcriptomic differences in Macaca mulatta and Macaca fascicularis with Plasmodium knowlesi infections

Plasmodium knowlesi, a model malaria parasite, is responsible for a significant portion of zoonotic malaria cases in Southeast Asia and must be controlled to avoid disease severity and fatalities. However, little is known about the host-parasite interactions and molecular mechanisms in play during the course of P. knowlesi malaria infections, which also may be relevant across Plasmodium species. Here we contrast P. knowlesi sporozoite-initiated infections in Macaca mulatta and Macaca fascicularis using whole blood RNA-sequencing and transcriptomic analysis. These macaque hosts are evolutionarily close, yet malaria-naïve M. mulatta will succumb to blood-stage infection without treatment, whereas malaria-naïve M. fascicularis controls parasitemia without treatment. This comparative analysis reveals transcriptomic differences as early as the liver phase of infection, in the form of signaling pathways that are activated in M. fascicularis, but not M. mulatta. Additionally, while most immune responses are initially similar during the acute stage of the blood infection, significant differences arise subsequently. The observed differences point to prolonged inflammation and anti-inflammatory effects of IL10 in M. mulatta, while M. fascicularis undergoes a transcriptional makeover towards cell proliferation, consistent with its recovery. Together, these findings suggest that timely detection of P. knowlesi in M. fascicularis, coupled with control of inflammation while initiating the replenishment of key cell populations, helps contain the infection. Overall, this study points to specific genes and pathways that could be investigated as a basis for new drug targets that support recovery from acute malaria.

Metabolic Dynamics in Escherichia coli-Based Cell-Free Systems

The field of metabolic engineering has yielded remarkable accomplishments in using cells to produce valuable molecules, and cell-free expression (CFE) systems have the potential to push the field even further. However, CFE systems still face some outstanding challenges, including endogenous metabolic activity that is poorly understood yet has a significant impact on CFE productivity. Here, we use metabolomics to characterize the temporal metabolic changes in CFE systems and their constituent components, including significant metabolic activity in central carbon and amino acid metabolism. We find that while changing the reaction starting state via lysate preincubation impacts protein production, it has a comparatively small impact on metabolic state. We also demonstrate that changes to lysate preparation have a larger effect on protein yield and temporal metabolic profiles, though general metabolic trends are conserved. Finally, while we improve protein production through targeted supplementation of metabolic enzymes, we show that the endogenous metabolic activity is fairly resilient to these enzymatic perturbations. Overall, this work highlights the robust nature of CFE reaction metabolism as well as the importance of understanding the complex interdependence of metabolites and proteins in CFE systems to guide optimization efforts.

Effective Use of Linear DNA in Cell-Free Expression Systems

Cell-free expression systems (CFEs) are cutting-edge research tools used in the investigation of biological phenomena and the engineering of novel biotechnologies. While CFEs have many benefits over in vivo protein synthesis, one particularly significant advantage is that CFEs allow for gene expression from both plasmid DNA and linear expression templates (LETs). This is an important and impactful advantage because functional LETs can be efficiently synthesized in vitro in a few hours without transformation and cloning, thus expediting genetic circuit prototyping and allowing expression of toxic genes that would be difficult to clone through standard approaches. However, native nucleases present in the crude bacterial lysate (the basis for the most affordable form of CFEs) quickly degrade LETs and limit expression yield. Motivated by the significant benefits of using LETs in lieu of plasmid templates, numerous methods to enhance their stability in lysate-based CFEs have been developed. This review describes approaches to LET stabilization used in CFEs, summarizes the advancements that have come from using LETs with these methods, and identifies future applications and development goals that are likely to be impactful to the field. Collectively, continued improvement of LET-based expression and other linear DNA tools in CFEs will help drive scientific discovery and enable a wide range of applications, from diagnostics to synthetic biology research tools.

SCOUR: a stepwise machine learning framework for predicting metabolite-dependent regulatory interactions

BACKGROUND

The topology of metabolic networks is both well-studied and remarkably well-conserved across many species. The regulation of these networks, however, is much more poorly characterized, though it is known to be divergent across organisms-two characteristics that make it difficult to model metabolic networks accurately. While many computational methods have been built to unravel transcriptional regulation, there have been few approaches developed for systems-scale analysis and study of metabolic regulation. Here, we present a stepwise machine learning framework that applies established algorithms to identify regulatory interactions in metabolic systems based on metabolic data: stepwise classification of unknown regulation, or SCOUR.

RESULTS

We evaluated our framework on both noiseless and noisy data, using several models of varying sizes and topologies to show that our approach is generalizable. We found that, when testing on data under the most realistic conditions (low sampling frequency and high noise), SCOUR could identify reaction fluxes controlled only by the concentration of a single metabolite (its primary substrate) with high accuracy. The positive predictive value (PPV) for identifying reactions controlled by the concentration of two metabolites ranged from 32 to 88% for noiseless data, 9.2 to 49% for either low sampling frequency/low noise or high sampling frequency/high noise data, and 6.6-27% for low sampling frequency/high noise data, with results typically sufficiently high for lab validation to be a practical endeavor. While the PPVs for reactions controlled by three metabolites were lower, they were still in most cases significantly better than random classification.

CONCLUSIONS

SCOUR uses a novel approach to synthetically generate the training data needed to identify regulators of reaction fluxes in a given metabolic system, enabling metabolomics and fluxomics data to be leveraged for regulatory structure inference. By identifying and triaging the most likely candidate regulatory interactions, SCOUR can drastically reduce the amount of time needed to identify and experimentally validate metabolic regulatory interactions. As high-throughput experimental methods for testing these interactions are further developed, SCOUR will provide critical impact in the development of predictive metabolic models in new organisms and pathways.

LK-DFBA: a linear programming-based modeling strategy for capturing dynamics and metabolite-dependent regulation in metabolism

BACKGROUND

The systems-scale analysis of cellular metabolites, "metabolomics," provides data ideal for applications in metabolic engineering. However, many of the computational tools for strain design are built around Flux Balance Analysis (FBA), which makes assumptions that preclude direct integration of metabolomics data into the underlying models. Finding a way to retain the advantages of FBA's linear structure while relaxing some of its assumptions could allow us to account for metabolite levels and metabolite-dependent regulation in strain design tools built from FBA, improving the accuracy of predictions made by these tools. We designed, implemented, and characterized a modeling strategy based on Dynamic FBA (DFBA), called Linear Kinetics-Dynamic Flux Balance Analysis (LK-DFBA), to satisfy these specifications. Our strategy adds constraints describing the dynamics and regulation of metabolism that are strictly linear. We evaluated LK-DFBA against alternative modeling frameworks using simulated noisy data from a small in silico model and a larger model of central carbon metabolism in E. coli, and compared each framework's ability to recapitulate the original system.

RESULTS

In the smaller model, we found that we could use regression from a dynamic flux estimation (DFE) with an optional non-linear parameter optimization to reproduce metabolite concentration dynamic trends more effectively than an ordinary differential equation model with generalized mass action rate laws when tested under realistic data sampling frequency and noise levels. We observed detrimental effects across all tested modeling approaches when metabolite time course data were missing, but found these effects to be smaller for LK-DFBA in most cases. With the E. coli model, we produced qualitatively reasonable results with similar properties to the smaller model and explored two different parameterization structures that yield trade-offs in computation time and accuracy.

CONCLUSIONS

LK-DFBA allows for calculation of metabolite concentrations and considers metabolite-dependent regulation while still retaining many computational advantages of FBA. This provides the proof-of-principle for a new metabolic modeling framework with the potential to create genome-scale dynamic models and the potential to be applied in strain engineering tools that currently use FBA.

Active Analyte Import Improves the Dynamic Range and Sensitivity of a Vitamin B12 Biosensor

Cell-free systems provide a versatile platform for the development of low-cost, easy-to-use sensors for diverse analytes. However, sensor affinity dictates response sensitivity, and improving binding affinity can be challenging. Here, we describe efforts to address this problem while developing a biosensor for vitamin B₁₂, a critical micronutrient. We first use a B₁₂-responsive transcription factor to enable B₁₂-dependent output in a cell-free reaction, but the resulting sensor responds to B₁₂ far above clinically relevant concentrations. Surprisingly, when expressed in cells, the same sensor mediates a much more sensitive response to B₁₂. The sensitivity difference is partly due to regulated import that accumulates cytoplasmic B₁₂. Overexpression of importers further improves sensitivity, demonstrating an inherent advantage of whole-cell sensors. The resulting cells can respond to B₁₂ in serum, can be lyophilized, and are functional in a minimal-equipment environment, showing the potential utility of whole-cell sensors as sensitive, field-deployable diagnostics.

The Cancer Microbiome: Distinguishing Direct and Indirect Effects Requires a Systemic View

The collection of microbes that live in and on the human body - the human microbiome - can impact on cancer initiation, progression, and response to therapy, including cancer immunotherapy. The mechanisms by which microbiomes impact on cancers can yield new diagnostics and treatments, but much remains unknown. The interactions between microbes, diet, host factors, drugs, and cell-cell interactions within the cancer itself likely involve intricate feedbacks, and no single component can explain all the behavior of the system. Understanding the role of host-associated microbial communities in cancer systems will require a multidisciplinary approach combining microbial ecology, immunology, cancer cell biology, and computational biology - a systems biology approach.

Metabolic Profiling of Escherichia coli-based Cell-Free Expression Systems for Process Optimization

Biotechnology has transformed the production of various chemicals and pharmaceuticals due to its efficient and selective processes, but it is inherently limited by its use of live cells as "biocatalysts." Cell-free expression (CFE) systems, which use a protein lysate isolated from whole cells, have the potential to overcome these challenges and broaden the scope of biomanufacturing. Implementation of CFE systems at scale will require determining clear markers of lysate activity and developing supplementation approaches that compensate for potential variability across batches and experimental protocols. Towards this goal, we use metabolomics to relate lysate preparation and performance to metabolic activity. We show that lysate processing affects the metabolite makeup of lysates, and that lysate metabolite levels change over the course of a CFE reaction regardless of whether a target compound is produced. Finally, we use this information to develop ways to standardize lysate activity and to design an improved CFE system.

Dynamic and tunable metabolite control for robust minimal-equipment assessment of serum zinc

Bacterial biosensors can enable programmable, selective chemical production, but difficulties incorporating metabolic pathways into complex sensor circuits have limited their development and applications. Here we overcome these challenges and present the development of fast-responding, tunable sensor cells that produce different pigmented metabolites based on extracellular concentrations of zinc (a critical micronutrient). We create a library of dual-input synthetic promoters that decouple cell growth from zinc-specific metabolite production, enabling visible cell coloration within 4 h. Using additional transcriptional and metabolic control methods, we shift the response thresholds by an order of magnitude to measure clinically relevant zinc concentrations. The resulting sensor cells report zinc concentrations in individual donor serum samples; we demonstrate that they can provide results in a minimal-equipment fashion, serving as the basis for a field-deployable assay for zinc deficiency. The presented advances are likely generalizable to the creation of other types of sensors and diagnostics.

Tightly controlling cell output is challenging, which has limited development and applications of bacterial sensors. Here the authors develop tunable, fast-responding sensors to control production of metabolic pigments and use them to assess zinc deficiency in a low-cost, minimal equipment fashion.

Point-of-care biomarker quantification enabled by sample-specific calibration

Easy-to-perform, relatively inexpensive blood diagnostics have transformed at-home healthcare for some patients, but they require analytical equipment and are not easily adapted to measuring other biomarkers. The requirement for reliable quantification in complex sample types (such as blood) has been a critical roadblock in developing and deploying inexpensive, minimal-equipment diagnostics. Here, we developed a platform for inexpensive, easy-to-use diagnostics that uses cell-free expression to generate colored readouts that are visible to the naked eye, yet quantitative and robust to the interference effects seen in complex samples. We achieved this via a parallelized calibration scheme that uses the patient sample to generate custom reference curves. We used this approach to quantify a clinically relevant micronutrient and to quantify nucleic acids, demonstrating a generalizable platform for low-cost quantitative diagnostics.

The Metabolomics Society-Current State of the Membership and Future Directions

Background: In 2017, the Metabolomics Society conducted a survey among its members to assess the degree of its current success, define opportunities for improving its service to the community and make plans to establish future goals and direction of the Society. Methods: A 32-question online survey was sent via e-mail to all Metabolomics Society members as of 19 June 2017 (n = 644). In addition to the direct e-mails, the link to access the survey was made available through social media. The survey was open until 10 August 2017. Question-specific data were reported using the summary data generated by SurveyMonkey and additional stratified analyses performed using Stata 15. Results: The number of respondents was 394 (61%) with 348 (88%) completing the multiple-choice questions in survey. Metabolomics Society annual meetings, networking and the opportunity to join the global metabolomics community were among the most important benefits expressed by the Metabolomics Society members. Conclusions: The survey collected the first data focusing on membership issues from Society members. The Society should focus on collecting and monitoring of demographic data during the membership registration process; continuing to support the early-career members of the Society; and developing initiatives that focus on member networking to retain and increase Society membership.

Gas Chromatography-Mass Spectrometry Microbial Metabolomics for Applications in Strain Optimization

Metabolomics is the systems-scale measurement of biochemical intermediates in biological systems; by virtue of its deep and broad study of metabolism, it has great potential for applications in metabolic engineering. While a number of the analytical techniques used widely in metabolomics are familiar to metabolic engineers performing post hoc analyses of product titers, the requirements for accurately capturing metabolism at a systems scale rather than just measuring a single secreted product are much more complicated. Nonetheless, metabolomics (which is still not widely available as an affordable consumer service like many molecular biology services) is within reach of many properly equipped metabolic engineering groups. To this end, we present a detailed metabolomics protocol with application to strain optimization. Specifically, we focus on characterizing metabolism in the yeast Saccharomyces cerevisiae using gas chromatography coupled to mass spectrometry. The measurement of metabolic intermediates that results from such approaches has the potential to enable more informed and rational efforts towards pathway engineering and strain optimization.

NS-kNN: a modified k-nearest neighbors approach for imputing metabolomics data

INTRODUCTION

A common problem in metabolomics data analysis is the existence of a substantial number of missing values, which can complicate, bias, or even prevent certain downstream analyses. One of the most widely-used solutions to this problem is imputation of missing values using a k-nearest neighbors (kNN) algorithm to estimate missing metabolite abundances. kNN implicitly assumes that missing values are uniformly distributed at random in the dataset, but this is typically not true in metabolomics, where many values are missing because they are below the limit of detection of the analytical instrumentation.

OBJECTIVES

Here, we explore the impact of nonuniformly distributed missing values (missing not at random, or MNAR) on imputation performance. We present a new model for generating synthetic missing data and a new algorithm, No-Skip kNN (NS-kNN), that accounts for MNAR values to provide more accurate imputations.

METHODS

We compare the imputation errors of the original kNN algorithm using two distance metrics, NS-kNN, and a recently developed algorithm KNN-TN, when applied to multiple experimental datasets with different types and levels of missing data.

RESULTS

Our results show that NS-kNN typically outperforms kNN when at least 20-30% of missing values in a dataset are MNAR. NS-kNN also has lower imputation errors than KNN-TN on realistic datasets when at least 50% of missing values are MNAR.

CONCLUSION

Accounting for the nonuniform distribution of missing values in metabolomics data can significantly improve the results of imputation algorithms. The NS-kNN method imputes missing metabolomics data more accurately than existing kNN-based approaches when used on realistic datasets.

Systems Biology-Based Investigation of Host-Plasmodium Interactions

Malaria is a serious, complex disease caused by parasites of the genus Plasmodium. Plasmodium parasites affect multiple tissues as they evade immune responses, replicate, sexually reproduce, and transmit between vertebrate and invertebrate hosts. The explosion of omics technologies has enabled large-scale collection of Plasmodium infection data, revealing systems-scale patterns, mechanisms of pathogenesis, and the ways that host and pathogen affect each other. Here, we provide an overview of recent efforts using systems biology approaches to study host-Plasmodium interactions and the biological themes that have emerged from these efforts. We discuss some of the challenges in using systems biology for this goal, key research efforts needed to address those issues, and promising future malaria applications of systems biology.

Metabolomics Analysis of the Toxic Effects of the Production of Lycopene and Its Precursors

Using cells as microbial factories enables highly specific production of chemicals with many advantages over chemical syntheses. A number of exciting new applications of this approach are in the area of precision metabolic engineering, which focuses on improving the specificity of target production. In recent work, we have used precision metabolic engineering to design lycopene-producing Escherichia coli for use as a low-cost diagnostic biosensor. To increase precursor availability and thus the rate of lycopene production, we heterologously expressed the mevalonate pathway. We found that simultaneous induction of these pathways increases lycopene production, but induction of the mevalonate pathway before induction of the lycopene pathway decreases both lycopene production and growth rate. Here, we aim to characterize the metabolic changes the cells may be undergoing during expression of either or both of these heterologous pathways. After establishing an improved method for quenching E. coli for metabolomics analysis, we used two-dimensional gas chromatography coupled to mass spectrometry (GCxGC-MS) to characterize the metabolomic profile of our lycopene-producing strains in growth conditions characteristic of our biosensor application. We found that the metabolic impacts of producing low, non-toxic levels of lycopene are of much smaller magnitude than the typical metabolic changes inherent to batch growth. We then used metabolomics to study differences in metabolism caused by the time of mevalonate pathway induction and the presence of the lycopene biosynthesis genes. We found that overnight induction of the mevalonate pathway was toxic to cells, but that the cells could recover if the lycopene pathway was not also heterologously expressed. The two pathways appeared to have an antagonistic metabolic effect that was clearly reflected in the cells’ metabolic profiles. The metabolites homocysteine and homoserine exhibited particularly interesting behaviors and may be linked to the growth inhibition seen when the mevalonate pathway is induced overnight, suggesting potential future work that may be useful in engineering increased lycopene biosynthesis.

Development of a Pigment-Based Whole-Cell Zinc Biosensor for Human Serum

Deficiencies in vitamins and minerals (micronutrients) are a critical global health concern, in part due to logistical difficulties in assessing population micronutrient status. Whole-cell biosensors offer a unique opportunity to address this issue, with the potential to move sample analysis from centralized, resource-intensive clinics to minimal-resource, on-site measurement. Here, we present a proof-of-concept whole-cell biosensor in Escherichia coli for detecting zinc, a micronutrient for which deficiencies are a significant public health burden. Importantly, the whole-cell biosensor produces readouts (pigments) that are visible to the naked eye, mitigating the need for measurement equipment and thus increasing feasibility for sensor field-friendliness and affordability at a global scale. Two zinc-responsive promoter/transcription factor systems are used to differentially control production of three distinctly colored pigments in response to zinc levels in culture. We demonstrate strategies for tuning each zinc-responsive system to turn production of the different pigments on and off at different zinc levels, and we demonstrate production of three distinct color regimes over a concentration range relevant to human health. We also demonstrate the ability of the sensor cells to grow and produce pigment when cultured in human serum, the ultimate target matrix for assessing zinc nutritional status. Specifically, we present approaches to overcome innate immune responses that would otherwise hinder bacterial sensor survival, and we demonstrate production of multiple pigment regimes in human serum with different zinc levels. This work provides proof of principle for the development of low-cost, minimal-equipment, field-deployable biosensors for nutritional epidemiology applications.

Metabolic modeling helps interpret transcriptomic changes during malaria

Disease represents a specific case of malfunctioning within a complex system. Whereas it is often feasible to observe and possibly treat the symptoms of a disease, it is much more challenging to identify and characterize its molecular root causes. Even in infectious diseases that are caused by a known parasite, it is often impossible to pinpoint exactly which molecular profiles of components or processes are directly or indirectly altered. However, a deep understanding of such profiles is a prerequisite for rational, efficacious treatments. Modern omics methodologies are permitting large-scale scans of some molecular profiles, but these scans often yield results that are not intuitive and difficult to interpret. For instance, the comparison of healthy and diseased transcriptome profiles may point to certain sets of involved genes, but a host of post-transcriptional processes and regulatory mechanisms renders predictions regarding metabolic or physiological consequences of the observed changes in gene expression unreliable. Here we present proof of concept that dynamic models of metabolic pathway systems may offer a tool for interpreting transcriptomic profiles measured during disease. We illustrate this strategy with the interpretation of expression data of genes coding for enzymes associated with purine metabolism. These data were obtained during infections of rhesus macaques (Macaca mulatta) with the malaria parasite Plasmodium cynomolgi or P. coatneyi. The model-based interpretation reveals clear patterns of flux redistribution within the purine pathway that are consistent between the two malaria pathogens and are even reflected in data from humans infected with P. falciparum. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.

Small molecule signaling, regulation, and potential applications in cellular therapeutics

Small molecules have many important roles across the tree of life: they regulate processes from metabolism to transcription, they enable signaling within and between species, and they serve as the biochemical building blocks for cells. They also represent valuable phenotypic endpoints that are promising for use as biomarkers of disease states. In the context of engineering cell-based therapeutics, they hold particularly great promise for enabling finer control over the therapeutic cells and allowing them to be responsive to extracellular cues. The natural signaling and regulatory functions of small molecules can be harnessed and rewired to control cell activity and delivery of therapeutic payloads, potentially increasing efficacy while decreasing toxicity. To that end, this review considers small molecule-mediated regulation and signaling in bacteria. We first discuss some of the most prominent applications and aspirations for responsive cell-based therapeutics. We then describe the transport, signaling, and regulation associated with three classes of molecules that may be exploited in the engineering of therapeutic bacteria: amino acids, fatty acids, and quorum-sensing signaling molecules. We also present examples of existing engineering efforts to generate cells that sense and respond to levels of different small molecules. Finally, we discuss future directions for how small molecule-mediated regulation could be harnessed for therapeutic applications, as well as some critical considerations for the ultimate success of such endeavors. WIREs Syst Biol Med 2018, 10:e1405. doi: 10.1002/wsbm.1405 This article is categorized under: Biological Mechanisms > Cell Signaling Biological Mechanisms > Metabolism Translational, Genomic, and Systems Medicine > Therapeutic Methods.

Integrative analysis associates monocytes with insufficient erythropoiesis during acute Plasmodium cynomolgi malaria in rhesus macaques

Background

Mild to severe anaemia is a common complication of malaria that is caused in part by insufficient erythropoiesis in the bone marrow. This study used systems biology to evaluate the transcriptional and alterations in cell populations in the bone marrow during Plasmodium cynomolgi infection of rhesus macaques (a model of Plasmodium vivax malaria) that may affect erythropoiesis.

Results

An appropriate erythropoietic response did not occur to compensate for anaemia during acute cynomolgi malaria despite an increase in erythropoietin levels. During this period, there were significant perturbations in the bone marrow transcriptome. In contrast, relapses did not induce anaemia and minimal changes in the bone marrow transcriptome were detected. The differentially expressed genes during acute infection were primarily related to ongoing inflammatory responses with significant contributions from Type I and Type II Interferon transcriptional signatures. These were associated with increased frequency of intermediate and non-classical monocytes. Recruitment and/or expansion of these populations was correlated with a decrease in the erythroid progenitor population during acute infection, suggesting that monocyte-associated inflammation may have contributed to anaemia. The decrease in erythroid progenitors was associated with downregulation of genes regulated by GATA1 and GATA2, two master regulators of erythropoiesis, providing a potential molecular basis for these findings.

Conclusions

These data suggest the possibility that malarial anaemia may be driven by monocyte-associated disruption of GATA1/GATA2 function in erythroid progenitors resulting in insufficient erythropoiesis during acute infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-2029-z) contains supplementary material, which is available to authorized users.

Precise control of lycopene production to enable a fast-responding, minimal-equipment biosensor

Pigmented metabolites have great potential for use in biosensors that target low-resource areas, since sensor output can be interpreted without any equipment. However, full repression of pigment production when undesired is challenging, as even small amounts of enzyme can catalyze the production of large, visible amounts of pigment. The red pigment lycopene could be particularly useful because of its position in the multi-pigment carotenoid pathway, but commonly used inducible promoter systems cannot repress lycopene production. In this paper, we designed a system that could fully repress lycopene production in the absence of an inducer and produce visible lycopene within two hours of induction. We engineered Lac, Ara, and T7 systems to be up to 10 times more repressible, but these improved systems could still not fully repress lycopene. Translational modifications proved much more effective in controlling lycopene. By decreasing the strength of the ribosomal binding sites on the crtEBI genes, we enabled full repression of lycopene and production of visible lycopene in 3-4h of induction. Finally, we added the mevalonate pathway enzymes to increase the rate of lycopene production upon induction and demonstrated that supplementation of metabolic precursors could decrease the time to coloration to about 1.5h. In total, this represents over an order of magnitude reduction in response time compared to the previously reported strategy. The approaches used here demonstrate the disconnect between fluorescent and metabolite reporters, help enable the use of lycopene as a reporter, and are likely generalizable to other systems that require precise control of metabolite production.

Metabolomics identifies the intersection of phosphoethanolamine with menaquinone-triggered apoptosis in an in vitro model of leukemia

Altered metabolism is increasingly acknowledged as an important aspect of cancer, and thus serves as a potentially fertile area for the identification of therapeutic targets or leads. Our recent work using transcriptional data to predict metabolite levels in cancer cells led to preliminary evidence of the antiproliferative role of menaquinone (vitamin K2) in the Jurkat cell line model of acute lymphoblastic leukemia. However, nothing is known about the direct metabolic impacts of menaquinone in cancer, which could provide insights into its mechanism of action. Here, we used metabolomics to investigate the process by which menaquinone exerts antiproliferative activity on Jurkat cells. We first validated the dose-dependent, semi-selective, pro-apoptotic activity of menaquinone treatment on Jurkat cells relative to non-cancerous lymphoblasts. We then used mass spectrometry-based metabolomics to identify systems-scale changes in metabolic dynamics that are distinct from changes induced in non-cancerous cells or by other chemotherapeutics. One of the most significantly affected metabolites was phosphoethanolamine, which exhibited a two-fold increase in menaquinone-treated Jurkat cells compared to vehicle-treated cells at 24 h, growing to a five-fold increase at 72 h. Phosphoethanolamine elevation was observed prior to the induction of apoptosis, and was not observed in menaquinone-treated lymphoblasts or chemotherapeutic-treated Jurkat cells. We also validated the link between menaquinone and phosphoethanolamine in an ovarian cancer cell line, suggesting potentially broad applicability of their relationship. This metabolomics-based work is the first detailed characterization of the metabolic impacts of menaquinone treatment and the first identified link between phosphoethanolamine and menaquinone-induced apoptosis.

A Multi-Modality CMOS Sensor Array for Cell-Based Assay and Drug Screening

In this paper, we present a fully integrated multi-modality CMOS cellular sensor array with four sensing modalities to characterize different cell physiological responses, including extracellular voltage recording, cellular impedance mapping, optical detection with shadow imaging and bioluminescence sensing, and thermal monitoring. The sensor array consists of nine parallel pixel groups and nine corresponding signal conditioning blocks. Each pixel group comprises one temperature sensor and 16 tri-modality sensor pixels, while each tri-modality sensor pixel can be independently configured for extracellular voltage recording, cellular impedance measurement (voltage excitation/current sensing), and optical detection. This sensor array supports multi-modality cellular sensing at the pixel level, which enables holistic cell characterization and joint-modality physiological monitoring on the same cellular sample with a pixel resolution of 80 μm × 100 μm. Comprehensive biological experiments with different living cell samples demonstrate the functionality and benefit of the proposed multi-modality sensing in cell-based assay and drug screening.

A tree-like Bayesian structure learning algorithm for small-sample datasets from complex biological model systems

Background

There are increasing efforts to bring high-throughput systems biology techniques to bear on complex animal model systems, often with a goal of learning about underlying regulatory network structures (e.g., gene regulatory networks). However, complex animal model systems typically have significant limitations on cohort sizes, number of samples, and the ability to perform follow-up and validation experiments. These constraints are particularly problematic for many current network learning approaches, which require large numbers of samples and may predict many more regulatory relationships than actually exist.

Results

Here, we test the idea that by leveraging the accuracy and efficiency of classifiers, we can construct high-quality networks that capture important interactions between variables in datasets with few samples. We start from a previously-developed tree-like Bayesian classifier and generalize its network learning approach to allow for arbitrary depth and complexity of tree-like networks. Using four diverse sample networks, we demonstrate that this approach performs consistently better at low sample sizes than the Sparse Candidate Algorithm, a representative approach for comparison because it is known to generate Bayesian networks with high positive predictive value. We develop and demonstrate a resampling-based approach to enable the identification of a viable root for the learned tree-like network, important for cases where the root of a network is not known a priori. We also develop and demonstrate an integrated resampling-based approach to the reduction of variable space for the learning of the network. Finally, we demonstrate the utility of this approach via the analysis of a transcriptional dataset of a malaria challenge in a non-human primate model system, Macaca mulatta, suggesting the potential to capture indicators of the earliest stages of cellular differentiation during leukopoiesis.

Conclusions

We demonstrate that by starting from effective and efficient approaches for creating classifiers, we can identify interesting tree-like network structures with significant ability to capture the relationships in the training data. This approach represents a promising strategy for inferring networks with high positive predictive value under the constraint of small numbers of samples, meeting a need that will only continue to grow as more high-throughput studies are applied to complex model systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0194-7) contains supplementary material, which is available to authorized users.

Precision metabolic engineering: The design of responsive, selective, and controllable metabolic systems

Metabolic engineering is generally focused on static optimization of cells to maximize production of a desired product, though recently dynamic metabolic engineering has explored how metabolic programs can be varied over time to improve titer. However, these are not the only types of applications where metabolic engineering could make a significant impact. Here, we discuss a new conceptual framework, termed "precision metabolic engineering," involving the design and engineering of systems that make different products in response to different signals. Rather than focusing on maximizing titer, these types of applications typically have three hallmarks: sensing signals that determine the desired metabolic target, completely directing metabolic flux in response to those signals, and producing sharp responses at specific signal thresholds. In this review, we will first discuss and provide examples of precision metabolic engineering. We will then discuss each of these hallmarks and identify which existing metabolic engineering methods can be applied to accomplish those tasks, as well as some of their shortcomings. Ultimately, precise control of metabolic systems has the potential to enable a host of new metabolic engineering and synthetic biology applications for any problem where flexibility of response to an external signal could be useful.

Improved metabolite profile smoothing for flux estimation

As genome-scale metabolic models become more sophisticated and dynamic, one significant challenge in using these models is to effectively integrate increasingly prevalent systems-scale metabolite profiling data into them. One common data processing step when integrating metabolite data is to smooth experimental time course measurements: the smoothed profiles can be used to estimate metabolite accumulation (derivatives), and thus the flux distribution of the metabolic model. However, this smoothing step is susceptible to the (often significant) noise in experimental measurements, limiting the accuracy of downstream model predictions. Here, we present several improvements to current approaches for smoothing metabolite time course data using defined functions. First, we use a biologically-inspired mathematical model function taken from transcriptional profiling and clustering literature that captures the dynamics of many biologically relevant transient processes. We demonstrate that it is competitive with, and often superior to, previously described fitting schemas, and may serve as an effective single option for data smoothing in metabolic flux applications. We also implement a resampling-based approach to buffer out sensitivity to specific data sets and allow for more accurate fitting of noisy data. We found that this method, as well as the addition of parameter space constraints, yielded improved estimates of concentrations and derivatives (fluxes) in previously described fitting functions. These methods have the potential to improve the accuracy of existing and future dynamic metabolic models by allowing for the more effective integration of metabolite profiling data.

Precise metabolic engineering of carotenoid biosynthesis in Escherichia coli towards a low-cost biosensor

Micronutrient deficiencies, including zinc deficiency, are responsible for hundreds of thousands of deaths annually. A key obstacle to allocating scarce treatment resources is the ability to measure population blood micronutrient status inexpensively and quickly enough to identify those who most need treatment. This paper develops a metabolically engineered strain of Escherichia coli to produce different colored pigments (violacein, lycopene, and β-carotene) in response to different extracellular zinc levels, for eventual use in an inexpensive blood zinc diagnostic test. However, obtaining discrete color states in the carotenoid pathway required precise engineering of metabolism to prevent reaction at low zinc concentrations but allow complete reaction at higher concentrations, and all under the constraints of natural regulator limitations. Hence, the metabolic engineering challenge was not to improve titer, but to enable precise control of pathway state. A combination of gene dosage, post-transcriptional, and post-translational regulation was necessary to allow visible color change over physiologically relevant ranges representing a small fraction of the regulator's dynamic response range, with further tuning possible by modulation of precursor availability. As metabolic engineering expands its applications and develops more complex systems, tight control of system components will likely become increasingly necessary, and the approach presented here can be generalized to other natural sensing systems for precise control of pathway state.

From genome-scale data to models of infectious disease: A Bayesian network-based strategy to drive model development

High-throughput, genome-scale data present a unique opportunity to link host to pathogen on a molecular level. Forging such connections will help drive the development of mathematical models to better understand and predict both pathogen behavior and the epidemiology of infectious diseases, including malaria. However, the datasets that can aid in identifying these links and models are vast and not amenable to simple, reductionist, and univariate analyses. These datasets require data mining in order to identify the truly important measurements that best describe clinical and molecular observations. Moreover, these datasets typically have relatively few samples due to experimental limitations (particularly for human studies or in vivo animal experiments), making data mining extremely difficult. Here, after first providing a brief overview of common strategies for data reduction and identification of relationships between variables for inclusion in mathematical models, we present a new generalized strategy for performing these data reduction and relationship inference tasks. Our approach emphasizes the importance of robustness when using data to drive model development, particularly when using genome-scale, small-sample in vivo data. We identify the use of appropriate feature reduction combined with data permutations and subsampling strategies as being critical to enable increasingly robust results from network inference using high-dimensional, low-observation data.

Differential metabolite levels in response to spawning-induced inappetence in Atlantic salmon Salmo salar

Atlantic salmon Salmo salar undergo months-long inappetence during spawning, but it is not known whether this inappetence is a pathological state or one for which the fish are adapted. Recent work has shown that inappetent whale sharks can exhibit circulating metabolite profiles similar to ketosis known to occur in humans during starvation. In this work, metabolite profiling was used to explore differences in analyte profiles between a cohort of inappetent spawning run Atlantic salmon and captively reared animals that were fed up to and through the time of sampling. The two classes of animals were easily distinguished by their metabolite profiles. The sea-run fish had elevated ɷ-9 fatty acids relative to the domestic feeding animals, while other fatty acid concentrations were reduced. Sugar alcohols were generally elevated in inappetent animals, suggesting potentially novel metabolic responses or pathways in fish that feature these compounds. Compounds expected to indicate a pathological catabolic state were not more abundant in the sea-run fish, suggesting that the animals, while inappetent, were not stressed in an unnatural way. These findings demonstrate the power of discovery-based metabolomics for exploring biochemistry in poorly understood animal models.

Distinct metabolic responses of an ovarian cancer stem cell line

BACKGROUND

Cancer metabolism is emerging as an important focus area in cancer research. However, the in vitro cell culture conditions under which much cellular metabolism research is performed differ drastically from in vivo tumor conditions, which are characterized by variations in the levels of oxygen, nutrients like glucose, and other molecules like chemotherapeutics. Moreover, it is important to know how the diverse cell types in a tumor, including cancer stem cells that are believed to be a major cause of cancer recurrence, respond to these variations. Here, in vitro environmental perturbations designed to mimic different aspects of the in vivo environment were used to characterize how an ovarian cancer cell line and its derived, isogenic cancer stem cells metabolically respond to environmental cues.

RESULTS

Mass spectrometry was used to profile metabolite levels in response to in vitro environmental perturbations. Docetaxel, the chemotherapeutic used for this experiment, caused significant metabolic changes in amino acid and carbohydrate metabolism in ovarian cancer cells, but had virtually no metabolic effect on isogenic ovarian cancer stem cells. Glucose deprivation, hypoxia, and the combination thereof altered ovarian cancer cell and cancer stem cell metabolism to varying extents for the two cell types. Hypoxia had a much larger effect on ovarian cancer cell metabolism, while glucose deprivation had a greater effect on ovarian cancer stem cell metabolism. Core metabolites and pathways affected by these perturbations were identified, along with pathways that were unique to cell types or perturbations.

CONCLUSIONS

The metabolic responses of an ovarian cancer cell line and its derived isogenic cancer stem cells differ greatly under most conditions, suggesting that these two cell types may behave quite differently in an in vivo tumor microenvironment. While cancer metabolism and cancer stem cells are each promising potential therapeutic targets, such varied behaviors in vivo would need to be considered in the design and early testing of such treatments.

Comparative transcriptomics and metabolomics in a rhesus macaque drug administration study

We describe a multi-omic approach to understanding the effects that the anti-malarial drug pyrimethamine has on immune physiology in rhesus macaques (Macaca mulatta). Whole blood and bone marrow (BM) RNA-Seq and plasma metabolome profiles (each with over 15,000 features) have been generated for five naïve individuals at up to seven timepoints before, during and after three rounds of drug administration. Linear modeling and Bayesian network analyses are both considered, alongside investigations of the impact of statistical modeling strategies on biological inference. Individual macaques were found to be a major source of variance for both omic data types, and factoring individuals into subsequent modeling increases power to detect temporal effects. A major component of the whole blood transcriptome follows the BM with a time-delay, while other components of variation are unique to each compartment. We demonstrate that pyrimethamine administration does impact both compartments throughout the experiment, but very limited perturbation of transcript or metabolite abundance was observed following each round of drug exposure. New insights into the mode of action of the drug are presented in the context of pyrimethamine's predicted effect on suppression of cell division and metabolism in the immune system.

Applications of metabolomics in cancer research

The first discovery of metabolic changes in cancer occurred almost a century ago. While the genetic underpinnings of cancer have dominated its study since then, altered metabolism has recently been acknowledged as a key hallmark of cancer and metabolism-focused research has received renewed attention. The emerging field of metabolomics – which attempts to profile all metabolites within a cell or biological system – is now being used to analyze cancer metabolism on a system-wide scale, painting a broad picture of the altered pathways and their interactions with each other. While a large fraction of cancer metabolomics research is focused on finding diagnostic biomarkers, metabolomics is also being used to obtain more fundamental mechanistic insight into cancer and carcinogenesis. Applications of metabolomics are also emerging in areas such as tumor staging and assessment of treatment efficacy. This review summarizes contributions that metabolomics has made in cancer research and presents the current challenges and potential future directions within the field.