- 1Department of Ophthalmology of Shanghai Tenth’s People Hospital and Laboratory of Clinical Visual Science of Tongji Eye Institute, Tongji University, Shanghai, China
- 2Department of Biochemistry and Molecular Biology, School of Medicine, Tongji University, Shanghai, China
- 3Department of General Surgery of Guizhou Provincial People’s Hospital, Guizhou, China
- 4Tongji University School of Medicine, Shanghai Tenth People’s Hospital of Tongji University, Shanghai, China
- 5Department of Pharmacology, School of Medicine, Tongji University, Shanghai, China
- 6Department of Ophthalmology, Shanghai First People’s Hospital, Shanghai Jiao Tong University, Shanghai, China
Hepatocellular carcinoma (HCC) is one of the most common types of cancer. The novel sensitive biomarkers and therapeutic targets are urgently needed for the early diagnosis of HCC and improvement of clinical outcomes. Glia maturation factor-β (GMFB) is a growth and differentiation factor for both glia and neurons and has been found to be tightly involved in inflammation and neurodegeneration conditions. In our study, the expression level of GMFB was significantly up-regulated in patients with HCC and positively co-expression with tumor node metastases (TNM) stage and histopathological grade of HCC. The high expression level of GMFB was remarkably associated with poor overall survival, which mainly occurred in males rather than females. Multivariate analysis revealed GMFB to be an independent prognostic factor for overall survival in patients with HCC. Results of Gene Ontology (GO) and KEGG pathways analysis showed that down-regulation of pathways related to protein translation and mitochondria function were enriched. Protein-protein interaction analysis revealed the central role of mitochondria protein in HCC. The downregulation of genes involved in glycolysis and gluconeogenesis was observed among the co-expression genes of GMFB. Knockdown of GMFB in Hep3B significantly inhibited proliferation, migration, and invasion of Hep3B cells, and also downregulated the expression levels of some of metal matrix proteinase (MMP), increased mtDNA copy number and loss of mitochondrial transmembrane potential. GMFB influences the malignancy rate of HCC possibly through regulation of the expression of MMPs, mtDNA function and glycolysis. We proposed that GMFB was a promising HCC diagnostic and prognostic biomarker and therapeutic target in HCC.
Introduction
Hepatocellular carcinoma (HCC) is the main component of primary liver cancer, accounting for 75%-85%, and remains the second leading cause of the death rate of cancers worldwide (1, 2). The major risk factors of HCC showed regional differences. Among the most important risk factors, chronic infection with hepatitis B virus (HBV) or hepatitis C virus (HCV), heavy alcohol intake, obesity and type 2 diabetes were highlighted (3). Although HCC can be prevented by the vaccine against HBV, the incidence and mortality of HCC have been increasing rapidly for the past several years worldwide and represents a considerable public health burden (4–8).
Previous studies have identified biomarkers, such as Alpha-fetoprotein (AFP), apyrimidinic endodeoxyribonuclease 1 (APEX1), glypican 3 (GPC3), Golgi protein-73 (GP73), and Dickkopf-1 (DKK1), for HCC diagnosis (9–13). Among them, AFP is considered the gold standard diagnostic marker for HCC. However, its sensitivity and specificity are low and its expression can be influenced by several non-HCC related factors (14). Nevertheless, the early and specific diagnosis of HCC remains challenging.
Molecularly targeted therapeutics are important methods for treating advanced HCC (15). To date, a large number of biomarkers have been served as therapeutic targets for HCC. For instance, as a multi-targeted tyrosine kinase receptor inhibitor Sorafenib exerts its function by targeting CRAF、BRAF、VEGFRl-3、PDGFR-p、cKIT、FLT-3 and RET (16). Although remarkable progress has been achieved, there are no robust biomarkers that could be applied in the early diagnosis and treatment of HCC (17). However, a limited effect with a five-year survival rate of 18% in HCC patients was obtained (18, 19). Therefore, there is an urgent need to develop a novel, highly sensitive and specific biomarker for the diagnosis, prognosis and therapeutic target for HCC (20). There is also a remarkable gender disparity of HCC incidence, primarily males (21). The current study provides evidence that androgen in males could play a critical role in hepatocarcinogenesis and development, and estrogen could protect women from hepatocarcinogenesis (22, 23). Therefore, the gender-related differences should be given special consideration in individual HCC treatment from the perspective of precision medicine.
GMFB is a growth and differentiation factor expressed predominantly in the central nervous system (CNS) and testis (24, 25). As for liver, the expression of GMFB was confirmed by previous proteomics analysis (26, 27). GMFB is involved in growth and differentiation in the vertebrate brain, neutrophil chemotaxis, migration of monocyte, migration, and adherence of T lymphocytes (28–31). GMFB was found to be upregulated in several neuroinflammation and neurodegeneration conditions (32). Increased GMFB expression was found highly associated with multiple types of cancer, including glioma and ovarian cancer, and GMFB overexpression was reported to be co-expression with poor prognosis in ovarian cancer and glioma (33, 34). To date, there have been no studies of the effects of GMFB on HCC, therefore its roles in HCC progression remain unclear.
This study aimed to characterize the association between GMFB and HCC, and reveal a potential underlying mechanism of GMFB of HCC. To verify analysis results, a series of experiments were carried out. The expression level of GMFB affects mtDNA copy number, mitochondrial membrane potential and MMPs expression levels in Hep3B cell, and impairs cell migration, invasion and adhesion. To our knowledge, this is the first study to report the diagnostic and prognostic values of GMFB in HCC and analyze the gender-based correlation between GMFB and overall survival in patients with HCC. Our study provided a new potential diagnostic and prognostic marker, and a novel bio-target for the treatment of HCC.
Materials and Methods
Ethics Statement
Our study was approved by the Academic Committee of Tongji University and conducted in accordance with the Declaration of Helsinki. We obtained all the data from the online databases, therefore the informed consent for data collection had already been obtained.
Oncomine Database Analysis
ONCOMINE (http://www.oncomine.org) is an integrated translational bioinformatics platform composed of datasets and analyses (35). Datasets include samples represented as microarray data measuring either mRNA expression or DNA copy number, which is used to set up analyses on groups of interest like cancer versus normal. Oncomine analyses result from computations that are performed on samples within a dataset. To investigate the transcription levels of GMFB in different types of cancers the Oncomine database was used. Student´s t-tests were performed, with results filtered by the cut-off of p-value <0.05, fold change of >1.5 and gene rank in the top 10%.
UALCAN Database Analysis
UALCAN (http://ualcan.path.uab.edu) is a comprehensive web resource based on 3 RNA-seq databases and the clinical data of 31 cancer types from the The Cancer Genome Atlas (TCGA) database (project ID, TCGA-LIHC) (36). To analyze the levels of the ADF family protein in liver hepatocellular carcinoma. We analyze the transcriptional expression patterns in normal liver tissue and liver hepatocellular carcinoma samples through UALCAN.
Human Protein Atlas Analysis
HPA (https://www.Proteinatlas.org/) is a valuable tool provided immunostaining on tissues and cell lines as well as differential expression analysis of proteins in normal and tumor tissues (37). In this study, we checked the protein expression of GMFB in the HPA database and analyzed the immunohistochemical results of GMFB in tumor tissues and normal tissues (Antibody: HPA053669).
cBioPortal
cBioportal (http://www.cbioportal.org/) is an open access powerful tool offered to explore, visualize and analyze multidimensional cancer genomics data. The genomic data include somatic mutations, DNA copy-number alterations (CNAs), mRNA and microRNA (miRNA) expression, DNA methylation, protein abundance and phosphoprotein abundance (38). The Oncoprint module in cBioportal was used to provide a graphic summary of major genetic alterations and changes in gene expression. Increases or decreases in mRNA level were based on a Z-score threshold of 2.0 or more standard deviations from the mean of the reference population. The reference population was the samples that are scored diploid for each gene. In our study, the relationship between GMFB in HCC patients and alteration frequency were analyzed through cBioportal. Venn diagram was created using online Venn software (http://bioinformatics.psb.ugent.be/webtools/Venn/).
LinkedOmics and GSE
The LinkedOmics database (http://www.linkedomics.org) is an open-access online biometrics platform that contains 32 cancer types comes from 11,158 patients from TCGA database (project ID, TCGA-LIHC) (39). In this study, we determined the GMFB associated gene enrichment using the “LinkInterpreter” module, which performs enrichment analysis based on Gene Ontology, biological pathways, network modules, among other functional categories.
Metascape, DAVID and Network Analyst Analysis
Metascape and DAVID are online software for gene annotation and gene set enrichment analysis. In this study, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of GMFB co-expression genes were performed using Metascape and DAVID. And the PPI networks were conducted by Metascape. In Metascape, the min overlap was set as 3, the P-value cutoff was 0.01, and the minimum enrichment factor was 1.5. The cutoff value was set as FDR of<0.05. NetworkAnalyst (http://www.networkanalyst.ca) was used to identify the differential expression of GMFB between HCC and normal samples from TCGA database (project ID, TCGA-LIHC) (40).
Construction of Prognostic Models and Survival Analysis
The Kaplan–Meier Plotter database (http://kmplot.com/analysis/) is an online public database, to draw survival plots using publicly available data (41). In this work, our survival analysis was carried on the Kaplan–Meier Plotter database contains survival information for 364 patients with HCC, 250 male patients and 121 female patients with the “Auto select best cutoff” option. The expression levels of GMFB and clinical features were merged to find independent prognostic factors through univariate and multivariate Cox regression analysis. The receiver operating characteristics curve (ROC) analysis was performed using the ‘survival ROC’ package in R. The area under the ROC curve (AUC) was measured for the prediction of GMFB.
Cell Lines
Hep3B cells purchased from ATCC were used in this study. All cells were maintained in High Glucose Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% FBS (Gibco, Thermo Fisher Scientific, Spain) at 37°C and 5% CO2. and cultured.
Cell Transfection
Hep3B cells were transfected using Lipofectamine 2000 (Invitrogen, USA), according to the manufacturer’s instructions. SiRNA and Lipofectamine 2000 were diluted in FBS free DMEM(Gibco) medium and incubated for 5 min, separately. Then mixed together, and incubated for 10-15 min at room temperature (RT) to form the DNA-Lipofectamine complexes. Finally added to the culture plates.
Western Blot Analysis
The samples were lysed with cell lysis buffer (RIPA Lysis Buffer; CAT: P0013B; Beyotime Institute of Biotechnology) to extract whole-cell protein. Protein concentration was quantified using a BCA kit (CAT: P0012; Beyotime Institute of Biotechnology Jiangsu, China), and 20–50 μg of each protein were separated by SDS-PAGE using a 12.5% SDS-PAGE gel. After transfer to a polyvinylidene fluoride membrane (Millipore, USA), membranes were blocked with 5% skim milk at room temperature for 1 hr and incubated overnight at 4°C with primary antibody against GMFB (1:2000; Cat no:10690-1-AP; Proteintech) and ACTB (1:5000; Cat no: 20536-1-AP; Proteintech). Secondary antibodies used were rabbit anti-goat Ig-HRP (1:5000; Proteintech). The membranes were then incubated with secondary antibodies for 2 hrs. Proteins were visualized using ECL (Millipore). The density of the bands was determined using ImageJ software (USA).
Wound Healing Assay
Hep3B cells were seeded in 6-well plates and grown to 90% confluence. The confluent cell monolayers were scratched using a sterile 200µl pipette tip and rinsed with PBS to remove scratched cells. Wound closure was observed for 36 hours. The wound closure rate was calculated as follows: wound closure (%) = (area of initial wound-area of final wound)/area of initial wound ×100.
Transwell Migration and Invasion Assays
For cell invasion analysis, transwell chambers were coated with Matrigel, and the cell migration (without matrigel) assay was also performed. Cell suspensions in FBS free medium were added to upper transwell chambers (24-well, 8-μm pores; BD Labware, USA), While 600 µL DMEM medium supplemented with 10% FBS was added into lower transwell chambers. Cell counting was performed using the cell counter plugin of ImageJ. The migratory capacity of Hep3B cells under various treatments was evaluated through scratch assay.
Measurement of mtDNA Copy Number
Total DNA (mtDNA and gDNA) was isolated using TIANamp Genomic DNA Kit from TIANGEN (TIANGEN, Beijing, China) according to the manual. Relative mtDNA copy number was quantified by quantitative PCR (qPCR) (42). The quantitative PCR was performed by using SuperReal PreMix Plus (SYBR Green) (TIANGEN, Beijing, China) on a BIO-RAD CFX96 Touch Real-Time PCR Detection System. The mtDNA (D-loop, MT-TL1 and MT-ND1) was normalized to nuclear DNA (NCOA3). The primers and primer sequences of mtDNA-specific primers used were as previously reported (43). The relative fold changes were calculated by 2−ΔΔCt.
Quantitative Real-Time PCR
Total RNA was isolated using the TRIzol reagent (Takara, Dalian, China) strictly according to the instructions. Then RNA was reverse transcribed to cDNA using PrimeScript RT polymerase (Takara, Dalian, China). SYBR Green Master Mix (Tiangen Biotech, China) on a LightCycler 96 Detection System (Roche) was used for RT-qPCR. Data were analyzed using the 2-ΔΔCt method. Oligos were synthesized by Sangon. Primer sequences were obtained from Origene.
Mitochondrial Membrane Potential
Hep3B cells were seeded onto coverslip placed in 24-well plates and cultured at 37°C in a 5% CO2 incubator. After 24 h, cells were transfected with scramble control siRNA or GMFB siRNA with Lipo2000 (Invitrogen). Scramble siRNA was used as a control. JC1 staining was performed 48 h after transfection. Cells were washed three times with PBS, and then stained with the 2 μM JC-1 dye (Beyotime, Shanghai, China) in the dark for 20 minutes. Cell visualization was performed by a fluorescence microscope.
Statistical Analysis
Image analysis was performed using Image J. Student t-test and one-way ANOVA were used in the statistical analysis. Data were expressed as means ± SEM. The p-value less than 0.05 was considered statistically significant.
Results
Increased Expression of GMFB in Liver HCC
We initially used the Oncomine database to analyze the expression profiles, and found that the GMFB transcription levels were markedly higher in breast cancer, head and neck cancer and liver cancer compared to the normal tissues (Figure 1A). To study the changes of GMFB mRNA expression levels between HCC and normal liver tissues, we then used the Oncomine and UALCAN database to analyze the differential. The mining of publicly available databases from Roessler Liver 2 Statistics showed that GMFB was up-regulated in HCC (Figure 1B). Roessler Liver 2 Statistics was containing 12,624 measured mRNAs from 445 samples. Samples were hybridized to Affymetrix Human Genome HT U133A Array. GEO sample IDs were given in Supplementary Table 1. Then, we explored the expression of GMFB in HCC using the UALCAN based on the data resources of The Cancer Genome Atlas database. The results were nearly identical (Figure 1C). Furthermore, the analysis of TCGA datasets revealed that the expression of GMFB in HCC was associated with TNM stage 1-3 (Figure 1D) and pathological grade 1-3 (Figure 1E). Considering the small number of samples of stage 4 and grade 4, the statistical conclusion may not be accurate, and future larger cohort validation is needed. The analysis of TCGA datasets revealed that GMFB was highly expressed in both male and female HCC patients using the UALCAN. And there was no difference was found between male and female HCC patients (Figure 1F). To solve the imbalance among normal, male and female data in Figure 1F, we downloaded normal (liver), male (HCC) and female (HCC) gene expression data from TCGA further validation. Each sample was assigned a number and 30 samples of each group were randomly selected (Supplementary Table 2). The randomization was computer-based, generated in R (version R 3.4.4). The result indicated that there was a significantly higher expression of GMFB in HCC patients in both genders, and no significant gender differences were found in normal (liver) or HCC tissues (Figure 1G). Furthermore, Detection of GMFB protein expression in normal and HCC liver tissue of females and males by immunohistochemistry retrieved from the Human Protein Atlas also indicated higher protein expression patterns of GMFB in HCC tissues when compared to normal samples (Figure 2A). Furthermore, we also performed GMFB copy number variation (CNV) and mRNA levels among normal and HCC tumor samples using Oncoprint algorithm from cBioPortal, and the result showed that 5% of the HCC cases had undergone genetic changes (Figure 2B).
Figure 1 GMFB expression levels in tumors based on sample types, tumor stage, tumor grade. (A) The transcription levels of GMFB in different types of cancers from the Oncomine database. The schematic in each cell reveals the numbers of datasets with statistically significant. The cell color is determined by the best gene rank percentile, red indicating high expression and blue indicating low expression. (B) The expression of GMFB in HCC and normal liver tissues from the Oncomine database. (C) The transcription levels of GMFB in HCC and normal tissues were analyzed by the UALCAN cancer database. The expression of GMFB in HCC is based on tumor stage (D), tumor grade (E), and gender (F) from the UALCAN database. (G) The transcription levels of GMFB in male normal, female normal, male HCC and female HCC tissues were analyzed by the NetworkAnalyst. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, ns, no significance.
Figure 2 Increased expression of GMFB in liver hepatocellular carcinoma. (A) Immunohistochemistry analysis of GMFB protein expression in normal liver tissues and hepatocellular carcinoma tissues of males and females in HPA. (B) OncoPrint of GMFB genetic alterations in HCC shows the frequency of genetic alterations included amplification, mRNA upregulation, and mRNA downregulation. The aberrant expression threshold was defined as z−score ±2.0.
The Significant Correlation Between Increased GMFB Expression and Poor Overall Survival in HCC
To clarify the associations between the expression levels of GMFB and patients’ clinical outcomes, Kaplan-Meier (KM)Plotter survival analysis was performed. The survival curves indicated that high GMFB expression was significantly associated with poor overall survival (OS) of HCC patients (P=0.0042) (Figure 3A). Next, we analyzed the effects of gender differences on the association between GMFB and OS. A significant difference was observed in male patients (Figure 3B), but not in females (P=0.36) (Figure 3C).
Figure 3 Kaplan-Meier curves showed high mRNA expression levels of GMFB was associated with unfavorable overall survival (A) High expression of GMFB was associated with poor survive in HCC patients (KM plot). Kaplan–Meier survival estimation for male (B) and female (C) HCC patients (KM plot). OS, Overall survival. (D) The receiver operating characteristics (ROC) curves for prediction of GMFB expression in patients with hepatocellular carcinomas vs. controls. ROC curves of the radiomics signature, clinical model, and combined model in HCC cohorts (E) Univariate and multivariate Cox analysis of the overall survival (OS) by proportional hazard analysis for cancer-specific survival in HCC patients. The red circle indicated the HR, and the transverse lines indicate the 95% CI.
GMFB Being an Independent Risk Factor for HCC
The AUC values with 95% CI of GMFB was 0.852, CI in a range from 0.809 to 0.895, indicating a significantly diagnostic accuracy for HCC (Figure 3D). Univariate analysis by Cox proportional hazard model showed that GMFB, pathologic stage and tumor status were responsible for OS in HCC. Furthermore, the multivariate Cox regression analysis confirmed that GMFB (P=0.039), pathologic stage(P<0.001) and tumor status(P=0.003) were independent factors for predicting the prognosis of HCC patients (Figure 3E). The clinical characteristics of the HCC patients were shown in Supplementary Table 3.
The Analysis of Co-Expression Genes With GMFB in HCC Patients
To deep mine the underlying mechanisms of the role of GMFB in HCC, LinkedOmics was utilized to analyze mRNA sequencing information from HCC patients in the TCGA. The plot showed the distribution of significant positive correlation with GMFB (red dots) or significant negative correlation (blue dots) in HCC patients (Figure 4A). Heat maps of the top 50 most positively and negatively co-regulated with GMFB in HCC patients were shown in Figures 4B, C. There were 3389 up-regulated genes and 1397 down-regulated genes in high expression level of GMFB in HCC patients (FDR (Benjamini & Hochberg (BH)-corrected p-value) <0.0001, P<0.0001) (Supplementary Table 4). GO (BP, MF, and CC) and KEGG enrichment analyses of co-expression genes were performed. As shown in Figure 4D, we found that the down-regulated co-expression genes of GMFB were associated with the mitochondrial electron transport chain and translation initiation, and the up-regulated co-expression genes of GMFB were relevant to transcription, DNA-template (BP), nucleoplasm (CC), protein binding (MF) and pathway in cancer (KEGG). The detailed information for the co-expression genes was obtained in Supplementary Tables 5, 6. Furthermore, results of Kaplan Meier survival analysis (males and females were shown separately) and reports associated with HCC of TOP100 DEGs in overlap, male-specific and female-specific groups were summarized in Supplementary Table 7.
Figure 4 Significant GO terms and KEGG pathways correlated with GMFB in HCC. (A) A Pearson test was used to identify the correlations between GMFB differentially expressed in HCC cohort (LinkedOmics). Red dots represent positively significantly correlated genes with GMFB and green dots represent negatively significantly correlated genes with GMFB. Heat maps showing genes positively (B) and negatively (C) correlated with GMFB in all HCC (TOP 50). (D) Enriched GO terms (BP, CC and MF) and KEGG pathways of upregulated DEGs and downregulated DEGs. Red indicated up-regulation, and green indicated down-regulated. Black circular nodes(count) represent the number of DEGs enriched in each term.
Functional Enrichment Analysis of GMFB in Patients With HCC Based on Gender
To elucidate the possible mechanism of underling GMFB expression level on the clinical outcome in HCC patients with different gender, male and female common, male-specific and female-specific GMFB co-expression genes were further analyzed. Venn diagram of GMFB co-expression genes (FDR (BH-corrected p-value) <0.01, P<0.01) in male and female HCC groups was illustrated that there were a larger number of regulated genes in males (6328 DEGs) than in females (2899 DEGs), and showing a substantial overlap (2513 DEGs) between males and females (Figure 5A). DAVID analysis software was utilized to identified functional enrichment terms, including GO (BP, MF, and CC) and KEGG of GMFB co-expression genes in overlap (Figures 5B, C), male-specific (Figures 5D, E) and female-specific (Figures 5F, G), respectively. The enrichment analysis disclosed the distinct role of GMFB co-expression genes in female and male HCC patients, despite increased expression of GMFB was observed in both male and female HCC patients. Results also showed a high expression level of GMFB related to mitochondrial dynamics. To further capture the relationships between the terms, PPI network and MCODE component analysis of top 500 DEGs of overlap (Figures 6A, B), male-specific (Figures 6C, D) and female-only (Figures 6E, F) were performed.
Figure 5 Sex differences in significant GO terms and KEGG pathways analysis. (A) Venn diagrams show the number of distinct and overlapping DEGs in male and female HCC patient group. (B) Enriched BP, CC and MF GO terms of 2513 overlap co-expression genes. Vertical axis represents GO terms. Fold enrichment was marked on the top of column. (C). Significant KEGG pathways of 2513 overlap co-expression genes, P values are marked in different colors. (D) Enriched BP, CC and MF GO terms in 3815 male-specific co-expression genes. (E) Significant KEGG pathways of 3815 male-specific co-expression genes. (F) Enriched BP, CC and MF GO terms in 386 female-specific co-expression genes. (G) Significant KEGG pathways of 386 female-specific co-expression genes.
Figure 6 Sex differences in the enrichment analysis of GMFB in HCC (Metascape). Network of enriched terms of top 500 overlap co-expression genes (A), male-specific co-expression genes (C) and female-specific co-expression genes (E), and colored by cluster ID. PPI network of proteins of top 500 overlap co-expression genes (B), male-specific co-expression genes (D) and female-specific co-expression genes (F).
Regardless of gender, co-expression genes of GMFB were mainly involved in oxidative phosphorylation and mitochondria function (Figures 6A, B). In male patients, peroxisome, cell cycle and inosine monophosphate biosynthesis, PRPP + glutamine => IMP were enriched from male-specific co-expression genes of GMFB (Figures 6C, D) while in female patients, KEGG pathways such as ribosome and ubiquitin mediated proteolysis were enriched (Figures 6E, F). Detailed information was available in Supplementary Table 8. Likewise, regardless of gender, PPI networks showed that GO annotation mainly related to MCODE 1 were mitochondrial complex and transferase complex. In males, GO annotation is mainly related to MCODE 1 of mRNA splicing, via spliceosome, RNA splicing, via transesterification reactions with bulged adenosine as nucleophile and RNA splicing, via transesterification reactions, while in females, MCODE 1 was anchoring translation, peptide biosynthetic process and ribosomal subunit. Supplementary Material reports detailed data (Supplementary Table 9).
The Increased mtDNA Copy Number and Loss of Mitochondrial Membrane Potential in Hep3B by Knocking Down GMFB
Based on the bioinformatic analysis above, mitochondria were identified to be a major target of GMFB and GMFB co-expression genes. Hence, the expression levels of mitochondrial expansion program and homeostasis-associated genes were analyzed. Results revealed that 52 genes required for mitochondrial ribosome function were down-regulated (Table 1). Mitochondrial dynamics-associated genes (MTP18, FIS1), mitochondrial protein import-associated genes (TOMM40, TOMM7), mitophagy-associated genes (PINK1) and glucose metabolism-associated genes (FBP1, TPI1 and PCK1) were significantly down-regulated in HCC patients with high GMFB expression (Table 1). Finally, we designed siRNA oligos against GMFB to verify the finding from the bioinformatic analysis. After the determination of the knockdown efficiency of siRNA against GMFB (Figures 7A, B), mtDNA copy number was measured by qPCR. We found that the expression level of GMFB negatively co-expression to mtDNA copy number in Hep3B cells (Figures 7C–E). Furthermore, we also detected mitochondrial membrane potential by JC-1 assay. As shown in Figure 7F, the transition of red to green fluorescent signal was observed in GMFB knockdown HCC cells, which indicated the loss of mitochondrial membrane potential is associated with low expression levels of GMFB in Hep3B cells.
Table 1 Relative expression of mitochondrial and glycolysis-associated genes in patients with high GMFB expression.
Figure 7 Knock down GMFB decreased mtDNA copy number and inhibits proliferation and invasion of Hep3B cell. (A) SiRNA efficiently knock down GMFB mRNA levels by qRT-PCR in Hep3B cells. (B) SiRNA efficiently knock down GMFB protein levels by western blot in Hep3B cells. (C–E) Relative mtDNA copy number (mtDNA/nDNA) was detected by qPCR and normalized by nuclear receptor coactivator 3(NCOA3) gene in knockdown GMFB. GMFB knockdown. (F) Fluorescence microscopic imaging of Hep3B cells stained with JC-1. Scale bar, 50 μm. (G) Image of the scratch assay. Scale bar, 100 μm. (H) Quantification of the scratch assay. Error bars represent standard error of means (SEM). (I) Transwell migration and invasion assays were used to examine the migration and invasion ability of HCC cells. Scale bar, 100 μm. Quantification of migration (J) and invasion assays (K). Error bars represent standard error of means (SEM). (L) Heat map of relative qPCR results of MMPs. **P < 0.01, ***P < 0.001, ****P < 0.0001.
The Inhibition of Hep3B Cell Proliferation, Migration and Invasion by GMFB Knockdown
To further verify the relationship between GMFB and the malignant behavior of the HCC, we also determined cell proliferation, invasive and migratory abilities by wound-healing and Transwell assays. Post-transfection 48 hours of siGMFB oligos in Hep3B cells, we found that decreased expression of GMFB inhibited wound closure significantly (Figures 7G, H). Transwell migration assay showed less motile with crystal violet staining in GMFB knockdown Hep3B cells (Figures 7I, J). Transwell-matrigel invasion assays showed that GMFB knockdown inhibited Hep3B invasion (Figures 7I, K). As shown in Table 2, MMPs expression patterns showed a gender disparity. There were more up-regulated MMPs appeared in the list of co-expression genes of GMFB in male HCC patients than in females. In females, the expression of MMP10, MMP12, MMP16, MMP23A were up-regulated while MMP19 and MMP24 were downregulated. In males, MMP1, MMP2, MMP3, MMP7, MMP8, MMP9, MMP10, MMP11, MMP12, MMP14, MMP16, MMP27 were up-regulated while MMP19 and MMP24 were downregulated. A validation experiment was conducted to examine MMPs expression after knockdown GMFB in Hep3B. qRT-PCR result showed that the expression levels of MMPs (except for MMP7, MMP8, MMP14, MMP17and MMP24) were decreased significantly after GMFB knockdown in Hep3B cells. Compared with MMPs expression patterns in male HCC patients, only MMP7, MMP8 and MMP14 expression levels were not consistent with bioinformatic analysis. Compared with MMPs expression patterns in female HCC patients, MMP12, MMP16, MMP23A and MMP24 expression levels were consistent with bioinformatic analysis. Taken together, our results supported that knockdown of GMFB suppressed the proliferation, invasion and migration of Hep3B cells, probably due to modulating MMPs expression.
Discussion
In the present study, we performed a comprehensive bioinformatic analysis on the association between GMFB and HCC. We found that GMFB was highly expressed in several types of cancer including HCC (Figures 1B, C, 2A), and confirmed the expression level of GMFB was closely associated with the TNM stage, clinical stage, and poor survival rates in HCC (Figures 1D, E). GMFB was considered to be a novel potential biomarker and therapeutic target HCC (Figure 3). Moreover, knockdown of GMFB with siRNA effectively inhibited the proliferation, migration and invasion in Hep3B cells and down-regulated the expression level of some of the MMPs. Targeting GMFB may represent a promising therapeutic strategy for HCC. However, the potential mechanisms underlying GMFB contributing to HCC remain unclear.
To clarify the potential mechanism, we found a large number of DEGs were co-expressed with the high level of GMFB, and most of them were highly associated with mitochondrial functions (Figure 4D). Thus, we compared transcriptional profiles of DEGs related to mitochondrial protein synthesis (MRPLs and, MRPSs), Mitochondrial dynamics (MTP18 and FIS1), Mitochondrial protein import (TOMM40 and TOMM7), mitophagy (PINK1) and glucose metabolism (FBP1, TPI1 and PCK1) (44). The result indicated that a high level of GMFB down-regulated 53 genes required for mitochondrial ribosome function (Table 1), suggesting that the high level of GMFB may suppress mtDNA replication, mitophagy, and energy metabolism in HCC.
Increased expression of GMFB rather than AFP was an independent risk factor for HCC (Figure 3E). Moreover, knockdown of GMFB with siRNA effectively inhibited the proliferation, migration and invasion in Hep3B cells and down-regulated the expression level of some of the MMPs.
The mtDNA copy number gained more and more attention in cancer research. Previous studies showed decreased copy numbers in cancer samples in HCC, bladder cancer, breast cancer, kidney clear cell carcinoma and myeloproliferative neoplasm, and increased copy number detected in chronic lymphocytic leukemia, lung squamous cell carcinoma and pancreatic adenocarcinoma (45–50). Furthermore, the decrease in mtDNA copy number was more significant in female HCC patients than in males (51). In our study, we found that the decreased expression of GMFB leads to up-regulation of mtDNA copy number, which means the higher level of GMFB negatively regulates the amount of mtDNA copy number (Figures 3C–E). Moreover, our result also suggested that low expression levels of GMFB induced loss of mitochondrial membrane potential in Hep3b cells (Figure 7F). Previous works showed that mitochondria is a central regulator of the decision between cellular survival and demise (52). Mitochondrial membrane potential is an indicator of mitochondrial function, and loss of mitochondrial membrane potential resulting in the intrinsic apoptotic pathway (53). To date, there is no study has explored the role of GMFB in mitochondrial dysfunction. which might partially contribute to the malignancy of HCC. Therefore, these results suggested that GMFB may regulate mtDNA function and consequently be involved in HCC tumorigenesis and progression.
It is well known that cancer cells are commonly exposed to nutritional deficiency and hypoxia, and these factors adversely impact cancer cell metastasis (54). Cancer cells use aerobic glycolysis to support proliferation instead of mitochondrial oxidative phosphorylation (55). Therefore, inhibition glycolysis in cancer cells represents a kind of therapeutic strategy (56). In our study, GMFB down-regulated the expression level of the negative regulator of glycolysis, fructose-1,6-bisphosphatase (FBP1), Triosephosphate isomerase (TPI1) and phosphoenolpyruvate carboxykinase 1 (PCK1) (Table 1). This suggests that high expression levels of GMFB co-expression with the promotion of glycolysis and suppression of gluconeogenesis. Furthermore, previous works showed that suppressed FBP1, a negative regulator of aerobic glycolysis, led to promoting HCC growth and metastasis (57). Knockout PCK1, a step limiting enzyme of gluconeogenesis, markedly enhanced the global O-GlcNAcylation levels, which is an emerging hallmark of HCC (58). TPI1 functioned as a tumor suppressor in HCC (59). In conjunction with down-regulated KEGG pathways from co-expression genes of GMFB, we favored that targeting GMFB/mtDNA/glycolysis may represent a novel therapeutic strategy.
According to epidemiological and clinical characteristics of HCC, there is a remarkable gender disparity of HCC incidence, primarily males. Liver cancer ranks fifth in terms of global cases and second in terms of deaths for males. Women appeared to have better survival rates than men in HCC (60). In 2018, the distribution of cases in both sexes, males, and females were 6th, 5th and 9th, and for deaths is 4th, 2nd and 6th respectively (1). When we did KM Plotter survival analysis, we found that high GMFB expression was significantly associated with poor OS of HCC patients (P=0.0042) (Figure 3A). Unexpectedly, there was a significant OS difference between male patients with low and high expression of GMFB (P=0.00014) (Figure 3B). However, it did not occur in females with HCC (P=0.36) (Figure 3C).
It is widely recognized that the liver is one of the main responsible organs for estrogens (61, 62). Previous studies indicate that estrogen and estrogen receptors played a role in the development and progression of HCC (62–64). But, the role of estrogen and estrogen receptors of HCC is disputed (65). It seems that increased estrogen synthesis and variants estrogen receptor in liver leads to an increased risk of HCC (66, 67). On the contrary, estrogen and estrogen receptors can decrease the malignancy of HCC by arresting cell cycle progression and promoting apoptosis (68). In addition, distinct gender differences were observed in our study. GO and KEGG pathways enrichment analyses of DEGs (Figures 5, 6). Further analysis showed that more DEGs co-expressed with GMFB were regulated in the males than in the females (Figure 5A). Results of survival analysis and current research of top 100 DEGs were shown in Supplementary Table 7. These distinct results above indicated that sex hormones might play a key role in regulating the gene expression profiling of co-expression genes of GMFB. Among top 100DEGs, some of them with statistical significance of survival analysis, such as Dodecenoyl-CoA isomerase (DCI, p value=4.3e-5 in males), G2/M phase-specific E3 ubiquitin-protein ligase (G2E3, p value=9.5e-6 in males), LSM domain-containing protein 1 (LSMD1, p value=0.0001 in males), Neural precursor cell expressed developmentally down-regulated protein 1 (NEDD1, p value=0.0001 in males) and Ras-related GTP-binding protein C (RRAGC, p value=0.0001 in females) have not been reported in HCC yet. These genes merit further investigation.
It is well known that the MMPs play a vital role in cancer invasion and metastasis (69). We also found that expression patterns of co-expression MMPs genes of GMFB showed a gender disparity (Table 2), which may explain the significant correlation between OS and GMFB in male HCC, not in females. Knockdown of GMFB effectively inhibited the expression levels of MMPs. Our results supported that the expression level of GMFB in male HCC modulated a subset of MMPs expression, finally contributing to male OS.
Certain limitations exist in our study. The first limitation is that we lack validation of our novel findings with human HCC samples. Since we could not obtain the levels of sex hormones, the expression level of sex hormone receptors in HCC patients used in this study, the relationship between levels of sex hormones, sex hormones reporters and GMFB expression merits further investigation. The second limitation is that mechanistic study is preliminary and should be deeply explored.
In summary, this work will shed deep light onto GMFB function in HCC. Targeting GMFB may represent a promising therapeutic strategy for HCC patients. GMFB may be an invaluable diagnostic and prognostic biomarker for HCC from the perspective of precision medicine.
Data Availability Statement
Publicly available datasets analyzed in this study can be found here: http://www.oncomine.org; http://ualcan.path.uab.edu; https://www.Proteinatlas.org; http://www.cbioportal.org/; http://www.linkedomics.org; http://kmplot.com/analysis/. The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author Contributions
Among the authors in the list, WS designed the study, performed experiments and wrote the manuscript. CH performed part of data analysis and interpreted the experimental data. TW performed part of data analysis. JW, JPZ, FG, QO, HT, CJ, JX, and JFZ discussed the results. LL and G-TX encouraged us to investigate and supervised the findings of this work. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants obtained from the Ministry of Science and Technology of China (2017YFA0104100), and the National Natural Science Foundation (81670867, 81372071, 81770942).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank all members of our laboratory for many helpful discussions.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.744331/full#supplementary-material
References
1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2018) 68(6):394–424. doi: 10.3322/caac.21492
2. Stavraka C, Rush H, Ross P. Combined Hepatocellular Cholangiocarcinoma (cHCC-CC): An Update of Genetics, Molecular Biology, and Therapeutic Interventions. J Hepatocell Carcinoma (2019) 6:11–21. doi: 10.2147/JHC.S159805
3. Wang XY, Huang JM, Lu XM, Harrison TJ, Liu HB, Jia HH, et al. Changing Risk Factors for Hepatocellular Carcinoma in Hyperendemic Regions in the Era of Universal Hepatitis B Vaccination. Cancer Epidemiol (2020) 67:101775. doi: 10.1016/j.canep.2020.101775
4. Parkin DM, Bray F, Ferlay J, Pisani P. Global Cancer Statistics, 2002. CA Cancer J Clin (2005) 55(2):74–108. doi: 10.3322/canjclin.55.2.74
5. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global Cancer Statistics, 2012. CA Cancer J Clin (2015) 65(2):87–108. doi: 10.3322/caac.21262
6. Parkin DM. Global Cancer Statistics in the Year 2000. Lancet Oncol (2001) 2(9):533–43. doi: 10.1016/S1470-2045(01)00486-7
7. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global Cancer Statistics. CA Cancer J Clin (2011) 61(2):69–90. doi: 10.3322/caac.20107
8. Byass P. The Global Burden of Liver Disease: A Challenge for Methods and for Public Health. BMC Med (2014) 12:159. doi: 10.1186/s12916-014-0159-5
9. Cao L, Cheng H, Jiang Q, Li H, Wu Z. APEX1 is a Novel Diagnostic and Prognostic Biomarker for Hepatocellular Carcinoma. Aging (Albany NY) (2020) 12(5):4573–91. doi: 10.18632/aging.102913
10. Muscari F, Maulat C. Preoperative Alpha-Fetoprotein (AFP) in Hepatocellular Carcinoma (HCC): Is This 50-Year Biomarker Still Up-to-Date? Transl Gastroenterol Hepatol (2020) 5:46. doi: 10.21037/tgh.2019.12.09
11. Zhou F, Shang W, Yu X, Tian J. Glypican-3: A Promising Biomarker for Hepatocellular Carcinoma Diagnosis and Treatment. Med Res Rev (2018) 38(2):741–67. doi: 10.1002/med.21455
12. ElZefzafy WM, Hussien M, Mohmmed ZAZ, Abd Elbaky NM. The Diagnostic Value of Golgi Protien-73 and DICKKOPF-1 in Hepatocellular Carcinoma. J Immunoassay Immunochem (2021) 42(2):174–87. doi: 10.1080/15321819.2020.1844750
13. Ali OM, El Amin HA, Sharkawy YL, Mohamed Ali AA, Kholef EFM, Elsewify WAE. Golgi Protein 73 Versus Alpha-Fetoprotein as a New Biomarker in Early Diagnosis of Hepatocellular Carcinoma. Int J Gen Med (2020) 13:193–200. doi: 10.2147/IJGM.S253622
14. Gopal P, Yopp AC, Waljee AK, Chiang J, Nehra M, Kandunoori P, et al. Factors That Affect Accuracy of Alpha-Fetoprotein Test in Detection of Hepatocellular Carcinoma in Patients With Cirrhosis. Clin Gastroenterol Hepatol (2014) 12(5):870–7. doi: 10.1016/j.cgh.2013.09.053
15. Scott JM, Lakoski S, Mackey JR, Douglas PS, Haykowsky MJ, Jones LW. The Potential Role of Aerobic Exercise to Modulate Cardiotoxicity of Molecularly Targeted Cancer Therapeutics. Oncologist (2013) 18(2):221–31. doi: 10.1634/theoncologist.2012-0226
16. Llovet JM, Ricci S, Mazzaferro V, Hilgard P, Gane E, Blanc JF, et al. Sorafenib in Advanced Hepatocellular Carcinoma. N Engl J Med (2008) 359(4):378–90. doi: 10.1056/NEJMoa0708857
17. Dimitroulis D, Damaskos C, Valsami S, Davakis S, Garmpis N, Spartalis E, et al. From Diagnosis to Treatment of Hepatocellular Carcinoma: An Epidemic Problem for Both Developed and Developing World. World J Gastroenterol (2017) 23(29):5282–94. doi: 10.3748/wjg.v23.i29.5282
18. Jemal A, Ward EM, Johnson CJ, Cronin KA, Ma J, Ryerson B, et al. Annual Report to the Nation on the Status of Cancer, 1975-2014, Featuring Survival. J Natl Cancer Inst (2017) 109(9):djx030. doi: 10.1093/jnci/djx030
19. El-Serag HB, Rudolph KL. Hepatocellular Carcinoma: Epidemiology and Molecular Carcinogenesis. Gastroenterology (2007) 132(7):2557–76. doi: 10.1053/j.gastro.2007.04.061
20. Hartke J, Johnson M, Ghabril M. The Diagnosis and Treatment of Hepatocellular Carcinoma. Semin Diagn Pathol (2017) 34(2):153–9. doi: 10.1053/j.semdp.2016.12.011
21. Perz JF, Armstrong GL, Farrington LA, Hutin YJ, Bell BP. The Contributions of Hepatitis B Virus and Hepatitis C Virus Infections to Cirrhosis and Primary Liver Cancer Worldwide. J Hepatol (2006) 45(4):529–38. doi: 10.1016/j.jhep.2006.05.013
22. Yu MW, Chang HC, Chang SC, Liaw YF, Lin SM, Liu CJ, et al. Role of Reproductive Factors in Hepatocellular Carcinoma: Impact on Hepatitis B- and C-Related Risk. Hepatology (2003) 38(6):1393–400. doi: 10.1016/j.hep.2003.09.041
23. Kanda T, Jiang X, Yokosuka O. Androgen Receptor Signaling in Hepatocellular Carcinoma and Pancreatic Cancers. World J Gastroenterol (2014) 20(28):9229–36. doi: 10.3748/wjg.v20.i28.9229
24. Bourgeois F, Guimiot F, Mas C, Bulfone A, Levacher B, Moalic JM, et al. Identification and Isolation of a Full-Length Clone of Mouse GMFB (Gmfb), a Putative Intracellular Kinase Regulator, Differentially Expressed in Telencephalon. Cytogenet Cell Genet (2001) 92(3-4):304–9. doi: 10.1159/000056919
25. Kaplan R, Zaheer A, Jaye M, Lim R. Molecular Cloning and Expression of Biologically Active Human Glia Maturation Factor-Beta. J Neurochem (1991) 57(2):483–90. doi: 10.1111/j.1471-4159.1991.tb03777.x
26. Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-Based Map of the Human Proteome. Science (2015) 347(6220):1260419. doi: 10.1126/science.1260419
27. Kampf C, Mardinoglu A, Fagerberg L, Hallstrom BM, Edlund K, Lundberg E, et al. The Human Liver-Specific Proteome Defined by Transcriptomics and Antibody-Based Profiling. FASEB J (2014) 28(7):2901–14. doi: 10.1096/fj.14-250555
28. Zaheer A, Fink BD, Lim R. Expression of Glia Maturation Factor Beta mRNA and Protein in Rat Organs and Cells. J Neurochem (1993) 60(3):914–20. doi: 10.1111/j.1471-4159.1993.tb03237.x
29. Aerbajinai W, Liu L, Zhu J, Kumkhaek C, Chin K, Rodgers GP. Glia Maturation Factor-Gamma Regulates Monocyte Migration Through Modulation of Beta1-Integrin. J Biol Chem (2016) 291(16):8549–64. doi: 10.1074/jbc.M115.674200
30. Aerbajinai W, Liu L, Chin K, Zhu J, Parent CA, Rodgers GP. Glia Maturation Factor-Gamma Mediates Neutrophil Chemotaxis. J Leukoc Biol (2011) 90(3):529–38. doi: 10.1189/jlb.0710424
31. Lippert DN, Wilkins JA. Glia Maturation Factor Gamma Regulates the Migration and Adherence of Human T Lymphocytes. BMC Immunol (2012) 13:21. doi: 10.1186/1471-2172-13-21
32. Fan J, Fong T, Chen X, Chen C, Luo P, Xie H. Glia Maturation Factor-Beta: A Potential Therapeutic Target in Neurodegeneration and Neuroinflammation. Neuropsychiatr Dis Treat (2018) 14:495–504. doi: 10.2147/NDT.S157099
33. Kuang XY, Jiang XF, Chen C, Su XR, Shi Y, Wu JR, et al. Expressions of Glia Maturation Factor-Beta by Tumor Cells and Endothelia Correlate With Neovascularization and Poor Prognosis in Human Glioma. Oncotarget (2016) 7(52):85750–63. doi: 10.18632/oncotarget.5509
34. Li YL, Ye F, Cheng XD, Hu Y, Zhou CY, Lu WG, et al. Identification of Glia Maturation Factor Beta as an Independent Prognostic Predictor for Serous Ovarian Cancer. Eur J Cancer (2010) 46(11):2104–18. doi: 10.1016/j.ejca.2010.04.015
35. Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, et al. ONCOMINE: A Cancer Microarray Database and Integrated Data-Mining Platform. Neoplasia (2004) 6(1):1–6. doi: 10.1016/s1476-5586(04)80047-2
36. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia (2017) 19(8):649–58. doi: 10.1016/j.neo.2017.05.002
37. Thul PJ, Lindskog C. The Human Protein Atlas: A Spatial Map of the Human Proteome. Protein Sci (2018) 27(1):233–44. doi: 10.1002/pro.3307
38. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the Cbioportal. Sci Signal (2013) 6(269):pl1. doi: 10.1126/scisignal.2004088
39. Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: Analyzing Multi-Omics Data Within and Across 32 Cancer Types. Nucleic Acids Res (2018) 46(D1):D956–63. doi: 10.1093/nar/gkx1090
40. Xia J, Gill EE, Hancock RE. NetworkAnalyst for Statistical, Visual and Network-Based Meta-Analysis of Gene Expression Data. Nat Protoc (2015) 10(6):823–44. doi: 10.1038/nprot.2015.052
41. Menyhart O, Nagy A, Gyorffy B. Determining Consistent Prognostic Biomarkers of Overall Survival and Vascular Invasion in Hepatocellular Carcinoma. R Soc Open Sci (2018) 5(12):181006. doi: 10.1098/rsos.181006
42. Yang SC, Yu EJ, Park JK, Kim TH, Eum JH, Paek SK, et al. The Ratio of Mitochondrial DNA to Genomic DNA Copy Number in Cumulus Cell May Serve as a Biomarker of Embryo Quality in IVF Cycles. Reprod Sci (2021) 28(9):2495–502. doi: 10.1007/s43032-021-00532-3
43. O'Hara R, Tedone E, Ludlow A, Huang E, Arosio B, Mari D, et al. Quantitative Mitochondrial DNA Copy Number Determination Using Droplet Digital PCR With Single-Cell Resolution. Genome Res (2019) 29(11):1878–88. doi: 10.1101/gr.250480.119
44. Shpilka T, Du Y, Yang Q, Melber A, Uma Naresh N, Lavelle J, et al. UPR(mt) Scales Mitochondrial Network Expansion With Protein Synthesis via Mitochondrial Import in Caenorhabditis Elegans. Nat Commun (2021) 12(1):479. doi: 10.1038/s41467-020-20784-y
45. Yuan Y, Ju YS, Kim Y, Li J, Wang Y, Yoon CJ, et al. Comprehensive Molecular Characterization of Mitochondrial Genomes in Human Cancers. Nat Genet (2020) 52(3):342–52. doi: 10.1038/s41588-019-0557-x
46. Li Y, Guo X, Guo S, Wang Y, Chen L, Liu Y, et al. Next Generation Sequencing-Based Analysis of Mitochondrial DNA Characteristics in Plasma Extracellular Vesicles of Patients With Hepatocellular Carcinoma. Oncol Lett (2020) 20(3):2820–8. doi: 10.3892/ol.2020.11831
47. Yu C, Wang X, Huang L, Tong Y, Chen L, Wu H, et al. Deciphering the Spectrum of Mitochondrial DNA Mutations in Hepatocellular Carcinoma Using High-Throughput Sequencing. Gene Expr (2018) 18(2):125–34. doi: 10.3727/105221618X15185539348147
48. Qiao L, Ru G, Mao Z, Wang C, Nie Z, Li Q, et al. Mitochondrial DNA Depletion, Mitochondrial Mutations and High TFAM Expression in Hepatocellular Carcinoma. Oncotarget (2017) 8(48):84373–83. doi: 10.18632/oncotarget.21033
49. Hsu CC, Lee HC, Wei YH. Mitochondrial DNA Alterations and Mitochondrial Dysfunction in the Progression of Hepatocellular Carcinoma. World J Gastroenterol (2013) 19(47):8880–6. doi: 10.3748/wjg.v19.i47.8880
50. Yamada S, Nomoto S, Fujii T, Kaneko T, Takeda S, Inoue S, et al. Correlation Between Copy Number of Mitochondrial DNA and Clinico-Pathologic Parameters of Hepatocellular Carcinoma. Eur J Surg Oncol (2006) 32(3):303–7. doi: 10.1016/j.ejso.2006.01.002
51. Yin PH, Lee HC, Chau GY, Wu YT, Li SH, Lui WY, et al. Alteration of the Copy Number and Deletion of Mitochondrial DNA in Human Hepatocellular Carcinoma. Br J Cancer (2004) 90(12):2390–6. doi: 10.1038/sj.bjc.6601838
52. Bras M, Queenan B, Susin SA. Programmed Cell Death via Mitochondria: Different Modes of Dying. Biochem (Mosc) (2005) 70(2):231–9. doi: 10.1007/s10541-005-0105-4
53. Jannus F, Medina-O'Donnell M, Rivas F, Diaz-Ruiz L, Rufino-Palomares EE, Lupianez JA, et al. A Diamine-PEGylated Oleanolic Acid Derivative Induced Efficient Apoptosis through a Death Receptor and Mitochondrial Apoptotic Pathway in HepG2 Human Hepatoma Cells. Biomolecules (2020) 10(19):1375. doi: 10.3390/biom10101375
54. Liu C, Li Y, Xing Y, Cao B, Yang F, Yang T, et al. The Interaction Between Cancer Stem Cell Marker CD133 and Src Protein Promotes Focal Adhesion Kinase (FAK) Phosphorylation and Cell Migration. J Biol Chem (2016) 291(30):15540–50. doi: 10.1074/jbc.M115.712976
55. Zhang B, Wu J, Cai Y, Luo M, Wang B, Gu Y. TCF7L1 Indicates Prognosis and Promotes Proliferation Through Activation of Keap1/NRF2 in Gastric Cancer. Acta Biochim Biophys Sin (Shanghai) (2019) 51(4):375–85. doi: 10.1093/abbs/gmz015
56. Sun Y, Long J, Zhou Y. Angiopoietin-Like 4 Promotes Melanoma Cell Invasion and Survival Through Aldolase a. Oncol Lett (2014) 8(1):211–7. doi: 10.3892/ol.2014.2071
57. Zou J, Zhu X, Xiang D, Zhang Y, Li J, Su Z, et al. LIX1-Like Protein Promotes Liver Cancer Progression via miR-21-3p-Mediated Inhibition of Fructose-1,6-Bisphosphatase. Acta Pharm Sin B (2021) 11(6):1578–91. doi: 10.1016/j.apsb.2021.02.005
58. Xiang J, Chen C, Liu R, Gou D, Chang L, Deng H, et al. Gluconeogenic Enzyme PCK1 Deficiency Promotes CHK2 O-GlcNAcylation and Hepatocellular Carcinoma Growth Upon Glucose Deprivation. J Clin Invest (2021) 131(8):e144703. doi: 10.1172/JCI144703
59. Jiang H, Ma N, Shang Y, Zhou W, Chen T, Guan D, et al. Triosephosphate Isomerase 1 Suppresses Growth, Migration and Invasion of Hepatocellular Carcinoma Cells. Biochem Biophys Res Commun (2017) 482(4):1048–53. doi: 10.1016/j.bbrc.2016.11.156
60. Zhang W, Jin K, Wang F, Zhangyuan G, Yu W, Liu Y, et al. Differences in the Prognostic Value of Tumor Size on Hepatocellular Cancer-Specific Survival Stratified by Gender in a SEER Population-Based Study. United Eur Gastroenterol J (2019) 7(7):933–41. doi: 10.1177/2050640619845602
61. Mauvais-Jarvis F, Clegg DJ, Hevener AL. The Role of Estrogens in Control of Energy Balance and Glucose Homeostasis. Endocr Rev (2013) 34(3):309–38. doi: 10.1210/er.2012-1055
62. Nagasue N, Kohno H. Hepatocellular Carcinoma and Sex Hormones. HPB Surg (1992) 6(1):1–6. doi: 10.1155/1992/72761
63. Carruba G. Estrogens in Hepatocellular Carcinoma: Friends or Foes? Cancers (Basel) (2021) 13(9):2085. doi: 10.3390/cancers13092085
64. Sun RF, Zhao CY, Chen S, Yu W, Zhou MM, Gao CR. Androgen Receptor Stimulates Hexokinase 2 and Induces Glycolysis by PKA/CREB Signaling in Hepatocellular Carcinoma. Dig Dis Sci (2021) 66(3):802–13. doi: 10.1007/s10620-020-06229-y
65. Kalra M, Mayes J, Assefa S, Kaul AK, Kaul R. Role of Sex Steroid Receptors in Pathobiology of Hepatocellular Carcinoma. World J Gastroenterol (2008) 14(39):5945–61. doi: 10.3748/wjg.14.5945
66. Murakami K, Hata S, Miki Y, Sasano H. Aromatase in Normal and Diseased Liver. Horm Mol Biol Clin Investig (2018) 41(1):20170081. doi: 10.1515/hmbci-2017-0081
67. Vizoso FJ, Rodriguez M, Altadill A, Gonzalez-Dieguez ML, Linares A, Gonzalez LO, et al. Liver Expression of Steroid Hormones and Apolipoprotein D Receptors in Hepatocellular Carcinoma. World J Gastroenterol (2007) 13(23):3221–7. doi: 10.3748/wjg.v13.i23.3221
68. Xu H, Wei Y, Zhang Y, Xu Y, Li F, Liu J, et al. Oestrogen Attenuates Tumour Progression in Hepatocellular Carcinoma. J Pathol (2012) 228(2):216–29. doi: 10.1002/path.4009
Keywords: hepatocellular carcinoma, GMFB, prognosis, bioinformatics, mitochondria function
Citation: Sun W, Hu C, Wang T, Wang J, Zhang J, Gao F, Ou Q, Tian H, Jin C, Xu J, Zhang J, Xu G-T and Lu L (2021) Glia Maturation Factor Beta as a Novel Biomarker and Therapeutic Target for Hepatocellular Carcinoma. Front. Oncol. 11:744331. doi: 10.3389/fonc.2021.744331
Received: 20 July 2021; Accepted: 07 October 2021;
Published: 02 November 2021.
Edited by:
Jianhui Zhu, University of Michigan, United StatesReviewed by:
Apurva Patel, Gujarat Cancer & Research Institute, IndiaDevanjan Sinha, Banaras Hindu University, India
Copyright © 2021 Sun, Hu, Wang, Wang, Zhang, Gao, Ou, Tian, Jin, Xu, Zhang, Xu and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guo-Tong Xu, gtxu@tongji.edu.cn; Lixia Lu, lulixia@tongji.edu.cn
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work