fig4

Figure 4. RNA-seq analysis of metabolic gene expression alteration between the MDA-MB-231HM.LNm5 and parental MDA-MB-231 cell lines. Expression level of all mitochondrial genes (MitoCarta 2.0) were compared (A), as well as genes involved in key processes of energy metabolism, including glycolysis (B) (glycolytic process: GO: 0061621 & 0006096; positive regulator of glycolytic process: GO: 0045821; negative regulator of glycolytic process: GO: 0045820; regulation of glycolytic process: GO: 0006110), tricarboxylic acid (TCA) cycle (C) (GO: 0006099), and the electron transport chain (D) (mitochondrial respiratory chain complexes: HGNC family ID: 639 & mitochondrial respiratory chain complex assembly factors HGNC family ID: 645). The log₂FC (y-axis) is derived from counts per million (CPM) values for MDA-MB-231HM.LNm5 divided by CPM values for MDA-MB-231, where a positive FC value represents up-regulation in the MDA-MB-231HM.LNm5 cells and a negative value represents down-regulation. Genes with a CPM value of < 1 across both samples were not included. FC: fold change; ENO3: enolase 3; HKDC1: hexokinase domain containing 1; MLXIPL: MLX interacting protein-like; FBP1: fructose-1,6-bisphosphatase-1; PRKAG2: protein kinase, AMP-activated, gamma 2 non-catalytic subunit; IDH2: isocitrate dehydrogenase-2; BCS1L: ubiquinol-cytochrome C reductase complex III chaperone