Thyroid hormone (T3) and its receptor (TR) are involved in cancer progression. While deregulation of long non-coding RNA (lncRNA) expression has been detected in many tumor types, the mechanisms underlying specific involvement of lncRNAs in tumorigenicity remain unclear. Experiments from the current study revealed negative regulation of BC200 expression by T3/TR. BC200 was highly expressed in hepatocellular carcinoma (HCC) and effective as an independent prognostic marker. BC200 promoted cell growth and tumor sphere formation, which was mediated via regulation of cell cycle-related genes and stemness markers. Moreover, BC200 protected cyclin E2 mRNA from degradation. Cell growth ability was repressed by T3, but partially enhanced upon BC200 overexpression. Mechanistically, BC200 directly interacted with cyclin E2 and promoted CDK2–cyclin E2 complex formation. Upregulation of cell cycle-related genes in hepatoma samples was positively correlated with BC200 expression. Our collective findings support the utility of a potential therapeutic strategy involving targeting of BC200 for the treatment of HCC.
Thyroid hormone (T3) and its receptor (TR) are involved in cell growth, metabolism, autophagy and cancer progression (Wu et al. 2013, Chi et al. 2016). Two TR genes, TRα1 and TRβ1, have been identified on chromosomes 17 and 3, respectively. The receptors function as ligand-dependent transcriptional factors through binding to specific regions known as thyroid hormone response elements (TREs) that are usually located in the promoter region of target genes. T3/TR has been shown to reduce tumor formation in various cancer types, including hepatocellular carcinoma (HCC) (Chi et al. 2013) and breast cancer (Martinez-Iglesias et al. 2009), suggestive of a tumor suppressor function.
According to the tumor-initiating cell (TIC) concept, a subset of cancer cells possesses stem cell features that are indispensable for tumor formation (Plaks et al. 2015). TICs are generally characterized by their capacity for self-renewal and generate progeny to create the tumor bulk. Accumulating evidence supports the involvement of TICs in perpetuation of various cancers, including liver (Zhu et al. 2016), breast (Nadal et al. 2013), prostate (Vander Griend et al. 2008), brain (Brescia et al. 2013) and colon (Jing et al. 2015). In liver cancers, CD133/Prominin-1, a transmembrane hematopoietic stem cell antigen, has been identified as a putative marker of TICs. Notably, CD133 promotes tumorigenic capacity in cancer stem cells through activation of PI3K/AKT/mTOR (Xia & Xu 2015) or β-catenin signaling pathways (Mak et al. 2012).
Long non-coding RNAs (lncRNAs) are a class of non-protein coding transcripts longer than 200 nucleotides that regulate complex cellular functions, such as cell growth, differentiation, metabolism and metastasis (Rinn & Chang 2012). Recently, dysregulation of many HCC-related lncRNAs such as lncRNA-ATB (Yuan et al. 2014), NEAT1 (Mang et al. 2017), HOTAIR (Gao et al. 2016) and MALAT1 (Malakar et al. 2017) have been identified. Brain cytoplasmic RNA 1 (BCYRN1 or BC200), hereafter referred to as BC200, is overexpressed in several tumor types, including esophageal squamous cell carcinoma (Zhao et al. 2016), breast (Iacoangeli et al. 2004) and lung cancer (Hu & Lu 2015). However, the mechanisms underlying functional impairment and specific involvement of BC200 in HCC remain to be established.
Here, we aimed to elucidate the involvement of specific deregulated lncRNAs and target genes mediated by T3/TR in tumor formation. Experiments from the current study revealed BC200 as a target gene downregulated by T3/TR. BC200 promoted cell growth and oncogenic sphere formation of hepatocarcinoma cells through upregulation of cell cycle-related genes. Furthermore, T3/TR signaling inhibited cell growth through suppression of BC200. Our results collectively demonstrate novel associations among T3/TR, BC200, cyclin E2 and cyclin-dependent kinases (CDK) 2 that are involved in regulation of the tumor formation.
Materials and methods
The human hepatoma cell lines, HepG2, Hep3B, SK-Hep1 (obtained from American Type Culture Collection), Huh7 (gift from Dr. T.Y Hsieh, Tri-Service General Hospital, Taiwan) (Shiu et al. 2013) and J7 (gift from Dr. C S Yang, National Taiwan University, Taiwan) (Chen et al. 2002) were cultured in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% (v/v) fetal bovine serum (FBS). Cells were grown at 37°C in a humidified atmosphere of 95% air and 5% CO2. The cell lines were authenticated using the Promega StemElite ID System, which is a short tandem repeat-based assay (Supplementary Fig. 1, see section on supplementary data given at the end of this article). All experiments were conducted with cells from passage numbers 5–20. Serum was depleted of T3, as described previously (Chen et al. 2008).
Human Disease-Related lncRNA Profiler (System Biosciences, Mountain View, CA, USA) consisting of 83 lncRNAs was used. Briefly, total RNA was isolated from HepG2-TRα1 cells treated with/without T3 for 24 h and two paired HCC specimens. The detection of lncRNA was performed according to the manufactory’s instructions. Values are expressed as log 2-transformed relative fold increase or decrease in lncRNA expression, relative to that in without T3 treatment or adjacent non-tumorous tissues after normalization to the internal control. A positive log 2-transformed fold-change indicates higher expression in T3 treatment and tumor specimens, whereas a negative value signifies relatively decreased expression. A Student’s t-test was performed to identify significantly and differentially expressed lncRNAs with a fold-change ≥2.0 or ≤0.5, P < 0.05.
Human hepatoma specimens
Hepatoma samples from Taiwan Liver Cancer Network (TLCN) were selected for study and subjected to qRT-PCR and western blot analyses (Chang et al. 2016). These analyses were carried out under informed consent. The protocol was approved by the Medical Ethics and Human Clinical Trial Committee at Chang Gung Memorial Hospital (Institutional Review Board, No: 103-4866B).
Quantitative reverse transcription-PCR (qRT-PCR)
To quantify lncRNA transcripts, total RNA was extracted from cells with TRIzol reagent kit (Life Technologies Inc.) and converted into cDNA using reverse transcriptase (Life Technologies). qRT-PCR was conducted in a total reaction volume of 15 µL containing forward and reverse primers and 1X SYBR Green mix (Applied Biosystems). The ABI Prism 7500 Fast Real-Time PCR system (Life Technologies) was employed for qRT-PCR analysis. The primer sequences were listed in Supplementary Table 1.
RNA in situ hybridization (RNA-ISH)
In situ detection of BC200 expression was performed on formalin-fixed paraffin-embedded HCC samples. BC200 anti-sense RNA probes were labeled with digoxigenin (DIG) using a DIG RNA labeling kit (Roche). Before the labeling reaction, BC200 plasmid was digested with restriction enzymes for linearization. Further, in vitro transcription from plasmid was carried out using SP6 RNA polymerase. After de-paraffinization and rehydration, samples were treated with proteinase K (20 μg/mL) to digest tissues before hybridization, which was conducted at 37°C for 30 min. Sections were prehybridized in buffer containing 4xSSC and 50% deionized formamide for at least 10 min at 37°C and hybridized overnight with DIG-labeled probe at 62°C in solution containing 40% deionized formamide, 10% dextran sulfate, 1× Denhardt’s solution, 4× SSC, 10 mM DTT and 1 µg/mL yeast t-RNA. After stringent washing, signals were detected using anti-DIG-AP antibodies (1:800 dilution, Roche) and NBT/BCIP substrate (Roche).
The immunoblot procedure was performed as described previously (Wu et al. 2011). Antibodies specific for cyclin E1, cyclin E2 and CDK2 (Santa Cruz Biotechnology Inc.), p21 (Thermo Fisher Scientific Inc.), p27 (Sigma-Aldrich), CD133 (Miltenyi Biotec, Auburn, CA, USA), CD44, Sox2, Nanog (GeneTex, Inc., Irvine, CA, USA) and GAPDH (Merck Millipore) were used. Band intensities were calculated using Image Gauge software (Fujifilm, Tokyo, Japan). The PVDF membranes were reprobed for different antibodies after using TOOLStripping Buffer (BIOTOOLS CO., LTD. Taiwan). The signal intensities of expression of target genes were normalized to those of GAPDH.
Establishment of overexpression and knockdown stable cell lines
For ectopic expression of BC200 ncRNA, the sequences were amplified via PCR and cloned into pLKO-TRbib1 vector. BC200 shRNAs (shBC200#1 and shBC200#2) were designed and cloned into pLKO-TRbib1 vector. The shRNA sequences were shown below:
Single shRNA plasmid and virus package plasmids (pCMV-ΔR8.91 and pMD.G) were co-transfected into 293FT, and the virus harvested after 72 h of transfection. A pool of stably infected cells was selected in medium containing puromycin.
Cells (1 × 105) were grown on 6 cm dishes. At the indicated time-points, cell growth rates were determined with trypan blue exclusion and quantitated using the LUNA Automated Cell Counter.
Soft agar assay
Stable cells (10,000/well) were seeded in 12-well plates for layer agar cultures. Cells were resuspended in 0.33% agar in DMEM containing 10% FBS, and the upper layer replaced once a week in complete medium. After 3 weeks, colonies were stained with 0.01% (w/v) crystal violet. Images of all plates were obtained under a microscope (Olympus IX71), and colony numbers scanned and counted with Image J.
RNA immunoprecipitation assay
The RNA immunoprecipitation (RIP) assay was performed as described previously (Hsieh et al. 2014). Antibodies for RIP assay against cyclin E2 (Santa Cruz) were used.
Cell sorting and flow cytometry
Cell sorting via flow cytometry was performed on HCC cells using PE-conjugated monoclonal mouse anti-human CD133/1 (Miltenyi Biotec). Isotype control mouse IgG1k-PE (eBioscience, San Diego, CA, USA) served as the negative control.
Sphere formation assay
For the sphere formation assay, cells (1 × 103) were plated to the ultralow attachment 6-well plate (Corning Inc., Corning, NY, USA) with DMEM/F12 medium containing 20 ng/mL basic fibroblast growth factor, 20 ng/mL epidermal growth factor and B27. Two weeks later, the whole image of each well was taken and the number of spheres larger than 100 μm was manually counted using ImageJ software.
Model I: Five-week-old male nude mice were subcutaneously injected with BC200-overexpressing Hep3B cells (5 × 106, n = 3 per group) and BC200-depleted SK-Hep1 cells (2 × 106, n = 4 per group) to assess the effects of the lncRNA on tumor formation ability. Tumor volumes (mm3) were measured using the formula: (W2 × L)/2 (W, smallest diameter; L, longest diameter). All animals were killed at the indicated time-points after tumor inoculation. Model II: Hyper- and eu-thyroid mice (n = 8 per group) were generated and samples preparations were described in our previous report (Chi et al. 2016). Briefly, the effects of T3 repressed diethylnitrosamine (DEN)-induced liver injury were determined in C57BL/6 male mice. Mice were treated with T3 (10 μg/100 g body weight) before intraperitoneal injection of DEN (100 mg/kg). All groups of mice received continued T3 after DEN injection. Animal experiments were approved by the committee and performed according to the guidelines of the Chang Gang Institutional Animal Care and Use Committee Guide for the Care and Use of Laboratory animals (CGU13-106).
Results are presented as means ± s.d. of three independent experiments. Statistical analysis was performed with SPSS version 15 software (SPSS Inc., Chicago, IL, USA) using the Mann–Whitney test for comparison of two groups and one-way ANOVA, followed by Tukey’s post hoc test, for two or more groups. Kaplan–Meier survival curves were employed to analyze survival outcomes. Overall survival (OS) with death as an event was analyzed using the log-rank test. P values <0.05 were considered significant.
To identify T3/TR-related lncRNAs that are differentially expressed in HCC, lncRNA expression profile screening was performed using the SYBR Green-based qRT-PCR array in TR-overexpressing HepG2 cell lines and HCC specimens. Candidate lncRNAs that were simultaneously downregulated by T3/TR and upregulated in HCC were selected for further study, leading to the identification of BC200 (Supplementary Fig. 2A). Expression levels of BC200 decreased by 50–80% after incubation of HepG2-TRα1 and HepG2-TRβ1 with T3 for 24–72 h (Fig. 1A). In contrast, BC200 expression levels were only marginally regulated by T3 in HepG2-neo cells (Fig. 1A). In other hepatoma cell lines, including J7-TRα1 and Huh7 expressing exogenous and endogenous TR protein, T3 suppressed BC200 expression (Fig. 1B), similar to results obtained with HepG2-TR cells. These findings indicate that T3/TR signaling negatively regulates BC200 expression in a dose- and time-dependent manner in vitro. Promoter activity analysis was further performed to determine whether BC200 is regulated by T3 at the transcriptional level. Putative negative TREs (nTRE) in the BC200 promoter region (I-VI), based on the sequences from nTRE-thyroid-stimulating hormone (TSH) (Nakano et al. 2004) or nTRE-Nm23-H1 (Lin et al. 2000), were identified with bioinformatics tools (Supplementary Fig. 2B). In the presence of T3, luciferase activity was decreased by 30–40% in constructs I, II, III and IV. Repression of luciferase activity was alleviated upon deletion of region d (−305/−289) in constructs V and VI (Supplementary Fig. 2B), indicating that potential nTRE of the BC200 promoter is located within this region. TR binding to promoter fragment VI of BC200 in vivo was validated using the ChIP assay. Specific binding of TR to the BC200 promoter fragment was confirmed based on pulldown with the TR antibody, but not mopc21 (immunoglobulin) antibody (Supplementary Fig. 2C).
BC200 expression is upregulated in HCC
qRT-PCR and ISH were performed to validate the profiling results. Notably, BC200 expression was significantly upregulated in HCC tissues compared with benign and adjacent normal tissues (Fig. 1C). Consistently, the ISH assay revealed high expression of BC200 in HCC (Fig. 1D). Furthermore, BC200 expression in HCC was positively correlated with tumor type, tumor size, vascular invasion and pathological stage (Table 1). Kaplan–Meier survival analysis showed association of high BC200 expression with poor OS rate in HCC patients (Fig. 1E). Multivariate Cox proportional hazard regression analysis revealed that BC200 is an independent prognostic factor associated with survival (Table 2). These data clearly support a critical role of BC200 in hepatocarcinogenesis.
Clinicopathological correlations of BC200 in HCC specimens.
|Parameters||Cases (n = 240)||Mean ± s.e.||Pa|
|<65||148||7.971 ± 1.838||0.4403|
|≥65||92||5.122 ± 0.6504|
|Male||131||4.340 ± 0.4438||0.1600|
|Female||109||3.721 ± 0.3843|
|No||142||5.891 ± 1.106||0.3239|
|Yes||98||8.310 ± 2.354|
|Low||107||5.059 ± 1.400||0.1404|
|Medium||49||5.485 ± 1.438|
|High||84||10.1 ± 2.654|
|NBNC||32||5.206 ± 1.452||0.9081|
|HBV||122||7.812 ± 1.890|
|HCV||81||6.335 ± 1.859|
|HBV & HCV||5||3.629 ± 1.146|
|Solitary||181||5.460 ± 0.9396||0.0322|
|Multiple||59||11.23 ± 3.715|
|<5 cm||137||5.172 ± 1.192||0.0343|
|≥5 cm||103||9.150 ± 2.183|
|No||128||4.090 ± 0.6392||0.01|
|Yes||112||10.07 ± 2.351|
|I||116||4.185 ± 0.6977||0.0376|
|II||71||7.686 ± 2.102|
|III||53||11.69 ± 4.126|
|1||11||2.173 ± 0.5413||0.2778|
|2||166||5.898 ± 1.047|
|3||63||10.28 ± 3.444|
Univariate and multivariate Cox regression analyses for overall survival.
|Variable||Univariate analysis||P||Multivariate analysis||P|
|HR (95% CI)||HR (95% CI)|
|BC200 (high vs low)||1.838 (1.156–2.922)||0.01||1.618 (1.011–2.589)||0.045|
|Age (≥50 vs <50)||1.154 (0.767–1.736)||0.493||–||–|
|Gender (male vs female)||0.912 (0.608–1.366)||0.654||–||–|
|Cirrhosis (yes vs no)||1.206 (0.804–1808)||0.366||–||–|
|AFP (high vs medium vs low)||1.310 (1.046–1.640)||0.019||1.080 (0.858–1.360)||0.51|
|Tumor size (≥5 vs <5 cm)||2.717 (1.796–4.108)||<0.001||1.705 (1.035–2.810)||0.036|
|Tumor type (multiple vs solitary)||2.135 (1.415–3.223)||<0.001||1.170 (0.627–2.181)||0.622|
|Vascular invasion (present vs absent)||1.824 (1.519–2.191)||<0.001||1.561 (1.124–2.166)||0.008|
|Pathological stage (III vs II vs I)||2.176 (1.695–2.794)||<0.001||1.043 (0.549–1.980)||0.898|
|Grading (3 vs 2 vs 1)||1.147 (0.777–1.694)||0.49||–||–|
BC200 promotes cell growth and transformation in vitro and in vivo
The precise functions of the majority of lncRNAs are still unknown. The co-expression network analysis has been used to predict function of unknown non-coding RNA since genes regulated by the specific regulator or a set of genes with the same function would be co-expressed (Lee et al. 2004, Gong et al. 2016, Shang et al. 2016). This type of approach may be effectively used for exploring the roles of lncRNAs in cancer progression (Lee et al. 2004). The Oncomine database (
BC200 is required for self-renewal maintenance of liver cancer stem cells
Our findings suggest that BC200 is responsible for conferring growth advantages at an early stage of tumor development. As specified earlier, TICs have the capacity to self-renew and regenerate new tumors. To determine the association between TIC behavior and BC200, flow cytometry of the CD133+ TIC population in HCC cell lines was performed. The data revealed the presence of the CD133+ population in HepG2 and Huh7 cells, with expression ranging from 8 to 68% (Supplementary Fig. 4A). HepG2 and Huh7 cell lines were further sorted into CD133- and CD133+ populations and confirmed via western blot (Supplementary Fig. 4B). BC200 was upregulated in CD133+ HepG2 cells (Supplementary Fig. 4C). To determine whether BC200 functions in liver TIC self-renewal, the sphere formation assay was performed. Knockdown of BC200 dramatically reduced the tumor sphere formation frequencies of CD133+ subsets of HepG2 and Huh7 cells, compared with sh-luc cells (Fig. 2D). To ascertain the role of T3/TR in liver TIC self-renewal, CD133+ HepG2 cells were infected with Adenovirus-TRα1 (Ad-TRα1) or Adenovirus-TRβ1 (Ad-TRβ1) and tumor sphere formation with/without T3 examined. T3 clearly suppressed tumor sphere formation in CD133+ HepG2 cells overexpressing Ad-TRα1 or Ad-TRβ1 relative to non-treated cells (Fig. 2E). In view of these results, we propose that BC200 is required for maintaining stemness ability in HCC.
CDKs and cyclin-dependent kinase inhibitors play important roles in control of cell growth (Abbas & Dutta 2009, Malumbres 2014). Deregulation of these proteins is a hallmark of several diseases, including cancers. To assess the contributory role of BC200 in cell growth, its effects on the mRNA levels of cell cycle-related genes were examined. The p21 transcript levels were increased in SK-Hep1 cells after knockdown of BC200 expression while mRNA levels of CDK2, cyclin E1 and cyclin E2 were reduced (Fig. 3A, left panel). Conversely, overexpression of BC200 led to opposite effects (Fig. 3A, right panel). As expected, levels of the corresponding proteins were consistent with mRNA results, both in vitro and in vivo (Fig. 3B and C). In addition, BC200 depletion in CD133+ HepG2 cells led to reduced expression of pluripotent transcription factors (Nanog and Sox2), CD44 and cell cycle-related genes, compared with control (sh-luc) cells (Fig. 3D). To confirm the involvement of BC200 in T3/TR-regulated cell growth, HepG2-TRβ1 cells were transfected with BC200-expressing or control (vc) plasmid and proliferation examined in the presence or absence of T3. T3 clearly repressed HepG2-TRβ1 cell growth. Interestingly, this effect was partially rescued by BC200 in the presence of T3 (Fig. 3E). Notably, protein levels of cyclin E2 and CDK2 were clearly increased in BC200-overexpressing HepG2-TRβ1 cells treated with T3, compared with their T3-treated control counterparts (Fig. 3E, lanes 2 vs 4). Our results collectively demonstrate that T3/TR suppresses BC200 expression, leading to modulation of cell cycle-related genes, and ultimately, inhibition of cell proliferation.
BC200 protects cyclin E2 mRNA from degradation
To clarify the regulatory mechanism of cyclin E2, its expression was measured using qRT-PCR after treatment with actinomycin D (ActD, a transcriptional inhibitor). Knockdown of BC200 promoted cyclin E2 mRNA degradation in SK-Hep1 cells (Fig. 4A), suggesting a protective effect of BC200 against degradation. To determine the functional sequence of BC200, we established stable expression of BC200-wt (1–200), BC200-F1 (1–119), BC200-F2 (122–200) and BC200-F3 (51–119) in BC200 knockdown cell lines. Overexpression of BC200-wt and BC200-F1 partially delayed cyclin E2 mRNA degradation, compared to the control group, while overexpression of BC200-F2 and BC200-F3 had no effect on the stability of these mRNAs (Fig. 4B). These findings suggest that a functional sequence of 1–50 nucleotides participates in regulation of cyclin E2 mRNA stability. Previously, cyclin E-CDK2 complex has been considered an essential and master regulator of progression through G1 phase of the cell cycle (Koff et al. 1992). To test this hypothesis, co-immunoprecipitation assay was performed in BC200 stable cell lines. We found BC200 promoted association with CDK2-cyclin E2, but not CDK2-cyclin E1 (Fig. 4C). These findings indicated that the critical role of BC200 in cyclin E2-CDK2 activity in HCC cells. Next, the physical interaction of BC200 and cyclin E2 was examined using RIP assay. The results showed a specific enrichment of BC200 (but not GAPDH mRNA) coprecipitated within the cyclin E2 (Fig. 4D), indicating BC200 associated with cyclin E2 through direct binding and enhanced CDK2–cyclin E2 complex formation. Furthermore, levels of the p27 and p21 were decreased in cyclin E2-overexpressing cells compared to control cells (Fig. 4E).
Cyclin E2 and CDK2 are highly expressed in HCC and positively correlated with BC200
T3 protects hepatocyte from DEN-induced HCC-like liver tumors in mice model (Chi et al. 2016). To further determined T3 represses this phenotype through BC200/cyclin E2/CDK2 cascade, the expression levels of these genes were determined in DEN/hyper- and DEN/euthyroid mice. The expression levels of BC200 and cell cycle-related genes were repressed in DEN/hyperthyroid group, compared with DEN/euthyroid group (Fig. 5A and B). Our results suggest that T3 suppresses DEN-induced hepatocarcinogenesis through BC200/cyclin E2/CDK2 cascade. To determine whether cyclin E2 and CDK2 are associated with tumor formation, mRNA levels of these molecules were determined in HCC specimens using qRT-PCR. The markers were significantly expressed in HCC compared with the corresponding adjacent non-tumor counterparts (Fig. 5C and Supplementary Fig. 5A), consistent with high expression in HCC reported earlier from two publicly available datasets (Supplementary Fig. 5B and C). Furthermore, both TRα1 and TRβ1 were downregulated in HCC relative to the normal control group (Supplementary Fig. 5B and C). Notably, cyclin E2 and CDK2 expression were significantly correlated with OS of HCC patients (Supplementary Fig. 6A), similar to findings from the public dataset (Supplementary Fig. 6B). Spearman’s correlation coefficient analysis revealed significant positive correlation of BC200 with cyclin E2 and CDK2 (Supplementary Table 2). Considering the positive correlation between these markers in HCC tumors, their combined influence on patient outcomes were evaluated. Patients were classified into three groups (I–III) using the median value as cutoff. Patients in group III (higher expression of both genes) showed significantly poorer OS than those in group I (lower expression of both genes) (Fig. 5D). Our findings clearly support an important regulatory function of the BC200/cyclin E2/CDK2 axis in hepatoma (Fig. 5E).
Previous studies have confirmed that lncRNAs and transcriptional factors form feedback or feed-forward loops, which play critical roles in biological processes. For example, the MYC-induced long non-coding RNA (MINCR) is activated by MYC and modulates the MYC transcription network (Doose et al. 2015). LincRNA-p21 is upregulated by HIF1α and essential for hypoxia-enhanced glycolysis (Yang et al. 2014). In view of these associations, it is proposed that deregulation of lncRNAs by transcriptional factors could cause disease or cancer progression. Therefore, transcriptional factor-mediated lncRNA regulation is a significant focus of cancer biology research. In the current study, BC200 was identified as a novel gene downregulated by T3/TR that modulates thyroid hormone-mediated functions.
BC200 is reported to be expressed at high levels in invasive carcinoma of the breast (Iacoangeli et al. 2004). Earlier receiver-operating characteristics analysis of sensitivity and specificity confirmed the diagnostic power of BC200 RNA as a molecular marker of invasive breast cancer. However, the functional and clinical significance of BC200 in HCC remains to be established. Here, we showed that BC200 promotes cancer cell growth, both in vitro and in vivo and is positively correlated with tumor type, tumor size, vascular invasion and pathological stage, supporting its utility as an independent prognostic factor associated with OS. Additionally, BC200 is known to act as a translational regulator in neurons (Iacoangeli & Tiedge 2013). Data from the current study support a novel role of BC200 in regulation of cyclin E2 mRNA stability. Functional sequences of BC200 1–50 nucleotides in length were responsible for regulating cyclin E2 mRNA stability. On the other hand, BC200 expression appeared positively correlated with the cell proliferation rate. These results were similar to previous findings (Booy et al. 2017). Interestingly, BC200 was upregulated in CD133+ HepG2 cells and stimulated stemness marker expression (Sox2 and Nanog), in turn, promoting sphere-formation capacity. Accordingly, we concluded that BC200 regulates early-stage HCC formation.
In our experiments, several genes that are positively or negatively regulated by T3 were identified in the HepG2-TRα1 cell line, including endoglin (Lin et al. 2013), pituitary tumor-transforming 1 (PTTG1) (Chen et al. 2008), Dickkopf 4 (DKK4) (Liao et al. 2012), Ubiquitin-like with PHD and ring finger domains 1 (UHRF1) (Wu et al. 2015) and Death-associated protein kinase 2 (DAPK2) (Chi et al. 2016), proposed to suppress tumor formation. Previously, hypothyroidism was shown to be associated with significantly elevated risk for HCC (Hassan et al. 2009). A recent study by our group suggests that the T3/PTEN-induced kinase 1/Parkin pathway plays an important role in protecting hepatocytes from HBx-induced HCC (Chi et al. 2017). Catalano and co-workers demonstrated that T3 reduces the tumorigenic potential of colorectal cancer-CSC to a significant extent through regulating Wnt and BMP signaling (Catalano et al. 2016). In contrast, the group of Wang showed that the thyroid hormone (T4) promotes cell self-renewal in HCC cells through TRα, but not TRβ (Wang et al. 2016). Kress and colleagues reported that TRα1 overexpression and Wnt pathway activation enhance colorectal cancer formation (Kress et al. 2010). Specifically, the group demonstrated that TRα1 promotes tumorigenesis in APC+/1638N mice with a Wnt-activated genetic background but TRα1 overexpression alone is unable to trigger tumor formation. Thus, T3/TRs may perform a dual function as an oncogene or tumor suppressor in different genetic backgrounds.
In conclusion, we have made conceptual progress toward advancing the theory that T3/TR suppresses tumor growth and tumor sphere formation via reduction of BC200. Mechanistically, BC200 promotes cyclin E2 expression via regulation of mRNA stability. The identification of a novel pathway interlinking T3/TR, BC200, cyclin E2 and CDK2, which regulates proliferation and tumor sphere formation of hepatoma cells, presents potential therapeutic strategies involving targeting of BC200 and associated molecules for treatment of HCC.
This is linked to the online version of the paper at
Declaration of interest
The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.
This work was supported by grants from Chang Gung Memorial Hospital, Taoyuan, Taiwan (CMRPD1D0383; CMRPD1G0421, CMRPD1G0422, CRRPD1F0011, CRRPD1F0012, CRRPD1F0013, NMRPD1G0951 to K-H L) and from the Ministry of Science and Technology of the Republic of China (MOST 103-2320-B-182-017-MY3; 106-2320-B-182-031-MY3, 106-2320-B-182-032-MY3 to K-H L; MOST 106-2321-B-182-003-MY3 to Y-H L).
Author contribution statement
Conception and design of the study: Yang-Hsiang Lin, Meng-Han Wu, Kwang-Huei Lin; Acquisition of data: Yang-Hsiang Lin, Meng-Han Wu, Ya-Hui Huang, Hsiang-Cheng Chi, Wen-Yu Chuang, Chia-Jung Yu, Ching-Ying Chen; Analysis and interpretation of data: Yang-Hsiang Lin, Ya-Hui Huang, Chau-Ting Yeh, Chung-Ying Tsai, I-Hsiao Chung; Drafting of the manuscript: Yang-Hsiang Lin, Meng-Han Wu, Kwang-Huei Lin; Statistical analysis: Yang-Hsiang Lin, Ya-Hui Huang, Chau-Ting Yeh; Study supervision: Kwang-Huei Lin.
The authors would like to thank Taiwan Liver Cancer Network (TLCN) for providing the hepatoma tissue samples and related clinical data (all are anonymous) for our research work. They thank Chuen Hsueh for helping to improve RNA-ISH detection.
CatalanoVDenticeMAmbrosioRLuongoCCarolloRBenfanteATodaroMStassiGSalvatoreD 2016 Activated thyroid hormone promotes differentiation and chemotherapeutic sensitization of colorectal cancer stem cells by regulating Wnt and BMP4 signaling. Cancer Research 76 1237–1244. (https://doi.org/10.1158/0008-5472.CAN-15-1542)
ChangICHuangSFChenPJChenCLChenCLWuCCTsaiCCLeePHChenMFLeeCM, 2016 The hepatitis viral status in patients with hepatocellular carcinoma: a study of 3843 patients from Taiwan Liver Cancer Network. Medicine 95 e3284. (https://doi.org/10.1097/MD.0000000000003284)
ChenCJTsaiNMLiuYCHoLIHsiehHFYenCYHarnHJ 2002 Telomerase activity in human hepatocellular carcinoma: parallel correlation with human telomerase reverse transcriptase (hTERT) mRNA isoform expression but not with cell cycle modulators or c-Myc expression. European Journal of Surgical Oncology 28 225–234. (https://doi.org/10.1053/ejso.2001.1237)
DooseGHaakeABernhartSHLopezCDuggimpudiSWojciechFBergmannAKBorkhardtABurkhardtBClaviezA, 2015 MINCR is a MYC-induced lncRNA able to modulate MYC’s transcriptional network in Burkitt lymphoma cells. PNAS 112 E5261–E5270. (https://doi.org/10.1073/pnas.1505753112)
KressESkahSSirakovMNadjarJGadotNScoazecJYSamarutJPlaterotiM 2010 Cooperation between the thyroid hormone receptor TRalpha1 and the WNT pathway in the induction of intestinal tumorigenesis. Gastroenterology 138 1863–1874. (https://doi.org/10.1053/j.gastro.2010.01.041)
MakABNixonAMKittanakomSStewartJMChenGICurakJGingrasACMazitschekRNeelBGStagljarI, 2012 Regulation of CD133 by HDAC6 promotes beta-catenin signaling to suppress cancer cell differentiation. Cell Reports 2 951–963. (https://doi.org/10.1016/j.celrep.2012.09.016)
MalakarPShiloAMogilevskyASteinIPikarskyENevoYBenyaminiHElgavishSZongXPrasanthKV, 2017 Long noncoding RNA MALAT1 promotes hepatocellular carcinoma development by SRSF1 upregulation and mTOR activation. Cancer Research 77 1155–1167. (https://doi.org/10.1158/0008-5472.CAN-16-1508)
MangYLiLRanJZhangSLiuJLiLChenYLiuJGaoYRenG 2017 Long noncoding RNA NEAT1 promotes cell proliferation and invasion by regulating hnRNP A2 expression in hepatocellular carcinoma cells. OncoTargets and Therapy 10 1003–1016. (https://doi.org/10.2147/OTT.S116319)
Martinez-IglesiasOGarcia-SilvaSTenbaumSPRegaderaJLarcherFParamioJMVennstromBArandaA 2009 Thyroid hormone receptor beta1 acts as a potent suppressor of tumor invasiveness and metastasis. Cancer Research 69 501–509. (https://doi.org/10.1158/0008-5472.CAN-08-2198)
NadalROrtegaFGSalidoMLorenteJARodriguez-RiveraMDelgado-RodriguezMMaciaMFernandezACorominasJMGarcia-PucheJL, 2013 CD133 expression in circulating tumor cells from breast cancer patients: potential role in resistance to chemotherapy. International Journal of Cancer 133 2398–2407. (https://doi.org/10.1002/ijc.28263)
NakanoKMatsushitaASasakiSMisawaHNishiyamaKKashiwabaraYNakamuraH 2004 Thyroid-hormone-dependent negative regulation of thyrotropin beta gene by thyroid hormone receptors: study with a new experimental system using CV1 cells. Biochemical Journal 378 549–557. (https://doi.org/10.1042/bj20031592)
RinnJLChangHY 2012 Genome regulation by long noncoding RNAs. Annual Review of Biochemistry 81 145–166. (https://doi.org/10.1146/annurev-biochem-051410-092902)
ShiuTYHuangTYHuangSMShihYLChuHCChangWKHsiehTY 2013 Nuclear factor kappaB down-regulates human UDP-glucuronosyltransferase 1A1: a novel mechanism involved in inflammation-associated hyperbilirubinaemia. Biochemical Journal 449 761–770. (https://doi.org/10.1042/BJ20121055)
WuSMChengWLLiaoCJChiHCLinYHTsengYHTsaiCYChenCYLinSLChenWJ, 2015 Negative modulation of the epigenetic regulator, UHRF1, by thyroid hormone receptors suppresses liver cancer cell growth. International Journal of Cancer 137 37–49. (https://doi.org/10.1002/ijc.29368)