Metabolic reprogramming is a hallmark of cancer. The main aim of this paper is to integrate a genome-scale metabolic description of tumor cells into a tumor growth model that accounts for the spatiotemporally heterogeneous tumor microenvironment, in order to study the effects of microscopic characteristics on tumor evolution. A lactate maximization metabolic strategy that allows near-optimal growth solution, while maximizing lactate secretion, is assumed. The proposed sub-cellular metabolic model is then incorporated into a hybrid discrete-continuous model of tumor growth. We produced several phenotypes by applying different constraints and optimization criteria in the metabolic model and explored the tumor evolution of the various phenotypes in different vasculature conditions and extracellular matrix densities. At first, we showed that the metabolic capabilities of phenotypes depending on resource availability can vary in a counter-intuitive manner. We then showed that: first, tumor population, morphology, and spread are affected differently in different conditions, allowing thus phenotypes to be superior than others in different conditions; and second, polyclonal tumors consisting of different phenotypes can exploit their different metabolic capabilities to enhance further tumor evolution. The proposed framework comprises a proof-of-concept demonstration showing the importance of considering the metabolic capabilities of phenotypes on predicting tumor evolution. The proposed framework allows the incorporation of context-specific and patient-specific data for the study of personalized tumor evolution and therapy efficacy, linking genome to metabolic capabilities and tumor dynamics.