JARID2 promotes invasion and metastasis of hepatocellular carcinoma by facilitating epithelial-mesenchymal transition through PTEN/AKT signaling

JARID2 is crucial for maintenance of pluripotency and differentiation of embryonic stem cells. However, little is known about the role of JARID2 in metastasis of hepatocellular carcinoma (HCC). This study found that JARID2 expression was significantly higher in HCC tissues than that in adjacent non-tumor liver tissues (ANLTs), and its expression level correlated with HCC metastasis. High JARID2 expression was significantly correlated with multiple tumor nodules, high Edmondson-Steiner grade, microvascular invasion, advanced TNM stage and advanced BCLC stage (all P < 0.05) and indicated poor prognosis of HCC in training and validation cohorts (all P < 0.05) totaling 182 patients. High JARID2 expression was an independent and significant risk factor for disease-free survival (DFS; P = 0.017) and overall survival (OS; P = 0.041) after curative liver resection in training cohort, and also validated as an independent and significant risk factor for DFS (P = 0.033) and OS (P = 0.031) in validation cohort. Moreover, down-regulation of JARID2 dramatically inhibited HCC cell migration, invasion, proliferation in vitro and metastasis in vivo, whereas overexpression of JARID2 significantly increased migration, invasion, proliferation in vitro and metastasis in vivo. Mechanistically, the data showed that JARID2 exerted its function by repressing PTEN expression through increasing H3K27 trimethylation (H3K27me3) at PTEN promoter region, which subsequently resulted in activation of protein kinase B (AKT) and enhanced epithelial-mesenchymal transition (EMT). In conclusion, this study revealed that JARID2 promotes invasion and metastasis of HCC by facilitating EMT through PTEN/AKT signaling.


INtrODUctION
Hepatocellular carcinoma (HCC) is one of the most common cancers and ranks the second leading cause of cancer death in men and sixth in women over the world [1]. Liver resection is still considered as the prime choice and the most effective ''curative'' treatment for HCC, especially for solitary large HCC [2]. Despite the improvement in treatment of HCC during recent decades, the overall survival of patients with HCC remains unsatisfactory due to high rate of recurrence and metastasis after liver resection [2,3]. Thus, it is essential to explore effective biomarkers for recurrence and metastasis to improve treatment strategy for better clinical outcome. Tumor metastasis is a complex process involving a complicate succession of invasion-metastasis steps and associated with numerous molecular genetic mechanisms. Though several molecules associated with metastasis of HCC had been found by our research group and others, [4][5][6][7][8] the mechanisms underlying HCC metastasis are still elusive. Recently, some researches have illustrated that the epithelial-mesenchymal transition (EMT) is an important mechanism promoting tumor metastasis, including HCC [9,10]. EMT, a normal embryological process, is

Research Paper
Oncotarget 40267 www.impactjournals.com/oncotarget characterized by the loss of epithelial markers and the gain of mesenchymal markers [9]. HCC cells undergoing EMT acquire enhanced migratory and invasive properties, thereby resulting in HCC metastasis [11]. Therefore, it is necessary to elucidate the new regulatory mechanisms of EMT in HCC invasion and metastasis.
JARID2, one member of the Jumonji protein family, is crucial for maintenance of pluripotency and differentiation of embryonic stem cells [12,13]. Its deletion was previously shown to result in severe defects in liver development [14]. Abnormal JARID2 expression had been reported in rhabdomyosarcomas and leukemia, [15,16] and contributed to metastatic behavior of cancer cells by promoting epithelial-mesenchymal transition (EMT) [17]. JARID2 is also associated with leukemic transformation of chronic myeloid malignancies [16,18]. These data suggest that JARID2 functions to the malignant character in tumors and also led us to hypothesize that JARID2 may play roles in HCC pathogenesis, especially in HCC metastasis. However, little is known about the function of JARID2 in the development of HCC. So we attempted to study the roles of JARID2 in the invasion and metastasis of human HCC.
Phosphatase and tension homolog (PTEN), a wellstudied tumor suppressor phosphatase, plays a central role in regulating tumor growth and metastasis [19,20]. It has been reported that PTEN mutation or deletion frequently occurs in various human cancers and is associated with cancer progression [20,21]. However, increasing studies showed that downregulation of PTEN in HCC, rather than PTEN mutations or deletions, is believed to contribute to HCC metastasis, [11,22] suggesting other mechanisms in HCC may be responsible for the downregulation of PTEN. Recently, it was shown downregulation of PTEN was also associated with trimethylation of H3 lysine 27 (H3K27me3) induced by polycomb repressive complex 2 (PRC2) [23]. JARID2 is essential for the binding of PRC2 to the targeted genes, leading to silence of gene expression by H3K27me3, [13] However, whether JARID2 could interact with the PRC2 complex to regulate PTEN expression in HCC remains unknown. Therefore, we also examined molecular mechanisms whether JARID2 exerted functions by regulating PTEN expression in HCC.
In this study, the data indicated that JARID2 was highly expressed in HCC tissues relative to the adjacent non-tumor liver tissues (ANLTs), which correlated with HCC metastasis and predicted poor prognosis of HCC patients. JARID2 promoted epithelial-mesenchymal transition (EMT)-mediated metastasis of HCC by suppressing PTEN and activating of AKT signaling.

JARID2 is significantly upregulated in HCC
First, quantitative real-time polymerase chain reaction (qRT-PCR) was used to detect JARID2 expression in 30 paired fresh HCC samples. The results showed that, compared with matched adjacent non-tumor liver tissues (ANLTs), the expression level of JARID2 in HCC tissues (T) was significantly up-regulated (fold change (T/ANLT) > 2) in 73.3% cases (22/30) (Figure 1A). RT-PCR ( Figure 1B) and western blotting ( Figure 1C) showed the similar results as qRT-PCR in these matched specimens. Interestingly, among three different subtypes of HCC including small hepatocellular carcinoma (SHCC , the diameter of HCC ≤ 5 cm), solitary large hepatocellular carcinoma (SLHCC, only one nodule, and diameter > 5 cm, and grows expansively within an intact capsule or pseudocapsule), nodular hepatocellular carcinoma (NHCC, has more than 2 nodules), [2,24] JARID2 expression level in NHCC with the greatest metastatic potential and poorest clinical outcome was significant higher than that in SLHCC and SHCC with relatively better clinical outcome [2,24] (Figure 1D1). But the median of JARID2 expression in SHCC was similar to that in SLHCC ( Figure 1D1). JARID2 was significantly higher in tumors with microvascular invasion (MVI) than tumors without MVI ( Figure 1D2). Patients with metastasis and/or recurrences of HCC also exhibited higher JARID2 mRNA expression than who without metastasis and/or recurrences ( Figure 1D3). To confirm the PCR and western blot results, immunohistochemical staining was performed to detect the JARID2 in paraffin-embedded paired HCC samples from training cohort and validation cohort (Supplementary Figure 1 and Supplementary Table 1). The results showed JARID2 protein was primarily localized to the nucleus in HCC ( Figure 1E1, Figure 2A) and was highly expressed in 59.4% and 57.6% (69/116; 38/66; Table 1) of HCC tissues, as compared with 12.9% and 3.0% (15/116; 2/66; Supplementary Table 2) in corresponding ANLTs. Its expression was also found to be significantly upregulated in HCC tissues and correlated with HCC subtypes in both training cohort ( Figure 1E2) and validation cohort ( Figure 1E3). Above of all, these data show that JARID2 is overexpressed in HCC and may contribute to HCC metastasis.

High JARID2 expression associates with the clinicopathologic features for HCC
Next, we sought to explore the association of JARID2 expression with the clinicopathologic features for HCC in two independent cohorts-training and validation cohort from two research centers (Supplementary Figure 1  and Supplementary Table 1). In training cohort, high JARID2 expression in HCC tissues positively correlated with tumor number (P = 0.002), microvascular invasion (P = 0.001), Edmondson-Steiner grade (P = 0.004), HCC subtype (P = 0.007), TNM stage (P = 0.042) and BCLC stage (P = 0.049) (  ANLTs, the expression level of JARID2 was significantly up-regulated in HCC tissues. The value of JARID2 mRNA in ANLTs was used to normalize (value set to 1) the value of JARID2 mRNA in HCC tissues. Red line indicates fold change of JARID2 equal to 2. T, HCC tissue; ANLT, adjacent nontumorous liver tissue. (b) The representative images of semiquantitative RT-PCR showed that JARID2 mRNA in HCC tissues was significantly higher than those in ANLTs. (c) Western blot results showed that JARID2 protein in HCC tissues was significantly higher than those in ANLTs. (D) JARID2 mRNA expression level was detected by qRT-PCR in NHCC, SLHCC and SHCC. The expression of JARID2 in NHCC was significantly higher than in SLHCC or SHCC, but there was no significant difference between SLHCC and SHCC (D1). JARID2 mRNA expression level was detected by qRT-PCR in HCC tissues with MVI and those without MVI (D2) and in HCC tissues with metastasis and/or recurrence and those without metastasis and/or recurrence (D3). SHCC, small hepatocellular carcinoma; SLHCC, solitary large hepatocellular carcinoma; NHCC, nodular hepatocellular carcinoma; MVI, with microvascular invasion; NMVI, without microvascular invasion; MR, with metastasis and/or recurrence; NMR, without metastasis and/or recurrence. (E) JARID2 expression in training cohort (n = 116) and in validation cohort (n = 66). The JARID2 expression was determined by immunohistochemistry (IHC) and calculated by percent of positive cells. Representative IHC images of normal liver (NL), ANLT, HCC tissues (T) (E1). The expression of JARID2 in ANLTs, HCC tissues (T) and in HCC subtypes (SHCC, SLHCC and NHCC) in training cohort (E2) and validation cohort (E3). *P < 0.05; **P < 0.01; ***P < 0.001.

JARID2 promotes invasion and metastasis of HCC Cells in vitro and in vivo
The difference of JARID2 expression was examined in a normal liver cell line (L02 cells) and five HCC cell lines with various metastatic potential by qRT-PCR (Supplementary Figure 2A1 Figure 3B2). Conversely, JARID2 overexpression significantly enhanced motility of HepG2 cells by wound healing assay ( Figure 3A3) and increased invasive ability of HepG2 cells by transwell assays ( Figure  3B3). Immunofluorescence staining of F-actin was used to analyze cytoskeleton. The results showed that the stress fiber-like structures disappeared in HCCLM3 shJARID2 cells compared to HCCLM3 control cells ( Figure Figure 4E, 4F). These findings implicate that JARID2 functions as an oncogene and is involved in HCC cell invasion.
Next, the effect of JARID2 on HCC invasion and metastasis was analyzed by in vivo metastatic assays. The mouse xenograft orthotopic liver cancer model assays showed that the growth of tumors (in terms of mean tumor volume) originally formed from HCCLM3 shJARID2 cells were significantly smaller than that from HCCLM3 Control cells ( Figure 3D1), and the mean tumor volume of the MHCC97-H shJARID2 group was also significantly smaller than that of the MHCC97-H Control group ( Figure 3D2), but the mean tumor volume of the HepG2 JARID2 group was also significantly larger than that of the HepG2 Control group ( Figure 3D3). Consistently, the number of metastatic nodules in the liver ( Figure 3E1) or lung ( Figure 3E2) in mice implanted with JARID2-knockdown cells was significantly smaller than that mice implanted with cells transfected with the control vector ( Figure 3F1, 3F2). However, the number of metastatic nodules in the liver (Supplementary Figure 6A) or lung (Supplementary Figure 6B) in mice implanted with JARID2-overexpressed cells was significantly larger than that mice implanted with cells transfected with the control vector ( Figure 3F3, 3F4). These data suggest that JARID2 promotes tumor formation and metastasis.

JARID2 promotes HCC invasion and metastasis by inhibiting PTEN expression through H3K27 trimethylation
To define the mechanisms by which JARID2 promotes HCC invasion and metastasis, we tried to identify potential targeted genes regulated by JARID2. It was reported that JARID2 physically associated with PRC2 and was functionally involved in trimethylation of histone H3 lysine 27 (H3K27me3) [12,13]. Likewise, the interaction between JARID2 and PRC2 was further validated by coimmunoprecipitation (co-IP) in HCC cells (Supplementary Figure 7A-7C). To further explore JARID2 is functionally together with PRC2 in HCC, we next investigated whether JARID2 knockdown or overexpression would affect H3K27me3 expression in HCC cells. Western blot showed JARID2 knockdown decreased the level of H3K27me3 expression, whereas JARID2 overexpression increased the level of H3K27me3 expression ( Figure 4A). However, H3K4me3 was not affected ( Figure 4A).
PRC2 could inhibit the transcription of PTEN, [23] and JARID2 is functionally together with PRC2 in HCC. Notably, PTEN plays a critical role in HCC proliferation, invasion and metastasis, [20] supporting the function of JARID2 in regulating proliferation, invasion and metastasis. These also suggest that JARID2 possibly coordinates with PRC2 to exert functions by inhibiting PTEN expression in HCC. As expected, JARID2 knockdown in HCCLM3 or MHCC97-H cells led to increased expression of PTEN at both the transcriptional level ( Figure 4B1, 4B2) and the protein level ( Figure 4B3). On the contrary, JARID2 overexpression in HepG2 cells resulted in decreased expression of PTEN at both the transcriptional level ( Figure 4B1, 4B2) and the protein level ( Figure 4B3). Furthermore, the expression level of JARID2 showed an inverse trend with that of PTEN in HCC samples, as measured by qRT-PCR ( Figure 4C1) and western blotting ( Figure 4C2). These data support the observation that PTEN are suppressed by JARID2 at the transcriptional level. Notably, quantitative chromatin immunoprecipitation (qChIP) assays showed the inhibition of JARID2 expression led to a strong reduction of JARID2 binding to PTEN promoter region and, consistent with this, to a reduction of H3K27me3 levels at PTEN promoter region ( Figure 4D1), whereas JARID2 overexpression led to a dramatic increase of JARID2 binding to PTEN promoter region and, similarly, to an increase of H3K27me3 levels at PTEN promoter region ( Figure 4D2). These findings indicate that JARID2 represses the expression of PTEN through increasing H3K27me3 at PTEN promoter region.
To determine whether PTEN is a mediator for JARID2-promoted invasion and metastasis in HCC, a complementary approach of gain-and loss-off function of PTEN was used. PTEN was overexpressed in HCCLM3 shcontrol (Supplementary Figure 7D

JARID2 promotes metastasis of HCC via facilitating EMT through activating PTEN/AKT Signaling Pathway
HCCLM3 cells with high metastatic capability exhibit high expression of EMT markers [27]. Interestingly, it was observed that HCCLM3 cells exhibited high JARID2 expression (Supplementary Figure 2) and knockdown of JARID2 resulted in the cell appearance with more cobble-like, epithelial morphology in HCCLM3 ( Figure 3C1) and MHCC97-H cells ( Figure 3C2), whereas overexpression of JARID2 changed the cell appearance to more spindle-like, fibroblastic morphology in HepG2 cells ( Figure 3C3). Thus, these indicated that JARID2 might promote metastasis of HCC via EMT.
To further confirm whether JARID2 promotes EMT, immunofluorescence (IF) analyses showed knockdown of JARID2 in HCCLM3 and MHCC97-H significantly reduced the expression level of mesenchymal marker vimentin, but significantly increased the expression level of epithelial marker E-cadherin ( Figure 5A). Conversely, overexpression of JARID2 in HepG2 cells significantly increased the expression level of mesenchymal marker vimentin, but dramatically reduced the expression level of epithelial marker E-cadherin ( Figure 5A). Similar with IF results, western blotting also showed that JARID2 knockdown in HCCLM3 and MHCC97-H decreased the expression of vimentin and increased the expression of E-cadherin ( Figure 5B1). However, JARID2 overexpression in HepG2 increased the expression of vimentin and decreased the expression of E-cadherin ( Figure 5B2). Moreover, the changes of these EMT markers and JARID2 in protein expression were also observed in transplanted orthotropic tumors of mice (Supplementary Figure 9A and 9B). Furthermore, IHC showed high JARID2 expression was co-location with low E-cadherin and high vimentin expression in HCC samples ( Figure 5C). A similar correlation of the expression of E-cadherin, vimentin with the expression of JARID2 was also observed in the HCC tissues (Supplementary Table 3). These results suggest that JARID2-promoted HCC invasion and metastasis is associated with EMT.
Previous studies have demonstrated that PTEN negatively regulated AKT signaling pathway, which is important for EMT [28,29]. Western blotting was performed to examine whether JARID2 suppressed PTEN expression and subsequently activated AKT signaling to facilitate EMT. Our data showed that inhibition of JARID2 reduced phosphorylation level of AKT ( Figure 5B1), and was associated with decreased vimentin expression and increased E-cadherin expression ( Figure 5B1). In contrast, overexpression of JARID2 up-regulated phosphorylation level of AKT ( Figure 5B2), and was associated with increased vimentin expression and decreased E-cadherin expression ( Figure 5B2). Furthermore, western blotting showed suppression of PTEN in JARID2-knockdown cells rescued AKT activities of HCC cells ( Figure 5B1), whereas up-regulation of PTEN abolished the AKT activities induced by JARID2 ectopic expression ( Figure  5B2). Moreover, IHC showed that HCC tissues with higher JARID2 expression, weaker PTEN expression exhibited stronger staining of p-AKT ( Figure 5C). The similar results were also observed in transplanted orthotropic tumors of mice (Supplementary Figure 9A abd 9B). The correlation analysis revealed JARID2 expression was negatively correlative with PTEN expression, but positively correlative with p-AKT level in HCC sample (Supplementary Table 3). Furthermore, PTEN expressions reversely correlated with p-AKT level, vimentin expression and positively with E-cadherin in HCC sample (Supplementary Table 4). Taken together, overexpression and their corresponding control by qPCR (B1), semiquantitative RT-PCR (b2), western blot (b3). (c) The expression level of JARID2 exhibited an inverse trend with that of PTEN in HCC samples. The expression of JARID2 and PTEN was detected in 30 paired HCC tissues and ANLTs. (C1) The relative level of JARID2 mRNA expression was plotted against the relative level of PTEN mRNA expression. (c2) Western blot detected the expression of JARID2 and PTEN protein in the same HCC tissues and ANLTs. (D) JARID2 was associated with increased H3K27me3 at PTEN promoter region. Quantitative chromatin immunoprecipitation (qChIP) was used to assess H3K27me3 at PTEN promoters region in HCCLM3, MHCC97-H cells with shJARID2 knockdown, HepG2 cells with JARID2 overexpression and their corresponding control vector. qChIP was performed using the indicated antibodies. Purified rabbit IgG was a negative control. The percentage fold difference of shJARID2/control, Vector/JARID2 was calculated as the difference in percentage of the qChIP in HCCLM3, MHCC97-H cells transfected with shJARID2 relative to cells transfected with control vector, in HepG2 cells transfected with control vector relative to cells transfected with JARID2 vector respectively. Each bar represents the mean ± SEM for triplicate experiments. (E, F) Gain-and loss-off function studied with PTEN expression vector, shPTENs or its control showed PTEN is a critical downstream effector in JARID2-promoted invasion and metastasis in HCC. HCC cells were transfected with shRNAs for PTEN or PTEN expression vector to inhibit or restore the PTEN expression. Wound-healing assays (E) and transwell assays (F) were determined in HCCLM3, MHCC97-H, HepG2 cells. **, P < 0.01; ***, P < 0.001.
Oncotarget 40278 www.impactjournals.com/oncotarget these data indicate that JARID2 promotes EMT by suppressing PTEN expression and subsequently activating AKT signaling.

DIscUssION
Recurrence and metastasis remains the main reason for poor prognosis of HCC after curative liver resection, [3,30] and the molecular mechanisms underlying HCC recurrence and metastasis have not been thoroughly understood. Although JARID2 was studied in the malignancies such as rhabdomyosarcoma and leukemia, [15,17] it remains unclear whether and how JARID2 contributes to HCC progression. In this study, we found that JARID2 was significantly upregulated in human HCC tissues compared with ANLTs. qRT-PCR showed that JARID2 mRNA was significantly higher in primary tumors with microvascular invasion (MVI) than those without MVI. Patients with metastasis and/or recurrence had higher JARID2 mRNA expression than those without metastasis and/or recurrence. Interestingly, NHCC with the greatest metastatic potential and poorest clinical outcome had the highest JARID2 expression in comparison to SLHCC and SHCC, which had relatively better clinical outcome based on our previous clinical studies [2,24]. The differential expression of JARID2 among three subtypes is an another important evidence to the diversity of molecular biological characters and consolidate our evidence that differentially expressed gene profiles between SLHCC and NHCC [31]. IHC also demonstrated that high Oncotarget 40279 www.impactjournals.com/oncotarget JARID2 protein expression was significantly correlated with multiple tumor numbers, MVI, high Edmondson-Steiner grade, advanced TNM stage and BCLC stage. The presence of MVI in surgical specimens is known to be an important cause of HCC recurrence after resection, and generally have a poor prognosis [32]. Kaplan-Meier analysis showed HCC patients with high JARID2 expression had worse clinical outcomes than those with low JARID2 expression. A multivariate analysis revealed that high JARID2 expression in HCC tissues was also found to be a significant and independent prognosis factor for DFS and OS. Moreover, high JARID2 expression in HCC tissues was tested to be a good predictor for HCC early recurrence (within 2 years), which is typically caused by dissemination of metastatic HCC cells [33]. This indicates that we may identify the subgroup of HCC patients who are at high risk of developing metastatic disease in advance according to JARID2 expression in HCC tissues. The predictive value of prognosis based on JARID2 expression level further validated in validation cohort, indicating the anticipated availability as a novel molecular biomarker. These data also suggest that JARID2 contributes to HCC invasion and metastasis.
A recent study reported that JARID2 promoted invasion of A549 lung cancer cell line and HT29 colon cancer cell line by the in vitro assays [17]. Here we showed that JARID2 overexpression was associated with metastatic potential of HCC cells. Downregulation of endogenous JARID2 expression in HCC cells suppressed cell migration and invasion, while ectopic JARID2 expression significantly promoted cell migration and invasion in vitro. Moreover, the in vivo metastatic assays showed that downregulation of JARID2 expression significantly decreased intrahepatic and lung metastasis, whereas JARID2 overexpression obviously increased intrahepatic and lung metastasis in vivo. Above of all, our data comprehensively shows that JARID2 promotes HCC invasion and metastasis.
The underlying mechanism of JARID2 in promoting HCC invasion and metastasis need to be further elucidated. In this study, we showed that PTEN was the key downstream effector of JARID2. PTEN is often down-regulated in cancers and is associated with advanced stages of cancers or metastasis, [21,34] including HCC. [22,35] Our data uncovered that JARID2 was associated with increasing H3K27me3, not H3K4me3, at the PTEN promoter region in HCC, and this resulted in down-regulation of PTEN expression, which consequently promotes HCC cell invasion and metastasis. By gain-and loss-off function assays, we also exhibited that overexpression of PTEN in high JARID2 expression HCC cells attenuated JARID2-promoted invasion and metastasis, while depletion of PTEN in JARID2knockdown HCC cells could mimick the roles of JARID2. Thus, these results show that JARID2 transcriptionally represses PTEN expression through increasing H3K27me3 at the gene promoter region, and resulting in increase in HCC invasion and metastasis.
An increase in migration and invasion ability is an important character of EMT, which is essential for tumor cells to disseminate to adjacent or distant tissues, and finally causes metastasis-related recurrence [36]. It has been shown that PTEN plays a critical role in the EMT by negative activation of downstream signaling AKT [11,22]. In HCC, activation of AKT signaling is frequently observed due to loss of PTEN, [20,35] which triggers EMT program and endows cancer cells with enhanced metastatic potential [28,29]. Therefore, it is not difficult to understand that JARID2-mediated suppression of PTEN can activate downstream AKT signaling, and enhance cancer cell invasion ability by inducing EMT in HCC cells. Moreover, this study also proved that JARID2 downregulated PTEN expression, which increased the activated phosphorylation of AKT, consequently enhanced vimentin expression and decreased E-cadherin expression in HCC cells. Notably, PTEN knockdown or overexpression abolished the effects of JARID2 knockdown or overexpression on EMT. Further correlation analysis in HCC samples by IHC showed JARID2 expression was positively associated with p-AKT, vimentin and negatively associated with PTEN, E-cadherin. Furthermore, HCC tissues with low PTEN expression, exhibited higher p-AKT level, higher vimentin expression, but lower E-cadherin expression. All together, these observations support that JARID2 promotes EMTinduced metastasis by suppressing PTEN and activating AKT signaling.
In conclusion, our study shows that JARID2 is significantly upregulated in HCC and high JARID2 expression is associated with poor prognosis, and EMT phenotype in HCC tissue and HCC cell lines. JARID2 promotes HCC invasion and metastasis by facilitating EMT through PTEN/AKT signal pathway. Our data suggest JARID2 functions as a potential oncogene, supporting the pursuit of JARID2 as a prospective therapeutic target for HCC.

HCC Samples
The study was approved by the Ethics Committee of Xiangya School of Medicine, CSU. From January 2009 to December 2010, 30 pairs of frozen fresh tumor liver tissues and corresponding ANLTs were collected after radical surgical resection at Department of surgery, Xiangya Hospital of Central South University. These tissues were used to detect the mRNA and protein expression of JARID2 and PTEN. In addition, normal liver tissues were obtained from 5 patients with hepatic hemangioma for hepatic resection. Another two independent cohorts of subjects including training cohort www.impactjournals.com/oncotarget (n = 116) and validation cohort (n = 66) from 2 different centers were used for prognostic study according to REMARK guidelines for reporting prognostic biomarkers in cancer [37]. In training cohort, paraffin-embedded paired HCC samples (including HCC tissues and ANLTs) were obtained from HCC patients undergoing radical surgical resection without any preoperative treatment at Department of Surgery, the Xiangya Hospital of Central South University (CSU) from January 2005 to December 2008. In validation cohort, paraffin-embedded paired HCC samples (including HCC tissues and ANLTs) were obtained from HCC patients who underwent radical surgical resection without any preoperative treatment at Department of Abdominal Surgical Oncology, Affiliated Cancer Hospital of Xiangya School of Medicine of CSU from January 2003 to December 2008. All samples randomly selected were histopathologically diagnosed and had completed clinicopathologic and follow-up data. Clinicopathologic characteristics of the patients in training cohort and validation cohort used for prognostic study were shown in Supplementary Table 1 and Supplementary  Figure 1. HCC samples were divided into three subgroups as we previously defined: [2,24] small hepatocellular carcinoma (SHCC, the diameter of HCC ≤ 5 cm), solitary large hepatocellular carcinoma (SLHCC, only one nodule, and diameter > 5 cm, and grows expansively within an intact capsule or pseudocapsule), nodular hepatocellular carcinoma (NHCC, has more than 2 nodules).

Prognostic study
All HCC patients were regularly followed-up by the same experienced surgical team. The follow-up period was defined as the interval between the date of operation and that of the patient's death or the last followup. Deaths from other causes were treated as censored cases. The recurrence and metastasis was surveillance by clinical examination, serial monitoring of alphafetoprotein levels and ultrasonography or CT scan or MRI at a 3-4 months' interval. Recurrence and metastasis were diagnosed by clinical examination, serial AFP level, and ultrasonography or computed tomography (CT) scan. Disease-free survival was defined as the length of time after liver resection during which a patient survived without sign of HCC. Data of conventional clinical and pathological variables were also collected for analysis, including age, gender, hepatitis B status, liver cirrhosis, Edmondson-Steiner grade, capsular formation, size of the tumor, number of tumor nodes, microvascular invasion, Child-Pugh classification, TNM stage and BCLC stage. Microvascular invasion was defined as tumor cells forming a thrombus in peritumoral vessels, can only be assessed after careful histological assessment of the whole surgical specimen [38,39]. The follow-up status and any recurrence were regularly updated in the database for each patient. The study was approved by the Ethics Committee of Xiangya School of Medicine, CSU.

Cell lines
HepG2 cells were purchased from the American Type Culture Collection (ATCC, Rockville, MD). L02 and SMMC-7721 cells were gifted by the Tumor Institute of Central South University (Changsha, China). MHCC97-L, MHCC97-H and HCCLM3 cells were gifted from the Liver Cancer Institute of Fudan University (Shanghai, China). Short tandem repeat (STR) DNA fingerprinting was used to authenticate all cell lines before experiments. All cell lines were routinely cultured with the high glucose DMEM supplemented with 10% fetal bovine serum, and maintained in 5% CO 2 humidified incubator at 37°C.

Vector construction and transfection
The plasmids carrying short hairpin RNAs (shRNAs) for JARID2 knockdown, JARID2 expression vector inserted with JARID2 coding sequences (CDS) for JARID2 overexpression and corresponding controls were purchased from GeneChem Company (Shanghai, China). The sequences of five shRNAs for JARID2 knockdown and its control were shown in Supplementary Table 5. HCCLM3 and MHCC97-H cells were transfected with the shRNA plasmid, and HepG2 cells were transfected with the JARID2 expression plasmid. Cells transfected with control plasmids were used as controls. Transfected cells were selected with 3 µg/mL puromycin. Down-regulated expression or overexpression of JARID2 was confirmed by western blot (Supplementary Figure 2C, 2D, 2E). The inhibitory efficiency of five shRNAs was validated and the JARID2-shRNA-Seq3, JARID2-shRNA-Seq5 (named shJARID2-1, shJARID2-2 respectively in the figures) were adopted for subsequent study (Supplementary Figure 2D and 2E) because of highly effective inhibition of JARID2 expression in HCCLM3 and MHCC97-H cells.
PTEN ectopic expression and knockdown lentivirus as well as their negative control lentivirus were purchased from GeneChem (Shanghai, China). The sequences of shRNA for PTEN were seen in Supplementary Table 5. The PTEN expression vector was constructed by inserting its CDS into the vector under the regulation of a CMV promoter as this might counteract the effects of JARID2mediated silencing. The lentivirus was transfected into the HCCLM3, MHCC97-H and HepG2 cells with an optimal multiplicity of infection (MOI) of 50 TU/mL. To increase the number of infected cells, the number of lentivirus copies should increase in each HCC cells. Down-regulated expression (at least 80%) or overexpression of PTEN was confirmed by western blot (Supplementary Figure 7D, 7E, 7F). The inhibitory efficiency of three shRNAs for PTEN was validated and the shPTEN-1, shPTEN-3 were adopted for subsequent study.

Western blot analysis
Total proteins were extracted and separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and then transferred onto PVDF membranes (Millipore, Bedford, MA). The blotted membranes were incubated with the primary antibodies and then an appropriate HRP-conjugated secondary antibody (KPL, Gaithersburg, MD) in order. Band was detected with enhanced chemiluminescence regents (Thermo Scientific, Rockford, IL). Beta-actin protein was also determined by using the specific antibody (Sigma, St Louis, MO) as a loading control. Protein expression were quantified by BandScan software (BioRad, Hercules, CA) and defined as the ratio of target protein relative to Beta-actin. Antibodies for JARID2, p-AKT, AKT, vimentin, E-cadherin, PTEN and corresponding secondary antibodies were purchased from Santa Cruz Biotechnology (Santa Cruz Biotechnology, Santa Cruz, CA).

Immunohistochemistry
Paraffin-embedded tissues were sectioned and microwave-pretreated in EDTA buffer (1mM, pH 8.0) for 10 minutes for antigen retrieval. Then, formalin-fixed paraffin sections were stained for JARID2 (Santa Cruz Biotechnology) using the streptavidin-peroxidase system (Zhong-shan Goldenbridge Biotechnology, Beijing, China). To support the validity of the immunohistochemistry, HCCLM3 shcontrol, HCCLM3 shJARID2-1 cells were chosen used to perform cell immunohistochemistry to demonstrate the effectiveness of the antibodies JARID2 used, as PCR and Western blot showed HCCLM3 cells exhibited high JARID2 expression, and HCCLM3 transfected with the shJARID2-1 plasmid showed dramatically reduced JARID2 expression. Positive, negative control slides were probed with JARID2 antibody, goat serum respectively, and followed by the secondary antibody under the same conditions (Supplementary Figure 10A and 10B). The results proved that JARID2 antibody used was useful and specific for JARID2 (Supplementary Figure 10A 10B). The expression levels of JARID2 were scored using a four-point scale according to the percentage of positive hepatocytes: [40,41] 0, ≤ 10% positive; 1+, 11%-25% positive; 2+, 26%-50% positive; 3+, ≥ 51% positive. The protein expression of JARID2 was thus considered negative if scored 0, and 1+, 2+, and 3+ as positive. According to the score of JARID2 expression, HCC specimens was also divided into a low expression group (0 or 1+) and a high expression group (2+ or 3+). For PTEN, p-AKT, vimentin, E-cadherin, the procedures as JARID2 and immunostaining score were adopt as described elsewhere [42]. Antibodies for PTEN, p-AKT, vimentin, E-cadherin were all purchased from Santa Cruz Biotechnology (Santa Cruz, CA).

Cell proliferation and colony formation assays
For cell proliferation assay, 3 × 10 4 HCC cells expressing a short hairpin targeting JARID2 or the control hairpin were seeded into a 35mm dishes (Corning Costar Corp, Corning, NY) in triplicate and cultured in 5% CO 2 at 37°C. The cells were trypsinized and the number of cells was counted by using cell counter for seven consecutive days. For colony formation assays, 500 cells were seeded into 35 mm dishes (Corning Costar Corp) and cultured in 5% CO 2 for 2 weeks at 37°C. The number of colonies per dish was counted after staining with crystal violet. Only positive colonies (diameter > 40 um) in the dishes were counted and compared [29,43]. These experiments were performed in triplicate.

Wound healing assays
Wound healing assays were performed and analyzed as described elsewhere [44]. Cells were seeded into sixwell culture plates (Corning Costar Corp) coated with fibronectin. After cells reached 100% confluence, to suppress cell proliferation which could confound the analysis of cell migration, cells were preincubated with mitomycin (Sigma, St. Louis, MO, 10 μg/ml) for 1h at 37°C. Wound healing assays were performed with a sterile pipette tip to make a scratch through the confluent monolayer. Medium was changed and cells were cultured for 24 hours. The percent wound closure was calculated for three randomly chosen fields.

Transwell assays
For the transwell assays, about 1 × 10 5 cells in serum-free medium were placed into the upper chamber of the insert with matrigel plated (BD Biosciences, Franklin Lakes, NJ). After 24 hours of incubation in 5% CO 2 at 37°C, the cells in upper chamber were removed with cotton swabs, following fixed by 20% methanol, and then stained with a solution containing 0.1% crystal violet (Beyotime Institute of Biotechnology, Beijing, China). The number of cells that adhered to the lower membrane of the inserts was counted. For each experimental group, the assays were performed in triplicates, and five random fields were chosen for analysis. www.impactjournals.com/oncotarget Immunofluorescence Immunofluorescence was performed as described previously [45]. Cells were seeded into the 6-well culture plate (Corning Costar Corp) to prepare for performing cell immunofluorescence (IF). After incubating with primary antibodies, cells then incubated with corresponding fluorescence labeled secondary antibody. The slides were photographed using the inverted fluorescence microscope TE-2000S (Nikon, Tokyo, Japan). Primary antibodies for E-cadherin, vimentin were purchased from Santa Cruz Biotechnology (Santa Cruz, CA). Rhodamine-conjugated phalloidin, DAPI and fluorescence labeled secondary antibody were obtained from Beyotime Institute of Biotechnology (Shanghai, China).

Co-immunoprecipitation
Co-immunoprecipitation was performed as described previously [46]. For Co-immunoprecipitation assays, the cells were lysed with cold lysis buffer at 4°C for 30 min. Supernatants of cell lysates were then collected by centrifugation at 50,000 rpm for 10min, followed by incubation with an appropriate primary antibody overnight at 4°C. The immunocomplexes were next precipitated by protein A/G sepharose beads. After washes, the protein complexes were eluted out by sample buffer and then subjected to western blot analysis by probing with antibodies against EZH2 (Santa Cruz Biotechnology), SUZ12 (BD Biosciences), RbAp46/48 (BD Biosciences), EED (BD Biosciences) and JARID2 (Santa Cruz Biotechnology).

Chromatin immunoprecipitation (ChIP) assays
ChIP was performed using a ChIP assay kit (Millipore, Billerica, MA) according to manufacturer's instructions. Briefly, cells were cross-linked with 1% formaldehyde. The chromatin was sonicated into fragments ranging between 200 and 1000 bp and then was pulled down by antibody for real-time quantitative PCR amplification. The antibodies used were anti-JARID2 (Santa Cruz), anti-H3K27me3 (Santa Cruz). The primers for amplifying the fragments of the PTEN promoter are as follows: 5ʹ-CCGTGCATTTCCCTCTACAC-3ʹ (sense) and 5ʹ-GAGGCGAGGATAACGAGCTA-3ʹ (antisense).

In vivo assays for metastasis
For the in vivo metastasis assays, the HCC metastatic model in mice was constructed as previously described [29]. Briefly, 5 × 10 6 HCC cells were injected subcutaneously into the left upper flank regions of nude mouse (4 weeks of age, male, BALB/c). After 35 days, the subcutaneous tumors were removed and divided into commensurate fragments of approximately 1 mm 3 , then implanted into the liver of nude mouse (3 mice in each group). After 6 weeks, mice were killed, and all livers and lungs were harvested. The tumor size was calculated as follows: tumor volume (mm 3 ) = (L × W 2 )/2, [47] where L = long axis and W = short axis. All livers and lungs were fixed with 10% phosphate-buffered neutral formalin, sectioned serially and stained with hematoxylin and eosin (H&E) for histological examination. The expression of JARID2, PTEN, p-AKT, vimentin, E-cadherin in orthotropic tumor tissues was also determined by immunohistochemistry. Mice were performed and housed in the Animal Institute of CSU according to the protocols approved by the Medical Experimental Animal Care Commission.

Statistical analysis
All data were analyzed using the statistical software SPSS 18.0 for Windows (SPSS Inc., Chicago, IL). The differences between groups were analyzed by Student's t test between two groups or by one-way analysis of variance (ANOVA) in more than two groups when the variance is homogeneous. If the variance is not homogeneous, the differences between groups were analyzed by Mann-Whitney U test between two groups or by Kruskal-Wallis H test in more than two groups. χ 2 analysis was used to analyze the correlation between JARID2 expression and clinicopathologic features. Linear regression was used to analyze the relationship between JARID2 and PTEN mRNA expression. Spearman's rank analysis was used to analyze the correlations between different protein expressions level. Survival curves were constructed using the Kaplan-Meier method and evaluated using the log-rank test. The Cox proportional hazards regression model was established to identify factors which were independently associated with the overall survival (OS) and disease-free survival (DFS) of HCC patients. All the tests were two-tailed and P < 0.05 was considered statistically significant.

cONFLIct OF INtErEst
There are no financial conflicts of interest from any authors.