Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer

*Na Li, Xiaohan Zhan and Xianquan Zhan*

## **Abstract**

Energy metabolism heterogeneity is a hallmark in ovarian cancer; namely, the Warburg and reverse Warburg effects coexist in ovarian cancer. Exploration of energy metabolism heterogeneity benefits the discovery of the effective biomarkers for ovarian cancers. The integrative analysis of transcriptomics (20,115 genes in 419 ovarian cancer samples), proteomics (205 differentially expressed proteins), and mitochondrial proteomics (1198 mitochondrial differentially expressed proteins) revealed (i) the upregulations of rate-limiting enzymes PKM2 in glycolysis, IDH2 in Krebs cycle, and UQCRH in oxidative phosphorylation (OXPHOS) pathways, (ii) the upregulation of PDHB that converts pyruvate from glycolysis into acetyl-CoA in Krebs cycle, and (iii) that miRNA (hsa-miR-186-5p) and RNA-binding protein (EIF4AIII) had target sites in those key proteins in energy metabolism pathways. Furthermore, lncRNA SNHG3 interacted with miRNA (hsa-miR-186-5p) and RNA-binding protein (EIF4AIII). Those results were confirmed in the ovarian cancer cell model and tissues. It clearly concluded that lncRNA SNHG3 regulates energy metabolism through miRNA (hsa-miR-186-5p) and RNA-binding protein (EIF4AIII) to regulate the key proteins in the energy metabolism pathways. SNHG3 inhibitor might interfere with the energy metabolism to treat ovarian cancers. These findings provide more accurate understanding of molecular mechanisms of ovarian cancers and discovery of effective energy-metabolism-heterogeneity therapeutic drug for ovarian cancers.

### **Highlights**

**Keywords:** ovarian cancer, iTRAQ, mitochondrial proteomics, TCGA, energy metabolism, SNHG3

### **1. Introduction**

Ovarian cancer is a common gynecologic cancer with high mortality [1]. Despite chemotherapy, radiotherapy, surgery, and target therapy has previously been developed in ovarian cancers [2], the 5-year overall survival rate for patients who diagnosed with late stage III−IV disease is still very poor (about 30%). Because of the site of the ovaries and the certain clinical characteristics of epithelial cancers, it is a challenge to make early diagnosis [3]. Women with high-risk factors (e.g., family history, or BRCA mutations) plan for a follow-up visit with cancer antigen 125 (CA-125) monitoring and ultrasound, however, prospective validation of these physical examination or lab tests remain elusive [4]. The changes of energy metabolism are common in cancer cells, which might be potential biomarkers and therapeutics targets [5]. During the last decade, a great attention has been paid to metabolic reprogramming of cancer. However, cancer basic studies fail to reach a consistent conclusion on mitochondrial function in cancer energy metabolism [6]. The traditional view of Warburg was that cancer cells undergo aerobic glycolysis, which refers to the fermentation of glucose to lactate in the presence of oxygen as opposed to the complete oxidation of glucose, thus brought attention to the role of mitochondria in tumorigenesis [7]. A previous study found that the glycolysis enzyme PKM2 is important for cancer metabolism and tumor growth, which can improve activity and expression of PKM2 [8]. On the contrary, mitochondria were observed dysfunction, including the decreased effectiveness of Krebs cycle and electron transfer chain (ETC) complexes decoupling [9]. However, a novel 'reverse Warburg effect', was put forward in 2009 and impacted previous perceptions on cancer metabolism [10]. In this model of reverse Warburg chain, cancer cells and the cancer-associated fibroblasts (CAFs) become metabolically coupled. Interactions between cancer cells and tumor-microenvironment (TME) highly affect proliferation, energy metabolism, metastasis, and relapse of carcinoma [11]. Cancer cells secrete a large amount of ROS into microenvironment to enhance oxidative stress in CAFs. If the inflammatory reaction, autophagy, loss of stromal caveolin-1 (Cav-1), and nitric oxide synthase (NOS) are increased in CAFs, there is a good chance for progression of aerobic glycolysis [12]. Consequently, CAFs secrete plenty of energy-rich fuels to TME, including ketone bodies, lactate, pyruvate, and fatty acids. In turn, the nourishment 'feed' mitochondrial oxidative phosphorylation and ATP supplements [13]. In this process, mono-carboxylate transporters (MCTs) were highly expressed in both cancer cells and CAFs to be involved in some regulations. Immunochemistry result demonstrates that MCT4 was distributed specifically in CAFs in human breast cancers, which implicated in lactate efflux progress; while MCT1 participated in lactate uptake, and significantly upregulated specifically in kinds of cancer cells [14]. Thus evidence indicates limitations of 'the Warburg effect'. However, some studies demonstrated that aerobic glycolysis was not the dominant energy metabolism approach for many human cancer cell lines. In the past decades, studies on Warburg and reverse Warburg effects in cancers have formed a new frontier regarding additional roles of mitochondria in a cancer, and multiple functions of mitochondria have been identified in tumorigenesis [15].

High-throughput proteomics approach provides a scientific evaluation of protein expression. Functional proteomics offers more subtle clues, due to a greater attention paid to subcellular proteome research [16]. However, the subcellular proteomics of ovarian cancer mitochondrial proteins has not been elucidated. Mitochondria are the center of energy metabolism in eukaryotic cells, and also involved in other functions, such as cell signaling, cellular differentiation, cell death, and maintaining control of the cell cycle and oxidative stress regulation [17]. Those mitochondria-mediated biological processes are so closely associated

**73**

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

with tumor relapse or metastasis. Thus, cancer therapeutics should urgently find a way to explore molecular mechanisms of mitochondrion during tumorigenesis and tumor progression [18]. Inside the cancer cell appeared structural and morphological alterations of the mitochondria, and variations of morphology and performance are presumably associated with mitochondrial differentially expressed proteins (mtDEPs) [19]. A slight increase in research on ovarian cancer has occurred in recent years, quantitative proteomic analysis of mitochondria from human ovarian cancer cells and their paclitaxel-resistant sublines proved that the chemoresistance mechanisms were partly related to the mitochondria [20]. Mitochondria similarly impart considerable flexibility for tumor cell growth and survival in otherwise harsh environments such as during nutrient depletion, hypoxia and cancer treatments, and are therefore key players in tumorigenesis [15]. The subcellular proteomics of ovarian cancer mitochondrial proteins may offer new insights into aspect

Regeneration of energy metabolism plays crucial roles in the pathogenesis and development of cancer since it accelerates cancer cell growth, cell cycle, proliferation and metastasis [21]. The impact of non-coding RNAs (ncRNAs) has profoundly touched the fields of human cancers, cell biology, functional genomics, and drug therapy. Long non-coding RNAs (lncRNAs) (>200 nucleotides) and microRNAs (20–24 nucleotides) have attracted much attention, which acted as key regulators in the cellular biological processes, gene expression, gene regulation, basic biological functions of eukaryotic genomes, and post-transcriptional regulation of mRNA [22]. Recent studies demonstrated that lncRNAs were widely used as biomarkers for the diagnosis and prognosis of malignant tumors [23], and some lncRNAs can even act as the new therapeutic targets [24]. More and more researchers have turned attention to the mechanism between non-coding RNAs and malignant tumors. LncRNAs affects on energy metabolism-related signaling pathways induced epigenetic regulation [25]. MicroRNAs can silence gene expression by binding to 3′ untranslated region (3'UTR) sequences in their target messenger RNAs (mRNAs), resulting in the inhibition of translation or mRNA degradation, but the interaction of lncRNAs with microRNAs can hamper this effect [26]. The present results revealed that lncRNA FOXD2-AS1 acted as a tumor promoters partly through EphB3 inhibition by directly interacting with lysine (K)-specific demethylase 1A (LSD1)

and zeste homolog 2 (EZH2), which indicates that lncRNA-target gene-

reveal energy heterogeneity and its molecular mechanisms.

carcinogenesis axis for cancers does exist [27]. Here, it emphasizes important scientific associations of lncRNAs with energy metabolism in cancer cells. Increasing evidence indicates that lncRNAs play significant roles in cancer metabolism, and explore the potential mechanisms that could help elucidate regulation axis or network and provide a new direction for clinical management of different malignant phenotypes [28]. In our previous research, iTRAQ-based quantitative proteomics identified 1198 mitochondrial differentially expressed proteins (mtDEPs) between mitochondria samples isolated from human ovarian cancer and control tissues [29] and 205 differentially expressed proteins (DEPs) between human ovarian cancers and controls tissues [39]. The TCGA database includes 20,115 genes in 419 ovarian cancer samples. The conjoint analysis of 1198 mtDEPs, 205 DEPs, and 20,115 TCGA data in ovarian cancers investigated the biological pathways and molecular mechanisms of SNHG3-downstream genes-energy metabolism axis. LncRNA SNHG3 was associated with survival for ovarian cancers, and further gene set enrichment analysis proved the roles of SNHG3 in the energy metabolism through miRNAs and RNA binding protein EIF4AIII to target genes, including PKM, PDHB, IDH2, and UQCRH [29]. **Figure 1** showed the experimental flow-chart of integrative analysis of 1198 mtDEPs [29], 205 DEPs [39], and 20,115 TCGA data in ovarian cancers [29] to

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

of tumor development.

#### *Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

with tumor relapse or metastasis. Thus, cancer therapeutics should urgently find a way to explore molecular mechanisms of mitochondrion during tumorigenesis and tumor progression [18]. Inside the cancer cell appeared structural and morphological alterations of the mitochondria, and variations of morphology and performance are presumably associated with mitochondrial differentially expressed proteins (mtDEPs) [19]. A slight increase in research on ovarian cancer has occurred in recent years, quantitative proteomic analysis of mitochondria from human ovarian cancer cells and their paclitaxel-resistant sublines proved that the chemoresistance mechanisms were partly related to the mitochondria [20]. Mitochondria similarly impart considerable flexibility for tumor cell growth and survival in otherwise harsh environments such as during nutrient depletion, hypoxia and cancer treatments, and are therefore key players in tumorigenesis [15]. The subcellular proteomics of ovarian cancer mitochondrial proteins may offer new insights into aspect of tumor development.

Regeneration of energy metabolism plays crucial roles in the pathogenesis and development of cancer since it accelerates cancer cell growth, cell cycle, proliferation and metastasis [21]. The impact of non-coding RNAs (ncRNAs) has profoundly touched the fields of human cancers, cell biology, functional genomics, and drug therapy. Long non-coding RNAs (lncRNAs) (>200 nucleotides) and microRNAs (20–24 nucleotides) have attracted much attention, which acted as key regulators in the cellular biological processes, gene expression, gene regulation, basic biological functions of eukaryotic genomes, and post-transcriptional regulation of mRNA [22]. Recent studies demonstrated that lncRNAs were widely used as biomarkers for the diagnosis and prognosis of malignant tumors [23], and some lncRNAs can even act as the new therapeutic targets [24]. More and more researchers have turned attention to the mechanism between non-coding RNAs and malignant tumors. LncRNAs affects on energy metabolism-related signaling pathways induced epigenetic regulation [25]. MicroRNAs can silence gene expression by binding to 3′ untranslated region (3'UTR) sequences in their target messenger RNAs (mRNAs), resulting in the inhibition of translation or mRNA degradation, but the interaction of lncRNAs with microRNAs can hamper this effect [26]. The present results revealed that lncRNA FOXD2-AS1 acted as a tumor promoters partly through EphB3 inhibition by directly interacting with lysine (K)-specific demethylase 1A (LSD1) and zeste homolog 2 (EZH2), which indicates that lncRNA-target genecarcinogenesis axis for cancers does exist [27]. Here, it emphasizes important scientific associations of lncRNAs with energy metabolism in cancer cells. Increasing evidence indicates that lncRNAs play significant roles in cancer metabolism, and explore the potential mechanisms that could help elucidate regulation axis or network and provide a new direction for clinical management of different malignant phenotypes [28]. In our previous research, iTRAQ-based quantitative proteomics identified 1198 mitochondrial differentially expressed proteins (mtDEPs) between mitochondria samples isolated from human ovarian cancer and control tissues [29] and 205 differentially expressed proteins (DEPs) between human ovarian cancers and controls tissues [39]. The TCGA database includes 20,115 genes in 419 ovarian cancer samples. The conjoint analysis of 1198 mtDEPs, 205 DEPs, and 20,115 TCGA data in ovarian cancers investigated the biological pathways and molecular mechanisms of SNHG3-downstream genes-energy metabolism axis. LncRNA SNHG3 was associated with survival for ovarian cancers, and further gene set enrichment analysis proved the roles of SNHG3 in the energy metabolism through miRNAs and RNA binding protein EIF4AIII to target genes, including PKM, PDHB, IDH2, and UQCRH [29].

**Figure 1** showed the experimental flow-chart of integrative analysis of 1198 mtDEPs [29], 205 DEPs [39], and 20,115 TCGA data in ovarian cancers [29] to reveal energy heterogeneity and its molecular mechanisms.

*Molecular Medicine*

**1. Introduction**

Ovarian cancer is a common gynecologic cancer with high mortality [1]. Despite

chemotherapy, radiotherapy, surgery, and target therapy has previously been developed in ovarian cancers [2], the 5-year overall survival rate for patients who diagnosed with late stage III−IV disease is still very poor (about 30%). Because of the site of the ovaries and the certain clinical characteristics of epithelial cancers, it is a challenge to make early diagnosis [3]. Women with high-risk factors (e.g., family history, or BRCA mutations) plan for a follow-up visit with cancer antigen 125 (CA-125) monitoring and ultrasound, however, prospective validation of these physical examination or lab tests remain elusive [4]. The changes of energy metabolism are common in cancer cells, which might be potential biomarkers and therapeutics targets [5]. During the last decade, a great attention has been paid to metabolic reprogramming of cancer. However, cancer basic studies fail to reach a consistent conclusion on mitochondrial function in cancer energy metabolism [6]. The traditional view of Warburg was that cancer cells undergo aerobic glycolysis, which refers to the fermentation of glucose to lactate in the presence of oxygen as opposed to the complete oxidation of glucose, thus brought attention to the role of mitochondria in tumorigenesis [7]. A previous study found that the glycolysis enzyme PKM2 is important for cancer metabolism and tumor growth, which can improve activity and expression of PKM2 [8]. On the contrary, mitochondria were observed dysfunction, including the decreased effectiveness of Krebs cycle and electron transfer chain (ETC) complexes decoupling [9]. However, a novel 'reverse Warburg effect', was put forward in 2009 and impacted previous perceptions on cancer metabolism [10]. In this model of reverse Warburg chain, cancer cells and the cancer-associated fibroblasts (CAFs) become metabolically coupled. Interactions between cancer cells and tumor-microenvironment (TME) highly affect proliferation, energy metabolism, metastasis, and relapse of carcinoma [11]. Cancer cells secrete a large amount of ROS into microenvironment to enhance oxidative stress in CAFs. If the inflammatory reaction, autophagy, loss of stromal caveolin-1 (Cav-1), and nitric oxide synthase (NOS) are increased in CAFs, there is a good chance for progression of aerobic glycolysis [12]. Consequently, CAFs secrete plenty of energy-rich fuels to TME, including ketone bodies, lactate, pyruvate, and fatty acids. In turn, the nourishment 'feed' mitochondrial oxidative phosphorylation and ATP supplements [13]. In this process, mono-carboxylate transporters (MCTs) were highly expressed in both cancer cells and CAFs to be involved in some regulations. Immunochemistry result demonstrates that MCT4 was distributed specifically in CAFs in human breast cancers, which implicated in lactate efflux progress; while MCT1 participated in lactate uptake, and significantly upregulated specifically in kinds of cancer cells [14]. Thus evidence indicates limitations of 'the Warburg effect'. However, some studies demonstrated that aerobic glycolysis was not the dominant energy metabolism approach for many human cancer cell lines. In the past decades, studies on Warburg and reverse Warburg effects in cancers have formed a new frontier regarding additional roles of mitochondria in a cancer, and multiple functions of mitochondria have been identified in tumorigenesis [15]. High-throughput proteomics approach provides a scientific evaluation of protein expression. Functional proteomics offers more subtle clues, due to a greater attention paid to subcellular proteome research [16]. However, the subcellular proteomics of ovarian cancer mitochondrial proteins has not been elucidated. Mitochondria are the center of energy metabolism in eukaryotic cells, and also involved in other functions, such as cell signaling, cellular differentiation, cell death, and maintaining control of the cell cycle and oxidative stress regulation [17]. Those mitochondria-mediated biological processes are so closely associated

**72**

**Figure 1.**

*The experimental flow-chart to study energy metabolic heterogeneity and its molecular mechanisms.*

### **2. Methods**

### **2.1 Ovarian cancer mitochondrial DEP data and bioinformatic analysis**

Mitochondria were separated from 7 ovarian cancer tissues (high-degrade, poorly or moderately differentiated carcinoma cells) (cancer group) and 11 control ovaries with benign gynecologic diseases (fibroids, adenomyosis, ovary serous cystadenoma, cervical intraepithelial neoplasia, atypical hyperplasia of endometrium, and pelvic organ prolapse) (control group), respectively [29]. The separated mitochondria were validated with electron microscopy and Western blotting. The extracted proteins from the prepared mitochondrial samples were used for iTRAQ-quantitative proteomics analysis. The extracted mitochondrial proteins from ovarian cancers and controls were analyzed with 6-plex iTRAQ labeling, SCX fraction, and LC-MS/ MS. MS/MS data were used to determine proteins, and the intensities of iTRAQ reporter ions were used to determine each mitochondrial DEP. The mitochondrial DEPs were further analyzed by bioinformatics including GO functional enrichment and KEGG pathway enrichment with DAVID Bioinformatics Resources 6.7.

### **2.2 Ovarian cancer DEP data and bioinformatic analysis**

Proteins were extracted from ovarian cancer and control tissues. The extracted proteins from ovarian cancers and controls were analyzed with 6-plex iTRAQ labeling, SCX fraction, and LC-MS/MS. MS/MS data were used to determine proteins, and the intensities of iTRAQ reporter ions were used to determine each mitochondrial DEP [39]. The mitochondrial DEPs were further analyzed by bioinformatics including GO functional enrichment and KEGG pathway enrichment with DAVID Bioinformatics Resources 6.7.

### **2.3 TCGA data of ovarian cancer patients and bioinformatic analysis**

TCGA (http://cancergenome.nih.gov/) includes 20,115 genes of 419 ovarian cancer patients, in the level of transcriptome. Those genes were classified as coding/

**75**

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

non-coding RNAs (mRNAs/ncRNAs) provided by the GENCODE/ENSEMBL pipeline. lncRNA genes were considered as a type of genes that exclusively produce transcripts of the 'antisense'. The lncRNA survival analysis was performed by TANRIC (http://ibl.mdanderson.org/tanric/\_design/basic/index.html). The Kaplan-Meier method was used to calculate overall survival. According to median value (3.39) of SNHG3 RNA expressions, 419 ovarian cancer patients were divided into SNHG3 high (>3.39; n = 210) vs. low (<3.39; n = 209) expression groups. TCGA data of two groups were analyzed with GSEA enrichment analysis. Moreover, the lncRNA expressions from Cancer Cell Line Encyclopedia (https://portals.broadinstitute. org/ccle), and chemosensitivity of tamoxifen from Genomics of Drug Sensitivity in Cancer (http://www.cancerrxgene.org/) were obtained for ovarian cancer cell lines. GraphPad Prism v6.0 (GraphPad Software, San Diego, CA, USA) was used to

**2.4 Integrative analysis of mitochondrial DEPs, tissue DEPs, and TCGA data**

The integrated miRNA-lncRNA SNHG3, miRNA-target gene, RNA binding protein-lncRNA SNHG3, RNA binding protein-mRNA, and protein-protein signatures were identified. STRING 10.0 was used predict interactions of chemicals and proteins. Chemicals were linked to other chemicals and proteins by evidence derived from experiments, databases and literature (http://string-db.org/cgi/input.pl). The large-scale CLIP-Seq data by starBasev 2.0 (http://starbase.sysu.edu.cn/mirCircRNA.php) was used to construct SNHG3-miRNA, protein-miRNA, SNHG3-RNA binding protein, mRNA-RNA binding protein, and mRNA-microRNA-lncRNA interaction networks. The mitochondrial DEPs in ovarian cancers were input into STRING for protein–protein interaction analysis. Network visualizations were performed with Cytoscape 3.4.0 (http://www.cytoscape.org/). The binding sites of 3′UTR region of targeted genes were predicted with three publicly available databases (TargetScan, NCBI, and RNAhybrid), sequences of microRNA (>hsa-miR-186-5p MIMAT0000456 CAAAGAAUUCUCCUUUUGGGCU) and PDHB 3′UTR region. MicroRNA binding sites with PDHB were predicted with

Three ovarian cancer cell lines (TOV-21G, SK-OV3, and OVCAR-3), and one normal control cell line (IOSE80) from Keibai Academy of Science (Nanjing, China) were used. RPMI-1640 medium were used to culture TOV-21G and OVCAR-3 cells in 5% CO2 atmosphere at 37°C. DMEM medium (Corning, NY, USA) were used to culture IOSE80 and SK-OV3 in 5% CO2 atmosphere at 37°C, with supplementation of 10% fetal bovine serum (FBS, GIBCO, South America, NY, USA). (i) Transient transfection was performed with Lipofectamine 3000 reagents according to the manufacturer's instructions (Invitrogen, USA). SK-OV3, OVCAR-3, and TOV-21G were seeded in 6-well plates at 30–50% density. Cells were collected at 24–48 h after transfection, for next-step experiments. (ii) RNA extraction and quantitative real-time PCR (qRT-PCR) analyses. TRizol® Reagent (Invitrogen, CA, USA) was used to extract total RNAs. total RNAs were reversely transcribed into cDNAs and then used to perform qRT-PCR analysis to detect SNHG3 and its target genes, with β-actin as an internal control. (iii) 1D-SDS-PAGE and Western blotting was used to detect PKM, PFKM, PDHB, IDH2, CS, OGDHL, and UQCRH against the corresponding antibodies, with β-actin as internal control. (iv) Data were expressed as the mean ± SD of triplicates. Each experiment was repeated at least three times. In all cases, P < 0.05 was considered as statistical significance.

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

construct histogram.

RNAhybrid database.

**2.5 Experimental validation in cell models**

**with bioinformatics**

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

non-coding RNAs (mRNAs/ncRNAs) provided by the GENCODE/ENSEMBL pipeline. lncRNA genes were considered as a type of genes that exclusively produce transcripts of the 'antisense'. The lncRNA survival analysis was performed by TANRIC (http://ibl.mdanderson.org/tanric/\_design/basic/index.html). The Kaplan-Meier method was used to calculate overall survival. According to median value (3.39) of SNHG3 RNA expressions, 419 ovarian cancer patients were divided into SNHG3 high (>3.39; n = 210) vs. low (<3.39; n = 209) expression groups. TCGA data of two groups were analyzed with GSEA enrichment analysis. Moreover, the lncRNA expressions from Cancer Cell Line Encyclopedia (https://portals.broadinstitute. org/ccle), and chemosensitivity of tamoxifen from Genomics of Drug Sensitivity in Cancer (http://www.cancerrxgene.org/) were obtained for ovarian cancer cell lines. GraphPad Prism v6.0 (GraphPad Software, San Diego, CA, USA) was used to construct histogram.

### **2.4 Integrative analysis of mitochondrial DEPs, tissue DEPs, and TCGA data with bioinformatics**

The integrated miRNA-lncRNA SNHG3, miRNA-target gene, RNA binding protein-lncRNA SNHG3, RNA binding protein-mRNA, and protein-protein signatures were identified. STRING 10.0 was used predict interactions of chemicals and proteins. Chemicals were linked to other chemicals and proteins by evidence derived from experiments, databases and literature (http://string-db.org/cgi/input.pl). The large-scale CLIP-Seq data by starBasev 2.0 (http://starbase.sysu.edu.cn/mirCircRNA.php) was used to construct SNHG3-miRNA, protein-miRNA, SNHG3-RNA binding protein, mRNA-RNA binding protein, and mRNA-microRNA-lncRNA interaction networks. The mitochondrial DEPs in ovarian cancers were input into STRING for protein–protein interaction analysis. Network visualizations were performed with Cytoscape 3.4.0 (http://www.cytoscape.org/). The binding sites of 3′UTR region of targeted genes were predicted with three publicly available databases (TargetScan, NCBI, and RNAhybrid), sequences of microRNA (>hsa-miR-186-5p MIMAT0000456 CAAAGAAUUCUCCUUUUGGGCU) and PDHB 3′UTR region. MicroRNA binding sites with PDHB were predicted with RNAhybrid database.

#### **2.5 Experimental validation in cell models**

Three ovarian cancer cell lines (TOV-21G, SK-OV3, and OVCAR-3), and one normal control cell line (IOSE80) from Keibai Academy of Science (Nanjing, China) were used. RPMI-1640 medium were used to culture TOV-21G and OVCAR-3 cells in 5% CO2 atmosphere at 37°C. DMEM medium (Corning, NY, USA) were used to culture IOSE80 and SK-OV3 in 5% CO2 atmosphere at 37°C, with supplementation of 10% fetal bovine serum (FBS, GIBCO, South America, NY, USA). (i) Transient transfection was performed with Lipofectamine 3000 reagents according to the manufacturer's instructions (Invitrogen, USA). SK-OV3, OVCAR-3, and TOV-21G were seeded in 6-well plates at 30–50% density. Cells were collected at 24–48 h after transfection, for next-step experiments. (ii) RNA extraction and quantitative real-time PCR (qRT-PCR) analyses. TRizol® Reagent (Invitrogen, CA, USA) was used to extract total RNAs. total RNAs were reversely transcribed into cDNAs and then used to perform qRT-PCR analysis to detect SNHG3 and its target genes, with β-actin as an internal control. (iii) 1D-SDS-PAGE and Western blotting was used to detect PKM, PFKM, PDHB, IDH2, CS, OGDHL, and UQCRH against the corresponding antibodies, with β-actin as internal control. (iv) Data were expressed as the mean ± SD of triplicates. Each experiment was repeated at least three times. In all cases, P < 0.05 was considered as statistical significance.

*Molecular Medicine*

**2. Methods**

**Figure 1.**

**2.1 Ovarian cancer mitochondrial DEP data and bioinformatic analysis**

*The experimental flow-chart to study energy metabolic heterogeneity and its molecular mechanisms.*

and KEGG pathway enrichment with DAVID Bioinformatics Resources 6.7.

**2.3 TCGA data of ovarian cancer patients and bioinformatic analysis**

TCGA (http://cancergenome.nih.gov/) includes 20,115 genes of 419 ovarian cancer patients, in the level of transcriptome. Those genes were classified as coding/

Proteins were extracted from ovarian cancer and control tissues. The extracted proteins from ovarian cancers and controls were analyzed with 6-plex iTRAQ labeling, SCX fraction, and LC-MS/MS. MS/MS data were used to determine proteins, and the intensities of iTRAQ reporter ions were used to determine each mitochondrial DEP [39]. The mitochondrial DEPs were further analyzed by bioinformatics including GO functional enrichment and KEGG pathway enrichment with DAVID

**2.2 Ovarian cancer DEP data and bioinformatic analysis**

Mitochondria were separated from 7 ovarian cancer tissues (high-degrade, poorly or moderately differentiated carcinoma cells) (cancer group) and 11 control ovaries with benign gynecologic diseases (fibroids, adenomyosis, ovary serous cystadenoma, cervical intraepithelial neoplasia, atypical hyperplasia of endometrium, and pelvic organ prolapse) (control group), respectively [29]. The separated mitochondria were validated with electron microscopy and Western blotting. The extracted proteins from the prepared mitochondrial samples were used for iTRAQ-quantitative proteomics analysis. The extracted mitochondrial proteins from ovarian cancers and controls were analyzed with 6-plex iTRAQ labeling, SCX fraction, and LC-MS/ MS. MS/MS data were used to determine proteins, and the intensities of iTRAQ reporter ions were used to determine each mitochondrial DEP. The mitochondrial DEPs were further analyzed by bioinformatics including GO functional enrichment

**74**

Bioinformatics Resources 6.7.

### **3. Results and discussion**

### **3.1 The changes of key proteins in the energy metabolism signaling pathways**

The iTRAQ-based quantitative proteomics identified 1198 DEPs between mitochondria samples isolated from ovarian cancer and control tissues [29]. The statistically significant KEGG pathways were mined with DAVID Bioinformatics Resources from those mitochondrial DEPs between EOCs and controls, among which those DEPs were significantly enriched in the processes of Krebs cycle, and oxidative phosphorylation (OXPHOS) pathways. The key proteins (PDHB, IDH2, and UQCRH) were associated with aerobic oxidation to supply in the Krebs cycle, and oxidative phosphorylation was upregulated (**Figure 2**). Interestingly, those results were coincided with the reverse Warburg effect proposed in 2009 [10].

The iTRAQ-based quantitative proteomics identified 205 DEPs between ovarian cancer and control tissues [13], which revealed the upregulation of the key enzyme PKM2 in glycolysis pathway I in ovarian cancers. It was coincided with the Otto Warburg effect proposed in 1956 [30]. Warburg discovered that cancer cells tend to produce ATP by aerobic glycolysis, even though it is a less efficient pathway contrasted with OXPHOS. This popular system, called 'Warburg effect', has been the dominant mechanism of tumors for energy generations, while its relationship with tumorigenesis remains still unclear.

The research of the 'Warburg effect' mechanism of a cancer cell has never interrupted at home and abroad. PKM2, a splice isoform of the pyruvate kinase, serves as a major metabolic reprogramming regulator with an adjustable activity subjected to numerous allosteric effectors and post-translational modifications [31]. One observed that PKM2 modification was associated with an enhanced glucose consumption, level of lipid and DNA synthesis, and lactate productions, indicating

#### **Figure 2.**

*The changes of key proteins in oxidative phosphorylation, Krebs cycle, and glycolysis pathways. A. The network derived from identified genes and pathways. The square box was pathways. The circular box was DEPs. The edge was the gene–gene and gene-pathway interactions. The node color from green to red was the gene expression level from low to high. The edge color was the correlation level from low to high (gray = 0.70–0.90, black = 0.90–1.00). B. Oxidative phosphorylation, Krebs cycle, and glycolysis pathways were altered in ovarian cancers. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

**77**

**Figure 3.**

still unclear.

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

that PKM2 transformation promotes the Warburg effect [32]. Then, a novel series of inhibitors were developed to anti-Warburg-effect drugs for cancer treatment. For example, erastin-like anti-Warburg agents prevent mitochondrial depolarization induced by free tubulin and decrease lactate formation in cancer cells [33]. However, Warburg effect also has some limitations, because it completely ignored these facts that cancer cells had a great interaction with tumor microenvironment. In the 2009, a new model for understanding the Warburg effect was proposed in tumor energy metabolism. The hypothesis is that cancer cells induce the aerobic glycolysis in neighboring stromal fibroblasts. These cancer-associated fibroblasts (CAFs) secrete energy-rich substances, including lactate and pyruvate, to tumor microenvironment. These energy-rich metabolites were eaten up by adjacent cancer cell and used by mitochondrial TCA cycle, resulting in a higher energy producing capacity. It termed this new idea as the "reverse Warburg effect" [10]. Taken all together, the reverse Warburg effect is a new energy metabolic pattern identified between cancer cells and CAFs, but this novel pattern does not deny Warburg effect status and still cannot replace it. Actually, the reverse Warburg effect extends energy metabolism content, which explained the nature of the heterogeneity and plasticity of cancer metabolism [34]. Although it's validated that the 'reverse Warburg effect' can be initiated by oxidative stress in two compartment metabolic coupling and change of cellular electromagnetic filed, detailed mechanisms remain

In order to verify the above of views, each EOC cell line (SK-OV3, TOV-21G, and OVCAR-3) showed high expression of energy metabolism-related genes relative to control cells IOSE80 by qRTPCR, such as PKM, PDHB, IDH3A, IDH3B, ND5, ND2,

*The gene expression changes of key proteins in glycolysis, Krebs cycle, and oxidative phosphorylation pathways confirmed by qRT-PCR analysis in ovarian cancer cells. \*p < 0.05; \*\*p < 0.01; \*\*\*p < 0.001. N.S., nonsignificance. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

and CYB in EOC cell lines relative to IOSE80 (**Figure 3**).

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

that PKM2 transformation promotes the Warburg effect [32]. Then, a novel series of inhibitors were developed to anti-Warburg-effect drugs for cancer treatment. For example, erastin-like anti-Warburg agents prevent mitochondrial depolarization induced by free tubulin and decrease lactate formation in cancer cells [33]. However, Warburg effect also has some limitations, because it completely ignored these facts that cancer cells had a great interaction with tumor microenvironment.

In the 2009, a new model for understanding the Warburg effect was proposed in tumor energy metabolism. The hypothesis is that cancer cells induce the aerobic glycolysis in neighboring stromal fibroblasts. These cancer-associated fibroblasts (CAFs) secrete energy-rich substances, including lactate and pyruvate, to tumor microenvironment. These energy-rich metabolites were eaten up by adjacent cancer cell and used by mitochondrial TCA cycle, resulting in a higher energy producing capacity. It termed this new idea as the "reverse Warburg effect" [10]. Taken all together, the reverse Warburg effect is a new energy metabolic pattern identified between cancer cells and CAFs, but this novel pattern does not deny Warburg effect status and still cannot replace it. Actually, the reverse Warburg effect extends energy metabolism content, which explained the nature of the heterogeneity and plasticity of cancer metabolism [34]. Although it's validated that the 'reverse Warburg effect' can be initiated by oxidative stress in two compartment metabolic coupling and change of cellular electromagnetic filed, detailed mechanisms remain still unclear.

In order to verify the above of views, each EOC cell line (SK-OV3, TOV-21G, and OVCAR-3) showed high expression of energy metabolism-related genes relative to control cells IOSE80 by qRTPCR, such as PKM, PDHB, IDH3A, IDH3B, ND5, ND2, and CYB in EOC cell lines relative to IOSE80 (**Figure 3**).

#### **Figure 3.**

*The gene expression changes of key proteins in glycolysis, Krebs cycle, and oxidative phosphorylation pathways confirmed by qRT-PCR analysis in ovarian cancer cells. \*p < 0.05; \*\*p < 0.01; \*\*\*p < 0.001. N.S., nonsignificance. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

*Molecular Medicine*

**3. Results and discussion**

tumorigenesis remains still unclear.

**3.1 The changes of key proteins in the energy metabolism signaling pathways**

The iTRAQ-based quantitative proteomics identified 1198 DEPs between mitochondria samples isolated from ovarian cancer and control tissues [29]. The statistically significant KEGG pathways were mined with DAVID Bioinformatics Resources from those mitochondrial DEPs between EOCs and controls, among which those DEPs were significantly enriched in the processes of Krebs cycle, and oxidative phosphorylation (OXPHOS) pathways. The key proteins (PDHB, IDH2, and UQCRH) were associated with aerobic oxidation to supply in the Krebs cycle, and oxidative phosphorylation was upregulated (**Figure 2**). Interestingly, those results were coincided with the reverse Warburg effect proposed in 2009 [10].

The iTRAQ-based quantitative proteomics identified 205 DEPs between ovarian cancer and control tissues [13], which revealed the upregulation of the key enzyme PKM2 in glycolysis pathway I in ovarian cancers. It was coincided with the Otto Warburg effect proposed in 1956 [30]. Warburg discovered that cancer cells tend to produce ATP by aerobic glycolysis, even though it is a less efficient pathway contrasted with OXPHOS. This popular system, called 'Warburg effect', has been the dominant mechanism of tumors for energy generations, while its relationship with

The research of the 'Warburg effect' mechanism of a cancer cell has never interrupted at home and abroad. PKM2, a splice isoform of the pyruvate kinase, serves as a major metabolic reprogramming regulator with an adjustable activity subjected to numerous allosteric effectors and post-translational modifications [31]. One observed that PKM2 modification was associated with an enhanced glucose consumption, level of lipid and DNA synthesis, and lactate productions, indicating

**76**

**Figure 2.**

*The changes of key proteins in oxidative phosphorylation, Krebs cycle, and glycolysis pathways. A. The network derived from identified genes and pathways. The square box was pathways. The circular box was DEPs. The edge was the gene–gene and gene-pathway interactions. The node color from green to red was the gene expression level from low to high. The edge color was the correlation level from low to high (gray = 0.70–0.90, black = 0.90–1.00). B. Oxidative phosphorylation, Krebs cycle, and glycolysis pathways were altered in ovarian*

*cancers. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

### **3.2 SNHG3 was significantly related to EOC survival through the key molecules in the energy metabolism pathways by their RNA-binding proteins or miRNA**

Among identified 1198 mitochondrial DEPs and 205 tissue DEPs, PFKM, PKM, PDHB, CS, IDH2, IDH3A, IDH3B, OGDHL, ND2, ND5, CYB, and UQCRH were enriched in glycolysis, Krebs cycle, and oxidative phosphorylation pathways. Six RNA-binding proteins (EIF4AIII, IGF2BP2, C22ORF28, UPF1, SFRS1, and EWSR1) were the iTRAQ-identified proteins in ovarian cancers (**Figure 4A**) based on Starbase v2.0 database. However, only EIF4AIII was associated with energy metabolic pathway, when did protein-protein network by STRING 10.0 software (**Figure 4B**). Furthermore, overlapping analysis between survival-related lncRNAs of EOC and lncRNAs binding with EIF4AIII obtained 16 lncRNAs (LINC00517, SNHG3, LBX2-AS1, ZNRF3-AS1, LINC00565, AL109767.1, WWTR1-AS1, HCG15, LEMD1-AS1, PDCD4-AS1, KIF9-AS1, SOS1-IT1, STARD13-IT1, PLCH1-AS1, ZNF674-AS1, and HOXC-AS3) existed those two groups. Among those 16 overlapped lncRNAs, only lncRNA SNHG3 was associated with energy metabolic pathways by GSEA analysis (**Figure 4C** and **D** and **Figure 4A** and **B**). The expression levels of lncRNA SNHG3 in different ovarian cancer cell lines indicated that poorly differentiated cell lines existed high SNHG3 expression, such as TYK-nu ovarian cancer cell line (**Figure 4E**). Additionally, q-PCR data demonstrated that SNHG3 was upregulated in SKOV3, TOV21G, and OVCAR3 relative to control cell line (IOSE80) (**Figure 4F**).

#### **Figure 4.**

*SNHG3 was significantly related to ovarian cancer survival through the key molecules in the energy metabolism pathways by their RNA-binding proteins or miRNA. A. Overlapping analysis of identified proteins and RNA binding proteins. B. Target DEPs-based protein-protein interaction network (STRING 10.0). C. Overlapping analysis of EIF4AIII-binding lncRNAs and lncRNAs involved in ovarian cancer survival. D. Kaplan–Meier survival analysis based on ovarian cancer SNHG3. E and F. SNHG3 expressions in ovarian cell lines. \*\*p < 0.01; \*\*\*p < 0.001. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

**79**

**Figure 5.**

*Elsevier, copyright 2018.*

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

may lead to ways to improve drug sensitivity of tamoxifen in EOCs.

One should also notice that gene sets were also significantly enriched in pathway

of tamoxifen response (**Figure 5**), which indicated SNHG3 was associated with drug sensitivity or multidrug resistance. The comprehensive evaluation on SNHG3

Gene sets enrichment analysis showed that mRNA metabolism and 3′UTRmediated translational regulation (**Figure 6**). Overlap analysis of RNAs-RNAs interaction networks showed that SNHG3 may regulated PDHB through binding hsa-miR-186-5p or hsa-miR-590-3p (**Figure 7A**), especially, hsa-miR-186-5p obtained high stringency to target PDHB with Starbase 2.0 analysis. Meanwhile, two binding sites were predicted between putative hsa-miR-186-5p and PDHB 3'UTR with RNAhybrid database (**Figure 7B** and **C**). Here, it can be forecasted boldly that SNHG3 might regulate the EOC energy metabolism by binding EIF4AIII and hsa-miR-186-5p, functioned as efficient sponges to regulate energy metabolism

pathways though mitochondrial key molecules (**Figures 7D** and **8A** and **B**).

miRNA sponge, ceRNA (competing endogenous RNAs) [38].

tumor energy metabolism would be rather influential in EOCs.

**3.3 Potential therapeutic targets in metabolic symbiosis**

To further verify that SNHG3 can lead to the carcinogenesis in vivo, SKOV3 cells were transfected with either si-SNHG3 or a si-RNA negative control. Target genes, including PFKM, PKM, PDHB, CS, IDH2, IDH3A, IDH3B, OGDHL, ND5, ND2, CYB and UQCRH turned significant decrease expression (**Figure 9**). The results were further validated to a reasonable degree by Western blot (**Figure 10**). Non coding RNAs, as one of epigenetic regulation form, play an important role in activation and suppression in a tumor by altering cell energy metabolism or biological behaviors [35]. However, lncRNAs have been identified and reported to be related to many kinds of carcinomas, little is known about lncRNAs whole molecular mechanisms in tumor energy metabolism. Recently, discovery of novel biomarkers focuses on ncRNAs, such as miR-125a, MALAT1, let-7a, miR-196a, HOXA11-AS, and lncRNA FAL1 [36]. Some biomarkers have been verified consistency in both tissues and serum, which improved clinical application value to use in early diagnosis or monitoring patient prognosis [37]. A number of studies have shown that lncRNAs can play an important role in tumorigenesis and progression through a variety of mechanisms, such as binding transcription factor, acting as

Therefore, lncRNA as an effective screening and their potential mechanisms in

Tumor tissues were made up by parenchymal cells and stromal elements. Parenchymal cells probably showed metabolism heterogeneity. So some cancer cells were high glycolytic cancer cells consisting with "Warburg effect", while

*The results of SNHG3 by GSEA analysis revealed tamoxifen pathway. A. Genes were enriched in response to tamoxifen. B. IC50 of tamoxifen in EOC cell lines. Reproduced from Li et al. [29], with permission from*

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

#### *Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

One should also notice that gene sets were also significantly enriched in pathway of tamoxifen response (**Figure 5**), which indicated SNHG3 was associated with drug sensitivity or multidrug resistance. The comprehensive evaluation on SNHG3 may lead to ways to improve drug sensitivity of tamoxifen in EOCs.

Gene sets enrichment analysis showed that mRNA metabolism and 3′UTRmediated translational regulation (**Figure 6**). Overlap analysis of RNAs-RNAs interaction networks showed that SNHG3 may regulated PDHB through binding hsa-miR-186-5p or hsa-miR-590-3p (**Figure 7A**), especially, hsa-miR-186-5p obtained high stringency to target PDHB with Starbase 2.0 analysis. Meanwhile, two binding sites were predicted between putative hsa-miR-186-5p and PDHB 3'UTR with RNAhybrid database (**Figure 7B** and **C**). Here, it can be forecasted boldly that SNHG3 might regulate the EOC energy metabolism by binding EIF4AIII and hsa-miR-186-5p, functioned as efficient sponges to regulate energy metabolism pathways though mitochondrial key molecules (**Figures 7D** and **8A** and **B**).

To further verify that SNHG3 can lead to the carcinogenesis in vivo, SKOV3 cells were transfected with either si-SNHG3 or a si-RNA negative control. Target genes, including PFKM, PKM, PDHB, CS, IDH2, IDH3A, IDH3B, OGDHL, ND5, ND2, CYB and UQCRH turned significant decrease expression (**Figure 9**). The results were further validated to a reasonable degree by Western blot (**Figure 10**).

Non coding RNAs, as one of epigenetic regulation form, play an important role in activation and suppression in a tumor by altering cell energy metabolism or biological behaviors [35]. However, lncRNAs have been identified and reported to be related to many kinds of carcinomas, little is known about lncRNAs whole molecular mechanisms in tumor energy metabolism. Recently, discovery of novel biomarkers focuses on ncRNAs, such as miR-125a, MALAT1, let-7a, miR-196a, HOXA11-AS, and lncRNA FAL1 [36]. Some biomarkers have been verified consistency in both tissues and serum, which improved clinical application value to use in early diagnosis or monitoring patient prognosis [37]. A number of studies have shown that lncRNAs can play an important role in tumorigenesis and progression through a variety of mechanisms, such as binding transcription factor, acting as miRNA sponge, ceRNA (competing endogenous RNAs) [38].

Therefore, lncRNA as an effective screening and their potential mechanisms in tumor energy metabolism would be rather influential in EOCs.

#### **3.3 Potential therapeutic targets in metabolic symbiosis**

Tumor tissues were made up by parenchymal cells and stromal elements. Parenchymal cells probably showed metabolism heterogeneity. So some cancer cells were high glycolytic cancer cells consisting with "Warburg effect", while

#### **Figure 5.**

*The results of SNHG3 by GSEA analysis revealed tamoxifen pathway. A. Genes were enriched in response to tamoxifen. B. IC50 of tamoxifen in EOC cell lines. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

*Molecular Medicine*

**78**

**Figure 4.**

*SNHG3 was significantly related to ovarian cancer survival through the key molecules in the energy metabolism pathways by their RNA-binding proteins or miRNA. A. Overlapping analysis of identified proteins and RNA binding proteins. B. Target DEPs-based protein-protein interaction network (STRING 10.0). C. Overlapping analysis of EIF4AIII-binding lncRNAs and lncRNAs involved in ovarian cancer survival. D. Kaplan–Meier survival analysis based on ovarian cancer SNHG3. E and F. SNHG3 expressions in ovarian cell lines. \*\*p < 0.01; \*\*\*p < 0.001. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

**3.2 SNHG3 was significantly related to EOC survival through the key molecules in the energy metabolism pathways by their RNA-binding proteins or miRNA**

Among identified 1198 mitochondrial DEPs and 205 tissue DEPs, PFKM, PKM, PDHB, CS, IDH2, IDH3A, IDH3B, OGDHL, ND2, ND5, CYB, and UQCRH were enriched in glycolysis, Krebs cycle, and oxidative phosphorylation pathways. Six RNA-binding proteins (EIF4AIII, IGF2BP2, C22ORF28, UPF1, SFRS1, and EWSR1) were the iTRAQ-identified proteins in ovarian cancers (**Figure 4A**) based on Starbase v2.0 database. However, only EIF4AIII was associated with energy metabolic pathway,

when did protein-protein network by STRING 10.0 software (**Figure 4B**). Furthermore, overlapping analysis between survival-related lncRNAs of EOC and lncRNAs binding with EIF4AIII obtained 16 lncRNAs (LINC00517, SNHG3, LBX2-AS1, ZNRF3-AS1, LINC00565, AL109767.1, WWTR1-AS1, HCG15, LEMD1-AS1, PDCD4-AS1, KIF9-AS1, SOS1-IT1, STARD13-IT1, PLCH1-AS1, ZNF674-AS1, and HOXC-AS3) existed those two groups. Among those 16 overlapped lncRNAs, only lncRNA SNHG3 was associated with energy metabolic pathways by GSEA analysis (**Figure 4C** and **D** and **Figure 4A** and **B**). The expression levels of lncRNA SNHG3 in different ovarian cancer cell lines indicated that poorly differentiated cell lines existed

high SNHG3 expression, such as TYK-nu ovarian cancer cell line (**Figure 4E**). Additionally, q-PCR data demonstrated that SNHG3 was upregulated in SKOV3, TOV21G, and OVCAR3 relative to control cell line (IOSE80) (**Figure 4F**).

#### **Figure 6.**

*The results of SNHG3 by GSEA analysis revealed mRNA metabolism and reactome 3*′*UTR-mediated translational regulation. A. Genes were enriched in metabolism of mRNAs. B. Genes were mainly enriched in reactome 3*′*UTR-mediated translational regulation. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

#### **Figure 7.**

*Overlapping analysis of the microRNA/mRNA and microRNA/lncRNAs interactive network. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

other cancer cells were oxidative cancer cells consisting with "reverse Warburg effect". Cancer cells and stroma cells (especially CAFs) have metabolic symbiosis, so cancer cells induce oxidative stress of CAFs by secreting ROS to enhance aerobic glycolysis of CAFs. In turn, CAFs produced lots of nourishment to be 'eaten' up by cancer cells for producing ATP through Krebs cycle and oxidative phosphorylation [13].

MCT-1 and MCT-4 were overexpressed in EOC cells by qRT-PCR experiments, including SKOV3, TOV21G and OVCAR3 (**Figure 11**). Even though tumors were characterized by metabolic heterogeneity, MCT-1 and MCT-4 were just like lactate shuttle between cancer cells and stroma cells. The nanomaterial-siRNAs of SNHG3 might be promising for EOC patients to block the abnormal energy metabolism (**Figure 12**).

**81**

**Figure 10.**

**Figure 8.**

**Figure 9.**

*2018.*

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

*Schematic model of the potential signaling mechanisms between SNHG3 and energy metabolism in the EOC*

*The mRNA expression levels of target genes of SNHG3 in EOC cells were determined by qRT-PCR. \*p < 0.05, \*\* p < 0.01. N.S., non-significance. Reproduced from Li et al. [29], with permission from Elsevier, copyright*

*The protein expression levels of target genes of SNHG3 in EOC cells were determined by Western blot.*

*Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

*regulation. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

#### **Figure 8.**

*Molecular Medicine*

**Figure 6.**

*Elsevier, copyright 2018.*

*The results of SNHG3 by GSEA analysis revealed mRNA metabolism and reactome 3*′*UTR-mediated translational regulation. A. Genes were enriched in metabolism of mRNAs. B. Genes were mainly enriched in reactome 3*′*UTR-mediated translational regulation. Reproduced from Li et al. [29], with permission from*

other cancer cells were oxidative cancer cells consisting with "reverse Warburg effect". Cancer cells and stroma cells (especially CAFs) have metabolic symbiosis, so cancer cells induce oxidative stress of CAFs by secreting ROS to enhance aerobic glycolysis of CAFs. In turn, CAFs produced lots of nourishment to be 'eaten' up by cancer cells for producing ATP through Krebs cycle and oxidative

*Overlapping analysis of the microRNA/mRNA and microRNA/lncRNAs interactive network. Reproduced*

*from Li et al. [29], with permission from Elsevier, copyright 2018.*

MCT-1 and MCT-4 were overexpressed in EOC cells by qRT-PCR experiments, including SKOV3, TOV21G and OVCAR3 (**Figure 11**). Even though tumors were characterized by metabolic heterogeneity, MCT-1 and MCT-4 were just like lactate shuttle between cancer cells and stroma cells. The nanomaterial-siRNAs of SNHG3 might be promising for EOC patients to block the abnormal energy metabolism

**80**

(**Figure 12**).

**Figure 7.**

phosphorylation [13].

*Schematic model of the potential signaling mechanisms between SNHG3 and energy metabolism in the EOC regulation. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

#### **Figure 9.**

*The mRNA expression levels of target genes of SNHG3 in EOC cells were determined by qRT-PCR. \*p < 0.05, \*\* p < 0.01. N.S., non-significance. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

#### **Figure 10.**

*The protein expression levels of target genes of SNHG3 in EOC cells were determined by Western blot. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

#### **Figure 11.**

*The expression levels of MCT-1 and MCT-4 in EOC cells were determined by qRT-PCR. \*\*p < 0.01, \*\*\*p < 0.001. N.S., non-significance. Reproduced from Li et al. [29], with permission from Elsevier, copyright 2018.*

#### **Figure 12.**

*Energy metabolic heterogeneity-based potential therapeutic targets model. Parenchymal cells demonstrated energy metabolic heterogeneity. Some cancer cells showed the "Warburg effect" with highly glycolytic functions, and other cancer cells showed the "reverse Warburg effect" with oxidative cancer cells. The metabolic symbiosis existed between tumor cells and CAFs through MCTs. The RNA interference sequence of SNHG3 might be effective. Modified from Li et al. [29], with permission from Elsevier, copyright 2018.*

### **4. Conclusions**

The identified 1198 mitochondrial DEPs, 205 tissue DEPs, and TCGA data in ovarian cancers provide new insights into human ovarian cancers, particularly the energy metabolism heterogeneity that 'Warburg effect' and 'reverse Warburg effect' were coexisted in ovarian cancer tissues. It emphasizes the important scientific merit in identity of new useful biomarkers within EOC energy metabolism heterogeneity system for the diagnosis and prognosis of ovarian cancer, and discovery of some potential therapeutic targets in energy metabolic interactions. Moreover, SNHG3 was related to energy metabolism through regulating hsamiR-186-5p or RNA binding protein EIF4AIII, and those two molecules had target

**83**

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

sites with key proteins in TCA cycle and oxidative phosphorylation pathways (PDHB, IDH2, and UQCRH). Therefore, energy metabolism-based target treatments might be very promising for ovarian cancer patients to block both 'Warburg

The authors acknowledge the financial supports from the Xiangya Hospital Funds

for Talent Introduction (to X.Z.), the Hunan Provincial "Hundred Talent Plan" program (to X.Z.), the National Natural Science Foundation of China (Grant nos. 81572278 and 81272798 to X.Z.), China "863" Plan Project (Grant No. 2014AA020610-1 to X.Z.), and the Hunan Provincial Natural Science Foundation of China (Grant No.

We declare that we have no financial and personal relationships with other

N.L. prepared figures, and wrote manuscript draft of the book chapter. X.H.Z. participated in writing and revising manuscript. X.Z. conceived the concept, designed the book chapter, prepared figures, wrote and critically revised the book chapter, coordinated and was responsible for the correspondence work and financial support.

iTRAQ isobaric target for relative and absolute quantification

mtDEPs mitochondrial differentially expressed proteins

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

effect' and 'reverse Warburg effect'.

**Acknowledgements**

14JJ7008 to X.Z.).

**Conflict of interest**

people or organizations.

**Author's contributions**

**Acronyms and abbreviations**

CAFs cancer-associated fibroblasts DEPs differentially expressed proteins

MCTs mono-carboxylate transporters

MS/MS tandem mass spectrometry

ETC electron transfer chain

LC liquid chromatography lncRNAs long non-coding RNAs

mRNAs messenger RNAs

ncRNAs non-coding RNAs NOS nitric oxide synthase OXPHOS oxidative phosphorylation SCX strong cation exchange TME tumor-microenvironment

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

sites with key proteins in TCA cycle and oxidative phosphorylation pathways (PDHB, IDH2, and UQCRH). Therefore, energy metabolism-based target treatments might be very promising for ovarian cancer patients to block both 'Warburg effect' and 'reverse Warburg effect'.

### **Acknowledgements**

*Molecular Medicine*

**Figure 11.**

*2018.*

**82**

**4. Conclusions**

**Figure 12.**

The identified 1198 mitochondrial DEPs, 205 tissue DEPs, and TCGA data in ovarian cancers provide new insights into human ovarian cancers, particularly the energy metabolism heterogeneity that 'Warburg effect' and 'reverse Warburg effect' were coexisted in ovarian cancer tissues. It emphasizes the important scientific merit in identity of new useful biomarkers within EOC energy metabolism heterogeneity system for the diagnosis and prognosis of ovarian cancer, and discovery of some potential therapeutic targets in energy metabolic interactions. Moreover, SNHG3 was related to energy metabolism through regulating hsamiR-186-5p or RNA binding protein EIF4AIII, and those two molecules had target

*Energy metabolic heterogeneity-based potential therapeutic targets model. Parenchymal cells demonstrated energy metabolic heterogeneity. Some cancer cells showed the "Warburg effect" with highly glycolytic functions, and other cancer cells showed the "reverse Warburg effect" with oxidative cancer cells. The metabolic symbiosis existed between tumor cells and CAFs through MCTs. The RNA interference sequence of SNHG3 might be*

*effective. Modified from Li et al. [29], with permission from Elsevier, copyright 2018.*

*The expression levels of MCT-1 and MCT-4 in EOC cells were determined by qRT-PCR. \*\*p < 0.01, \*\*\*p < 0.001. N.S., non-significance. Reproduced from Li et al. [29], with permission from Elsevier, copyright*

The authors acknowledge the financial supports from the Xiangya Hospital Funds for Talent Introduction (to X.Z.), the Hunan Provincial "Hundred Talent Plan" program (to X.Z.), the National Natural Science Foundation of China (Grant nos. 81572278 and 81272798 to X.Z.), China "863" Plan Project (Grant No. 2014AA020610-1 to X.Z.), and the Hunan Provincial Natural Science Foundation of China (Grant No. 14JJ7008 to X.Z.).

### **Conflict of interest**

We declare that we have no financial and personal relationships with other people or organizations.

### **Author's contributions**

N.L. prepared figures, and wrote manuscript draft of the book chapter. X.H.Z. participated in writing and revising manuscript. X.Z. conceived the concept, designed the book chapter, prepared figures, wrote and critically revised the book chapter, coordinated and was responsible for the correspondence work and financial support.

### **Acronyms and abbreviations**

*Molecular Medicine*

### **Author details**

Na Li1,2, Xiaohan Zhan1,2 and Xianquan Zhan1,2\*

1 Key Laboratory of Cancer Proteomics of Chinese Ministry of Health, Xiangya Hospital, Central South University, Changsha, China

2 State Local Joint Engineering Laboratory for Anticancer Drugs, Xiangya Hospital, Central South University, Changsha, China

\*Address all correspondence to<EMAIL_ADDRESS>
© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

**85**

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

[8] Christofk HR, Vander Heiden MG, Harris MH, Ramanathan A, Gerszten RE, Wei R, et al. The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth.

Nature. 2008;**452**:230-233. DOI:

[9] Birsoy K, Wang T, Chen WW, Freinkman E, Abu-Remaileh M, Sabatini DM. An essential role of the mitochondrial electron transport chain in cell proliferation is to enable aspartate synthesis. Cell. 2015;**162**:540-551. DOI:

[10] Pavlides S, Whitaker-Menezes D, Castello-Cros R, Flomenberg N, Witkiewicz AK, Frank PG, et al. The reverse Warburg effect: Aerobic glycolysis in cancer associated

fibroblasts and the tumor stroma. Cell Cycle. 2009;**8**:3984-4001. DOI: 10.4161/

[11] Catteau X, Simon P, Noel JC. Stromal expression of matrix metalloproteinase 2 in cancer-associated fibroblasts is strongly related to human epidermal growth factor receptor 2 status in invasive breast carcinoma. Molecular and Clinical Oncology. 2016;**4**:375-378.

[12] Lee M, Yoon JH. Metabolic interplay between glycolysis and mitochondrial oxidation: The reverse Warburg effect and its therapeutic implication. World Journal of Biological Chemistry. 2015;**6**:148-161. DOI: 10.4331/wjbc.

[13] Bar-Or D, Carrick M, Tanner A 2nd, Lieser MJ, Rael LT, Brody E. Overcoming the Warburg effect: Is it the key to survival in sepsis? Journal of Critical Care. 2018;**43**:197-201. DOI:

[14] Whitaker-Menezes D, Martinez-Outschoorn UE, Lin Z, Ertel A,

10.1016/j.jcrc.2017.09.012

DOI: 10.3892/mco.2015.721

10.1038/nature06734

10.1016/j.cell.2015.07.016

cc.8.23.10238

v6.i3.148

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

[1] Torre LA, Trabert B, DeSantis CE. Ovarian cancer statistics. CA: A Cancer Journal for Clinicians. 2018;**68**:284-296.

[2] Kreienbring K, Franz A, Richter R, Dragun D, Heidecke H, Muller D, et al. The role of PAR1 autoantibodies in patients with primary epithelial ovarian cancer. Anticancer Research. 2018;**38**:3619-3625. DOI: 10.21873/

[3] Luo L, Dong LY, Yan QG, Cao SJ, Wen XT, Huang Y, et al. Research progress in applying proteomics technology to explore early diagnosis biomarkers of breast cancer, lung cancer and ovarian cancer. Asian Pacific Journal of Cancer Prevention.

[4] Urban RR, Pappas TC, Bullock RG, Munroe DG, Bonato V, Agnew K, et al. Combined symptom index and secondgeneration multivariate biomarker test for prediction of ovarian cancer in patients with an adnexal mass. Gynecologic Oncology. 2018. DOI: 10.1016/j.ygyno.2018.06.004

[5] Antonosante A, d'Angelo M, Castelli V, Catanesi M, Iannotta D, Giordano A, et al. The involvement of PPARs in the peculiar energetic metabolism of tumor cells. International

Journal of Molecular Sciences. 2018;**19**(7). pii: E1907. DOI: 10.3390/

[6] Sebastian C. Tracking down the origin of cancer: Metabolic

Oncogenesis. 2014;**19**:363-382

effect in tumor progression:

10.1016/j.canlet.2014.04.001

reprogramming as a driver of stemness and tumorigenesis. Critical Reviews in

[7] Lu J, Tan M, Cai Q. The Warburg

Mitochondrial oxidative metabolism as an anti-metastasis mechanism. Cancer Letters. 2015;**356**:156-164. DOI:

ijms19071907

DOI: 10.3322/caac.21456

anticanres.12636

2014;**15**:8529-8538

**References**

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

### **References**

*Molecular Medicine*

**84**

**Author details**

provided the original work is properly cited.

Na Li1,2, Xiaohan Zhan1,2 and Xianquan Zhan1,2\*

Central South University, Changsha, China

Hospital, Central South University, Changsha, China

\*Address all correspondence to<EMAIL_ADDRESS>
© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

1 Key Laboratory of Cancer Proteomics of Chinese Ministry of Health, Xiangya

2 State Local Joint Engineering Laboratory for Anticancer Drugs, Xiangya Hospital,

[1] Torre LA, Trabert B, DeSantis CE. Ovarian cancer statistics. CA: A Cancer Journal for Clinicians. 2018;**68**:284-296. DOI: 10.3322/caac.21456

[2] Kreienbring K, Franz A, Richter R, Dragun D, Heidecke H, Muller D, et al. The role of PAR1 autoantibodies in patients with primary epithelial ovarian cancer. Anticancer Research. 2018;**38**:3619-3625. DOI: 10.21873/ anticanres.12636

[3] Luo L, Dong LY, Yan QG, Cao SJ, Wen XT, Huang Y, et al. Research progress in applying proteomics technology to explore early diagnosis biomarkers of breast cancer, lung cancer and ovarian cancer. Asian Pacific Journal of Cancer Prevention. 2014;**15**:8529-8538

[4] Urban RR, Pappas TC, Bullock RG, Munroe DG, Bonato V, Agnew K, et al. Combined symptom index and secondgeneration multivariate biomarker test for prediction of ovarian cancer in patients with an adnexal mass. Gynecologic Oncology. 2018. DOI: 10.1016/j.ygyno.2018.06.004

[5] Antonosante A, d'Angelo M, Castelli V, Catanesi M, Iannotta D, Giordano A, et al. The involvement of PPARs in the peculiar energetic metabolism of tumor cells. International Journal of Molecular Sciences. 2018;**19**(7). pii: E1907. DOI: 10.3390/ ijms19071907

[6] Sebastian C. Tracking down the origin of cancer: Metabolic reprogramming as a driver of stemness and tumorigenesis. Critical Reviews in Oncogenesis. 2014;**19**:363-382

[7] Lu J, Tan M, Cai Q. The Warburg effect in tumor progression: Mitochondrial oxidative metabolism as an anti-metastasis mechanism. Cancer Letters. 2015;**356**:156-164. DOI: 10.1016/j.canlet.2014.04.001

[8] Christofk HR, Vander Heiden MG, Harris MH, Ramanathan A, Gerszten RE, Wei R, et al. The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth. Nature. 2008;**452**:230-233. DOI: 10.1038/nature06734

[9] Birsoy K, Wang T, Chen WW, Freinkman E, Abu-Remaileh M, Sabatini DM. An essential role of the mitochondrial electron transport chain in cell proliferation is to enable aspartate synthesis. Cell. 2015;**162**:540-551. DOI: 10.1016/j.cell.2015.07.016

[10] Pavlides S, Whitaker-Menezes D, Castello-Cros R, Flomenberg N, Witkiewicz AK, Frank PG, et al. The reverse Warburg effect: Aerobic glycolysis in cancer associated fibroblasts and the tumor stroma. Cell Cycle. 2009;**8**:3984-4001. DOI: 10.4161/ cc.8.23.10238

[11] Catteau X, Simon P, Noel JC. Stromal expression of matrix metalloproteinase 2 in cancer-associated fibroblasts is strongly related to human epidermal growth factor receptor 2 status in invasive breast carcinoma. Molecular and Clinical Oncology. 2016;**4**:375-378. DOI: 10.3892/mco.2015.721

[12] Lee M, Yoon JH. Metabolic interplay between glycolysis and mitochondrial oxidation: The reverse Warburg effect and its therapeutic implication. World Journal of Biological Chemistry. 2015;**6**:148-161. DOI: 10.4331/wjbc. v6.i3.148

[13] Bar-Or D, Carrick M, Tanner A 2nd, Lieser MJ, Rael LT, Brody E. Overcoming the Warburg effect: Is it the key to survival in sepsis? Journal of Critical Care. 2018;**43**:197-201. DOI: 10.1016/j.jcrc.2017.09.012

[14] Whitaker-Menezes D, Martinez-Outschoorn UE, Lin Z, Ertel A,

Flomenberg N, Witkiewicz AK, et al. Evidence for a stromal-epithelial "lactate shuttle" in human tumors: MCT4 is a marker of oxidative stress in cancer-associated fibroblasts. Cell Cycle. 2011;**10**:1772-1783. DOI: 10.4161/ cc.10.11.15659

[15] Vyas S, Zaganjor E, Haigis MC. Mitochondria and cancer. Cell. 2016;**166**:555-566. DOI: 10.1016/j. cell.2016.07.002

[16] Williams EG, Wu Y, Ryu D, Kim JY, Lan J, Hasan M, et al. Quantifying and localizing the mitochondrial proteome across five tissues in a mouse population. 2018;**17**:1766-1777. DOI: 10.1074/mcp.RA118.000554

[17] Sassano ML, van Vliet AR, Agostinis P. Mitochondria-associated membranes as networking platforms and regulators of cancer cell fate. Frontiers in Oncology. 2017;**7**:174. DOI: 10.3389/ fonc.2017.00174

[18] Strickertsson JAB, Desler C, Rasmussen LJ. Bacterial infection increases risk of carcinogenesis by targeting mitochondria. Seminars in Cancer Biology. 2017;**47**:95-100. DOI: 10.1016/j.semcancer.2017.07.003

[19] Mintz HA, Yawn DH, Safer B, Bresnick E, Liebelt AG, Blailock ZR, et al. Morphological and biochemical studies of isolated mitochondria from fetal, neonatal, and adult liver and from neoplastic tissues. The Journal of Cell Biology. 1967;**34**:513-523

[20] Chen M, Huang H, He H, Ying W, Liu X, Dai Z, et al. Quantitative proteomic analysis of mitochondria from human ovarian cancer cells and their paclitaxel-resistant sublines. Cancer Science. 2015;**106**:1075-1083. DOI: 10.1111/cas.12710

[21] Scalise M, Pochini L, Galluccio M, Console L, Indiveri C. Glutamine transport and mitochondrial

metabolism in cancer cell growth. Frontiers in Oncology. 2017;**7**:306. DOI: 10.3389/fonc.2017.00306

[22] Huarte M. The emerging role of lncRNAs in cancer. Nature Medicine. 2015;**21**:1253-1261. DOI: 10.1038/ nm.3981

[23] Mou Y, Wang D, Xing R, Nie H, Mou Y, Zhang Y, et al. Identification of long noncoding RNAs biomarkers in patients with hepatitis B virusassociated hepatocellular carcinoma. Cancer Biomarkers. 2018. DOI: 10.3233/ CBM-181424

[24] Rodriguez Bautista R, Ortega GA. Long non-coding RNAs: Implications in targeted diagnoses, prognosis, and improved therapeutic strategies in human non- and triple-negative breast cancer. 2018;**10**:88. DOI: 10.1186/ s13148-018-0514-z

[25] Chen B, Li Y, He Y, Xue C, Xu F. The emerging roles of long non-coding RNA in gallbladder cancer tumorigenesis. Cancer Biomarkers. 2018;**22**:359-366. DOI: 10.3233/cbm-170979

[26] Chen Z, Wu J, Huang W, Peng J, Ye J, Yang L, et al. Long non-coding RNA RP11-789C1.1 suppresses epithelial to mesenchymal transition in gastric cancer through the RP11- 789C1.1/MiR-5003/E-cadherin axis. Cellular Physiology and Biochemistry. 2018;**47**:2432-2444. DOI: 10.1159/000491617

[27] Xu TP, Wang WY, Ma P, Shuai Y, Zhao K, Wang YF, et al. Upregulation of the long noncoding RNA FOXD2-AS1 promotes carcinogenesis by epigenetically silencing EphB3 through EZH2 and LSD1, and predicts poor prognosis in gastric cancer. Oncogene. 2018. DOI: 10.1038/s41388-018-0308-y

**87**

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer*

Current Cancer Drug Targets. 2018. DOI: 10.2174/1568009618666180709163718

[36] Adamopoulos PG, Tsiakanikas P,

peptidases and associated microRNAs as promising prognostic biomarkers in gastrointestinal malignancies. Biological Chemistry. 2017;**399**:821-836. DOI:

[37] Wang Y, Zhao Z, Zhang S, Li Z, Li D, Yang S, et al. LncRNA FAL1 is a negative prognostic biomarker and exhibits pro-oncogenic function in osteosarcoma. 2018. DOI: 10.1002/

[38] Tian L, He Y, Zhang H, Wu Z, Li D, Zheng C. Comprehensive analysis of differentially expressed profiles of lncRNAs and mRNAs reveals ceRNA networks in the transformation of diffuse large B-cell lymphoma. Oncology Letters. 2018;**16**:882-890.

Scorilas A. Kallikrein-related

10.1515/hsz-2017-0342

DOI: 10.3892/ol.2018.8722

[39] Wang LN, Tong SW, Hu HD, et al. Quantitative proteome analysis of ovarian cancer tissues using a iTRAQ approach. Journal of Cellular Biochemistry. 2012;**113**(12):3762-3772

jcb.27074

*DOI: http://dx.doi.org/10.5772/intechopen.80622*

[29] Li N, Zhan X, Zhan X. The lncRNA SNHG3 regulates energy metabolism of ovarian cancer by an analysis of mitochondrial proteomes. Gynecologic

[30] Warburg O. On the origin of cancer

Oncology. 2018. DOI: 10.1016/j.

cells. Science. 1956;**123**:309-314

[31] Kobierzycki C, Piotrowska A, Latkowski K, Zabel M, Nowak-Markwitz E, Spaczynski M, et al. Correlation of pyruvate kinase M2 expression with clinicopathological data in ovarian cancer. Anticancer Research. 2018;**38**:295-300. DOI: 10.21873/

[32] Wang Y, Liu J, Jin X, Zhang D, Li D, Hao F, et al. O-GlcNAcylation destabilizes the active tetrameric PKM2 to promote the Warburg effect. Proceedings of the National Academy of Sciences of the United States of America. 2017;**114**:13732-13737. DOI:

10.1073/pnas.1704145115

[34] Yoshida GJ. Metabolic

10.1186/s13046-015-0221-y

[33] DeHart DN, Lemasters JJ, Maldonado EN. Erastin-like anti-Warburg agents prevent

mitochondrial depolarization induced by free tubulin and decrease lactate formation in cancer cells. SLAS Discovery. 2018;**23**:23-33. DOI: 10.1177/2472555217731556

reprogramming: The emerging concept and associated therapeutic strategies. Journal of Experimental & Clinical Cancer Research. 2015;**34**:111. DOI:

[35] Liu P, Tang H, Wu J, Qiu X, Kong Y, Zhang L, et al. Linc01638 promotes tumorigenesis in HER2+ breast cancer.

ygyno.2018.06.013

anticanres.12221

epithelial-mesenchymal transition and the Warburg effect in aggressive tumors. Critical Reviews in Oncology/ Hematology. 2018;**126**:112-120. DOI: 10.1016/j.critrevonc.2018.03.028

[28] Hua Q , Mi B, Huang G. The emerging co-regulatory role of long noncoding RNAs in

*Energy Metabolism Heterogeneity-Based Molecular Biomarkers for Ovarian Cancer DOI: http://dx.doi.org/10.5772/intechopen.80622*

epithelial-mesenchymal transition and the Warburg effect in aggressive tumors. Critical Reviews in Oncology/ Hematology. 2018;**126**:112-120. DOI: 10.1016/j.critrevonc.2018.03.028

*Molecular Medicine*

cc.10.11.15659

cell.2016.07.002

fonc.2017.00174

Flomenberg N, Witkiewicz AK, et al. Evidence for a stromal-epithelial "lactate shuttle" in human tumors: MCT4 is a marker of oxidative stress in cancer-associated fibroblasts. Cell Cycle. 2011;**10**:1772-1783. DOI: 10.4161/ metabolism in cancer cell growth. Frontiers in Oncology. 2017;**7**:306. DOI:

[22] Huarte M. The emerging role of lncRNAs in cancer. Nature Medicine. 2015;**21**:1253-1261. DOI: 10.1038/

[23] Mou Y, Wang D, Xing R, Nie H, Mou Y, Zhang Y, et al. Identification of long noncoding RNAs biomarkers in patients with hepatitis B virusassociated hepatocellular carcinoma. Cancer Biomarkers. 2018. DOI: 10.3233/

[24] Rodriguez Bautista R, Ortega GA. Long non-coding RNAs: Implications in targeted diagnoses, prognosis, and improved therapeutic strategies in human non- and triple-negative breast cancer. 2018;**10**:88. DOI: 10.1186/

[25] Chen B, Li Y, He Y, Xue C, Xu F. The emerging roles of long non-coding RNA in gallbladder cancer tumorigenesis. Cancer Biomarkers. 2018;**22**:359-366.

[26] Chen Z, Wu J, Huang W, Peng J, Ye J, Yang L, et al. Long non-coding RNA RP11-789C1.1 suppresses epithelial to mesenchymal transition in gastric cancer through the RP11- 789C1.1/MiR-5003/E-cadherin axis. Cellular Physiology and

Biochemistry. 2018;**47**:2432-2444. DOI:

[27] Xu TP, Wang WY, Ma P, Shuai Y, Zhao K, Wang YF, et al. Upregulation of the long noncoding RNA FOXD2-AS1

epigenetically silencing EphB3 through EZH2 and LSD1, and predicts poor prognosis in gastric cancer. Oncogene. 2018. DOI: 10.1038/s41388-018-0308-y

promotes carcinogenesis by

[28] Hua Q , Mi B, Huang G. The emerging co-regulatory role of long noncoding RNAs in

10.3389/fonc.2017.00306

nm.3981

CBM-181424

s13148-018-0514-z

10.1159/000491617

DOI: 10.3233/cbm-170979

[15] Vyas S, Zaganjor E, Haigis MC. Mitochondria and cancer. Cell. 2016;**166**:555-566. DOI: 10.1016/j.

[16] Williams EG, Wu Y, Ryu D, Kim JY, Lan J, Hasan M, et al. Quantifying and localizing the mitochondrial proteome across five tissues in a mouse population. 2018;**17**:1766-1777. DOI:

[17] Sassano ML, van Vliet AR, Agostinis P. Mitochondria-associated membranes as networking platforms and regulators

10.1074/mcp.RA118.000554

of cancer cell fate. Frontiers in Oncology. 2017;**7**:174. DOI: 10.3389/

[18] Strickertsson JAB, Desler C, Rasmussen LJ. Bacterial infection increases risk of carcinogenesis by targeting mitochondria. Seminars in Cancer Biology. 2017;**47**:95-100. DOI: 10.1016/j.semcancer.2017.07.003

[19] Mintz HA, Yawn DH, Safer B, Bresnick E, Liebelt AG, Blailock ZR, et al. Morphological and biochemical studies of isolated mitochondria from fetal, neonatal, and adult liver and from neoplastic tissues. The Journal of Cell

[20] Chen M, Huang H, He H, Ying W, Liu X, Dai Z, et al. Quantitative proteomic analysis of mitochondria from human ovarian cancer cells and their paclitaxel-resistant sublines. Cancer Science. 2015;**106**:1075-1083.

[21] Scalise M, Pochini L, Galluccio M, Console L, Indiveri C. Glutamine transport and mitochondrial

Biology. 1967;**34**:513-523

DOI: 10.1111/cas.12710

**86**

[29] Li N, Zhan X, Zhan X. The lncRNA SNHG3 regulates energy metabolism of ovarian cancer by an analysis of mitochondrial proteomes. Gynecologic Oncology. 2018. DOI: 10.1016/j. ygyno.2018.06.013

[30] Warburg O. On the origin of cancer cells. Science. 1956;**123**:309-314

[31] Kobierzycki C, Piotrowska A, Latkowski K, Zabel M, Nowak-Markwitz E, Spaczynski M, et al. Correlation of pyruvate kinase M2 expression with clinicopathological data in ovarian cancer. Anticancer Research. 2018;**38**:295-300. DOI: 10.21873/ anticanres.12221

[32] Wang Y, Liu J, Jin X, Zhang D, Li D, Hao F, et al. O-GlcNAcylation destabilizes the active tetrameric PKM2 to promote the Warburg effect. Proceedings of the National Academy of Sciences of the United States of America. 2017;**114**:13732-13737. DOI: 10.1073/pnas.1704145115

[33] DeHart DN, Lemasters JJ, Maldonado EN. Erastin-like anti-Warburg agents prevent mitochondrial depolarization induced by free tubulin and decrease lactate formation in cancer cells. SLAS Discovery. 2018;**23**:23-33. DOI: 10.1177/2472555217731556

[34] Yoshida GJ. Metabolic reprogramming: The emerging concept and associated therapeutic strategies. Journal of Experimental & Clinical Cancer Research. 2015;**34**:111. DOI: 10.1186/s13046-015-0221-y

[35] Liu P, Tang H, Wu J, Qiu X, Kong Y, Zhang L, et al. Linc01638 promotes tumorigenesis in HER2+ breast cancer.

Current Cancer Drug Targets. 2018. DOI: 10.2174/1568009618666180709163718

[36] Adamopoulos PG, Tsiakanikas P, Scorilas A. Kallikrein-related peptidases and associated microRNAs as promising prognostic biomarkers in gastrointestinal malignancies. Biological Chemistry. 2017;**399**:821-836. DOI: 10.1515/hsz-2017-0342

[37] Wang Y, Zhao Z, Zhang S, Li Z, Li D, Yang S, et al. LncRNA FAL1 is a negative prognostic biomarker and exhibits pro-oncogenic function in osteosarcoma. 2018. DOI: 10.1002/ jcb.27074

[38] Tian L, He Y, Zhang H, Wu Z, Li D, Zheng C. Comprehensive analysis of differentially expressed profiles of lncRNAs and mRNAs reveals ceRNA networks in the transformation of diffuse large B-cell lymphoma. Oncology Letters. 2018;**16**:882-890. DOI: 10.3892/ol.2018.8722

[39] Wang LN, Tong SW, Hu HD, et al. Quantitative proteome analysis of ovarian cancer tissues using a iTRAQ approach. Journal of Cellular Biochemistry. 2012;**113**(12):3762-3772

**89**

**1. Introduction**

**Chapter 6**

**Abstract**

A Novel P53/POMC/Gas/

Feedback Loop and Pathologic

*Ding'an Zhou, Jiawei Zeng, Xing Zeng, Yadong Li,*

*Zhixiong Wu, Xin Wan, Pingshen Hu and Xiaodong Su*

P53-regulated proteins in transcriptional level are associated with many signal transduction pathways and p53 plays a pivotal role in a number of positive and negative autoregulatory feedback loops. Although POMC/α-MSH productions induced by ultraviolet (UV) are directly mediated by p53, p53 is related to UV-independent pathological pigmentation. In the process of identifying the causative gene of dyschromatosis universalis hereditaria (DUH), three mutations encoding amino acid substitutions were found in the gene SAM and SH3 domain containing 1 (SASH1). SASH1 was identified to interact with guanine nucleotide-binding protein subunit-alpha isoforms short (Gαs). However, for about 90 years, the pathological gene and the pathological mechanism of DUH are unclear. Our study indicates that SASH1 is physiologically medicated by p53 upon UV stimulation and a reciprocal SASH-p53 inducement is existed physiologically and pathophysiologically. A novel p53/POMC/α-MSH/Gαs/SASH1 signal cascade regulates SASH1 to foster melanogenesis. SASH1 mutations control a novel p53/POMC/Gαs/SASH1 autoregulatory positive feedback loop to promote pathological hyperpigmentation phenotype in DUH-affected individuals. Our work illustrates a novel p53/POMC/Gαs/SASH1 autoregulatory positive feedback loop that is mediated by SASH1 mutations to foster

**Keywords:** SASH1 substitution mutations, p53, autoregulatory feedback loop, DUH,

The skin pigmentation is formed by the synthesis of melanin in the melanocytes. Melanocyte is a kind of epithelial cells mainly locating basal cell layers of epidermis, and a few number of melanocytes are located in mucosa. Pigment granules constituted with melanin can distribute and transport to neighboring keratinocytes [1]. Mutations in melanocortin-1-receptor (MC1R) are pivotal for human skin's tanning and pigmentation. MC1R belongs to a G-protein-coupled receptor (GPCR) that is expressed in epidermal melanocytes in a preferential

SASH1 Autoregulatory

Hyperpigmentation

pathological hyperpigmentation phenotype.

pathological hyperpigmentation phenotype

### **Chapter 6**