Development of a genome-scale metabolic model for the lager hybrid yeast S. pastorianus to understand the evolution of metabolic pathways in industrial settings

Abstract

In silico tools such as genome-scale metabolic models have shown to be powerful for metabolic engineering of microorganisms. Saccharomyces pastorianus is a complex aneuploid hybrid between the mesophilic Saccharomyces cerevisiae and the cold-tolerant Saccharomyces eubayanus . This species is of biotechnological importance because it is the primary yeast used in lager beer fermentation and is also a key model for studying the evolution of hybrid genomes, including expression pattern of ortholog genes, composition of protein complexes, and phenotypic plasticity. Here, we created the iSP_1513 GSMM for S. pastorianus CBS1513 to allow top-down computational approaches to predict the evolution of metabolic pathways and to aid strain optimization in production processes. The iSP_1513 comprises 4,062 reactions, 1,808 alleles, and 2,747 metabolites, and takes into account the functional redundancy in the gene-protein-reaction rule caused by the presence of orthologous genes. Moreover, a universal algorithm to constrain GSMM reactions using transcriptome data was developed as a python library and enabled the integration of temperature as parameter. Essentiality data sets, growth data on various carbohydrates and volatile metabolites secretion were used to validate the model and showed the potential of media engineering to improve specific flavor compounds. The iSP_1513 also highlighted the different contributions of the parental sub-genomes to the oxidative and non-oxidative parts of the pentose phosphate pathway. Overall, the iSP_1513 GSMM represent an important step toward understanding the metabolic capabilities, evolutionary trajectories, and adaptation potential of S. pastorianus in different industrial settings. IMPORTANCE Genome-scale metabolic models (GSMM) have been successfully applied to predict cellular behavior and design cell factories in several model organisms, but no models to date are currently available for hybrid species due to their more complex genetics and general lack of molecular data. In this study, we generated a bespoke GSMM, iSP_1513, for this industrial aneuploid hybrid Saccharomyces pastorianus , which takes into account the aneuploidy and functional redundancy from orthologous parental alleles. This model will (i) help understand the metabolic capabilities and adaptive potential of S. pastorianus (domestication processes), (ii) aid top-down predictions for strain development (industrial biotechnology), and (iii) allow predictions of evolutionary trajectories of metabolic pathways in aneuploid hybrids (evolutionary genetics).