Integrating gene expression data into a genome-scale metabolic model to identify reprogramming during adaptive evolution

The development of a method for identifying latent reprogramming in gene expression data resulting from adaptive laboratory evolution (ALE) in response to genetic or environmental perturbations has been a challenge. In this study, a method called Metabolic Reprogramming Identifier (MRI), based on the integration of expression data to a genome-scale metabolic model has been developed. To identify key genes playing the main role in reprogramming, a MILP problem is presented and maximization of an adaptation score as a criterion indicating a pattern of using metabolism with maximum utilization of gene expression resources is defined as an objective function. Then, genes with complete expression usage and significant expression differences between wild-type and evolved strains were selected as key genes for reprogramming. This score is also applied to evaluate the compatibility of expression patterns with maximal use of key genes. The method was implemented to investigate the reprogramming of Escherichia coli during adaptive evolution caused by changing carbon sources. cyoC and cydB responsible for establishing proton gradient across the inner membrane were identified to be vital in the E. coli reprogramming when switching from glucose to lactate. These results indicate the importance of the inner membrane in reprogramming of E. coli to adapt to the new environment. The method predicts no reprogramming occurs during the evolution for growth on glycerol.