Abstract

The BolA gene family member (BOLA1–3) plays an important role in regulating normal and pathological biological processes including liver tumorigenesis. However, their expression patterns as prognostic factors in hepatocellular carcinoma (HCC) patients have not to be elucidated. We examined the transcriptional expressions and survival data of BolA family member in patients with HCC from online databases including ONCOMINE, TCGA, UALCAN, Gene Expression Profiling Interactive Analysis (GEPIA), Kaplan-Meier plotter, SurvExpress, cBioPortal, and Exobase. Network molecular interaction views of BolA family members and their neighborhoods were constructed by the IntAct web server. In our research, we had found that the expression levels of BolA /2/3 mRNA were higher in HCC tissue than in normal liver tissues from TGCA databases. Moreover, the BolA family gene expression level is significantly associated with distinct tumor pathological grade, TMN stage, and overall survival (OS). The BolA family can be considered as prognostic risk biomarkers of HCC. A small number of BolA gene-mutated samples were detected in the HCC tissue. IntAct analysis revealed that BolA1/2/3 was closely associated with the GLRX3 expression in HCC, which is implicated in the regulation of the cellular iron homeostasis and tumor growth. Furthermore, prognostic values of altered BolAs and their neighbor GLRX3 gene in HCC patients were validated by SurvExpress analysis. In conclusion, the membrane BolA family identified in this study provides very useful information for the mechanism of hepatic tumorigenesis.

1. Introduction

Hepatocellular carcinoma (HCC) has very aggressive neoplasms and describes as a major health problem worldwide [1]. Genetic and epigenetic alterations, which lead to uncontrolled cellular proliferation and metastasis, are the characters of HCC development.

Recent research has revealed a critical role for cellular iron homeostasis in the clinical context of liver tumorigenesis [2, 3]. Although significant progress has been made in understanding the iron homeostasis disruption associated with HCC, the precise molecular signals that trigger initiation and progression of HCC remain to be identified.

The human BolA gene family consists of BOLA1, BOLA2, and BOLA3 [4]. It has been suggested that BolA family members serve as assembly factors for mitochondrial iron-sulfur (Fe/S) cluster proteins that has involvement in cancer cell biology [5, 6]; although, the functions of BOLA1 and BOLA3 are still undefined in cancer. Prior research has highlighted the importance role of BOLA3 in human endothelial metabolism and cardiovascular disease pathogenesis [7]. More specifically, evidence points out that BOLA2 has been shown to be highly correlated with hepatic iron homeostasis [8]. And yet, even the overexpression of BOLA2 is required to drive HCC tumor growth and tumor hemorrhage [9, 10], and high BOLA2 can promote tumor growth and predict the HCC prognosis [11]. BOLA1 plays a leading role in mitochondrial morphology by potential regulation and can induce diseases [12]. In the ovarian cancer, the BOLA2 and BOLA3 were higher in cancer tissues and may act as prognostic biomarkers [13], and in the lung adenocarcinoma, the BOLA3 was correlated with the immune cell infiltrates [14]. However, it has been poorly characterized whether the expression of BolA family members in HCC is correlated with clinical outcomes.

In our research, we analyzed the BolA family member mRNA level in HCC tissues and nontumor liver tissues by the public database. In addition, we investigated correlation between their expressions and clinical characteristics and performed SurvExpress analysis of prognostic risks for overall survival. The results showed that BOLA1\2\3 may be a promising biomarker for the prognosis in HCC.

2. Material and Methods

2.1. ONCOMINE Database Analysis

The difference mRNA expression level of the BolA family gene in human cancer was identified in the ONCOMINE online microarray database (http://www.oncomine.org). For each BolA family gene, the thresholds were set as the following values: value of 0.01, fold change of 2, and gene ranking of all. Analysis type was set as follows: cancer vs. normal analysis.

2.2. UALCAN Database Analysis

The UALCAN online database (http://ualcan.path.uab.edu) was used to calculate the BolA gene expression level and clinicopathologic parameters in the TCGA database on patient with LIHC (liver hepatocellular carcinoma) [15] .

2.3. cBioPortal and Exobase Database Analysis

The cBio Cancer Genomics Portal (http://www.cbioportal.org/) performed to estimate the cancer genomics data sets of BolA family gene using TCGA-LIHC data [16]. The exoRBase database (http://www.exoRBase.org) can analysis the human blood exosomes, including circRNA, lncRNA, and mRNA [17].

2.4. GEPIA Database Analysis

Gene Expression Profiling Interactive Analysis (GEPIA) web server (http://gepia.cancer-pku.cn/) was used to study the correlation mRNA expression of BolA family members and overall survival (OS) in LIHC [18]. A total of 331 LIHC patients were enrolled, and “median” was regarded as group cutoff value.

2.5. Kaplan-Meier Plotter Analysis

The Kaplan-Meier (KM) plotter database (http://kmplot.com/analysis) was used to calculate the survival time in LIHC patients [19]. Briefly, each BolA family member was individually analyzed to obtain KM plots. Group cutoff was set as “median.” Hazard ratios (HR) with 95% confidence intervals (CI) were extracted from the KM plotter webpage. Overall survival (OS) data from 364 patients with HCC were enrolled.

2.6. SurvExpress Database Analysis

SurvExpress (http://bioinformatica.mty.itesm.mx/SurvExpress) was used for obtaining survival data for the expression of BolA family members in patient with LIHC, for which information was not available on the GEPIA and KM plotter database [20]. Briefly, in the TCGA-LIHC datasets containing 381 samples, BOLA1, BOLA2, and BOLA3 were entered into the number-at-risk cases, median mRNA expression levels, HRs, 95% confidence interval (CI), and values that were displayed.

2.7. IntAct Database Analysis

IntAct (http://www.ebi.ac.uk/intact) was applied to identify densely connected network components and BolA family members, for which protein-protein interaction enrichment analysis data populated by either curated from the literature or from direct data depositions [21].

2.8. Western Blot Analysis

Western blot analysis was performed as previously described [11]. The antibody dilutions were 1 : 1,000 for BOLA1 polyclonal antibody (Cat. # 18017-1-AP, Proteintech), 1 : 1,000 for BOLA2 polyclonal antibody (Cat. # ab169481, Abcam), 1 : 1,000 for BOLA3 polyclonal antibody (Cat. # ab185339, Abcam), and 1 : 5,000 for the β-actin mouse monoclonal antibody (Sigma-Aldrich, Cat. # A1978).

3. Results

3.1. BolA Family Members Are Frequently Upregulated in HCC

In order to analysis the expression differences of the BolA family, we first performed an analysis using the ONCOMINE database to investigate differences in the mRNA levels of each BolA family in cancers. As shown in Figure 1, the number of the upregulation BOLA1\2\3 expression was found in tumors compared with normal tissues in various types of cancers. Significantly higher mRNA expressions of BOLA1\2 were found in multiple HCC tissues datasets. The BOLA1\2 overexpression was found in HCC tissues compared with normal tissues in Roessler Liver 2 dataset (1.51-fold increase, ; 2.67-fold increase, , respectively) [22], while were observed in Wurmbach liver dataset (1.65-fold increase, ; 2.21-fold increase, , respectively) [23]. Significant upregulation of BOLA1\2 was also found in Chen Liver dataset (1.57-fold increase, ; 1.56-fold increase, , respectively) [24]. The abovementioned observations suggest that the overexpression of BOLA family members is associated with cancer progression and might be of clinical importance.

3.2. BolA Family Member Expression Was Higher in HCC

To further validate the observations made in the ONCOMINE database, TCGA-LIHC cohort performed a retrospective study. As shown in Figure 2, the BolA family expression level in HCC was higher than in the normal liver tissues (). In order to confirm this, we investigated protein levels of BolA family members by the Human Protein Atlas database (http://www.proteinatlas.org/pathology) [25]. As shown in Figure 3, BOLA1\2\3 proteins had lower level in the normal liver, while medium and high level were observed in HCC. And we also found that the BOLA1\2\3 expression level was higher in HCC than in the nontumor using our HCC samples. Human BolA proteins (BOLA1\2\3) are novel nonclassical secreted proteins [4]. In addition, a very low mutation rate of BOLA1\2\3 was observed in HCC patients (Figure S1A), the BOLA1 mutation rate was 4%, and there was no mutation in BOLA3. Intriguingly, using the Exosomes web-accessible database (http://www.exoRBase.org) analysis, the increased expression of BOLA2 may be used as circulating biomarkers for HCC patients (Figure S1B). Taken together, BOLA2 may had the potential ability for HCC diagnose.

3.3. Association between BolA Family Member and Tumor Grades and Stages

Both the mRNA and protein expression of BolA family members were found to be overexpressed in HCC; we next analyzed the relationship between mRNA expressions of each BolA family members with clinicopathological parameters of HCC patients by UALCAN. As was shown in Figure 4(a), we found that the elevated level of BOLA1\2\3 mRNA had a higher proportion of high-grade tumors (G3/G4). The BOLA1\2\3 mRNA level had significantly correlated with tumor stage in HCCs, which means that the advanced stage HCCs can express higher BolA mRNA (Figure 4(b)). The reason why mRNA expressions of BOLA1\2\3 in stage 3 seemed to be higher than that in stage 4 may be due to the small sample size (only 6 HCC patients were at stage 4). These findings indicated that the BOLA1\2\3 may accelerate HCC growth and progression.

3.4. BolA Family Member Predicts the Prognosis in HCC Patients

We used GEPIA web server to analyze the prognostic values of BolAs in TCGA-LIHC patients. As were shown in Figure 5, upregulation of BOLA1, BOLA2, and BOLA3 were significantly associated with shorter OS (, ; , ; , , respectively, Figures 5(a)5(c)). The relationship between combinatory mRNA expressions of all 3 BolA family members and prognosis of liver cancer patients were further analyzed by SurvExpress. In our study, we also found that higher combinatory mRNA expressions of all 3 BolA family members were associated with poorer OS in LIHC patients (, 95% CI: 1.14-2.26, and , Figure 5(d)). And then, we anlayed the prognostic role of BolA family members in HCC patients.. As was shown in Figure 5(e), the higher mRNA expression of BOLA1 (), BOLA2 (), and BOLA3 () was significantly associated with shorter OS of LIHC patients. These results indicated that mRNA expressions of BOLA1\2\3 may be exploited as useful biomarkers for prediction of HCC patient’s survival.

3.5. Identification of Hub BolA Family Member and Their Clinical Value in HCC

After analyzing the genetic alterations in BolAs and their prognostic value in HCC patients, we further analyzed the protein-protein interaction network among BolAs using IntAct databases. The top hub genes were GLRX3, DDIT4L, BIRC7, HDX, C1orf94, XIAP, BIRC2, BCKDHA, PMPCA, PMPCB, GLRX5, and SDHAF3 (Figure 6). As was shown in Figure 7, Kaplan-Meier (KM) plotter survival analysis, based on clinical information from the TCGA liver cancer datasets, revealed that the low expression of BIRC2 (, 95% CI: 0.46-0.96, and , Figure 7(c)), BCKDHA (, 95% CI: 0.34-0.74, and , Figure 7(d)), PMPCB (, 95% CI: 0.49-0.99, and , Figure 7(e)), and GLRX5 (, 95% CI: 0.5-1, and , Figure 7(f)) significantly correlated with shorter OS of LIHC patients. GLRX3 (, 95% CI: 1.44-2.92, and , Figure 7(a)) and BIRC7 (, 95% CI: 1.09-2.18, and , Figure 7(b)) were quite the contrary. Notably, higher combinatory mRNA expressions of BOLA2 with GLRX3 were associated with poorer OS in HCC patients (, 95% CI: 1.1-2.22, and , respectively, Figure 8). Many studies have investigated the expression of GLRX3 imply in regulating HCC cell proliferation, growth, and microvascular invasion via disruption of iron homeostasis [26]. Thus, we could guess that BOLA2 has the ability to promote the development of HCC and maintains cancer cell growth in the condition of metabolic stress.

4. Discussion

HCC is one of the leading causes of lethal, and there is great interest in understanding the underlying differentially expressed genes involved in the development and progression of individual tumors. In this study, we investigated the relationship between BolA family members and HCC patients using comprehensive data mining. We found that BolA family members are uniquely overexpressed in HCCs. Moreover, the mRNA expression levels of BolA family genes are associated with distinct tumor grade, TMN stage, and OS. Thus, BOLA1, BOLA2, and BOLA3 can predict the prognosis of HCC patients and may serve as oncogenes that promote HCC growth.

It has been proved that HCC development is a multistep process, including cell proliferation, adhesion, and metabolism. Iron metabolism plays an important role in both normal and cancer cells. In the process of HCC development, more iron is required to maintain the cancer cell proliferation, growth, and self-renewal in stem cells [27]. BOLA1, a mitochondrial protein, makes balances the effect of L-buthionine-(S, R)-sulfoximine (BSO)-induced glutathione (GSH) depletion on the mitochondrial thiol redox potential [12]. BOLA3 plays an important role in form [2Fe-2S] cluster-bridged dimeric heterocomplexes with the human monothiol glutaredoxin GRX5 [28]. A recent study indicated that BOLA1 and BOLA3 are associated with clinical outcomes in many diseases [5]. However, a thoughtful description of the relationship between expression level and cancer prognosis has not been analyzed. Although the increased expression of BOLA1/3 was obverse in present study, a correlation was observed between BOLA1/3 expression and defined genes in LIHC, such as oncogenic activity of BIRC2 [29] and tumor suppressor PMPCB [30]. Therefore, we can speculate that the BOLA1/3 expression in HCCs contributes to uncontrolled cell cellular proliferation. Further studies will be needed to clarify its role in HCC.

BOLA2, a gene associated with iron homeostasis, has been described in its biological function by the animal model [10]. The mechanisms of BOLA2 regulation are as follows: (i) specific in-frame fusion transcript regulation [31], (ii) monothiol CGFS glutaredoxin binding partners, (iii) GRX3-dependent anamorsin maturation pathway [32], and (iv) as c-Myc-regulated gene in HCC [10]. In our study, BOLA2 and GLRX3 are frequently overexpressed in HCC tumors tissues. Interestingly, our study revealed that the upregulation of BOLA2 and GLRX3 was associated with worse OS in patients with HCC. Up to now, more and more novel biomarkers, such as circular RNAs (circRNAs) [33], circulating microRNAs [34], and serum extracellular vesicles [35], had appeared for diagnosing HCC and predicting clinical outcomes. Our study analyzed the relationship between BOLA2 and serum extracellular vesicles. Hence, we postulate that the BOLA2 may have the potential for predicting the prognosis in HCC patients. Due to the limitations in our study, the relationship between the BOLA2 protein expression was not be clearly assessed, and further researches were needed to elaborate.

5. Conclusion

In our study, we found that BolA gene family members (BOLA1-3) may serve as prognostic biomarkers of HCC. In addition, BolA family members and their neighborhood GLRX3 play a leading role in HCC stage and tumor grade. These interesting results have important implications that can identify novel therapeutic targets in HCC.

Abbreviations

NP:Normal person
CHD:Coronary heart disease
CRC:Colorectal cancer
HCC:Hepatocellular carcinoma
PAAD:Pancreatic adenocarcinoma
WhB:Whole blood.

Data Availability

All the data used to support the findings of this study are available online.

Conflicts of Interest

The authors declare that they have no potential competing interests in this work.

Authors’ Contributions

YiMing Tao designed the study and wrote the manuscript. Dong Wang analyzed the BOLA family members in all the database. Dong Wang and ZhiMing Wang performed all the figures in this study.

Acknowledgments

This work was supported by the Natural Science Foundation of Shandong Province (No. ZR202103020004), Clinical Medicine Plus X Project of Affiliated Hospital of Qingdao University (202102006), and National Nature Science Foundation of China under Grant (Nos. 81372630, 81372631, 81702956). We thank the people for their efforts for providing data and databases.

Supplementary Materials

Figure S1: genetic alterations of 3 BolA family members were shown in HCC patients (cBioPortal). Notes: (a) OncoPrint of 3 BolA family member alterations in LIHC. (b) Using Exosomes web-accessible database (http://www.exoRBase.org) analysis, the increased expression of BOLA2 may be used as circulating biomarkers for HCC patients. (Supplementary Materials)