Abstract
Hypoxia affects the tumor microenvironment and is considered important to metastasis progression and therapy resistance. Thus far, the majority of global analyses of tumor hypoxia responses have been limited to just a single omics level. Combining multiple omics data can broaden our understanding of tumor hypoxia. Here, we investigate the temporal change of the metabolite composition with gene expression data from literature to provide a more comprehensive insight into the system level in response to hypoxia. Nuclear magnetic resonance spectroscopy was used to perform metabolomic profiling on the MDA-MB-231 breast cancer cell line under hypoxic conditions. Multivariate statistical analysis revealed that the metabolic difference between hypoxia and normoxia was similar over 24 h, but became distinct over 48 h. Time dependent microarray data from the same cell line in the literature displayed different gene expressions under hypoxic and normoxic conditions mostly at 12 h or earlier. The direct metabolomic profiles show a large overlap with theoretical metabolic profiles deduced from previous transcriptomic studies. Consistent pathways are glycolysis/gluconeogenesis, pyruvate, purine and arginine and proline metabolism. Ten metabolic pathways revealed by metabolomics were not covered by the downstream of the known transcriptomic profiles, suggesting new metabolic phenotypes. These results confirm previous transcriptomics understanding and expand the knowledge from existing models on correlation and co-regulation between transcriptomic and metabolomics profiles, which demonstrates the power of integrated omics analysis.