Liver cancer is a malignant cancer with poor prognosis that is responsible for more than 780,000 deaths annually, making it the second most common cause of cancer-related mortality worldwide (Bray et al., 2018). Although liver cancer can be alleviated or cured through hepatectomy, orthotopic liver transplantation, and/or ablative procedures (Hanouneh, Alkhouri & Singal, 2019). More than 65% of patients fail to be cured and frequent recurrence contributes to poor survival. Predicting the overall 1-year survival rates for liver cancer have remained a challenge (Galle et al., 2018; Hanouneh, Alkhouri & Singal, 2019). New biomarkers that can predict liver cancer recurrence are urgently required to improve prognosis.
Liver cancer is influenced by post-transcriptional mechanisms that dynamically regulate protein expression (Wong, Tsang & Ng, 2018; Goldstrohm, Hall & McKenney, 2018). Cis-regulatory RNA elements and trans-acting factors (Gerstberger, Hafner & Tuschl, 2014; Moore, 2005) including RNA-binding proteins (RBPs) play an essential role in gene expression in cancer cells (Masuda & Kuwano, 2019). Recently, a group of RBPS termed MEX-3 RNA binding family member (MEX3) was identified in the nematode Caenorhabditis elegans (Ciosk, DePalma & Priess, 2006) and revealed one of the few RBPs with carcinogenic or tumor suppressor activity (Kim, Hur & Jeong, 2009; Pereira et al., 2013a).
MEX3 proteins are evolutionarily conserved RNA-binding proteins that consist of four homologous genes (MEX3A–D) (Buchet-Poyau et al., 2007; Courchet et al., 2008; Pereira et al., 2013a). They contained highly conservative one carboxy-terminal RING finger module and two K homology domains, the former mediating E3 ubiquitin ligase activity, the latter providing RNA-binding capacity (Buchet-Poyau et al., 2007). Available evidence implicates the MEX3 family in epithelial homeostasis, embryonic development, metabolism, immune responses and cancer, but the specific mechanisms of these effects require elucidation (Pereira et al., 2013a). MEX3A is a member of the MEX3 family (also known as RKHD4 or RNF162) that is expressed in endometrium tissue and the ovaries. MEX3A is a novel component of GW-182 or Dcp-containing bodies in mammals that represent cellular sites of mRNA degradation and the sequestration of non-translated transcripts (Cougot, Babajko & Seraphin, 2004; Eystathioy et al., 2002; Sheth & Parker, 2003).
MEX3A mRNA was recently shown to be overexpressed in Wilms tumors (Krepischi et al., 2016), gastric cancer (Jiang et al., 2012), bladder cancer (Huang et al., 2017), and bladder urothelial carcinoma (Shi & Huang, 2017). MEX3A promotes cell proliferation in bladder (Huang et al., 2017) and gastric (Jiang et al., 2012) cancer and shows potential as a biomarker to predict carcinogenesis (Pereira et al., 2013a). In this study, we analyzed the expression of MEX3A in liver cancer and assessed its clinicopathological potential. We further investigated the potential of MEX3A as an independent predictor of liver cancer prognosis.
Clinical and RNA-Seq analysis
We downloaded the all RNA-Seq expression matrix from the Cancer Genome Atlas (TCGA) database and obtained MEX3A mRNA expression data from liver cancer vs. normal liver tissue using the matrix. We further obtained corresponding clinical and pathological information from TCGA database. The basic clinical data included age, gender, histological grade, stage, T∕N∕M classification and vital status. MEX3A mRNA expression were estimated as log2(x+1) values and transformed RSEM normalized counts.
We retrospectively analyzed all data using R (version 3.5.1) (R Core Team, 2009). We used non-parametric rank sum tests to analyze MEX3A mRNA expression levels according to different variables and boxplots were visualized. Wilcoxon rank sum tests were used to compare the two subgroups, including disease, age, gender and vital status. Kruskal–Wallis tests were used for the comparison of three or more subgroups, including clinical stage, histologic grade and T∕N∕M classification. We used the pROC package to draw ROC curves for the evaluation of MEX3A diagnosis through the calculation of AUC values and the measurement of optimal cutoff point to divide samples into high and low MEX3A expression groups (Robin et al., 2011). Further, correlation between clinical features and MEX3A expression group were analyzed through chi-square tests with Fisher’s exact test.
To evaluate prognosis, Kaplan–Meier curves were used based on log-rank tests to compare differences in survival status, including overall survival (OS) and relapse-free survival (RFS) between the high and low MEX3A expression groups using the survival package in R (Therneau, 1994; Therneau & Grambsch, 2000). Univariate Cox analysis was used to select factors associated with prognosis, with calculations of hazard ratios (HRs) and 95% confidence intervals (95% CIs). Independent prognostic values of OS and RFS in the patients were determined through Multivariate Cox analysis. P-values < 0.05 were deemed statistically significant.
Clinical characteristics and RNA-Seq analysis
A total of 423 tissue samples with MEX3A mRNA expression data, including 373 liver cancer and 50 normal liver tissues were obtained from the TCGA. All patients were diagnosed with primary liver cancer. Corresponding patient demographic and clinical characteristics such as age, gender, histologic grade, TNM stage, vital status, and radiation therapy were obtained. All patient data is shown in Table 1.
|Characteristics||Number of pts(%)|
NA, not available.
MEX3A is highly expressed in liver cancer
In liver cancer, MEX3A mRNA expression level was significantly up-regulated compared to normal tissues (p = 1.1e−14; Fig. 1A) and increased with higher histological grades (p = 0.00016; Fig. 1C). MEX3A expression level was significantly associated with vital status (p = 0.032; Fig. 1I) and age (p = 0.0011; Fig. 1G).
MEX3A as a liver cancer diagnostic.
ROC curve analysis showed that MEX3A had moderate diagnostic ability in patients with liver cancer (AUC=0.837; Fig. 2A). The diagnostic ability of MEX3A was comparable in all clinical stages (stage I: AUC = 0.823, stage II: AUC = 0.844, stage III: AUC = 0.835, stage IV: AUC = 0.888; Figs. 2B–2E).
Relationship between MEX3A expression and clinical characteristics
Figure S1 showed that the cutoff point was 7.266 when samples were divided according to high and low MEX3A expression. Chi-square tests were used to analyze clinical variables between the two groups, in which high MEX3A expression was associated with cancer related mortality (p = 0.001; Table 2). High MEX3A expression was also associated with a deterioration in liver tumor histopathology (p < 0.001; Table 2).
|Clinical characteristics||MEX3A expression|
|Variable||No. of patients||High||%||Low||%||χ2||p-value|
|Histological type||Fibrolamellar Carcinoma||3||0||0||3||1.17||1.8006||0.547|
% represents the distribution of different clinical features in the single MEX3A expression group P-value in bold represent significant clinical significance (p < 0.05).
MEX3A is an independent prognostic to evaluate the survival of liver cancer patients
Kaplan–Meier curves showed that patients with high MEX3A expression were more likely to have a poor OS (p < 0.0001; Fig. 3A). Further subgroup analysis showed that high MEX3A expression was associated with poor OS for all variables: stage I/II (p = 0.0011; Fig. 3B), stage III/IV (p = 0.00022; Fig. 3C), stage G1/G2 (p < 0.0001; Fig. 3D), stage G3/G4 (p = 0.044; Fig. 3E), male (p < 0.0001; Fig. 3F), female (p = 0.0066; Fig. 3G), younger (p = 0.0026; Fig. 3H), older (p = 0.00022; Fig. 3I). Through Univariate and Multivariate Cox analysis of OS, MEX3A (HR = 2.26, 95% CI [1.58–3.23], p < 0.0001) was identified as an independent risk factor for the prognosis of liver cancer along with T stage (p < 0.0001) and residual tumors (p = 0.026; Table 3).
|Univariate analysis||Multivariate analysis|
|Parameters||Hazard ratio||95%CI (lower∼upper)||P value||Hazard ratio||95%CI (lower–upper)||P value|
P-value in bold represent significant clinical significance (p < 0.05)
Based on the OS, we further explored the connection between RFS and MEX3A expression, and found that high MEX3A expression was associated with poor RFS (p < 0.0001; Fig. 4A). RFS was related to the expression of MEX3A for some variables, including stage I/II (p = 0.013; Fig. 4B), stage III/IV (p = 0.00027; Fig. 4C), stage G1/G2 (p < 0.0001; Fig. 4D), male (p = 0.0023; Fig. 4F), female (p = 0.0049; Fig. 4G), younger (p = 0.0025; Fig. 4H) and older (p = 0.0047; Fig. 4I). Univariate and Multivariate Cox analysis suggested that MEX3A (HR = 2.19, 95% CI [1.54–3.12], p < 0.0001) was an independent risk factor to evaluate the RFS for liver cancer along with T stage (p < 0.0001) and residual tumor status (p = 0.024; Table 4).
Liver cancer is one of the deadliest tumors worldwide. The main risk factors for liver cancer include hepatitis B or C virus infections, the intake of aflatoxin, alcohol abuse and non-alcoholic fatty liver disease (NAFLD) (Zhang, Yang & Ericsson, 2019). Due to the lack of effective diagnostic and prognostic evaluation methods, the mortality rates of liver cancer patients have gradually increased. New molecular markers that can guide prognosis and improve the survival rates of liver cancer patients are urgently required. Our team has devoted to exploring diagnostic and prognostic biomarkers in various cancers (Hou et al., 2019; Jiao et al., 2018; Jiao et al., 2019a; Jiao et al., 2019b; Jiao et al., 2019c; Jiao et al., 2019d; Li et al., 2019; Sun et al., 2019). In this study, MEX3A mRNA was identified as overexpressed in liver cancer tissue and could effectively evaluate the prognosis of liver cancer patients as an independent predictor. A strong correlation between high MEX3A expression and liver malignancy was also observed.
MEX3A is known to be upregulated in Wilms renal cancer (Krepischi et al., 2016), gastric cancer (Jiang et al., 2012), bladder cancer (Huang et al., 2017) and bladder urothelial cancer (Shi & Huang, 2017). These results are consistent with our finding that MEX3A mRNA is overexpressed in liver cancer (p = 1.1e−14). Area under the ROC curves was 0.837, which provided evidence that MEX3A was a potential biomarker for liver cancer diagnosis. Interestingly, the expression of MEX3A was higher with increased histological grade (p < 0.0001), suggesting that MEX3A is related to tumor progression.
The molecular mechanism underlying the oncogenic effects of MEX3A remain poorly understood. Huang et al. (2017) found that MEX3A silencing significantly inhibits the proliferation of bladder cancer cells and promotes apoptosis. Jiang et al. similarly reported that MEX3A silencing delays the cell cycle progression of gastric cancer cells. MEX3A silencing significantly inhibited cell migration and anchorage-independent growth (Jiang et al., 2012). In colorectal cells, MEX3A is a stemness-related gene (Barriga et al., 2017; Chatterji & Rustgi, 2018; Fernandez-Barral et al., 2019) that acts as a repressive factor through controlling the expression of CDX2. CDX2 inhibits colorectal tumor cells growth, invasion, progression and migration and plays an essential regulatory role in intestinal homeostasis (Bonhomme et al., 2003; Brabletz et al., 2004; Gross et al., 2008; Pereira et al., 2013b). If MEX3A is overexpressed in colonic cell lines, cell polarity and differentiation become impaired leading to carcinogenesis (Pereira et al., 2013b). Combined with our findings, the role of MEX3A in cancer progression explains its clinical links to poor histological grade and poor patient prognosis in liver cancer. It is therefore necessary to explore the relationship between MEX3A and survival.
The MEX3 family shows promise as a biomarker for both cancer progression and prognosis (Pereira et al., 2013a). However, in studies by Huang et al. (Shi & Huang, 2017), no prognostic significance for cancer bladder urothelial carcinoma samples following MEX3A overexpression were observed. These results contrasted our findings and may highlight differential roles of MEX3A in cancer progression according to cancer-type. Of note, upon assessment of each sub-variable group, Kaplan–Meier curves revealed that the liver cancer patients with high MEX3A expression had a poor OS. RFS were evaluated according to MEX3A expression and showed a similar relationship with OS, excluding G3/G4 group. This highlights the unique superiority of MEX3A expression for the assessment of liver cancer survival. In particular, as shown in Tables 3 and 4, MEX3A may represent a prognostic marker for liver cancer survival (OS and RFS) under the strong confounding effects of clinicopathological features, providing useful references for clinicians to aid the development of individualized patient’s treatments.
In summary, this is the first study to report the association between MEX3A mRNA and the clinical characteristics and survival of liver cancer patients. MEX3A has great potential to predict the prognosis of liver cancer patients. Future studies should explore the mechanisms by which MEX3A promotes liver cancer in vivo and in vitro. Further clinicopathological information and corresponding clinical tissue samples should be obtained to further validate these findings and to establish MEX3A as a novel prognostic for patients with liver cancer.
|Univariate analysis||Multivariate analysis|
|Parameters||Hazard ratio||95%CI (lower∼upper)||P value||Hazard Ratio||95%CI (lower–upper)||P value|
P-value in bold represent significant clinical significance (p < 0.05).
This is the first study to investigate the expression of MEX3A mRNA in liver cancer, revealing its association with specific clinical features. Moreover, our results indicate that MEX3A plays a significant role in the prognosis of liver cancer and can be used as an independent factor to predict liver cancer progression.