PubMed Health⌕ Search

Biomedical subjects

Stephen S Fong

Publications and source records attributed to Stephen S Fong.

10 recordsLinked to original sources

Metabolic analysis of adaptive evolution for in silico-designed lactate-producing strains.

Experimental evolution is now frequently applied to many biological systems to achieve desired objectives. To obtain optimized performance for metabolite production, a successful strategy has been recently developed that couples metabolic engineering techniques with laboratory evolution of microorganisms. Previously, we reported the growth characteristics of three lactate-producing, adaptively evolved Escherichia coli mutant strains designed by the OptKnock computational algorithm. Here, we describe the use of (13)C-labeled experiments and mass distribution measurements to study the evolutionary effects on the fluxome of these differently designed strains. Metabolic flux ratios and intracellular flux distributions as well as physiological data were used to elucidate metabolic responses over the course of adaptive evolution and metabolic differences among strains. The study of 3 unevolved and 12 evolved engineered strains as well as a wild-type strain suggests that evolution resulted in remarkable improvements in both substrate utilization rate and the proportion of glycolytic flux to total glucose utilization flux. Among three strain designs, the most significant increases in the fraction of glucose catabolized through glycolysis (>50%) and the glycolytic fluxes (>twofold) were observed in phosphotransacetylase and phosphofructokinase 1 (PFK1) double deletion (pta- pfkA) strains, which were likely attributed to the dramatic evolutionary increase in gene expression and catalytic activity of the minor PFK encoded by pfkB. These fluxomic studies also revealed the important role of acetate synthetic pathway in anaerobic lactate production. Moreover, flux analysis suggested that independent of genetic background, optimal relative flux distributions in cells could be achieved faster than physiological parameters such as nutrient utilization rate.

Acetate-CoA Ligase↗

Systems approach to refining genome annotation.

Genome-scale models of Escherichia coli K-12 MG1655 metabolism have been able to predict growth phenotypes in most, but not all, defined growth environments. Here we introduce the use of an optimization-based algorithm that predicts the missing reactions that are required to reconcile computation and experiment when they disagree. The computer-generated hypotheses for missing reactions were verified experimentally in five cases, leading to the functional assignment of eight ORFs (yjjLMN, yeaTU, dctA, idnT, and putP) with two new enzymatic activities and four transport functions. This study thus demonstrates the use of systems analysis to discover metabolic and transport functions and their genetic basis by a combination of experimental and computational approaches.

Biological Transport↗

Identification of genome-scale metabolic network models using experimentally measured flux profiles.

Genome-scale metabolic network models can be reconstructed for well-characterized organisms using genomic annotation and literature information. However, there are many instances in which model predictions of metabolic fluxes are not entirely consistent with experimental data, indicating that the reactions in the model do not match the active reactions in the in vivo system. We introduce a method for determining the active reactions in a genome-scale metabolic network based on a limited number of experimentally measured fluxes. This method, called optimal metabolic network identification (OMNI), allows efficient identification of the set of reactions that results in the best agreement between in silico predicted and experimentally measured flux distributions. We applied the method to intracellular flux data for evolved Escherichia coli mutant strains with lower than predicted growth rates in order to identify reactions that act as flux bottlenecks in these strains. The expression of the genes corresponding to these bottleneck reactions was often found to be downregulated in the evolved strains relative to the wild-type strain. We also demonstrate the ability of the OMNI method to diagnose problems in E. coli strains engineered for metabolite overproduction that have not reached their predicted production potential. The OMNI method applied to flux data for evolved strains can be used to provide insights into mechanisms that limit the ability of microbial strains to evolve towards their predicted optimal growth phenotypes. When applied to industrial production strains, the OMNI method can also be used to suggest metabolic engineering strategies to improve byproduct secretion. In addition to these applications, the method should prove to be useful in general for reconstructing metabolic networks of ill-characterized microbial organisms based on limited amounts of experimental data.

Computational Biology↗

Latent pathway activation and increased pathway capacity enable Escherichia coli adaptation to loss of key metabolic enzymes.

The ability of biological systems to adapt to genetic and environmental perturbations is a fundamental but poorly understood process at the molecular level. By quantifying metabolic fluxes and global mRNA abundance, we investigated the genetic and metabolic mechanisms that underlie adaptive evolution of four metabolic gene deletion mutants of Escherichia coli (delta pgi, delta ppc, delta pta, and delta tpi) in parallel evolution experiments of each mutant. The initial response to the gene deletions was flux rerouting through local bypass reactions or normally latent pathways. The principal effect of evolution was improved capacity of already active pathways, whereas new flux distributions were not observed. Combinatorial changes in capacity and pathway activation, however, led to different intracellular flux states that enabled evolution in three of the four parallel cases tested. The molecular bases of the evolved phenotypes were then elucidated by global mRNA transcript analyses. Activation of latent pathways and flux changes in the tricarboxylic acid cycle were found to correlate well with molecular changes at the transcriptional level. Flux alterations in other central metabolic pathways, in contrast, were apparently not connected to changes in the transcriptional network. These results give new insight into the dynamics of the evolutionary process by demonstrating the flexibility of the metabolic network of E. coli to compensate for genetic perturbations and the utility of combining multiple high throughput data sets to differentiate between causal and noncausal mechanistic changes.

Animals↗

In silico design and adaptive evolution of Escherichia coli for production of lactic acid.

The development and validation of new methods to help direct rational strain design for metabolite overproduction remains an important problem in metabolic engineering. Here we show that computationally predicted E. coli strain designs, calculated from a genome-scale metabolic model, can lead to successful production strains and that adaptive evolution of the engineered strains can lead to improved production capabilities. Three strain designs for lactate production were implemented yielding a total of 11 evolved production strains that were used to demonstrate the utility of this integrated approach. Strains grown on 2 g/L glucose at 37 degrees C showed lactate titers ranging from 0.87 to 1.75 g/L and secretion rates that were directly coupled to growth rates.

Adaptation, Physiological↗

Parallel adaptive evolution cultures of Escherichia coli lead to convergent growth phenotypes with different gene expression states.

Laboratory evolution can be used to address fundamental questions about adaptation to selection pressures and, ultimately, the process of evolution. In this study, we investigated the reproducibility of growth phenotypes and global gene expression states during adaptive evolution. The results from parallel, replicate adaptive evolution experiments of Escherichia coli K-12 MG1655 grown on either lactate or glycerol minimal media showed that (1) growth phenotypes at the endpoint of evolution are convergent and reproducible; (2) endpoints of evolution have different underlying gene expression states; and (3) the evolutionary gene expression response involves a large number of compensatory expression changes and a smaller number of adaptively beneficial expression changes common across evolution strains. Gene expression changes initially showed a large number of differentially expressed genes in response to an environmental change followed by a return of most genes to a baseline expression level, leaving a relatively small set of differentially expressed genes at the endpoint that varied between evolved populations.

Biological Evolution↗

Metabolic gene-deletion strains of Escherichia coli evolve to computationally predicted growth phenotypes.

Genome-scale metabolic models have a promising ability to describe cellular phenotypes accurately. Here we show that strains of Escherichia coli carrying a deletion of a single metabolic gene increase their growth rates (by 87% on average) during adaptive evolution and that the endpoint growth rates can be predicted computationally in 39 of 50 (78%) strains tested. These results show that computational models can be used to predict the eventual effects of genetic modifications.

Directed Molecular Evolution↗

Description and interpretation of adaptive evolution of Escherichia coli K-12 MG1655 by using a genome-scale in silico metabolic model.

Genome-scale in silico metabolic networks of Escherichia coli have been reconstructed. By using a constraint-based in silico model of a reconstructed network, the range of phenotypes exhibited by E. coli under different growth conditions can be computed, and optimal growth phenotypes can be predicted. We hypothesized that the end point of adaptive evolution of E. coli could be accurately described a priori by our in silico model since adaptive evolution should lead to an optimal phenotype. Adaptive evolution of E. coli during prolonged exponential growth was performed with M9 minimal medium supplemented with 2 g of alpha-ketoglutarate per liter, 2 g of lactate per liter, or 2 g of pyruvate per liter at both 30 and 37 degrees C, which produced seven distinct strains. The growth rates, substrate uptake rates, oxygen uptake rates, by-product secretion patterns, and growth rates on alternative substrates were measured for each strain as a function of evolutionary time. Three major conclusions were drawn from the experimental results. First, adaptive evolution leads to a phenotype characterized by maximized growth rates that may not correspond to the highest biomass yield. Second, metabolic phenotypes resulting from adaptive evolution can be described and predicted computationally. Third, adaptive evolution on a single substrate leads to changes in growth characteristics on other substrates that could signify parallel or opposing growth objectives. Together, the results show that genome-scale in silico metabolic models can describe the end point of adaptive evolution a priori and can be used to gain insight into the adaptive evolutionary process for E. coli.

Adaptation, Physiological↗

Characterization and efficacy of PKH26 as a probe to study the replication history of the human hematopoietic KG1a progenitor cell line.

The PKH26 dye can, in principle, be used for the study of asymmetric cell divisions (ASDs). A requirement for the identification of ASDs based on fluorescence intensity is that the PKH26 dye is distributed equally between daughter cells at each division, but this has not been demonstrated at a single-cell level. The efficacy of PKH26 as a probe for the study of ASDs was examined using the human hematopoietic KG1a cell. An automated time-lapse fluorescent microscope system was used to determine changes in cell size and fluorescence intensity during culture, and track cell divisions. The images of daughter cells were analyzed using the Isee software to determine the distribution of PKH26 dye between daughter cells. Ratios of cell size, mean fluorescence intensity, and total fluorescence intensity were calculated by dividing the values for one daughter cell by the value of the other daughter cell. The ratios for cell size, mean intensity, and total intensity were 1.13 +/- 0.12, 1.08 +/- 0.07, and 1.15 +/- 0.14 (mean +/- SD), respectively. Thus, PKH26 is not distributed equally to both daughter cells upon cell division. However, the replication history of individual KG1a cells can be reliably deduced for up to three divisions based solely on the mean and total fluorescence intensity of the PKH26 dye, using PKH26 concentrations below the chemical and phototoxic limits (2 microM).

Cell Division↗