Inflammatory related gene IKKα, IKKβ, IKKγ cooperates to determine liver cancer stem cells progression by altering telomere via heterochromatin protein 1-HOTAIR axis

Cancer stem cells are associated with tumor recurrence. IKK is a protein kinase that is composed of IKKα, IKKβ, IKKγ. Herein, we demonstrate that IKKα plus IKKβ promoted and IKKγ inhibited liver cancer stem cell growth in vitro and in vivo. Mechanistically, IKKα plus IKKβ enhanced and IKKγ inhibited the interplay among HP1α, HP1β and HP1γ that competes for the interaction among HP1α, SUZ12, HEZ2. Therefore, IKKα plus IKKβ inhibited and IKKγ enhanced the activity of H3K27 methyltransferase SUZ12 and EZH2, which methylates H3K27 immediately sites on HOTAIR promoter region. Therefore, IKKα plus IKKβ increased and IKKγ decreased the HOTAIR expression. Strikingly, IKKα plus IKKβ decreases and IKKγ increases the HP1α interplays with DNA methyltransferase DNMT3b, which increases or decreases TERRA promoter DNA methylation. Thus IKKα plus IKKβ reduces and IKKγ increases to recruit TRF1 and RNA polymerase II deposition and elongation on the TERRA promoter locus, which increases or decreases TERRA expression. Furthermore, IKKα plus IKKβ decreases/increases and IKKγ increases/decreases the interplay between TERT and TRRRA/between TERT and TREC. Ultimately, IKKα plus IKKβ increases and IKKγ decreases the telomerase activity. On the other hand, at the telomere locus, IKKα plus IKKβ increases/drcreases and IKKγ decreases/increases TRF2, POT1, pPOT1, Exo1, pExo1, SNM1B, pSNM1B/CST-AAF binding, which keep active telomere regulatory genes and poised for telomere length. Strikingly, HOTAIR is required for IKKα plus IKKβ and IKKγ to control telomerase activity and telomere length. These observations suggest that HOTAIR operates the action of IKKα, IKKβ, IKKγ in liver cancer stem cells. This study provides a novel basis to elucidate the oncogenic action of IKKα, IKKβ, IKKγ and prompts that IKKα, IKKβ, IKKγ cooperate to HOTAR to be used as a novel therapeutic targets for liver cancer.


Research Paper
of an HP1-like protein regulates heterochromatin body assembly for DNA elimination [10]. Emerging evidence has shown that HP1α serves a biological role in cancer [11]. HP1 and the histone H3 lys9 methyltransferase Su(var)3-9 is requireed in stem cell self-renewal [12]. Loss of HP1 causes depletion of H3K27me3 and gain of H3K27me2 at constitutive heterochromatin [13]. Furthermore, HP1 regulates this gene's alternative splicing by recruiting splicing factors [14].
Accumulating evidence indicates that HOTAIR plays a critical role in cancer progression. HOTAIR could promote migration and invasion of hepatocellular carcinoma (HCC) cells by inhibiting RBM38 [15]. The Polycomb group (PcG) protein heterodimer EZH2-EED is necessary and sufficient for binding to the lncRNA HOTAIR [16]. Studies suggest that HOTAIR recruits chromatin-modifying complexes to specific target sequences and triggers stemness acquisition [17]. Moreover, Enforced expression of HOTAIR altered histone H3 lysine 27 methylation, gene expression, and increased cancer invasiveness in a manner dependent on Polycomb Repressive Complex 2 (PRC2) [18]. The PRC2 is composed of a trimeric core of SUZ12, EED and EZH1/2, and PRC2 catalyzes histone H3K27 trimethylation (H3K27me3), a hallmark of gene silencing [19]. H3K27 methylations contribute to the role of PRC2 in maintaining cellular growth [20].
TRF2 activity alteration is absorbing and of great concern. Our previous findings suggest that SET1A and CUDR increased TRF2 expression on the transcriptional and translational level through H3K4me3 [21]. The shelterin protein TRF2 is essential for chromosome-end protection. Dyefunction of TRF2 causes chromosome end-to-end fusions, initiating genomic instability. Elevated levels of TRF2 induce telomeric ultrafine anaphase bridges and rapid telomere deletions [22]. Accumulating evidence suggests that the human telomerase reverse transcriptase catalytic subunit (TERT) contributes to cell to elongate telomeres [23]. The telomeric long noncoding RNA Telomeric repeat-containing RNA (TERRA) is important for telomere regulation. TERRA G-quadruplex structure is critical for binding to telomeres [24]. Moreover, TERRA participates in the regulation of telomere length, telomerase activity and heterochromatinization [25,26]. Resaerches suggest the SNM1B/APOLLO DNA nuclease functions in resolution of replication stress and maintenance of common fragile site stability [27]. CST/ AAF, a DNA polα. primase accessory factor, binds POT1b and shortens the extended overhangs produced by Exo1 [28].
In this report, our findings suggest that IKKα plus IKKβ promoted and IKKγ inhibited liver cancer stem cell growth in vitro and in vivo. Mechanistically, IKKα plus IKKβ enhanced and IKKγ inhibited the interplay among HP1α, HP1β and HP1γ. Therefore, IKKα plus IKKβ increased and IKKγ decreased the TOTAIR expression.
Ultimately, IKKα plus IKKβ increases and IKKγ decreases the telomere length and telomerase activity. Of significance, HOTAIR is required for IKKα plus IKKβ and IKKγ to control telomere. This study provides a novel basis to elucidate the oncogenic action of IKKα, IKKβ, IKKγ.

IKKα, IKKβ, IKKγ regulate HOTAIR expression dependent on H3K27me3
To investigate whether IKKα, IKKβ, and IKKγ influence on HOTAIR expression, we first performed IP assay in liver cancer stem cells. As shown in Figure 4A, compared to the control, IKKα, IKKβ,

IKKβ+IKKγ,
IKKα+IKKγ, IKKα+IKKβ+IKKγ overexpression group, when IKKα and IKKβ were co-overexpressed, the interplay between HP1α and H3K27me1/2/3 was significantly decreased. Conversely, when single IKKγ was overexpressed, the interplay between HP1α and H3K27me3 was significantly decreased. Furthermore, the loadeing of H3K27me3 onto the HOTAIR promoter region was significantly decreased when IKKα and IKKβ were co-overexpressed. Conversely, when single IKKγ was overexpressed, the loading of H3K27me3 onto the HOTAIR promoter region was significantly decreased ( Figure 4B). However, the action was abolished when HP1α was knocked down in these liver cancer stem cells ( Figure  4C). As shown in Figure 4D, compared to the control, IKKα, IKKβ, IKKβ+IKKγ, IKKα+IKKγ, IKKα+IKKβ+IKKγ overexpression group, the HOTAIR promoter luciferase activity was significantly increased when IKKα and IKKβ were co-overexpressed. Conversely, the HOTAIR promoter luciferase activity was significantly decreased when single IKKγ was overexpressed. As shown in Figure 4E, when IKKα and IKKβ were co-overexpressed, the HOTAIR exptression was significantly increased. Conversely, when single IKKγ was overexpressed, the HOTAIR experssion was significantly decreased. Taken together, IKKα, IKKβ, and IKKγ regulate HOTAIR expression dependent on H3K27me3.

IKKα, IKKβ, and IKKγ control telomerase activity
To investigate whether IKKα, IKKβ, IKKγ influence on telomere, we first consider to analyse the with anti-HP1α followed by western blotting with, anti-H3K27me1, anti-H3K27me2, anti-H3K27me3 in liver cancer stem cells transfected with pcDNA3.1-IKKα, pcDNA3.1-IKKβ, pcDNA3.1-IKKγ respectively. IgG IP as negative control. Western blotting with anti-HP1 α as INPUT. B. Chromatin Immunoprecipitation (CHIP) with anti-H3K27me3 followed by PCR with HOTAIR promoter primers in liver cancer stem cells transfected with pcDNA3.1-IKKα, pcDNA3.1-IKKβ, pcDNA3.1-IKKγ respectively. IgG CHIP as negative control. PCR for HOTAIR promoter as INPUT. C. Chromatin Immunoprecipitation (CHIP) with anti-H3K27me3 followed by PCR with HOTAIR promoter primers in HP1 α depleted liver cancer stem cells transfected with pcDNA3.1-IKKα, pcDNA3.1-IKKβ, pcDNA3.1-IKKγ respectively. IgG CHIP as negative control. PCR for HOTAIR promoter as INPUT. D. HOTAIR promoter luciferase activity assay in liver cancer stem cells transfected with pcDNA3.1-IKKα, pcDNA3.1-IKKβ, pcDNA3.1-IKKγ, respectively. Data represent mean±SEM, n=3, **p<0.01, *p<0.05 represents difference significane. E. RT-PCR with HOTAIR primers in liver cancer stem cells transfected with pcDNA3.1-IKKα, pcDNA3.1-IKKβ, pcDNA3.1-IKKγ respectively. β-actin as internal control. alteration of telomerase activity. As shown in Figure  5A, compared to the control, IKKα, IKKβ, IKKβ+IKKγ, IKKα+IKKγ, IKKα+IKKβ+IKKγ overexpression group, the interplay between RNAPolIII and TRF1, DNMT3b and HP1α were significantly increased when IKKα and IKKβ were co-overexpressed,. Conversely, the interplay between RNAPolIII and TRF1, DNMT3b and HP1α were significantly decreased when single IKKγ was overexpressed. Moreover, the interplay between TERRA promoter DNA probe and DNMT3b was significantly increased, and the interplay between TERRA promoter DNA probe and RNAPolIII, TRF1 were significantly decreased when IKKα and IKKβ were co-overexpressed. Conversely, when single IKKγ was overexpressed, the interplay between TERRA promoter DNA probe and DNMT3b was significantly decreased, and the interplay between TERRA promoter DNA probe and RNAPolIII, TRF1 were significantly increased ( Figure  5B). Furthermore, the loading of DNMT3b on the TERRA promoter region was significantly increased, and the loading of RNAPolIII, TRF1 on the TERRA promoter region were significantly decreased when IKKα and IKKβ were co-overexpressed. Conversely, the loading of DNMT3b on the TERRA promoter region was significantly decreased, and the loading of RNAPolIII, TRF1 on the TERRA promoter region were significantly increased when only IKKγ was overexpressed ( Figure  5C). Thereby, the methylation of the TERRA promoter region was significantly increased when IKKα and IKKβ were co-overexpressed. Conversely, the methylation of the TERRA promoter region was significantly decreased when single IKKγ was overexpressed ( Figure 5D&5E). Furthermore, the expression of the TERRA promoter region was significantly decreased when IKKα and IKKβ were co-overexpressed. Conversely, the expression of the TERRA promoter region was significantly increased when single IKKγ was overexpressed ( Figure 5F). Moreover, compared to the control compared to the control, IKKα, IKKβ, IKKβ+IKKγ, IKKα+IKKγ, IKKα+IKKβ+IKKγ overexpression group, the interplay between TERRA probe and TERT was significantly decreased when IKKα and IKKβ were co-overexpressed,. Conversely, the interplay between TERRA probe and TERT was significantly increased when single IKKγ overexpressed ( Figure 6A). The interplay between TERRA and TERT was significantly decreased, and the interplay between TERC and TERT was significantly increased when IKKα and IKKβ were co-overexpressed. Conversely, the interplay between TERRA and TERT was significantly increased, and the interplay between TERC and TERT was significantly decreased when single IKKγ was overexpressed, ( Figure 6B). Significantly, the telomerase activity was significantly increased when IKKα and IKKβ were co-overexpressed. Conversely, the telomerase activity was significantly decreased when single IKKγ was overexpressed, (Figure 6C). Taken together, IKKα, IKKβ, and IKKγ control telomerase activity in liver cancer stem cells.

IKKα, IKKβ, and IKKγ alter telomere length
To explore whether IKKα, IKKβ, and IKKγ influence on telomere, we also detected the alteration of telomere length in these liver cancer stem cell. As shown in Figure 7A, compared to the control compared to the control, IKKα, IKKβ, IKKβ+IKKγ, IKKα+IKKγ, IKKα+IKKβ+IKKγ overexpression group, the interplay between the telomere DNA probe and POT1, Exo1, SNM1b, TRF2 were significantly increased, and the interplay between the telomere DNA probe and CST/AAF was significantly decreased when IKKα and IKKβ were co-overexpressed,. Conversely, the interplay between the telomere DNA probe and POT1, Exo1, SNM1b, TRF2 were significantly decreased, and the interplay between the telomere DNA probe and CST/AAF was significantly increased when single IKKγ was overexpressed,. Furthermore, the loading of POT1, Exo1, SNM1b, TRF2 on the telomere DNA were significantly increased when IKKα and IKKβ were co-overexpressed. Conversely, the loading of POT1, Exo1, SNM1b, TRF2 on the telomere DNA were significantly decreased when single IKKγ was overexpressed, (Figure 7B). The interplay between the telomere DNA probe and TRF2 were significantly increased when IKKα and IKKβ were co-overexpressed. Conversely, the interplay between the telomere DNA probe and TRF2 were significantly decreased when single IKKγ was overexpressed, (Figure 7C). Finally, compared to the control, IKKα, IKKβ, IKKβ+IKKγ, IKKα+IKKγ, IKKα+IKKβ+IKKγ overexpression group, the telomere length was significantly increased when IKKα and IKKβ were co-overexpressed. Conversely, the telomere length was significantly decreased when single IKKγ was overexpressed, (Figure 7D&7E). Collectively, these findings suggest that IKKα plus IKKβ increases and IKKγ decreases the telomere length in the liver cancer stem cells.

HOTAIR depletion blocks function of IKKα, IKKβ, and IKKγ
To determine whether HOTAOR regulated the IKKα, IKKβ, and IKKγ's function on telomeres, we performed rescued-test in liver cancer stem cells. As showed in the Figure 8A, compared to the control, the TERT expression was significantly increased, and the TERRA expression was significantly decreased when IKKα and IKKβ were co-overexpressed. However, the IKKα plus IKKβ function was fully abrogated when HOTAIR was knocked down. Conversely, when single IKKγ was overexpressed, the TERT expression was significantly decreased, and the TERRA expression was significantly dincreased. However, the IKKγ function was fully abrogated when HOTAIR was overexpressed. Furthermore, compared to the control group, when IKKα and IKKβ were cooverexpressed, the interaction between TERT and TERC was significantly increased, and the interaction between TRET andTERRA, DNMT3b and HP1α, RNAPolIII and TRF1 were significantly decreased. However, the IKKα plus IKKβ functions were fully abrogated when HOTAIR was knocked down. Conversely, when single IKKγ was overexpressed, the interaction between TERT and TERC was significantly decreased, and the interaction between TRET and TERRA, DNMT3b and HP1α, RNAPolIII and TRF1 were significantly increased. However, the IKKγ function was fully abrogated when HOTAIR was overexpressed ( Figure 8B). Intriguingly, HOTAIR depletion decreased the interplay between TERT and TERC, CST/AAF and telomere, increased the interaction between TERRA and TERT, TRF2 and telomere. However, compared to the IKKα, IKKβ, IKKβ+IKKγ, IKKα+IKKγ, IKKα+IKKβ+IKKγ overexpression group, the interplay between TERT and TERC, TERRA and TERT, TRF2 and telomere, CST/AAF and telomere were significantly not altered when IKKα plus IKKβ or IKKγ were overexpressed in HOTAIR knocked-down liver cancer stem cells ( Figure 8C). Ultimately, the telomerase activity and telomere length expression were significantly increased when IKKα and IKKβ were cooverexpressed. However, the IKKα plus IKKβ function was fully abrogated when HOTAIR was knocked down. Conversely, when single IKKγ was overexpressed, the telomerase activity and telomere length expression were significantly decreased. However, the IKKγ function was fully abrogated when HOTAIR was overexpressed ( Figure  8D&8E). Taken together, HOTAIR is required for IKKα plus IKKβ and IKKγ to control telomerase activity and telomere length in liver cancer stem cells.

DISCUSSION
Cancer stem cells (CSCs) are involved in tumorigenesis. The IKK is critical in cellular proliferation and in tumor-initiating [29]. However, the regulatory mechanism of IKKs has not been elucidated. To our knowledge, this is the first report demonstrating inflammatory related gene IKKα, IKKβ, and IKKγ cooperate to determine liver cancer stem cells progression by altering telomere via heterochromatin protein 1-HOTAIR axis. As shown in Figure 10, we demonstrate that IKKα plus IKKβ promoted and IKKγ inhibited liver cancer stem cell growth in vitro and in vivo. Mechanistically, IKKα plus IKKβ enhanced and IKKγ inhibited the interplay among HP1α, HP1β and HP1γ that competes for the interaction among HP1α, SUZ12, HEZ2. Therefore, IKKα plus IKKβ inhibited and IKKγ enhanced the activity of H3K27 methyltransferase SUZ12 and EZH2, which methylates H3K27 immediately sites on HOTAIR promoter region. Therefore, IKKα plus IKKβ increased and IKKγ decreased the TOTAIR expression. Strikingly, IKKα plus IKKβ decreases and IKKγ increases the HP1α interplays with DNA methyltransferase DNMT3b, which increases or decreases TERRA promoter DNA methylation. Thus IKKα plus IKKβ reduces and IKKγ increases to recruit TRF1 and RNA polymerase II deposition and elongation on the TERRA promoter locus, which increases or decreases TREEA expression.
Furthermore, IKKα plus IKKβ decreases/increases and IKKγ increases/decrease the interplay between TERT and TRRRA/between TERT and TREC. Ultimately, IKKα plus IKKβ increases and IKKγ decreases the telomerase activity. On the other hand, at the telomere locus, IKKα plus IKKβ increases/drcreases and IKKγ decreases/increases TRF2, POT1, pPOT1, Exo1, pExo1, SNM1B, pSNM1B/CST-AAF binding, which keep active telomere regulatory genes and poised for telomere length. Strikingly, HOTAIR is required for IKKα plus IKKβ and IKKγ to control telomerase activity and telomere length. HOTAIR depletion blocks the function of IKKα plus IKKβ, IKKγ on telomere. These observations suggest that HOTAIR operates the action of IKKα, IKKβ, IKKγ in liver cancer stem cells. This study provides a novel basis to elucidate the oncogenic action of IKKα, IKKβ, IKKγ and prompts that IKKα, IKKβ, IKKγ cooperate to HOTAR to be used as a new therapeutic targets for liver cancer.
It is worth mentioning that IKKα, IKKβ, IKKγ may play an important role in hepatocarcinogenesis. In this report, we focused mainly on the view how IKKα, IKKβ, IKKγ function during liver cancer stem cells malignant growth. To this date, accumulating evidence indicates that IKKα, IKKβ, IKKγ influence on cell proliferation. For examples, inhibition of IKK/NF-κB pathway controls stem cell proliferation [30]. IKKβ plays a role during intestinal tumorigenesis [31]. BRAF-induced tumorigenesis is dependent on IKKα [32]. IKKβ could regulates VEGF expression in ovarian cancer as an antiangiogenic target [33]. Our present findings are consistent with some reports.
It is worth noting that our findings in this study provide novel evidence for an active role of IKKα plus Figure 10: The schematic diagram illustrates a model that that IKKα plus IKKβ promoted and IKKγ inhibited liver cancer stem cell growth in vitro and in vivo. (Continued ) IKKβ promotion or IKKγ inhibition of liver cancer stem cell growth. Herein, the involvement of promotion or inhibition of liver cancer stem cells growth based on IKKα, IKKβ, IKKγ is supported by results from two parallel sets of experiments: (1) IKKα plus IKKβ promoted and IKKγ inhibited liver cancer stem cell growth in vitro; (2) IKKα plus IKKβ promoted and IKKγ inhibited liver cancer stem cell growth in vivo.
Strikingly, our observations suggest that IKKα plus IKKβ increased and IKKγ inhibited HOTAIR expression dependent on tri-methylation of Histone H3 on the twenty-seven lysine. This assertion is based on several observations in IKKα plus IKKβ or IKKγ overexpressed liver cancer stem cells: (1) IKKα plus IKKβ enhanced and IKKγ inhibited the interplay among HP1α, HP1β and HP1γ that competes for the interaction among HP1α, SUZ12, HEZ2. HP1α, HP1β and HP1γ that competes for the interaction among HP1α, SUZ12, HEZ2. Therefore, IKKα plus IKKβ inhibited and IKKγ enhanced the activity of H3K27 methyltransferase SUZ12 and EZH2, which methylates H3K27 immediately sites on HOTAIR promoter region. Therefore, IKKα plus IKKβ increased and IKKγ decreased the TOTAIR expression. Strikingly, IKKα plus IKKβ decreases and IKKγ increases the HP1α interplays with DNA methyltransferase DNMT3b, which increases or decreases TERRA promoter DNA methylation. So IKKα plus IKKβ reduces and IKKγ increases to recruit TRF1 and RNA polymerase II deposition and elongation on the TERRA promoter locus, which increases or decreases TREEA expression. Further on, IKKα plus IKKβ decreases/increases and IKKγ increases/ decrease the interplay between TERT and TRRRA/between TERT and TREC. Ultimately, IKKα plus IKKβ increases and IKKγ decreases the telomerase activity. On the other hand, at the telomere locus, IKKα plus IKKβ increases/drcreases and IKKγ decreases/increases TRF2, POT1, pPOT1, Exo1, pExo1, SNM1B, pSNM1B/CST-AAF binding, which keep active telomere regulatory genes and poised for telomere length. B. HOTAIR is required for IKKα plus IKKβ and IKKγ to control telomerase activity and telomere length. HOTAIR depletion blocks the function of IKKα plusIKKβ, IKKγ on telomere. κB through H3K27me3. (4) IKKα plus IKKβ increased and IKKγ decreased HOTAIR expression dependent on H3K27me3. Researches indicated heterochromatin causes epigenetic repression that control gene expression and function [34]. HP1α is an essential protein critical for heterochromatin assembly and regulation [35]. Strikingly, HP1 promotes tumor suppressor BRCA1 functions during the DNA damage response [36]. The trimethylation of histone H3 lysine 27 (H3K27me3) contributes to gene repression [37]. NF-κB is involved in inflammation and tumor growth [38].
On the other hand, we find that IKKα, IKKβ, and IKKγ control telomerase activity and telomere length. The stability of telomeres depends upon TRF2, which prevents inappropriate repair [39]. Upon telomere shortening or telomere uncapping induced by loss of TRF2, telomeres elicit a DNA-damage response [40]. Our previous study shows CUDR promotes liver cancer stem cell growth through upregulating TERT [41]. In addition, telomeres are protected from hyper-resection through the repression of the ATM and ATR kinases by TRF2 and TPP1-bound POT1a/b, respectively [42]. Shelterin can protect chromosome ends as a TRF2-tethered TIN2/TPP1/POT1 complex [43]. This assertion is based on several observations in TLR4 overexpression or knockdown liver cancer stem cells: (1) When IKKα and IKKβ were co-overexpressed, the telomerase activity and telomere length were significantly increased. (2) Conversely, when only IKKγ was overexpressed, the telomerase activity and telomere length were significantly decreased. (3) IKKα plus IKKβ decreased and IKKγ increased the HP1α interplay with DNA methyltransferase DNMT3b, which increased or decreased TERRA promoter DNA methylation. (2) IKKα plus IKKβ reduced and IKKγ increased to recruit TRF1 and RNA polymerase II deposition and elongation on the TERRA promoter locus, which increased or decreased TREEA expression. (3) IKKα plus IKKβ decreases/ increases and IKKγ increases/decrease the interplay between TERT and TRRRA/between TERT and TREC. (4) IKKα plus IKKβ increases/drcreases and IKKγ decreases/increases TRF2, POT1, pPOT1, Exo1, pExo1, SNM1B, pSNM1B/CST-AAF binding, which keep active telomere regulatory genes and poised for telomere length.
Significantly, HOTAIR is required for IKKα plus IKKβ and IKKγ to control telomerase activity, telomere length and tumorigenesis. HOTAIR is associated with a variety of human cancers, such as breast, liver and endometrial carcinoma [44]. HOTAIR induced colony formation and orthotopic tumor growth [45]. HOTAIR directly decreased WIF-1 expression by promoting its histone H3K27 methylation and then activated the Wnt/βcatenin signaling pathway [46]. Furthermore, HOTAIR promotes proliferation of ovarian cancer cells through regulating PIK3R3 [47]. Herein, the involvement of promotion or inhibition of telomere length or telomerase activity based on IKKα, IKKβ, IKKγ is supported by results from three parallel sets of experiments: (1) HOTAIR depletion abolished the function of IKKα plus IKKβ, IKKγ on telomere length. (2) HOTAIR depletion abolished the function of IKKα plus IKKβ, IKKγ on telomerase. (2) HOTAIR depletion abrogated IKKα plus IKKβ oncogenic function.
Especially, our results showed IKKα plus IKKβ enhanced the interplay among HP1α, HP1β and HP1γ by an unknown mechanism, and then inhibited the activity of H3K27 methyltransferase SUZ12 and EZH2 and increased the HOTAIR expression. However, single IKKα or IKKβ overexpression could not alter HOTAIR expression. Because of excessive HOTAIR, IKKα plus IKKβ increases the interplay between TERT and TREC, the telomerase activity and telomere length. However, single IKKα or IKKβ overexpression could not alter the telomerase activity, increases telomere length. Moreover, because of the increase of telomerase activity and telomere length, IKKα plus IKKβ simultaneous overexpression significantly increased liver cancer stem cells's growth, while single IKKα or IKKβ overexpression could not promote liver cancer stem cells's growth.
In addition, our results showed that HOTAIR depletion blocks IKKα plus IKKβ, IKKγ function on telomere. It suggest HOTAIR is required for IKKα plus IKKβ and IKKγ to control telomerase activity and telomere length. In this study, we found IKKα plus IKKβ enhanced and HP1γ decreased the interplay among HP1α, HP1β and HP1γ by an unknown mechanism, and then IKKα plus IKKβ inhibited and HP1γ enhanced the activity of H3K27 methyltransferase SUZ12 and EZH2. Ultimately, IKKα plus IKKβ increased and HP1γ decreased the HOTAIR expression. Because of HOTAIR, IKKα plus IKKβ increases and HP1γ decreased the interplay between TERT and TREC, the telomerase activity and telomere length. IKK is also a crucial protein kinase that activates NF-κB translocating from cytoplasm to nucleus for DNA binding. Especially, IKK phosphorylates NF-κB to modulate its promoter specificity and promote genes expression. NF-κB could regulates and controls telomere related genes. Although IKK could regulate NF-κB and HOTAIR, both abnormal NF-κB and HOTAIR could alter telomeres through different pathways.
We should explore the function of IKKα plus IKKβ, IKKγ in liver cancer stem cells. For example, what causes strong oncogenic action of IKKα plus IKKβ and triggers IKKγ's actions? How does IKKα plus IKKβ, IKKγ cooperates with HOTAIR? Does IKKα plus IKKβ, IKKγ regulates a series of molecular events liver stem cells malignant growth? Answering these questions will help understand the mechanism about liver stem cell malignant differentiation. In summary, our present data indicated that IKKα plus IKKβ, and IKKγ promotes/inhibits liver cancer stem cells malignant progression through altering telomere length and telomerase activity dependent on HOTAIR, with clinic implications. These observations provide insight into a novel link between IKKα, IKKβ, IKKγ and liver cancer stem malignant transformation.

Cell transfection and stable cell lines cell transfection and stable cell lines
Cells were transfected with DNA plasmids using transfast transfection reagent lipofectamine R 2000 (Invitrogen) according to manufacturer's instructions.

RT-PCR
Total RNA was purified using Trizol (Invitrogen) according to manufacturer's instructions. cDNA was prepared by using oligonucleotide (dT) [17][18] , random primers, and a SuperScript First-Strand Synthesis System (Invitrogen). PCR analysis was performed under the specical conditions. β-actin was used as an internal control.

Western blotting
The logarithmically growing cells were washed twice with ice-cold phosphate-buffered saline (PBS, Hyclone) and lysed in a RIPA lysis buffer. Cells lysates were centrifuged at 12,000g for 20 minutes at 4°C after sonication on ice, and the supernatant were separated. After being boiled for 5-10 minutes in the presence of 2-mercaptoethanol, samples containing cells proteins were separated on a 10% sodium dodecyl sulfatepolyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto a nitrocellulose membranes (Invitrogen, Carlsbad, CA, USA). Then blocked in 10% dry milk-TBST (20mM Tris-HCl [PH 7.6], 127mM NaCl, 0.1% Tween 20) for 1 h at 37°C. Following three washes in Tris-HCl pH 7.5 with 0.1% Tween 20, the blots were incubated with 0.2 μg/ml of antibody (appropriate dilution) overnight at 4°C. Following three washes, membranes were then incubated with secondary antibody for 60 min at 37°C or 4°C overnight in TBST. Signals were visualized by ODYSSEY infrared imaging system (LI-COR, Lincoln, Nebraska USA). IRDye 680LT secondary antibodies were purchased from LI-COR scientific company.

Co-immunoprecipitation (IP)
Cells were lysed in 1 ml of the whole-cell extract buffer A (50mM pH7.6 Tris-HCl, 150mMNaCl, 1%NP40, 0.1mMEDTA,1.0mM DTT,0.2mMPMSF, 0.1mM Pepstatine,0.1mM Leupeptine,0.1mM Aproine). Five-hundred-microliter cell lysates was used in immunoprecipitation with antibody. In brief, protein was pre-cleared with 30μl protein G/A-plus agarose beads (Santa Cruz, Biotechnology, Inc. CA) for 1 hour at 4°C and the supernatant was obtained after centrifugation (5,000rpm) at 4°C. Precleared homogenates (supernatant) were incubated with 2 μg of antibody and/ or normal mouse/rabbit IgG by rotation for 4 hours at 4°C, and then the immunoprecipitates were incubated with 30μl protein G/A-plus agarose beads by rotation overnight at 4°C, and then centrifuged at 5000rpm for 5 min at 4°C. The precipitates were washed five times×10min with beads wash solution (50 mM pH7.6 TrisCl,150mMNaCl,0.1%NP-40,1mM EDTA) and then resuspended in 60μl 2×SDS-PAGE sample loading buffer to incubate for 10 min at 100°C. Then Western blot was performed with a another related antibody indicated in Western blotting.

DNA pull down
Cells were lysed by sonication in HKMG buffer (10 mM HEPES, PH7.9, 100 mM KCl, 5 mM MgCl 2 , 100% glycerol, 1 mM DTT, and 0.5% NP40) containing protease and phosphatase inhibitors for the preparation of nuclear exact. Equal amount of cell nuclear extracts were precleared with Streptavidin-agarose Resin (Thermo) for 1 hours, and then were incubated with 1μg biotinylated double-stranded-oligonucleotides and together with 10μg poly(dI-dC) at 4°C for 24 hours. DNA-bound proteins were collected with the incubation with streptavidinagarose Resin for 1 hour with gently shaking to prevent precipitation in solution. Following 5 washings of the resin bound complex with 0.5-1.0 ml of binding buffer, the samples were boiled and subjected to SDS-PAGE and Western blot analysis.

Chromatin immunoprecipitation (CHIP) assay
Cells were cross-linked with 1% (v/v) formaldehyde (Sigma) for 10 min at room temperature and stopped with 125 mm glycine for 5 min. Crossed-linked cells were washed with phosphate-buffered saline, resuspended in lysis buffer, and sonicated for 8-10 min in a SONICS VibraCell to generate DNA fragments with an average size of 500 bp or so. Chromatin extracts were diluted 5-fold with dilution buffer, pre-cleared with Protein-A/G-Sepharose beads, and immunoprecipitated with specific antibody on Protein-A/G-Sepharose beads. After washing, elution and de-cross-linking, the ChIP DNA was detected by PCR.

Super-EMSA (Gel-shift)
Cells were washed and scraped in ice-cold PBS to prepare nuclei for electrophoretic gel mobility shift assay with the use of the gel shift assay system modified according to the manufacturer's instructions (Promega).

Quantitative telomerase detection
The telomerase activity was measured by using Quantitative Telomerase Detection Kit (MT3010) according to manufacturer's instructions (US Biomax, Inc).

Telomere length assay
Telomere length assay using Telo TAGGG PCR ELISApuls kit according to manufacturer's instructions (Roche). A standard curve is established by dilution of known quantities of a synthesised 84 mer oligonucleotide containing only TTAGGG repeats.

Cells proliferation CCK8 assay
Cells were synchronized in G0 phase by serum deprivation and then released from growth arrest by reexposure to serum, and then cells were grown in complete medium for assay. The cell proliferation reagent CCK8 is purchased from Roch and the operation according to the manufacturer instruction.
Soft agar colony formation capacity assay 2 × 10 2 cells were plated on a 6 well plate containing 0.5% (lower) and 0.35% (upper) double layer soft-agar. The 6 well plates were incubated at 37°C in humidified incubator for 14 days. The cells were fed 1-2 times per week with cell culture media (DMEM). Soft-agar colonies on the 6 well plates were stained with 0.5 ml of 0.05% Crystal Violet for more than 1 hour and the colonies were counted.

BrdU staining
70-80% confluent cells were cultured for 24 hour before treatment with 10μl BrdU (Roche) for 4 hours. Immunofluorescent staining with an anti-BrdU antibody was performed according to the manufacturer's instructions (Becton Dickinson). BrdU positive cells from ten random chosen fields of at least three independent samples were counted.

Xenograft transplantation in vivo
The Four-weeks athymic Balb/C mouse was injected the liver cancer cells at the armpit area subcutaneously. The mice were observed over 4 weeks, and then sacrificed to recover the tumors. The wet weight of each tumor was determined for each mouse. A portion of each tumor was fixed in 4% paraformaldehyde and embedded in paraffin for histological hematoxylineosin (HE) staining. The use of mice for this work was reviewed and approved by the institutional animal care and use committee in accordance with China national institutes of health guidelines. www.impactjournals.com/oncotarget

Statistical analysis
The significant differences between mean values obtained from at least three independent experiments. Each value was presented as mean±standard error of the mean (SEM) unless otherwise noted, with a minimum of three replicates. Student's t-test was used for comparisons, with P<0.05 considered significant.