Skip to main content
  • Research article
  • Open access
  • Published:

TUSC2P suppresses the tumor function of esophageal squamous cell carcinoma by regulating TUSC2 expression and correlates with disease prognosis



Pseudogenes are RNA transcripts with high homology with its parent protein-coding genes. Although pseudogenes lost the ability to produce protein, it still exert import biological function, and play important role in the pathogenesis of a wide varity of tumors; However, the role of pseudogenes in esophageal squamous cell carcinoma (ESCC) is poorly understood.


TUSC2P function in ESCC were explored using both in vitro and in vivo experiments cell proliferation, invasion and apoptosis assay was performed to evaluated the effect of TUSC2P on the tumor biology of ESCC. Expression of relative genes was assessed by quantitative real-time PCR (qRT-PCR) and western blotting in EC109 and TE-1 cell, as well as ESCC patients. 3’UTR luciferase assay was used to confirm the direct binding of miRNAs with TUSC2 and TUSC2P 3’UTR. Relation betweenTUSC2P, TUSC2 and ESCC prognosis was predicted by survival analysis (n = 56).


Pseudogene TUSC2P was down regulated in ESCC tissues compared with paired normal adjacent tissues, and the expression of TUSC2P was significantly correlated with survivalof ESCC patients. Over expression of TUSC2P in EC109 and TE-1 cells resulted in altered expression of TUSC2, thus inhibited proliferation, invasion and promoted apoptosis. Dual luciferase assay demonstrated that TUSC2P 3’UTR decoyed miR-17-5p, miR-520a-3p, miR-608, miR-661 from binding to TUSC2.


TUSC2P can suppresses the tumor function of esophageal squamous cell carcinoma by regulating TUSC2 expression and may also serve as a prognostic factor for ESCC patients.

Peer Review reports


Esophageal cancer is one of the most leading causes of cancer-related death worldwide, and esophageal squamous cell carcinoma (ESCC) is histologically the most frequent type of esophageal cancer [1]. Esophageal cancer is aggressive and has poor prognosis, and the overall 5-year survival rate of esophageal cancer is about 20% [2, 3]. Additionally, the incidence of ESCC is increasing rapidly [3]. ESCC is a multifactorial, multistage and complex pathological process, which involves the interaction of numerous oncogenes and tumor suppressor genes, as well as transcriptional and post-transcriptional levels of regulation. In the last two decades, non-coding RNA and the importance of its post transcriptional regulation has drawn extensive attentions. The non-coding RNAs family broadened as more new classes of non-coding RNAs have been identified and named, for example lncRNA, miRNA, circRNA, tiRNA etc. [4,5,6,7].

For the last few years, depending on the development of large scale sequencing technique and bioinformatics analysis, thousands of pseudogenes were identified [8]. Pseudogenes is a special form of RNA, which resembles mRNA but lost its protein coding function, thus belongs to the family of long non-coding RNA. In most of the cases, pseudogenes are duplicated DNA sequences of its corresponding parent coding gene, and can be transcribed along with other genes into RNA fragments. However the transcripts usually are incomplete or with mutations such as point or frameshift, which resulting in mutant coding-function [9,10,11]. Interestingly, the RNA fragments derived from pseudogenes can exert biological functions. They play important role in the post-transcriptional regulation at multiple levels including DNA, RNA and protein, and been proven involved in diverse physiological and pathological processes such as carcinogenesis [11, 12]. Though thousands of pseudogenes have been sequenced, very few of them have been functionally characterized so far.

Tumor suppressor candidate-2 (TUSC2) is a novel tumor suppressor gene that found played important role in the pathogenesis of cancer. TUSC2 was firstly reported in a study analyzing frequent deletions in the short arm of chromosome 3p. Frequent deletions in the short arm of chromosome 3p occurs in a wide variety of cancers [13, 14]. Exogenous expression of TUSC2 in non-small cell lung carcinoma cells significantly inhibited tumor cell growth by activating the apoptotic protease activating factor 1 (Apaf-1) [15]. Intravenous systemic delivery of TUSC2 to distant tumors, via intravenous cholesterol nanovesicles, suppressed tumor growth and progression in orthotopic human lung cancer xenograft models [16]. Numerous TUSC2-targeting miRNAs have been predicted by bioinformatics analysis and many targeting has already been validated, which indicating that TUSC2 abundance in cancer cells is largely dependent on post-transcriptional regulation.

Pseudogene of TUSC2 named TUSC2P, sequence of which shared 89% homology with the 3’-UTR of TUSC2. In addition, they shared binding sites of many miRNAs, including miRNA-17-5p, miRNA-608, miRNA-661, miRNA-520a-3p. In this study, we investigated the relationship between TUSC2 and its pseudogene TUSC2P and the potential role of common binding miRNAs in ESCC.


Tissue samples

Fifty-six human esophagus tumor samples were consecutively collected from Zhang Zhou Hospital, and The First Affiliated Hospital of Fujian Medical University between September 2012 and March 2013, and The Fujian Provincial Cancer Hospital between September 2014 and January 2015. Forty-nine oral cancer samples were consecutively collected from The First Affiliated Hospital of Fujian Medical University, within a period from January 2010 to December 2016.

Construct generation

The TUSC2 and TUSC2P 3’-UTR fragments were synthesized. The synthesized fragments and pcDNA3.1 vector were digested with NheI and ApaI. The digested fragments were then inserted into pcDNA3.1 vector to obtain 3’-UTR over-expression construct.

The fragment of 3’UTR of TUSC2 and TUSC2P which contained potential binding sites of miRNA-17-5p, miRNA-608, miRNA-661, miR-520a-3p were synthesized. The synthesized fragments and psiCHECKTM-2 luciferase vector were digested with XhoI and NotI. The digested fragments were then inserted into the opened psiCHECKTM-2 luciferase vector to obtain the luciferase constructs TUSC2–1, TUSC2–2, TUSC2–3 and TUSC2P-1, TUSC2P-2, TUSC2P-3 respectively.

To obtain a negative control, the corresponding fragments of 3’UTR of TUSC2 and TUSC2P 3’UTR of which the miRNAs binding sites were mutated was synthesized and cloned to psiCHECK™-2 luciferase vector to obtain corresponding mutant luciferase constructs.


EC109 cells or TE-1 cells were seeded at six-well plates for 3’-UTR over-expression construct transduction by using Lipofectamine 2000 (Invitrogen). Transfected ESCC cells were selected by G418 to gain ESCC cells with stable over-expression of TUSC2 and TUSC2P.

EC109 or TE-1 cells were seeded in 12-well dishes and cultured for 24 h, followed by transfection of TUSC2 siRNAs or miRNAs mimics/ inhibitors using Lipofectamine 2000 (Invitrogen). The concentration for siRNAs is 200 nmol/L, and concentrations for miRNAs mimics/ inhibitors is 400 nmol/L.

Proliferation assay

EC109 cells or TE-1 cells were transfected for 6 h, then were trypsinized, resuspended in 50 ml DMEM medium and then seeded in 5 sets of 12-well plate. Within 12-well plate, each treatment group were seeded triplicately. The cells were observed for 5 days. For each of the following day (d1, d2, d3, d4, d5), one set of plate was fixed using 10% formalin solution and kept in PBS at 4 °C. At the end of the observation, all the wells were stained with crystal violet, and then all the well were lysised with 10% acetic acid, optical density of each well was detected at 590 nm. The experiment were repeated for three times.

Cell invasion assay

Matrigel was diluted with DMEM medium in 1:10 and coated in cell transwell membrane inserts, 100 ml matrigel for each insert. Cell transwell membrane inserts were set in 24-well plates.

EC109 cells or TE-1 cells were suspended in serum free DMEM media and seeded in transwell membrane inserts, 2× 103 cell in 100 μl free DMEM media per inserts. 500 μl DMEM media with 10% FBS were added into the wells of the plate which at bottom of the inserts. Cells were incubated at 37 °C and allow cell migrate through the matrigel and the membrane in the inserts into bottom side of the membrane. After 24 h, the inserts were taken out, the cell contain DMEM media and matrigel were removed, while the cells on the bottom side of the membrane were fixed with 10% formalin solution and stained with crystal violet. Cells were photographed and counted.

Apoptosis analysis

Sub-confluent cells were treated with 5 μM H2O2 for 40 min to induce oxidative stress. Thereafter, cells were subjected to Apoptosis assays. An Annexin V-FITC apoptosis detection kit (Biovision Inc) was used to detect apoptotic activity. Cells (1 × 106) were collected and resuspended in binding buffer, and incubated with Annexin V-FITC and propidium iodide in the dark for 15 min. Annexin V-FITC binding was determined by flow cytometry (Ex = 488 nm; Em = 530 nm) using FITC signal detector (FL1) and propidium staining by the phycoerythrin emission signal detector (FL2).

Luciferase activity assays

Luciferase activity assays were performed using a dual-luciferase reporter system developed by Promega. In brief, EC109 cell were cultured in 24-well plate at a density of 3 × 10 4 cells per well. Cells were cultured at 37 °C for 24 h, TUSC2 and TUSC2P luciferase constructs or corresponding mutant constructs were co-transfected with miRNA mimics, or miRNAs controls, respectively using Lipofectamine 2000. After 12 h, cells were collected and lysed. Luciferase activity was measured using the Dual-Luciferase reporter assay system (Promega) according to the manufacturer’s instructions.

RNA analysis

Quantification of TUSC2 mRNA transcripts was performed by SYBR Green quantitative real-time PCR using the ABI Prism 7500 sequence detection system (Applied Biosystems) with normalization to the expression of GAPDH. MiRNAs were detected by specific stem-loop primers using miScript Reverse Transcription Kit (Qiagen) and miScriptSYBR Green PCR Kit (Qiagen) with normalization to the expression of human-U6RNA. All primer sequences are provided (see Additional file 1).

Western blot

EC109 cells or TE-1 cells were collected and lysed. Protein concentration was determined by BCA assay kit. Total protein (20 mg/well) was subjected to 12% SDS–PAGE and transferred to a nitrocellulose membrane, and hybridized with anti-TUSC2 (Abcam, Ab70182) at a dilution of 1:500 at 4 °C overnight. Secondary antibodies (Abcam, ab98488) were dialuted at 1:2,000 and incubated with the membranes at room temperature for 2 h. After secondary antibody incubation, the blot was washed and detected using ECL kit (Millipore) in autoradiography.

β-actin (Abcam, ab179467) was used to confirm equal sample loading.

Statistical analysis

Mann–Whitney test or Student t test were used for were used to compare quantitative data among groups using SPSS (version 16.0; SPSS Inc., Chicago, IL). TUSC2 or TUSC2P expression intensity in tissue sample were transformed into log10, and correlations between TUSC2 or TUSC2P expression in tissue were calculated using Spearman rank correlation coefficients. For survival analysis, TUSC2 and TUSC2P expression levels in human esophagus cancer were classified into low and high subgroups according to their median expression value. The survival rate was calculated using the Kaplan Meier method, and the log rank test was performed for significance test of TUSC2 and TUSC2P subgroups. All tests were considered significant at P value < 0.05.


TUSC2 and TUSC2P shared the same targeting miRNAs

The DNA sequence of pseudogene TUSC2P is highly homologous to its corresponding gene TUSC2, with only eighteen mismatches for the coding sequence. TUSC2P possesses a 3’UTR that is about 1.2 kilobase, which shared 89% homology with the 3’UTR of TUSC2 (see Additional file 2). When analyzing the sequence of the TUSC2P and TUSC2 3’UTR, miRNA-17-5p, miRNA-608, miRNA-661, miRNA-520a-3p were found to poss conserved binding sites for TUSC2P and TUSC2 (Fig. 1a). Among these potential targeting miRNAs, some of them exhibit more than one potential binding site.

Fig. 1
figure 1

TUSC2P is targeted by TUSC2-targeting miRNAs. a Binding of TUSC2 targeting miRNAs to TUSC2P. Seeds and seed matches, bold; canonical pairings, solid lines; b, c TUSC2-targeting miRNA mimics decrease while miRNA inhibitors increase TUSC2 and TUSC2P mRNA abundance in ESCC cell line Ec109 and TE-1. d TUSC2-targeting miRNA mimics decrease while miRNA inhibitors increase TUSC2 protein abundance. n = 3, **P < 0.01

To explore the potential targeting regulation of these miRNAs on both TUSC2 and TUSC2P in ESCC. In EC109 and TE-1 esophageal cancer cells, transient transfection of a pool of miRNA mimics suppressed both TUSC2 and TUSC2P mRNA abundance (Fig. 1b and c), as well as TUSC2 protein abundance (Fig. 1d). In these ESCC cells, a pool of miRNA inhibitors de-repressed both TUSC2 and TUSC2P transcript levels (Fig. 1b and c), as well as TUSC2 protein expression (Fig. 1d). Luciferase activity assays indicated that the miRNA–TUSC2 (Fig. 2ad) and miRNA–TUSC2P (Fig. 2eh) interaction was direct, which indicate that TUSC2 and TUSC2P are targeted by common miRNAs and in the post-transcriptional regulation manner.

Fig. 2
figure 2

Confirmation of TUSC2 and TUSC2P targeting by luciferase activity assays. a-d Fragments of TUSC2 harbouring different target sites of miR-661 (a), miR-17-5p/520a-3p (b), miR-608 (c) were cloned into the luciferase reporter vector psiCHECKTM-2(D). Ec109 cells were co-transfected with one of these miRNAs and the luciferase constructs or a mutant constructs, in which the miRNA target sites were mutated. Luciferase activity assays indicated that all four miRNAs repressed luciferase activities when it harboured the corresponding TUSC2 fragment, which was reversed when the potential miRNA target site was mutated; e-h Luciferase activity assays indicated that all four miRNAs repressed luciferase activities when it harboured the corresponding TUSC2P fragment, which was reversed when the potential miRNA target site was mutated; n = 3, **P < 0.01

Expression of TUSC2P suppresses tumour activity in ESCC

The TUSC2P 3’-UTR over-expression construct was generated to express TUSC2P 3’UTR, and was then transfected into EC109 cells or TE-1 cells to examined the ability of TUSC2P 3’ UTR as decoy of TUSC2 -targeting miRNAs.

Indeed, TUSC2P 3’UTR over expression resulted in a de-repression of both TUSC2 transcript and protein in EC109 and TE-1 esophageal cancer cells (Fig. 3a and b). Consistent with elevated TUSC2, TUSC2P 3’UTR over expression can induce growth (Fig. 3c) and invasion inhibition (Fig. 3d), and promote apoptosis (Fig. 3e).

Fig. 3
figure 3

TUSC2P 3’ UTR exerts a tumour suppressive function. a, b TUSC2P 3’UTR transfected Ec109 and TE-1 cell showed increased TUSC2 mRNA and protein levels. c-e TUSC2P 3’UTR-transfection decreased proliferation rate and invasion, while promoted apoptosis in Ec109 and TE-1 cell. n = 3, *P < 0.05,**P < 0.01

To examine the consequences of TUSC2 and TUSC2P down-regulation, siRNAs for TUSC2 (si-TUSC2/ TUSC2P) which can suppress both TUSC2 and TUSC2P expression was transfected into EC109 cells or TE-1 cells Transfection of si-TUSC2/ TUSC2P in EC109 and TE-1 cell reduced TUSC2 expression (Fig. 4a and b), and thus accelerated cell proliferation (Fig. 4c) and invasion (Fig. 4d), decelerated apoptosis (Fig. 4e), indicating that TUSC2 and its pseudogene have important roles in tumor cell biology.

Fig. 4
figure 4

Silencing of TUSC2 and TUSC2P exerts a tumour promotion function. siRNA pools for TUSC2 which can bind to common sequences in both TUSC2 and TUSC2P were transfected in Ec109 and TE-1 cell. a, b siRNAs transfection decreased TUSC2 mRNA and protein levels. c-e siRNAs transfection increased proliferation rate and invasion, while reduced apoptosis in Ec109 and TE-1 cell. *P < 0.05

The tumour suppressive function of TUSC2P is dependent on miRNA binding

DICER −/− EC109 cells were established to investigate the potential mechanism by which TUSC2P can regulate TUSC2 expression. DICER is a critic enzyme for miRNA maturation. In DICER −/− EC109 cells, expression of DICER was silenced, leading to decreased expression level of miRNAs compared to control EC109 cells (see Additional file 3). We observed that the de-repression of TUSC2 abundance by TUSC2P 3’ UTR over-expression was blunted under absent expression of miRNAs (Fig. 5a), which supported that regulation of TUSC2 by 3’UTR of TUSC2P requires mature miRNAs. Consistent with blunted TUSC2 de-repression, effects of TUSC2P 3’UTR over expression on growth (Fig. 5b), invasion (Fig. 5c) and apoptosis (Fig. 5d) were also blunted in DICER −/− EC109 cells.

Fig. 5
figure 5

The tumour suppressive function of TUSC2P is dependent on miRNA binding. TUSC2P 3’ UTR or empty vector was transfected in parental Ec109 or Ec109 DICER−/− cells. a De-repression of TUSC2 abundance by TUSC2P 3’ UTR over-expression was blunted in DICER −/− EC109 cells. b-d Effects of TUSC2P 3’UTR over expression on growth (b), invasion (c) and apoptosis (d) were also blunted in DICER −/− EC109 cells (n = 3)

Expression and losses of TUSC2 and TUSC2P in ESCC patients

TUSC2 and TUSC2P expression was detected in tumor samples and matched adjacent normal tissue from ESCC patient (n = 56), and then validated in oral cancer samples (n = 49), using qRT-PCR. Both TUSC2 and TUSC2P expression levels were decreased in esophagus tumour samples compared with matched adjacent normal tissue (Fig. 6a and b). Additionally, TUSC2 and TUSC2P expression showed highly correlation (r = 0.90, P < 0.001 and r = 0.85, P < 0.001, respectively) in both adjacent tissue and esophagus tumour samples, which suggests that they may be co-regulated in vivo (Fig. 6c and d). In further, higher level of TUSC2 and TUSC2P expression predicted better survival in esophagus cancer patients (Fig. 6e and f). Similar pattern was observed in oral cancer samples (see Additional file 4). This finding supports our molecular observations that TUSC2 and TUSC2P exert important role in tumor cell biology, and TUSC2P can regulate TUSC2 expression.

Fig. 6
figure 6

TUSC2 and TUSC2P were depressed in ESCC tissues, and decreased levels of TUSC2 and TUSC2P indicate worse ESCC prognosis. a, b The expression level of TUSC2 and TUSC2P is lower in ESCC tissues compared with adjacent normal tissues (n = 56). c, d TUSC2P expression is highly related to TUSC2 in both ESCC tissues and adjacent normal tissues (n = 56). e, f Kaplane-Meier survival curves according to the relative expression level of TUSC2 and TUSC2P in 56 ESCC patients. Decreased level of TUSC2 and TUSC2P indicate worse ESCC prognosis


Pseudogenes are fragments of DNA with variation in sequence relative to the parent coding gene. Pseudogenes have lost at least some functionality, relative to the complete gene, in cellular gene expression or protein-coding ability. Although not being fully functional, pseudogenes are similar to other noncoding RNA, which can perform regulatory functions. Many pseudogenes have important roles in normal physiologic and abnormal pathologic process [17, 18]. Similarly to other non-coding RNAs, pseudogenes can affect the function and/or expression of their parental genes or other coding genes through mechanisms as following:1) Producing antisense transcripts thus can completely bind to its parental gene and act as endogenous siRNAs [19, 20]; 2) Functioning a competing endogenous RNAs, which can regulate the stability and translation of other RNA transcripts by competing for shared microRNAs [21, 22]; 3) Encoding shortened proteins or peptides, which exert certain functions in cell [23, 24]. Pseudogene-mediated miRNA decoys is the most functional way since all transcribed pseudogenes can in principle compete with cognate genes for miRNA binding, only few pseudogenes undergo antisense transcription or produce functional protein. However, very few pseudogenes have been functionally well elucidated so far.

TUSC2P was a newly discovered and poorly elucidated pseudogene [25]. The only related study reported that TUSC2P may exert tumor suppression role in prostate tumor [26]. However, the involvement of TUSC2P in the development of ESCC have yet not been well characterized. In the present study, for the first time, we reported a functional role for TUSC2P is relevant to TUSC2 biology in esophageal carcinoma as minute changes in TUSC2 can have tumorigenic consequences. In our analysis, we found that TUSC2P derepress TUSC2 expression by competing with miRNA binding. Loss of TUSC2P and TUSC2 can lead to accelerated proliferation and invasion, decelerated apoptosis in esophageal cancer cell, suggesting TUSC2P is a bona fide tumor suppressor gene. Actually, TUSC2P is not the first functional pseudogene as miRNA decoy, which can regulated corresponding cognate genes. Pseudogene PTENP1 was mostly elucidated and was reported to regulate cellular levels of PTEN and exert a growth-suppressive role in prostate cancer cell [27]. Pseudogene PHBP1 promotes esophageal squamous cell carcinoma proliferation by increasing its cognate gene PHB expression [28]. Pseudogene GBAP1 regulates the glucocerebrosidase gene GBA by competively binding miR-22–3p [29].

Additionally, TUSC2P is repressed in ESCC tissues compared with adjacent normal tissues, and that expression level of TUSC2P and corresponding cognate gene TUSC2 is associated with survival outcome of ESCC.

Accumulating evidence showed that pseudogenes usually expressed in a cancer specific pattern [17, 27, 30]. They have shown differential expression profile between tumor species and normal control tissues, or sometimes their expression can be tumor tissue specific [31, 32]. To date, pseudogene detection has been neglected in most of the studies about their parental genes. However, more and more pseudogenes have been identified, and they exhibit cancer specific expression pattern and are proven to be related to cancer. Thus makes it as a potential diagnostic and prognostic indicator [33,34,35,36]. For example increased levels of OCT4 pseudogene POU5F1B has been reported to promote tumor growth and predict poor prognosis in stage IV gastric cancer patients [37]. Another OCT4 pseudogene OCT4-pg4 can competitively bind to miRNA-145 thus regulate OCT4 expression. Additionally OCT4-pg4 is over expressed in hepatocellular carcinoma and predict poor prognosis [38]. In light of this, pseudogenes can not only play structural and functional roles in the tumorigenesis, but also in disease diagnosis and prognosis. Therefore, a better understanding of pseudogenes expression change may provide important clues for both etiology and prognosis of cancer.


In summary, we provided a better understanding of the biology of ESCC carcinogenesis by TUSC2P. We investigate that the TUSC2P was significantly down regulated in ESCC tissues and increased expression of TUSC2P might play a tumor repressive role in ESCC carcinogenesis by acting like ‘endogenous competitors’ of miRNAs, consequently de-repressed TUSC2 mRNA. All of these findings illustrate the important roles of TUSC2P in ESCC carcinogenesis and the potential role of TUSC2P as a novel biomarker for ESCC.



Esophageal squamous cell carcinoma




Octamer-binding protein 4




POU domain class 5 transcription factor 1B


Phosphatase and tensin homolog


Tumor suppressor candidate-2


Tumor suppressor candidate-2 pseudogene


  1. Pennathur A, Gibson MK, Jobe BA, Luketich JD. Oesophageal carcinoma. Lancet. 2013;381:400–12.

    Article  Google Scholar 

  2. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87–108.

    Article  Google Scholar 

  3. Hur C, Miller M, Kong CY, Dowling EC, Nattinger KJ, Dunn M, Feuer EJ. Trends in esophageal adenocarcinoma incidence and mortality. Cancer. 2013;119:1149–58.

    Article  Google Scholar 

  4. Guil S, Esteller M. RNA-RNA interactions in gene regulation: the coding and noncoding players. Trends Biochem Sci. 2015;40:248–56.

    Article  CAS  Google Scholar 

  5. Schmitz SU, Grote P, Herrmann BG. Mechanisms of long noncoding RNA function in development and disease. Cell Mol Life Sci. 2016;73:2491–509.

    Article  CAS  Google Scholar 

  6. Di Leva G, Croce CM. Roles of small RNAs in tumor formation. Trends Mol Med. 2010;16:257–67.

    Article  CAS  Google Scholar 

  7. Garzon R, Calin GA, Croce CM. MicroRNAs in Cancer. Annu Rev Med. 2009;60:167–79.

    Article  CAS  Google Scholar 

  8. Shi X, Nie F, Wang Z, Sun M. Pseudogene-expressed RNAs: a new frontier in cancers. Tumour Biol. 2016;37:1471–8.

    Article  CAS  Google Scholar 

  9. Pink RC, Carter DR. Pseudogenes as regulators of biological function. Essays Biochem. 2013;54:103–12.

    Article  CAS  Google Scholar 

  10. Wang KC, Chang HY. Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011;43:904–14.

    Article  CAS  Google Scholar 

  11. Pavlicek A, Gentles AJ, Paces J, Paces V, Jurka J. Retroposition of processed pseudogenes: the impact of RNA stability and translational control. Trends Genet. 2006;22:69–73.

    Article  CAS  Google Scholar 

  12. Zheng LL, Zhou KR, Liu S, Zhang DY, Wang ZL, Chen ZR, Yang JH, Qu LH. dreamBase: DNA modification, RNA regulation and protein binding of expressed pseudogenes in human health and disease. Nucleic Acids Res. 2018;46:85–91.

    Article  Google Scholar 

  13. Lin J, Xu K, Gitanjali J, Roth JA, Ji L. Regulation of tumor suppressor gene FUS1 expression by the untranslated regions of mRNA in human lung cancer cells. Biochem Biophys Res Commun. 2011;410:235–41.

    Article  CAS  Google Scholar 

  14. Prudkin L, Behrens C, Liu DD, Zhou X, Ozburn NC, Bekele BN, Minna JD, Moran C, Roth JA, Ji L, Wistuba II. Loss and reduction of FUS1 protein expression is a frequent phenomenon in the pathogenesis of lung cancer. Clin Cancer Res. 2008;14:41–7.

    Article  CAS  Google Scholar 

  15. Yang Q, Yoshimura G, Mori I, Sakurai T, Kakudo K. Chromosome 3p and breast cancer. J Hum Genet. 2002;47:453–9.

    Article  CAS  Google Scholar 

  16. Ito I, Ji L, Tanaka F, Saito Y, Gopalan B, Branch CD, Xu K, Atkinson EN, Bekele BN, Stephens LC, Minna JD, Roth JA, Ramesh R. Liposomal vector mediated delivery of the 3p FUS1 gene demonstrates potent antitumor activity against human lung cancer in vivo. Cancer Gene Ther. 2004;11:733–9.

    Article  CAS  Google Scholar 

  17. Kalyana-Sundaram S, Kumar-Sinha C, Shankar S, Robinson DR, Wu YM, Cao X, Asangani IA, Kothari V, Prensner JR, Lonigro RJ, Iyer MK, Barrette T, Shanmugam A, Dhanasekaran SM, Palanisamy N, Chinnaiyan AM. Expressed pseudogenes in the transcriptional landscape of human cancers. Cell. 2012;149:1622–34.

    Article  CAS  Google Scholar 

  18. Poliseno L. Pseudogenes: newly discovered players in human cancer. Sci Signal. 2012;5:re5.

    Article  Google Scholar 

  19. Tam OH, Aravin AA, Stein P, Girard A, Murchison EP, Cheloufi S, Hodges E, Anger M, Sachidanandam R, Schultz RM, Hannon GJ. Pseudogene-derived small interfering RNAs regulate gene expression in mouse oocytes. Nature. 2008;453:534–8.

    Article  CAS  Google Scholar 

  20. Guo X, Zhang Z, Gerstein MB, Zheng D. Small RNAs originated from pseudogenes: cis- or trans-acting? PLoS Comput Biol. 2009;5:e1000449.

    Article  Google Scholar 

  21. Marques AC, Tan J, Lee S, Kong L, Heger A, Ponting CP. Evidence for conserved post-transcriptional roles of unitary pseudogenes and for frequent bifunctionality of mRNAs. Genome Biol. 2012;13:R102.

    Article  Google Scholar 

  22. Sumazin P, Yang X, Chiu HS, Chung WJ, Iyer A, Llobet-Navas D, Rajbhandari P, Bansal M, Guarnieri P, Silva J, Califano A. An extensive microRNA-mediated network of RNA-RNA interactions regulates established oncogenic pathways in glioblastoma. Cell. 2011;147:370–81.

    Article  CAS  Google Scholar 

  23. Kandouz M, Bier A, Carystinos GD, Alaoui-Jamali MA, Batist G. Connexin43 pseudogene is expressed in tumor cells and inhibits growth. Oncogene. 2004;23:4763–70.

    Article  CAS  Google Scholar 

  24. Betran E, Wang W, Jin L, Long M. Evolution of the phosphoglycerate mutase processed gene in human and chimpanzee revealing the origin of a new primate gene. Mol Biol Evol. 2002;19:654–63.

    Article  CAS  Google Scholar 

  25. Skaletsky H, Kuroda-Kawaguchi T, Minx PJ, Cordum HS, Hillier L, Brown LG, Repping S, Pyntikova T, Ali J, Bieri T, Chinwalla A, Delehaunty A, Delehaunty K, Du H, Fewell G, Fulton L, Fulton R, Graves T, Hou SF, Latrielle P, Leonard S, Mardis E, Maupin R, McPherson J, Miner T, Nash W, Nguyen C, Ozersky P, Pepin K, Rock S, Rohlfing T, Scott K, Schultz B, Strong C, Tin-Wollam A, Yang SP, Waterston RH, Wilson RK, Rozen S, Page DC. The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes. Nature. 2003;423:825–37.

    Article  CAS  Google Scholar 

  26. Rutnam ZJ, Du WW, Yang W, Yang X, Yang BB. The pseudogene TUSC2P promotes TUSC2 function by binding multiple microRNAs. Nat Commun. 2014;5:2914.

    Article  Google Scholar 

  27. Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010;465:1033–8.

    Article  CAS  Google Scholar 

  28. Feng F, Qiu B, Zang R, Song P, Gao S. Pseudogene PHBP1 promotes esophageal squamous cell carcinoma proliferation by increasing its cognate gene PHB expression. Oncotarget. 2017;8:29091–100.

    PubMed  PubMed Central  Google Scholar 

  29. Straniero L, Rimoldi V, Samarani M, Goldwurm S, Di Fonzo A, Kruger R, Deleidi M, Aureli M, Solda G, Duga S, Asselta R. The GBAP1 pseudogene acts as a ceRNA for the glucocerebrosidase gene GBA by sponging miR-22-3p. Sci Rep. 2017;7:12702.

    Article  Google Scholar 

  30. Kong Y, Zhang L, Huang Y, He T, Zhang L, Zhao X, Zhou X, Zhou D, Yan Y, Zhou J, Xie H, Zhou L, Zheng S, Wang W. Pseudogene PDIA3P1 promotes cell proliferation, migration and invasion, and suppresses apoptosis in hepatocellular carcinoma by regulating the p53 pathway. Cancer Lett. 2017;407:76–83.

    Article  CAS  Google Scholar 

  31. Hwang SL, Chang JH, Cheng CY, Howng SL, Sy WD, Lieu AS, Lin CL, Lee KS, Hong YR. The expression of rac1 pseudogene in human tissues and in human brain tumors. Eur Surg Res. 2005;37:100–4.

    Article  CAS  Google Scholar 

  32. Gong T, Zheng S, Huang S, Fu S, Zhang X, Pan S, Yang T, Sun Y, Wang Y, Hui B, Guo J, Zhang X. PTENP1 inhibits the growth of esophageal squamous cell carcinoma by regulating SOCS6 expression and correlates with disease prognosis. Mol Carcinog. 2017;56:2610–9.

    Article  CAS  Google Scholar 

  33. Wezel F, Pearson J, Kirkwood LA, Southgate J. Differential expression of Oct4 variants and pseudogenes in normal urothelium and urothelial cancer. Am J Pathol. 2013;183:1128–36.

    Article  CAS  Google Scholar 

  34. Lian Y, Xu Y, Xiao C, Xia R, Gong H, Yang P, Chen T, Wu D, Cai Z, Zhang J, Wang K. The pseudogene derived from long non-coding RNA DUXAP10 promotes colorectal cancer cell growth through epigenetically silencing of p21 and PTEN. Sci Rep. 2017;7:7312.

    Article  Google Scholar 

  35. Vitiello M, Tuccoli A, Poliseno L. Long non-coding RNAs in cancer: implications for personalized therapy. Cell Oncol (Dordr). 2015;38:17–28.

    Article  CAS  Google Scholar 

  36. Wei Y, Chang Z, Wu C, Zhu Y, Li K, Xu Y. Identification of potential cancer-related pseudogenes in lung adenocarcinoma based on ceRNA hypothesis. Oncotarget. 2017;8:59036–47.

    PubMed  PubMed Central  Google Scholar 

  37. Hayashi H, Arao T, Togashi Y, Kato H, Fujita Y, De Velasco MA, Kimura H, Matsumoto K, Tanaka K, Okamoto I, Ito A, Yamada Y, Nakagawa K, Nishio K. The OCT4 pseudogene POU5F1B is amplified and promotes an aggressive phenotype in gastric cancer. Oncogene. 2015;34:199–208.

    Article  CAS  Google Scholar 

  38. Wang L, Guo ZY, Zhang R, Xin B, Chen R, Zhao J, Wang T, Wen WH, Jia LT, Yao LB, Yang AG. Pseudogene OCT4-pg4 functions as a natural micro RNA sponge to regulate OCT4 expression by competing for miR-145 in hepatocellular carcinoma. Carcinogenesis. 2013;34:1773–81.

    Article  CAS  Google Scholar 

Download references


We are grateful to the subjects participated in the study. We thank for the staff members of The First Affiliated Hospital of Fujian Medical University and The Fujian Provincial Cancer Hospital for contributing to the collection of tumor tissue and data.


This work was supported by the Natural Science Foundation of Fujian Province [grant numbers2016J01356]; the Nature Science Foundation of Fujian Province (No.2015 J01305); and the National Key R&D Program of China (No.2017YFC0907100).

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its Additional files].

Author information

Authors and Affiliations



FQL and RJG performed all the experimental parts regarding to the EC109 and TE-1 cell, and FQL was a major contributor in writing the manuscript. BCH and FC analyzed and interpreted the patient data regarding the ESCC patients and oral cancer patients. ZJH provided all the patients sample and clinical information and proofed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Zhijian Hu.

Ethics declarations

Ethics approval and consent to participate

Experiments using freshly frozen clinical specimens were approved by the Ethics Committee of Fujian Medical University. The IRB approval information: No 95 (2014) from Ethics Committee of Fujian Medical University (Number 1, Xuefubei Road, Fuzhou, Fujian, China.). Written informed consent was obtained from all subjects.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Primers of mRNA and miRNA. Primers sequences of mRNA and miRNA used for PCR in the manuscript. (DOCX 16 kb)

Additional file 2:

Sequence Alignment of TUSC2 and TUCS2P. TUSC2P possesses a 3’UTR that is about 1.2 kilobase, which shared 89% homology with the 3’UTR of TUSC2. (DOCX 17 kb)

Additional file 3:

Silenced expression of miRNAs in DICER −/− EC109 cells. Expression of miR-17-5p, miR-520a-3p, miR-608, miR-661 in were silenced DICER −/− EC109 cells. (JPG 53 kb)

Additional file 4:

TUSC2 and TUSC2P were more lowly expressed in oral cancer tissues, and decreased levels of TUSC2 and TUSC2P indicate worse oral cancer prognosis. (A-B) The expression level of TUSC2 and TUSC2P is lower in oral cancer tissues compared with adjacent normal tissues (n = 49). (C-D) TUSC2P expression is related to TUSC2 in both oral cancer tissues and adjacent normal tissues (n = 49). (E-F) Kaplane-Meier survival curves according to the relative expression level of TUSC2 and TUSC2P in 49 oral cancer patients. Decreased level of TUSC2 and TUSC2P indicate worse oral cancer prognosis. (JPG 797 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, F., Gong, R., He, B. et al. TUSC2P suppresses the tumor function of esophageal squamous cell carcinoma by regulating TUSC2 expression and correlates with disease prognosis. BMC Cancer 18, 894 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: