TLE2 is associated with favorable prognosis and regulates cell growth and gemcitabine sensitivity in pancreatic cancer.

BACKGROUND: The transducin-like enhancer of split (TLE) proteins are a group of transcriptional corepressors. They play a crucial role in cellular homeostasis and are involved in various cancers. Compared with other TLE family members, little is known about the role and the underlying mechanism of TLE2 in human cancers. This study aimed to investigate the role of TLE2 in pancreatic ductal adenocarcinoma (PDAC) using in silico analysis and in vitro experiments.
METHODS: Data were obtained from the Cancer Genome Atlas (TCGA) database to evaluate the prognostic value of TLE2 in PDAC. The MiaPaCa-2 cell line was transfected with siRNA to inhibit endogenous TLE2 expression, and a PANC-1 cell line with stable TLE2 overexpression was constructed using lentiviral transfection, which were confirmed by real-time polymerase chain reaction and western blotting. MTT assay, transwell invasion assays, and flow cytometry were carried out to assess cell viability, invasion, and apoptosis, respectively. TLE2 expression in PDAC cells was altered to evaluate their sensitivity to gemcitabine. Gene set enrichment analysis (GSEA) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted to predict the biological role of TLE2.
RESULTS: High expression of TLEs was significantly associated with increased overall survival (OS) and disease-free survival (DFS) in patients with PDAC. Among the PDAC cell lines, TLE2 expression was lowest and highest in PANC-1 cells and MiaPaCa-2 cells, respectively. TLE2 overexpression impaired the proliferation ability of PANC-1 cells and downregulation of TLE2 promoted the proliferation of MiaPaCa-2 cells. Upregulation of TLE2 in PANC-1 cells induced S-phase accumulation and sensitivity to gemcitabine. In contrast, the downregulation of TLE2 in MiaPaCa-2 cells promoted resistance to gemcitabine. Moreover, bioinformatics analysis also revealed the potential tumor suppressor role of TLE2 and uncovered a close relationship between TLE2 expression and cell cycle regulation.
CONCLUSIONS: Our results suggest that TLE2 expression is correlated with prognosis in patients with PDAC and show that TLE2 plays a central role in the regulation of cell proliferation, the cell cycle, and gemcitabine sensitivity. This study provides new insights and evidence that TLE2 functions as a tumor suppressor gene and prognostic marker in PDAC. 2020 Annals of Translational Medicine. All rights reserved.

Introduction

Pancreatic ductal adenocarcinoma (PDAC) is the most common pancreatic malignancy and one of the most challenging cancers to treat (1). Despite the substantial therapeutic advancements, the prognosis for patients with PDAC is still dismal (2-4). Therefore, a better understanding of the molecular biology and tumor genetics of PDAC, and new therapeutic targets that can improve patient outcomes are urgently needed.

Transcriptional corepressors are a group of proteins that play vital roles in the transcriptional regulation of various genes (5,6). The Groucho (Gro)/transducin-like enhancer of split (TLE) proteins are a large family of gene regulators that are widely expressed in eukaryotes and influence the transcriptional activity of numerous genes (7). Instead of directly binding to DNA strands, TLE proteins bind to other transcription factors (8,9), usually resulting in transcriptional repression (10,11). TLE proteins have been observed to interact with transcriptional factors, including Hes, c-Myc, and Runx (12,13). Additionally, some studies have reported that TLE1 and TLE3 can participate in tumor development as tumor suppressor genes (14-17).

TLE2 is a member of the TLE gene family located on chromosome 19p13.3. Previous studies have revealed that TLE2 participates in the development of pancreatic tissue through its interactions with a series of transcription factors (18). Furthermore, TLE2 is involved in the ventral telencephalon from the early neural plate stage during neurogenesis (19). Many genes known to be involved in the regulation of histogenesis and tissue development, such as Robo and Wnt, may also perform pivotal roles in neoplastic conditions (20-22). Currently, the role TLE2 plays in most human cancers, especially pancreatic cancer, is a long way from being elucidated. Therefore, this study aimed to examine the prognostic value and cellular function of TLE2 in PDAC through in silico and in vitro study for the first time, in order to explore the role of TLE2 in PDAC and to elucidate the underlying mechanisms. We present the following article in accordance with the MDAR reporting checklist (available at http://dx.doi.org/10.21037/atm-20-5492).

Methods

Cells and reagents

PANC-1 and MiaPaCa-2 human pancreatic cancer cell lines were obtained from American Type Culture Collection (ATCC) and maintained in Dulbecco’s Modified Eagle Medium (DMEM; Invitrogen) with 10% fetal bovine serum (HyClone) at 37 °C with 5% CO2.

Transfection of siRNA

MiaPaCa-2 cells were transfected with siRNAs targeting the sequence 5'-CTCCAGAAATAACAAATACATCG-3' of TLE2. A scramble (5'-GGACGGAGCAGTACAA-3') without known targets was chosen as a non-sense control. Transfection was performed with Lipofectamine 2000 (Thermo Fisher Scientific, USA) in line with the manufacturer’s protocol. siRNA and control oligonucleotides were synthesized by Jennio Biological Technology (Guangzhou, China).

Lentiviral transfection for TLE2 overexpression

The TLE2 overexpression lentiviral vector was established by recombining the pGC-LV-GV287-GFP vector with the TLE2 gene. PANC-1 cells overexpressing TLE2 were constructed via lentiviral-mediated plasmid transfection (MOI =10). An empty vector served as a transfection control. Untransfected cells were used as a normal control. TLE2 expression was validated by quantitative reverse transcription polymerase chain reaction (qRT-PCT) and western blot analysis.

RNA extraction and qRT-PCR

Total RNA was extracted from cells using TRIZOL reagent (Invitrogen), and cDNA was reverse transcribed with the cDNA Reverse Transcription Kit (Takara Biotechnology). Then, qRT-PCR was performed on each cDNA template in 20-µL reaction systems using the SYBR Green Kit (Takara, Japan) on an ABI 7500HT Fast real-time PCR system. The following primers were used for the amplification of TLE2 fragments: forward 5'-TATTTCCGATTACGGCACTCG-3'; reverse 5'-CTACTCTCCATTCCGACCGC-3', GAPDH was used as internal control (forward 5'-TTGTCTCCTGCGACTTCAACAG-3'; reverse 5'-GGTCTGGGATGGAAATTGTGAG-3').

Western blot analysis

Pancreatic cancer cell lines were washed and resuspended with PBS twice. Total protein was isolated and quantified using RIPA extraction buffer and the BCA kit (Pierce). The protein samples were separated by 10% SDS-PAGE and transferred onto polyvinylidene fluoride (PVDF) membranes. The membranes were then blocked with buffer containing 5% non-fat milk and incubated with TLE2 antibody (rabbit polyclonal anti-human, Abcam, USA; 1:200, No. ab206147). Anti-rabbit IgG was used as the secondary antibody (Santa Cruz Biotechnology, No. sc-2357). The blots were visualized with an ELC kit (Millipore, USA).

Cell viability assays

For the evaluation of cell proliferation, cell viability assays were carried out. Briefly, different groups of PANC-1 cells and MiaPaCa-2 cells (control, scramble transfection, empty vector transfection, siTLE2, and TLE2 overexpression) were seeded in 96-well plates at a density of 1,000 cells/well and incubated for 72 hours at 37 °C with 5% CO2. For gemcitabine cytotoxicity experiments, groups of PANC-1 and MiaPaCa-2 cells, as described above, were seeded in 96-well plates at a density of 1,000 cells/well. The cells were treated with different concentrations of gemcitabine for 48 hours. For both assays, the MTT absorbance was finally evaluated at 570 nm to determine cell viability.

Invasion assays

PANC-1 and MiaPaCa-2 cells with the indicated treatments were cultured in 6-well plates and re-suspended in serum-free medium. Cell invasion was studied by performing a Transwell assay with the chamber precoated with Matrigel (8.0 µm pore size, Corning) and a cell density of 2.0×104 cells/well. The numbers of migrated cells on the lower side of the Transwell membrane were counted and compared.

Cell cycle analysis

Flow cytometry was carried out to detect the distribution of cell cycle. Briefly, cells were transfected with either siRNA targeting TLE2 or a plasmid overexpressing TLE2. Single cell suspensions were prepared 48 hours after transfection. Cell aliquots were fixed in 70% ethanol and stained with propidium iodide for flow cytometry analysis using a FACScan instrument (BD Biosciences, USA).

Survival analysis

The prognostic values of TLE1, TLE2, and TLE3 in cancers were analyzed according to the transcriptome data and clinical parameters obtained from The Cancer Genome Atlas (TCGA; http://cancergenome.nih.gov/), which were analyzed with transcripts per million (TPM) as log2 (TPM +1) in Gene Expression Profiling Interactive Analysis (GEPIA; http://gepia.cancer-pku.cn) portal (23).

Functional annotation analysis

The RNA-seq data (FPKM, fraction per kilobase per million reads) of PDAC were downloaded from the TCGA portal (TCGA-PAAD). The correlation between the expressions of TLE2 and other genes at the transcriptional level (FPKM) was estimated by Pearson’s correlation analysis. Genes exhibiting a strong positive or negative correlation with TLE2 expression were involved in the functional annotation analysis. Functional annotation analysis was performed with the clusterProfiler R package (version 3.7.0), gene set enrichment analysis (GSEA), and WEB-based GEne SeT AnaLysis Toolkit (WebGestalt; (http://bioinfo.vanderbilt.edu/webgestalt) (24).

Statistical analysis

In vitro experiments were performed in triplicate. Continuous variables were presented as mean ± SD (standard deviation). Student’s t-test was used for comparisons between two groups, and comparisons between more than three groups were conducted using one-way analysis of variance (ANOVA) with a post hoc least significant difference (LSD) test. Survival outcomes, including overall survival (OS) and disease-free survival (DFS), were compared by Kaplan–Meier curve analysis and the log-rank test. All statistical analyses were performed with SPSS 19.0 (SPSS Inc., Chicago, IL, USA) and R (version 3.7.0). A two-sided P value of <0.05 was considered as statistically significant.

Results

TLE2 was correlated with prognosis in PDAC patients

Survival analysis was performed to investigate the correlations between the expressions of TLE family members (TLE1, TLE2, TLE3, TLE4, and TLE5) and prognosis of human cancers. Our results revealed a certain degree of heterogeneity in the prognostic value of TLEs across different cancer types (). However, in PDAC patients, only increased TLE2 expression was significantly correlated with better overall survival (OS) (P<0.001) () and disease-free survival (DFS) (P=0.022) (). In contrast, no correlation was found between the expression of TLE1, TLE3, and TLE4, and the prognosis of patients with PDAC (). Increased expression of TLE5 was correlated with improved OS (P=0.019) (), but not with improved DFS (P=0.59) in PDAC patients (). Further, we investigated the relation between TLE2 expression and clinicopathological parameters, including TNM stage and differentiation. As shown in , TLE2 expression was significantly associated with tumor stage (P<0.001) and tumor differentiation (P<0.001).

Figure 1

Prognostic value of members of TLE family across cancers. (A) The relationship between TLEs and OS across multiple cancer types. (B) The relationship between TLEs and DFS across multiple cancer types. (C,D,E,F,G,H,I,J,K,L) Patients in the TCGA-PAAD cohort were stratified by the median values of TLE expressions, the prognostic value of each TLE family member was evaluated by Kaplan-Meier survival curves with the log-rank test. The figure was generated with GEPIA, as mentioned the methods above. (M) XXX. TLE, transducin-like enhancer of split; OS, overall survival; DFS, disease-free survival; TCGA-PAAD, the Cancer Genome Atlas portal; GEPIA, Gene Expression Profiling Interactive Analysis.

TLE2 expression affects the proliferation and cell cycle of PDAC cells

To further investigate the biological function of TLE2 in PDAC, the expression of TLE2 was investigated across various pancreatic cancer cell lines by qRT-PCR and western blotting. The MiaPaCa-2 cell line was found to have the highest expression of TLE2 among the six pancreatic cancer cell lines, while the PANC-1 cell line exhibited the lowest TLE2 expression level (). Consequently, TLE2 overexpression plasmid was constructed and introduced into PANC-1 cells. Also, siRNA oligonucleotides targeting TLE2 were transduced into MiaPaCa-2 cells to inhibit TLE2 expression, as validated by western blotting ().

Figure 2

Effect of TLE2 on the proliferation, invasion, and cell cycle of pancreatic cancer cells in vitro. (A,B) The expression of TLE2 in pancreatic cell lines was assessed by qRT-PCR (A) and western blotting (B), with GAPDH used as an internal control. (C) TLE2 protein expression was validated in pancreatic cancer cell lines with the indicated treatments by western blotting. (D,E) Cell growth in PANC-1 cells (D) and MiaPaCa-2 cells (E) was assessed by MTT assays at an absorption value of 570 nm (OD570). (F,G) Flow cytometry was used to detect cell cycle distribution in the PANC-1 cells (F) and MiaPaCa-2 cells (G). The proportions of cells in each phase are displayed in the histograms. (H) Cell invasion was evaluated using Transwell matrigel invasion assays. As showed in the histograms (I,J), the numbers of cells per HPF that passed through the Transwell chambers were counted, normalized, and compared. TLE, transducin-like enhancer of split; HPF, high power field.

The effect of TLE2 expression on cell proliferation, invasion, and cell cycle were evaluated by MTT assay, Transwell assay, and flow cytometry, respectively. Our data revealed that the exogenous overexpression of TLE2 significantly inhibited the proliferation ability of PANC-1 cells (). In contrast, downregulation of TLE2 significantly promoted the cell proliferation of MiaPaCa-2 cells (). Interestingly, TLE2 overexpression induced delayed M-phase entry and S-phase accumulation (). However, no obvious effect on cell invasion was observed upon alternation of TLE2 expression ().

Gene function and pathway enrichment analysis

To elucidate the potential functional role of TLE2, functional annotation and enrichment analysis was carried out. The correlation between TLE2 and other genes at the transcriptional level was determined by Pearson’s correlation analysis. Next, we selected the 200 genes most positively correlated with TLE2 expression, as well as the 200 genes most negatively correlated with TLE2 expression, for further functional analysis (https://cdn.amegroups.cn/static/application/15d2d5448ba8a81d254cae02e0738b3e/atm-20-5492.pdf). By using clusterProfiler package, 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified for the top 200 genes negatively correlated with TLE2 expression, including “Cell cycle”, “Oocyte meiosis”, “p53 signaling pathway”, “Proteasome”, “On carbon pool by folate”, and “Cellular senescence”, most of which are important pathways/mechanisms involved in the regulation of cell proliferation ().

Figure 3

Functional annotation and enrichment of key pathways associated with TLE2 expression in pancreatic cancer. (A) KEGG enrichment analysis of the genes most correlated with TLE2 was performed using the R package clusterProfiler. (B) KEGG pathway database enrichment analysis of the genes most correlated with TLE2 was performed using the WebGelstat tool. (C) The GSEA enrichment analysis results of the top KEGG pathways and oncogenic signatures in the TCGA pancreatic cancer samples according to TLE2 expression. The descriptions of the abbreviations of pathways are as follows: HINATA_NFKB_MATRIX, Matrix, adhesion or cytoskeleton genes induced by NF-kappaB in primary keratinocytes and fibroblasts; KEGG_Cell_Cycle, Cell cycle; KEGG_DNA_Replication, DNA replication; KEGG_Proteasome, Proteasome; MTOR_Up.V1_UP, Genes up-regulated by everolimus in prostate tissue; P53_DN.V1_UP, Genes up-regulated in NCI-60 panel of cell lines with mutated TP53; SINGH_KRAS_DEPENDENCY_SIGNATURE, genes defining the KRAS dependency signature. KEGG, Kyoto Encyclopedia of Genes and Genomes; GSEA, gene set enrichment analysis; TCGA, the Cancer Genome Atlas.

Furthermore, gene enrichment analysis using WEB-based GEne SeT AnaLysis Toolkit also demonstrated that the “Cell cycle” pathway was most significantly enriched in genes negatively correlated with TLE2, which indicated a close relationship between TLE2 and cell cycle (), and was consistent with our above in vitro finding. Moreover, GSEA analysis was performed using all of the gene sets in “oncogenic signatures” and “KEGG subset of Canonical pathways” obtained from the GSEA website (https://www.gsea-msigdb.org/gsea/msigdb). PDAC samples from the TCGA dataset were sorted according to expression of TLE2. As shown in , the enrichment results showed that several oncogenic pathways/mechanisms were enriched in samples with low TLE2 expression. The enrichment of “KEGG_Cell_Cycle” and “KEGG_DNA_Replication” pathways support the findings of our in vitro experiments and functional annotation analysis that TLE2 has an important role in the regulation of cell proliferation and cell cycle in pancreatic cancer cells. Moreover, the enrichment of Kras-related signaling pathways, NF-kb-related mechanisms, and mTOR-related gene sets in samples with low TLE2 expression further indicated the tumor suppressor role of TLE2.

Upregulation of TLE2 enhances the sensitivity to Gemcitabine in PDAC cells

The finding that the S-phase cell ratio was increased by TLE2 overexpression, as well as observations that the cell cycle pathway was enriched in PDAC samples with low TLE2 expression, impelled us to investigate the impact of TLE2 dysregulation on sensitivity to gemcitabine in pancreatic cancer cells. PANC-1 cells and MiaPaCa-2 cells were treated with different concentrations of gemcitabine. Through cell viability assays, we found that the cell viability of gemcitabine-treated PANC-1 cells was significantly inhibited after TLE2 overexpression (), suggesting that increased TLE2 expression was associated with gemcitabine sensitivity. In contrast, the proliferation capacity of MiaPaCa-2 cells treated with gemcitabine was increased after TLE2 silencing (). Overall, compared with the control group, TLE2 overexpression in PANC-1 cells significantly increased their sensitivity to gemcitabine, while knockdown of TLE2 substantially decreased the sensitivity.

Figure 4

TLE2 overexpression sensitizes pancreatic cells to gemcitabine. (A) Under treatment with gemcitabine at the indicated concentrations, cell viability was assessed by MTT assay, which showed that the proliferation capacity of PANC-1 cells was significantly decreased in TLE2-overexpressing cells in comparison with the control cells. (B) MiaPaCa-2 cells were treated with gemcitabine at the indicated concentrations, and the MTT assay showed that the proliferation capacity was significantly increased in TLE2-silenced cells compared with the control cells. *, P<0.05; ***, P<0.001.

Discussion

In the present study, the role of TLE2 in PAAD was investigated. In silico data revealed that decreased TLE2 expression was associated with a poorer survival outcome in PDAC patients. In vitro experiments showed that TLE2 overexpression impaired the proliferation ability of pancreatic cells, and downregulation of TLE2 improved the proliferation capacity. Moreover, upregulation of TLE2 induced S-phase accumulation and sensitized pancreatic cells to gemcitabine. Functional annotation revealed TLE2 to be closely correlated with several important mechanisms involved in tumor development, particularly cell cycle regulation, which also explained the in vitro findings. Thus, we concluded that TLE2 acts as a tumor suppressor gene in PDAC.

Members of the TLE protein family have structural redundancy but act as repressors for a broad range of gene targets (7,25). Experimental studies have demonstrated that TLE proteins are critical to cell fate determination and influence the differentiation and development of tissue dramatically (26). In neurogenesis, TLE1 has been recognized as a key player participating in multiple mechanisms that repress neuronal differentiation (27-29). Mutation of the TLE1 gene has been discovered in patients with neurodevelopmental disorders, including postnatal microcephaly (27). During bone development, a profound decrease in bone mineralization was displayed in TLE4-deficient, revealing the essential role of TLE4 in the intrinsic regulation of bone development (30). In embryonic development, TLE1 has been shown to regulate the differentiation of chicken embryo fibroblasts through directly binding to the transcription factor Qin (31).

The significant involvement of TLEs in development indicates that they play a role in the regulation of cellular homeostasis and cell fate. Therefore, despite the effect of TLEs on normal cells, some studies have focused on their roles in human cancers. TLE4 acts as a tumor suppressor in leukemia through regulating Wnt-mediated inflammation (32). A large amount of evidence suggests that TLE1 plays a central role in the progression and pathophysiology of many cancers (33). Meanwhile, in prostate cancer, breast cancer, and colorectal cancer, TLE3 acts as a tumor suppressor (17,34,35). Experiments have also revealed that TLE4 acts as a tumor suppressor in myeloid leukemia and colorectal cancer through regulation of the Wnt and JNK/c-Jun pathways, respectively (32,36). Compared with the other TLEs, only recently was an association between TLE2 and cancer revealed. He et al. found that TLE2 could inhibit replication and transcription in Kaposi’s sarcoma-associated herpesvirus (37). Wu et al. identified TLE2 as an independent prognostic factor in patients with muscle invasive bladder cancer (38). As far as we know, there is currently very limited evidence identifying the role of TLE2 in other cancers, including PDAC.

The TLEs have been shown to regulate multiple steps during the development of the pancreas including differentiation of progenitor cells and endocrine cell maturation (18). With respect to pancreatic cancer, we reported for the first time the tumor suppressor role of TLE2 in PDAC. Both in vitro and in silico results demonstrated a close relationship between TLE2 and cell proliferation and cell cycle, which mimics the findings that other TLEs showed a strong growth regulation function in other cell systems. Additionally, multidrug resistance of PDAC cells hinders the effectiveness of chemotherapy, and the relationship between TLE2 and cell cycle inspired us to explore the effect of TLE2 on gemcitabine sensitivity in pancreatic cancer cells. Here, we provided in vitro evidence that upregulating TLEs could enhance the sensitivity of pancreatic cancer cell lines to gemcitabine, which could partly be explained by the specificity of gemcitabine for S-phase cells.

Although our results revealed the important role of TLE2 as a tumor suppressor in PDAC, there were some deficiencies in this study. First, the experiments were in vitro, therefore further in vivo experiments were needed to compensate for the in vitro evidence. Second, we focused on the toughest malignancy of pancreas, the PDAC; considering that the TLEs were widely involved in the regulation of pancreas development and differentiation, the role of TLEs other pancreatic tumors deserves further investigation. Last, our bioinformatic analysis revealed several key pathways associated with TLE2 expression in PDAC, further confirmation of these mechanisms via more in-depth researches is still required.

In summary, high expression of TLE2 is associated with improved survival in patients with PDAC. We found a negative correlation between TLE2 expression and proliferation capacity of PDAC cells, probably through cell cycle regulation. Moreover, our present study revealed that TLE2 expression was positively correlated with gemcitabine sensitivity in PDAC cells. Our study revealed the tumor suppressor role of TLE2 in PDAC for the first time. Although our findings have significant implications for therapeutic development in PDAC, the mechanism that mediates the tumor suppressor effect of TLE2 in PDAC demands further investigation.

The article’s supplementary files as