Pathol. Oncol. Res., 01 April 2021

Identification of Potential Biomarkers From Hepatocellular Carcinoma With MT1 Deletion

www.frontiersin.orgRuohao Zhang1, www.frontiersin.orgMiao Huang2, www.frontiersin.orgHong Wang1, www.frontiersin.orgShengming Wu3, www.frontiersin.orgJiali Yao1, www.frontiersin.orgYingying Ge4, www.frontiersin.orgYufei Lu1* and www.frontiersin.orgQiping Hu1*
  • 1Department of Cell Biology and Genetics, School of Pre-clinical Medicine, Guangxi Medical University, Nanning, China
  • 2Radiology Department, Affiliated Tumor Hospital of Guangxi Medical University, Nanning, China
  • 3Department of Pathology, Affiliated Tumor Hospital of Guangxi Medical University, Nanning, China
  • 4Department of Histology and Embryology, School of Pre-clinical Medicine, Guangxi Medical University, Nanning, China

Background: Hepatocellular carcinoma (HCC) is one of the deadliest cancers worldwide. Metallothioneins (MTs) are metal-binding proteins involved in multiple biological processes such as metal homeostasis and detoxification, as well as in oncogenesis. Copy number variation (CNV) plays a vital role in pathogenesis and carcinogenesis. Nevertheless, there is no study on the role of MT1 CNV in HCC.

Methods: Array-based Comparative Genomic Hybridization (aCGH) analysis was performed to obtain the CNV data of 79 Guangxi HCC patients. The prognostic effect of MT1-deletion was analyzed by univariate and multivariate Cox regression analysis. The differentially expressed genes (DEGs) were screened based on The Gene Expression Omnibus database (GEO) and the Liver Hepatocellular Carcinoma of The Cancer Genome Atlas (TCGA-LIHC). Then function and pathway enrichment analysis, protein-protein interaction (PPI) and hub gene selection were applied on the DEGs. Lastly, the hub genes were validated by immunohistochemistry, tissue expression and prognostic analysis.

Results: The MT1-deletion was demonstrated to affect the prognosis of HCC and can act as an independent prognostic factor. 147 common DEGs were screened. The most significant cluster of DEGs identified by Molecular Complex Detection (MCODE) indicated that the expression of four MT1s were down-regulated. MT1X and other five hub genes (TTK, BUB1, CYP3A4, NR1I2, CYP8B1) were associated with the prognosis of HCC. TTK, could affect the prognosis of HCC with MT1-deletion and non-deletion. NR1I2, CYP8B1, and BUB1 were associated with the prognosis of HCC with MT1-deletion.

Conclusions: In the current study, we demonstrated that MT1-deletion can be an independent prognostic factor in HCC. We identified TTK, BUB1, NR1I2, CYP8B1 by processing microarray data, for the first time revealed the underlying function of MT1 deletion in HCC, MT1-deletion may influence the gene expression in HCC, which may be the potential biomarkers for HCC with MT1 deletion.


Globally, liver cancer is one of the leading causes of cancer death. HCC, a predominant type of liver cancer, which counts for approximately 80% of all the liver cancer cases [1]. Nevertheless, the underlying molecular mechanisms of HCC is currently unclear.

CNV is caused by the rearrangement of the genome. It generally refers to the duplication or deletion of large fragments of the genome with a length of more than 1 kB, can activate oncogenes and neoplastic pathways [2]. Some studies have shown that copy number variation in HCC may affect gene expression alterations, thus affecting the occurrence and development of HCC [3, 4]. Genomic copy number deletion at chromosome 14q31.1-32.13 was frequently observed in HCC, as well as zinc finger CCCH-type containing 14 (ZC3H14) is located on 14q31.1-32.13, which was demonstrated to be related with HCC [5].

Metallothioneins (MTs) are metal-binding proteins involved in diverse processes, including metal homeostasis and detoxification, oxidative stress response, and cell proliferation, which have four main isoforms (MT1, MT2, MT3, and MT4) [6, 7]. Emerging studies have demonstrated that the abnormal expression of MTs, such as MT1 are able to trigger the process of carcinogenesis in various types of human malignancies, including HCC [8]. However, currently, there is no research on the relationship between MT1 CNV and HCC.

Hence, in our present study, we performed an aCGH analysis to obtain the CNV data of Guangxi HCC patients, and conducted the bioinformatic analysis on HCC datasets from the online database (Figure 1), The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Four MT1s and five hub genes were screened to investigate the underlying connection of MT1 CNV and HCC. Then we investigated the prognostic significance of these genes for HCC outcomes. MT1X, TTK, BUB1, CYP3A4, NR1I2, and CYP8B1 were associated with the prognostic of HCC. Furthermore, the prognostic analysis of these genes in MT1-deletion and non-deletion cases were performed and TTK, BUB1, NR1I2, and CYP8B1 may be the potential biomarkers for HCC with MT1-deletion.


FIGURE 1. Flowchart of the identification and validation process.

Methods and Materials

Patient and Sample Collection

We selected 79 HCC patients (termed Guangxi HCC cohort) with pathologic diagnosis and clinical follow-up results in the Affiliated Tumor Hospital of Guangxi Medical University from 2014 to 2016. As of June 30, 2018, 31 of the 79 patients were alive, 38 died and 10 lost to follow-up. Formalin-fixed paraffin embedded (FFPE) tumor tissue blocks from these patients were collected from the tumor biobank of Guangxi Medical University. This research project was approved by the Medical Ethics Committee of Guangxi Medical University (20170228-31) and research applications of residual specimens followed university policies and laboratory standards [9].

Array-Based Comparative Genomic Hybridization Analysis

Genomic DNAs of Guangxi HCC cohort were extracted from fresh tissue using Puregene Kit following manufacturer’s instruction (Qiagen Inc., Valencia, CA). DNA concentration was measured using a NanoDrop spectrophotometer (ND-1000, Thermo Fisher Scientific Inc., Waltham, MA) and high molecular weight DNA quality was verified by agarose gel electrophoresis. For each aCGH analysis, 2.5 μg of test genomic DNA from the patient and 2.5 μg of control DNA from a sex-matched or -mismatched healthy individual were used following the manufacturer’s protocol for the Agilent Human Genome aCGH microarray 60 K (Agilent Technologies Inc., Santa Clara, CA). The differential labeling of test and control DNAs, comparative hybridization onto 8 × 60 K Agilent slides, posthybridization wash, slide scanning, image feature extraction were processed following manufacturer’s instruction. The data were analyzed using Agilent CytoGenomics (version with the built-in ADM-2 algorithm set at threshold value of 6, a cut off value of 0.2, and a filter of ten probes. All CNVs except the recognized copy number variants from the HCC cohort of Genomic Variants (http://projects.tcag.ca/variation/) were recorded. The base pair designations from the Agilent 60 K array are according to the Feb. 2009 (GRCh37/hg19) on the UCSC Human Genome browser (http://genome.ucsc.edu/).

Online Data Collection and Processing

The CNV data of MT1 were from Copy Number Variation category of LIHC cohort from The Cancer Genome Atlas (TCGA, https://www.cancer.gov/tcga) database (termed TCGA-LIHC). The microarray data of GSE112790 and GSE101685 of HCC including tumor and normal samples were from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/), the online data used in the differential expression gene analysis was shown in Table 1. To avoid errors caused by platform differences, the datasets from GEO were all based on the platform GPL570 ([HG-U133_ Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array).


TABLE 1. Online data used in the differential expression gene analysis.

DEGs Identification and Venn Screening

Limma package is a R package which is widely used in biology field [10]. We used Limma to determine the differentially expressed genes (DEGs) between MT1-deletion and non-deletion samples. |logFC| > 1 and FDR < 0.05 were considered statistically significant. Next, we performed the GEO2R, an official analysis tool to screen DEGs between tumor samples and normal samples of GEO datasets with |logFC| > 1 and adjust p Value < 0.05. The common genes among all the three datasets were investigated by VennDiagram package in R [11].

Protein-Protein Interaction Network Construction and Hub Gene Recognition

We used the STRING (https://string-db.org) database to analyze the up-regulated and down-regulated DEGs and analyzed the protein-protein interaction (PPI) [12]. Cytoscape is a software for visualizing interaction networks and biological pathways [13]. The MCODE plugin was used to find clusters in PPI networks with the degree cutoff, node score cutoff, k-core and max depth as 2, 0.2, 0.2, and 100 as threshold. Moreover, the cytoHubba plugin was used to identify hub genes of the network we imported by calculating the node scores [14]. To get a more reliable result, we analyzed the top 10 nodes with highest degree with all the 12 means. Then we ordered the number of occurrences of these genes, and the genes with the highest occurrence were the most significant hub genes.

GO and KEGG Pathway Analysis

The Gene Ontology (GO, http://www.geneontology.org/) knowledgebase is a source of information on the functions of genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource for understanding high-level functions and utilities of the biological system. To investigate the potential function among the DEGs, we applied GO and KEGG pathway enrichment on the Database for Annotation, Visualization and Integrated Discovery (DAVID, https://david.ncifcrf.gov/), a bioinformatic website that is designed to investigate the biological functions of a number of genes [15].

Survival and Expression Analysis of Hub Genes

Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/) is a web server for cancer and normal gene expression profiling and interactive analysis [16]. The clinical prognosis and expression of the hub genes we found were analyzed by GEPIA and Kaplan-Meier Plotter database (http://kmplot.com) based on TCGA-LIHC data. The prognostic effect of the hub genes in MT1-deletion and non-deletion cases were performed by Survival package in R. Then the Human Protein Atlas (HPA, https://www.proteinatlas.org), a database aims to map all the human proteins in cells, tissues and organs using integration of various omics technologies, were used to analyze the protein expression differentiation between tumor and normal tissue [17].


MT1 Deletion Has a High Incidence and Causes Poor Prognosis in Guangxi HCC Cohort

The aCGH analysis results show that MT1 deletion occurred in 35 of the 79 (44.3%) samples. The clinical information of 79 Guangxi HCC patients are displayed in Table 2. Then the prognostic analysis shows that the deletion of MT1 had a significant association with poor prognosis of HCC patients (Figure 2A). Then we performed univariate and multivariate Cox regression analysis in Guangxi cohort (Table 3). The results show that CNV (MT1 deletion), BCLC (Barcelona Clinic Liver Cancer staging) and postoperative recurrence can be independent prognostic factors in HCC. We used the MT1 deletion segment less than −0.2 as the screening threshold, and screened a total of 372 cancer tissue samples in TCGA-LIHC dataset, of which 127 were with MT1 deletion (34.14%). We further compared the MT1 deletion rate of cancer samples from Guangxi HCC cohort and TCGA-LIHC dataset, and found no significant difference between the two rates (p > 0.5, Figure 2B). Given the same MT1 deletion rate, we analyzed transcriptome data and clinical follow-up data from TCGA and GEO datasets to identify underlying molecules that may function in MT1 deletion in HCC.


TABLE 2. The clinical information of Guangxi cohort.


FIGURE 2. The CNV and prognostic analysis of MT1 deletion in Guangxi HCC cohort and TCGA-LIHC dataset. (A) prognostic analysis of MT1 deletion in Guangxi HCC cohort; (B) MT1 deletion rate comparing between Guangxi HCC cohort and TCGA-LIHC dataset.


TABLE 3. The univariate and multivariate Cox regression analysis in Guangxi cohort.

147 Common DEGs Were Screened Through the Intersection of Three Datasets

Based on the deletion or non-deletion of MT1, we got 869 DEGs from TCGA-LIHC datasets, including 461 up-regulated genes and 408 down-regulated genes. Based on the tumor and normal samples of GEO datasets, we got 1,915 and 1,664 differentially expressed genes in GSE112790 and GSE101685, respectively, of which 1,088 were up-regulated and 827 were down-regulated in GSE112790 and 733 were up-regulated and 931 were down-regulated in GSE101685 (Figures 3A–C). Taking the intersection of these datasets, we got 147 common genes, of which 71 were up-regulated and 76 were down-regulated (Figures 3D,E).


FIGURE 3. Volcano plot and Venn plot of differentially expressed genes of three datasets. (A) Differentially expressed genes of MT1-deletion and non-deletion HCC samples of TCGA-LIHC; (B–C) Differentially expressed genes of tumor and normal samples of GSE112790 and GSE101685 (D–E) Common differentially expressed genes of three datasets.

The Up-, and Down-Regulated DEGs May Be Involved in Cell Division and Metabolism Respectively

The 147 DEGs we identified were analyzed by DAVID further. As the figure shows, in up-regulated common genes, the major biological processes were “mitotic nuclear division” and “cell division,” the major cellular component were “cytosol” and “nucleoplasm” and the molecular function analysis indicated that these genes were mainly enriched in “protein binding” and “ATP binding” (Figure 4A). The KEGG pathway analysis showed these genes were mostly enriched in cell cycle pathway (Figure 4A). Likewise, in down-regulated common genes, biological process analysis indicated that these genes were mainly associated with “oxidation-reduction process” and “xenobiotic metabolic process,” the major cellular component was “mitochondrion” and the major molecular function was “zinc ion binding” (Figure 4B). Lastly, the KEGG analysis demonstrated that these genes were most significantly enriched in metabolic pathways (Figure 4B).


FIGURE 4. Functional annotation and pathway analysis of common DEGs. (A) GO and KEGG pathway enrichment of up-regulated common genes; (B) GO and KEGG pathway enrichment of down-regulated common genes. CYP3A4, NR1I2, CYP8B1, TTK, BUB1, and four MT1s were selected as high score hub genes from the 147 common DEGs.

The PPI network constructed by Cytoscape is shown in Figure 5. The MCODE plugin was applied to detect the cluster that has the highest score among the PPI network (Figures 6A–D). Finally, we got a cluster with 47 nodes and 1,047 edges in up-regulated genes and three clusters with 4 nodes 6 edges, 3 nodes 3 edges and 7 nodes 8 edges in down-regulated genes, respectively. The most significant cluster with the highest K-core contains four genes including MT1G, MT1E, MT1F, and MT1X, that were down-regulated and closely related to this study. Then the 10 most significant up-regulated hub genes that cytoHubba obtained were, TTK, BUB1, TOP2A, TPX2, DLGAP5, ASPM, UBE2C, CDC20, MKI67, BUB1B, of which TTK and BUB1 ranked first count 10 times were selected as up-regulated hub genes. Likewise, the 10 most significant down-regulated hub genes were CYP3A4, NR1I2, CYP8B1, CYP1A2, ESR1, TAT, PCK1, SLC10A1, AR, CTH, of which CYP3A4, NR1I2, and CYP8B1 ranked first were selected as down-regulated hub genes (Table 4).


FIGURE 5. PPI network construction. (A) Up-regulated genes; (B) Down-regulated genes.


FIGURE 6. Clusters of MCODE (K-score > 2). (A) Cluster of up-regulated genes (B–D) Clusters of down-regulated genes. The expression and prognostic effect of MT1X and other five hub genes are validated.


TABLE 4. Hub genes analysis of common DEGs by Cytoscape, MCODE and cytoHubba.

The expression profile of four MT1s were validated in GEPIA (Figures 7A–D), as well as the prognosis analysis was applied in Kaplan-Meier Plotter (Figures 8A–D). Except for MT1X which was related to the prognosis of HCC, the other three MT1s genes had no association with the prognosis of HCC patients, though all of the four genes were low expressed in HCC.


FIGURE 7. Expression analysis of four MT1s and five hub genes in HCC tissue and adjacent tissue based on TCGA-LIHC data. Orange and gray background represents tumor and normal tissue, respectively; *p < 0.05(A) MT1G; (B) MT1E; (C) MT1F; (D) MT1X; (E) TTK; (F) BUB1; (G) CYP3A4; (H) NR1I2; (I) CYP8B1.


FIGURE 8. The prognostic effect of four MT1s and five hub genes in Kaplan-Meier Plotter database. (A) MT1G; (B) MT1E; (C) MT1F; (D) MT1X; (E)TTK; (F) BUB1; (G) CYP3A4; (H) NR1I2; (I) CYP8B1.

Given that only the prognostic effect of MT1X was statistically significant, we explored the expression and prognostic effect of the other five hub genes. The results show that the expressions of two up-regulated hub genes (TTK, BUB1) and three down-regulated hub genes (CYP3A4, NR1I2, CYP8B1) in tumor and normal samples were consistent with their differential expression analysis: TTK and BUB1 were overexpressed, while CYP3A4, NR1I2, and CYP8B1 were low in expression in HCC respectively (Figures 7E–I). Similarly, both up-regulated genes were significantly associated with HCC patients’ prognosis, with higher expression and worse prognosis in HCC (p < 0.05) while all the three down-regulated genes were associated with significant worse outcomes (p < 0.05) (Figures 8E–I). In order to explore the prognostic effect of these genes in the MT1 deletion and non-deletion groups, we conducted survival analysis in these two groups separately. In MT1-deletion cases, only TTK is statistically significant; in MT1-non-deletion cases, NR1I2, CYP8B1, TTK, and BUB1 are statistically significant (Figures 9, 10). However, the p Values of NR1I2, CYP8B1, and BUB1 are all close to 0.05 in MT1-deletion cases. We think that the reason for this result is the MT1-deletion group has fewer samples than the non-deletion group (127 cases compared to 245 cases). At last, we checked the protein expression distribution of these genes in The Human Protein Atlas, TTK and CYP3A4 were found to be differentially expressed in tumor and normal tissues, which TTK was overexpressed in HCC and the expression of CYP3A4 was found decreased in HCC (Figure 11).


FIGURE 9. The prognostic effect of hub genes in MT1-deletion cases. (A) CYP3A4; (B) NR1I2; (C) CYP8B1; (D) TTK; (E) BUB1; (F) MT1X.


FIGURE 10. The prognostic effect of hub genes in MT1-non-deletion cases. (A) CYP3A4; (B) NR1I2; (C) CYP8B1; (D) TTK; (E) BUB1; (F) MT1X.


FIGURE 11. Validation of the protein expression in the Human Protein Atlas database. (A) TTK; (B) CYP3A4.


HCC is one of the most harmful diseases to human life throughout the world and its pathogenesis is extremely complicated. Multiple biological processes including genetic and environment factors are involved in the development of HCC. Thus, discovering potential biomarkers and molecular mechanisms of HCC genesis and development is vital for diagnosis and treatment. In recent years, as researchers have deepened their understanding of CNV, the application of CNV in HCC research has also increased. However, so far, there is no research on the relationship between the metallothionein family CNV and HCC.

In this study, we demonstrated that MT1 deletion had a high incidence in Guangxi HCC patients and was associated with poor HCC prognosis. Then we screened 147 common DEGs between MT1 deletion DEGs and HCC DEGs form TCGA and GEO datasets, including 71 up-regulated genes and 76 down-regulated genes. The up-regulated DEGs were mainly enriched in cell cycle and the p53 signaling pathway, while the down-regulated DEGs were mainly involved in several metabolic pathways. We have learnt that cancer is characterized by aberrant cell cycle activity, the aberrant activity of cell cycle protein is a hallmark of human cancer [18, 19]. p53 is well known as a tumor suppressor gene, and its abnormal expression is one of the main causes of cancer [20]. What’s more, the metabolic pathway in HCC is generally abnormal and it is essential to the development and growth of HCC [21].

Interestingly, through PPI analysis, we got a four-gene cluster with MT1G, MT1E, MT1F, and MT1X. It has been reported that the MT1 family are abnormally expressed in a variety of cancers including HCC, which have a diverse expression pattern in different type of cancers [2224]. In HCC, the expression of most MT1 family are significantly down-regulated, and will affect the occurrence and development of HCC [8]. With regard to the four MT1s we identified, studies have confirmed that the expression of MT1G could serve as tumor suppressor by inhibiting cell proliferation and enhancing apoptosis in HCC [25]; MT1F was proven to inhibit the HCC cell growth [26]; Liu et al. demonstrated that MT1X is probably a prognostic biomarker and reduces the progression and metastasis of HCC [27]. Likewise, we validated that four MT1s were all downregulated in HCC with MT1 deletion. Therefore, the four MT1s (MT1G, MT1E, MT1F and MT1X) are down regulated in tumor samples, at the same time, they are also down regulated in tumor samples with MT1 deletion, indicating that MT1 deletion may affect the expression of its genes in HCC. However, the prognosis results showed that the four MT1s had weak correlations with HCC patient prognosis, except MT1X. It may be due to differences between groups caused by different datasets, it requires further effort to investigate. Afterward, we turned our attention to other common DEGs to explore new biomarkers which possibly function in HCC.

Then five genes form PPI network were selected by comprehensive analysis as hub genes by cytoHubba, including TTK, BUB1, CYP3A4, NR1I2 and CYP8B1. TTK Protein kinase, also known as MPS1 is the core component of the spindle assembly checkpoint, which functions to ensure the chromosomes are properly separated [28]. A number of studies have shown that TTK has an important role in cancer, including HCC [2932]. Budding uninhibited by benzimidazoles 1 (BUB1) encodes serine/threonine protein kinases that are crucial in spindle assemble checkpoint [33], as well as investigations have demonstrated the role in oncogenesis of BUB1 mutation in several cancers [34]. Whereas, the role of BUB1 in carcinogenesis and development of HCC has not been well studied. In this study, we found BUB1 was significantly up-regulated and associated with poor prognosis in HCC with MT1 deletion, which may be involved in HCC carcinogenetic process.

CYP3A4 and CYP8B1 are the members of cytochrome P450 family which may rely on the function of metabolic carcinogens to play a vital role in chemoprevention, carcinogenesis, cancer therapy and metastasis [35]. CYP3A4 has been well studied for decades. Sandström et al. [36] have reported that the activation of CYP3A4 in colorectal cancer cells may affect tumor sensitivity to certain anti-colon cancer drugs. In breast cancer, Murray and colleagues found high expression of CYP3A4 was associated with a poor prognosis [37]. Nevertheless, Fujimura and colleagues verified that lower immunoreactivity score showed a poorer prognosis in prostate cancer patients [38]. Moreover, CYP3A4 has been elucidated to be low expressed and as a biomarker for predicting poor prognosis in HCC [39, 40]. Nuclear receptor subfamily 1 group I member 2 (NR1I2) is a member of nuclear receptor family that regulates the expression of metabolizing and detoxifying enzymes [41, 42]. It is reported that the expression of NR1I2 was often related to chemoresistance in various cancer cells including HCC [4346]. In addition, NR1I2 is also the main regulator of CYP3A4 [47]. The researches on CYP8B1 are currently scarce. Some studies demonstrated that CYP8B1 transcription was regulated by bile acids, which participated in lipid metabolisms and other liver diseases [48]. Whereas the role of CYP8B1 in HCC remains uncharted. Our results revealed that CYP8B1 was low in expression and had poor prognostic effect in HCC. Next, we used database analysis to verify the differential expression of these key genes in HCC tissue and adjacent normal tissue. Various results indicated that these genes were significantly differentially expressed in cancer and adjacent tissue, and significantly affected the prognosis of HCC. The above results indicated that MT1 deletion was related to the poor prognosis of HCC, which was possible mainly caused by the regulation of cell proliferation and metabolism through the differential expression of CYP3A4, TTK, NR1I2, CYP8B1, and BUB1.

This study also has some shortcomings, for example, we need to further verify the impact of these hub genes on HCC by experiment, and explore the specific mechanism between MT1 deletion and its expression in more depth.

In summary, this article demonstrated MT1-deletion can be an independent prognostic factor in HCC, screened 147 differential expression genes through aCGH and bioinformatics, and found multiple genes may be related to MT1 deletion in HCC. Here, we bridged the gap between the MT1 deletion and HCC, MT1 deletion may cause abnormal expression of genes (TTK, BUB1, NR1I2, CYP8B1) in HCC, thus affecting the prognosis of MT1 deleted HCC patients.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.cancer.gov/tcgahttps://www.ncbi.nlm.nih.gov/geo/.

Ethics Statement

This research project was approved by the Medical Ethics Committee of Guangxi Medical University (20170228-31).

Author Contributions

Study concepts and design: QH and YL; Sample collection and aCGH analysis: RZ and QH; Statistic analysis and data interpretation: RZ, MH, HW, SW, JY, YG, YL, and QH; Manuscript preparation: QH, RZ, and YL. All authors read and approved the final manuscript.


The present study was supported by the National Natural Science Foundation of China (No. 81760424).

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.


The results published here are in part based upon data generated by the TCGA Research Network (https://www.cancer.gov/tcga) and the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/).


1. Bray, F, Ferlay, J, Soerjomataram, I, Siegel, RL, Torre, LA, and Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: a Cancer J clinicians (2018). 68(6):394–424. doi:10.3322/caac.21492

CrossRef Full Text | Google Scholar

2. Shlien, A., and Malkin, D. Copy number variations and cancer. Genome Med (2009). 1(6):62. doi:10.1186/gm62

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Cao, P, Yang, A, Wang, R, Xia, X, Zhai, Y, Li, Y, et al. Germline duplication of SNORA18L5 increases risk for HBV-related hepatocellular carcinoma by altering localization of ribosomal proteins and decreasing levels of p53. Gastroenterology (2018). 155(2):542–56. doi:10.1053/j.gastro.2018.04.020

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Stranger, BE, Forrest, MS, Dunning, M, Ingle, CE, Beazley, C, Thorne, N, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science (2007). 315(5813):848–53. doi:10.1126/science.1136678

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zhang, C, Cao, P, Yang, A, Xia, X, Li, Y, Shi, M, et al. Downregulation of ZC3H14 driven by chromosome 14q31 deletion promotes hepatocellular carcinoma progression by activating integrin signaling. Carcinogenesis (2019). 40(3):474–86. doi:10.1093/carcin/bgy146

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Vasák, M. Advances in metallothionein structure and functions. J Trace Elem Med Biol (2005). 19(1):13–7. doi:10.1016/j.jtemb.2005.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Takahashi, S. Positive and negative regulators of the metallothionein gene (review). Mol Med Rep (2015). 12(1):795–9. doi:10.3892/mmr.2015.3459

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Si, M., and Lang, J. The roles of metallothioneins in carcinogenesis. J Hematol Oncol (2018). 11(1):107. doi:10.1186/s13045-018-0645-x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Grommisch, B, DiAdamo, AJ, Xu, ZY, Pan, XH, Ma, DQ, Xie, JS, et al. Biobanking of residual specimens from diagnostic genetic laboratories: standard operating procedures, ethical and legal considerations, and research applications. North Am J Med Sci (2013). 6(4):200–7. doi:10.7156/najms.2013.0604200

CrossRef Full Text | Google Scholar

10. Ritchie, ME, Phipson, B, Wu, D, Hu, Y, Law, CW, Shi, W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015). 43(7):e47. doi:10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Chen, H., and Boutros, PCJBB VennDiagram: a package for the generation of highly-customizable venn and euler diagrams in R. 12. BMC Bioinformatics (2011). 12:35. doi:10.1186/1471-2105-12-35

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Szklarczyk, D, Gable, AL, Lyon, D, Junge, A, Wyder, S, Huerta-Cepas, J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res (2019). 47(D1):D607–d613. doi:10.1093/nar/gky1131

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Shannon, P, Markiel, A, Ozier, O, Baliga, NS, Wang, JT, Ramage, D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res (2003). 13(11):2498–504. doi:10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chin, C-H, Chen, S-H, Wu, H-H, Ho, C-W, Ko, M-T, and Lin, C-Y cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol (2014). 8(Suppl 4):S11. doi:10.1186/1752-0509-8-s4-s11

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Dennis, G, Sherman, BT, Hosack, DA, Yang, J, Gao, W, Lane, HC, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol (2003). 4(5):P3. doi:10.1186/gb-2003-4-5-p3

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Tang, Z, Li, C, Kang, B, Gao, G, Li, C, and Zhang, Z GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res (2017). 45(W1):W98–w102. doi:10.1093/nar/gkx247

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Uhlén, M, Fagerberg, L, Hallström, BM, Lindskog, C, Oksvold, P, Mardinoglu, A, et al. Tissue-based map of the human proteome. Science (2015). 347(6220):1260419. doi:10.1126/science.1260419

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lapenna, S, and Giordano, A. Cell cycle kinases as therapeutic targets for cancer. Nat Rev Drug Discov (2009). 8(7):547–66. doi:10.1038/nrd2907

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Otto, T, and Sicinski, P. Cell cycle proteins as promising targets in cancer therapy. Nat Rev Cancer (2017). 17(2):93–115. doi:10.1038/nrc.2016.138

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Levine, AJ. p53: 800 million years of evolution and 40 years of discovery. Nat Rev Cancer (2020). 20(8):471–80. doi:10.1038/s41568-020-0262-1

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zarrinpar, A. Metabolic pathway inhibition in liver cancer. Slas Technology: Translating Life Sci Innovation (2017). 22(3):237–44. doi:10.1177/2472630317698683

CrossRef Full Text | Google Scholar

22. Wang, XT, Chen, CW, Zheng, XM, Wang, B, Zhang, SX, Yao, MH, et al. Expression and prognostic significance of melatonin receptor MT1 in patients with gastric adenocarcinoma. neo (2020). 67(2):415–20. doi:10.4149/neo_2019_190220N141

CrossRef Full Text | Google Scholar

23. Jagielski, L, Jelen, M, Kobierzycki, C, Jagielska, G, and Blok, R. Increase of nuclear expression of metallothionein I/II in neoplastic transformation of the endometrium. Ginekol Pol (2015). 86(3):182–7. doi:10.17772/gp/2060

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Wu, Y-J, Ko, B-S, Liang, S-M, Lu, Y-J, Jan, Y-J, Jiang, S-S, et al. ZNF479 downregulates metallothionein-1 expression by regulating ASH2L and DNMT1 in hepatocellular carcinoma. Cell Death Dis (2019). 10(6):408. doi:10.1038/s41419-019-1651-9

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wang, Y, Wang, G, Tan, X, Ke, K, Zhao, B, Cheng, N, et al. MT1G serves as a tumor suppressor in hepatocellular carcinoma by interacting with p53. Oncogenesis (2019). 8(12):67. doi:10.1038/s41389-019-0176-5

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lu, DD, Chen, YC, Zhang, XR, Cao, XR, Jiang, HY, and Yao, L. The relationship between metallothionein-1F (MT1F) gene and hepatocellular carcinoma. Yale J Biol Med (2003). 76(2):55–62.

PubMed AbstractGoogle Scholar

27. Liu, Z, Ye, Q, Wu, L, Gao, F, Xie, H, Zhou, L, et al. Metallothionein 1 family profiling identifies MT1X as a tumor suppressor involved in the progression and metastastatic capacity of hepatocellular carcinoma. Mol carcinogenesis (2018). 57(11):1435–44. doi:10.1002/mc.22846

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Musacchio, A. Spindle assembly checkpoint: the third decade. Phil Trans R Soc B (2011). 366(1584):3595–604. doi:10.1098/rstb.2011.0072

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Xie, Y, Wang, A, Lin, J, Wu, L, Zhang, H, Yang, X, et al. Mps1/TTK: a novel target and biomarker for cancer. J Drug Target (2017). 25(2):112–8. doi:10.1080/1061186x.2016.1258568

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Daniel, J, Coulter, J, Woo, J-H, Wilsbach, K, and Gabrielson, E. High levels of the Mps1 checkpoint protein are protective of aneuploidy in breast cancer cells. Proc Natl Acad Sci (2011). 108(13):5384–9. doi:10.1073/pnas.1007645108

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Maire, V, Baldeyron, C, Richardson, M, Tesson, B, Vincent-Salomon, A, Gravier, E, et al. TTK/hMPS1 is an attractive therapeutic target for triple-negative breast cancer. PloS one (2013). 8(5):e63712. doi:10.1371/journal.pone.0063712

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Jemal, A, Bray, F, Center, MM, Ferlay, J, Ward, E, and Forman, D. Global cancer statistics. CA: a Cancer J clinicians (2011). 61(2):69–90. doi:10.3322/caac.20107

CrossRef Full Text | Google Scholar

33. Manic, G, Corradi, F, Sistigu, A, Siteni, S, and Vitale, I. Molecular regulation of the spindle assembly checkpoint by kinases and phosphatases. Int Rev Cel Mol Biol (2017). 328:105–61. doi:10.1016/bs.ircmb.2016.08.004

CrossRef Full Text | Google Scholar

34. Bolanos-Garcia, VM, and Blundell, TL. BUB1 and BUBR1: multifaceted kinases of the cell cycle. Trends Biochemical Sciences (2011). 36(3):141–50. doi:10.1016/j.tibs.2010.08.004

CrossRef Full Text | Google Scholar

35. Alzahrani, AM, and Rajendran, P. The multifarious link between cytochrome P450s and cancer. Oxidative Med Cell longevity (2020). 2020:1. doi:10.1155/2020/3028387

CrossRef Full Text | Google Scholar

36. Sandström, R, Knutson, TW, Knutson, L, Jansson, B, and Lennernäs, H. The effect of ketoconazole on the jejunal permeability and CYP3A metabolism of (R/S)-verapamil in humans. Br J Clin Pharmacol (1999). 48(2):180–9. doi:10.1046/j.1365-2125.1999.00999.x

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Murray, GI, Patimalla, S, Stewart, KN, Miller, ID, and Heys, SD. Profiling the expression of cytochrome P450 in breast cancer. Histopathology (2010). 57(2):202–11. doi:10.1111/j.1365-2559.2010.03606.x

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Fujimura, T, Takahashi, S, Urano, T, Kumagai, J, Murata, T, Takayama, K, et al. Expression of cytochrome P450 3A4 and its clinical significance in human prostate cancer. Urology (2009). 74(2):391–7. doi:10.1016/j.urology.2009.02.033

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Ashida, R, Okamura, Y, Ohshima, K, Kakuda, Y, Uesaka, K, Sugiura, T, et al. Gene is a novel biomarker for predicting a poor prognosis in hepatocellular carcinoma. Cancer Genomics Proteomics (2017). 14(6):445–53. doi:10.21873/cgp.20054

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Yu, T, Wang, X, Zhu, G, Han, C, Su, H, Liao, X, et al. The prognostic value of differentially expressed CYP3A subfamily members for hepatocellular carcinoma. Cmar (2018). Vol. 10:1713–26. doi:10.2147/cmar.S159425

CrossRef Full Text | Google Scholar

41. Feng, F, Jiang, Q, Cao, S, Cao, Y, Li, R, Shen, L, et al. Pregnane X receptor mediates sorafenib resistance in advanced hepatocellular carcinoma. Biochim Biophys Acta (Bba)—Gen Subjects (2018). 1862(4):1017–30. doi:10.1016/j.bbagen.2018.01.011

CrossRef Full Text | Google Scholar

42. Harmsen, S, Meijerman, I, Maas-Bakker, RF, Beijnen, JH, and Schellens, JHM. PXR-mediated P-glycoprotein induction by small molecule tyrosine kinase inhibitors. Eur J Pharm Sci (2013). 48(4-5):644–9. doi:10.1016/j.ejps.2012.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Jiang, H, Chen, K, He, J, Pan, F, Li, J, Chen, J, et al. Association of pregnane X receptor with multidrug resistance-related protein 3 and its role in human colon cancer chemoresistance. J Gastrointest Surg (2009). 13(10):1831–8. doi:10.1007/s11605-009-0964-x

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Gupta, D, Venkatesh, M, Wang, H, Kim, S, Sinz, M, Goldberg, GL, et al. Expanding the roles for pregnane X receptor in cancer: proliferation and drug resistance in ovarian cancer. Clin Cancer Res (2008). 14(17):5332–40. doi:10.1158/1078-0432.Ccr-08-1033

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Chen, Y, Tang, Y, Wang, M-T, Zeng, S, and Nie, D. Human pregnane X receptor and resistance to chemotherapy in prostate cancer. Cancer Res (2007). 67(21):10361–7. doi:10.1158/0008-5472.Can-06-4758

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Yoshinari, K. Role of nuclear receptors PXR and CAR in xenobiotic-induced hepatocyte proliferation and chemical carcinogenesis. Biol Pharm Bull (2019). 42(8):1243–52. doi:10.1248/bpb.b19-00267

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Chen, H, Shen, ZY, Xu, W, Fan, TY, Li, J, Lu, YF, et al. Expression of P450 and nuclear receptors in normal and end-stage Chinese livers. Wjg (2014). 20(26):8681–90. doi:10.3748/wjg.v20.i26.8681

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Zhang, M, and Chiang, JYL. Transcriptional regulation of the human sterol 12α-hydroxylase gene (CYP8B1). J Biol Chem (2001). 276(45):41690–9. doi:10.1074/jbc.M105117200

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hepatocellular carcinoma, copy number variation, metallothionein, hub gene, deletion, screening, biomarker

Citation: Zhang R, Huang M, Wang H, Wu S, Yao J, Ge Y, Lu Y and Hu Q (2021) Identification of Potential Biomarkers From Hepatocellular Carcinoma With MT1 Deletion. Pathol. Oncol. Res. 27:597527. doi: 10.3389/pore.2021.597527

Received: 21 August 2020; Accepted: 01 February 2021;
Published: 01 April 2021.

Edited by:

József Tímár, Semmelweis University, Hungary

Copyright © 2021 Zhang, Huang, Wang, Wu, Yao, Ge, Lu and Hu. 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: Qiping Hu, huqiping@gxmu.edu.cn; Yufei Lu, lyff1002@126.com

These authors have contributed equally to this work