Upregulated long non‑coding RNA LincIN promotes tumor progression via the regulation of nuclear factor 90/microRNA‑7/HOXB13 in esophageal squamous cell carcinoma.

Long non‑coding RNA LincIN has been reported to be overexpressed and to be involved in the metastasis of breast cancer. However, the expression and role of LincIN in esophageal squamous cell carcinoma (ESCC) remain unsolved. In the present study, LincIN expression was examined in ESCC by RT‑qPCR, and the roles of LincIN in ESCC were determined using cell growth, migration and invasion assays. In addition, the effects of LincIN on nuclear factor 90 (NF90) and microRNA/miR (miR)‑7 were examined by RNA immunoprecipitation assay, RT‑qPCR, dual‑luciferase reporter assay and western blot analysis. The results revealed that LincIN expression was significantly increased in ESCC tissues and cell lines. The increased expression of LincIN was positively associated with invasion depth, lymph node metastasis, TNM stage and a poor prognosis. Functional assays revealed that the overexpression of LincIN promoted ESCC cell growth, migration and invasion. Mechanistic analysis revealed that LincIN physically bound to NF90, enhanced the binding between NF90 and primary miR‑7 (pri‑miR‑7), and further enhanced the inhibitory effects of NF90 on miR‑7 biogenesis. Therefore, LincIN downregulated miR‑7 expression in ESCC. The expression of miR‑7 inversely correlated with that of LincIN in ESCC tissues. By downregulating miR‑7, LincIN increased the expression of HOXB13, a target of miR‑7. The overexpression of miR‑7 or the depletion of HOXB13 both attenuated the tumor‑promoting roles of LincIN in ESCC cell growth, migration and invasion. On the whole, the findings of the present study suggest that LincIN is overexpressed and plays an oncogenic role in ESCC via the regulation of the NF90/miR‑7/HOXB13 axis. Thus, LincIN may prove to be a promising prognostic biomarker and therapeutic target for ESCC.

Introduction

Esophageal carcinoma is one of the most aggressive gastrointestinal malignancies and accounts for the sixth leading cause of cancer-related mortality worldwide; there were 572,000 estimated new cases and 509,000 estimated deaths in 2018 worldwide (1). Esophageal squamous cell carcinoma (ESCC) is the major histological type, which accounts for approximately 90% of all esophageal carcinoma cases (2). Due to the difficulty in the early diagnosis of ESCC, the majority of patients with ESCC are not suitable for radical surgical resection (3). The prognoses of patients with ESCC are extremely poor (4). Further revealing the critical molecular mechanisms underlying the initiation and development of ESCC would promote the findings of novel therapeutic targets and prognostic biomarkers for ESCC.

Genome and transcriptome sequencings have found that the majority of the human genome is transcribed; however, only approximately 2% of the human genome encodes for proteins (5). Non-coding RNAs comprise >90% of the human transcriptome (6). Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs or miRs) are two classes of important regulatory non-coding RNAs, which have been revealed to play critical roles in various pathophysiological statuses (6,7). lncRNAs are a class of RNAs with >200 nucleotides in length and limited protein-coding potential (8). Accumulating evidence has indicated that a number of lncRNAs are dysregulated and play important roles in various diseases, including cancers (9-17). lncRNAs have been revealed to regulate various biological functions of cancer cells, such as growth, cycle, apoptosis, migration, invasion, senescence, drug resistance and others (18-25). The mechanisms of action of lncRNAs are diverse. lncRNAs can directly bind proteins, mRNAs, miRNAs, or DNAs, and further modulate the location, expression and function of the interacted partners (26-29).

miRNAs are another class of regulatory non-coding RNAs with 19-25 nucleotides in length (30). Consistent with lncRNAs, a number of miRNAs have been revealed to be dysregulated and to play important roles in various diseases (31-35). miRNAs negatively regulate the expression of their targets by binding to the 3′-untranslated regions (3′UTRs) of target mRNAs, and further causing translational inhibition and/or target mRNAs degradation (36-39).

lncRNA LincIN (long intergenic non-coding RNA between ITGB1 and NRP1) was recently identified to be overexpressed in human breast tumors and to be associated with a poor prognosis of patients with breast cancer (40). Furthermore, LincIN was revealed to promote breast tumor cell migration and invasion by binding nuclear factor 90 (NF90) (40). However, the expression, clinical significance, biological role and mechanisms of action of LincIN in ESCC remain unknown.

In the present study, LincIN expression was measured in ESCC by RT-qPCR, and the roles of LincIN in ESCC were detected using cell growth, migration and invasion assays. In addition, and the mechanisms of action of LincIN were investigated by RNA immunoprecipitation assay, RT-qPCR, dual-luciferase reporter assay and western blot analysis.

Materials and methods

Clinical specimens

A total of 56 pairs of ESCC tissues and adjacent noncancerous tissue specimens were acquired from patients with ESCC who received surgical resection at Nanfang Hospital, Southern Medical University (Guangzhou, China) between February, 2016 and December, 2017. The clinical specimens were immediately frozen in liquid nitrogen and stored at −80°C after surgery. All clinical specimens were diagnosed by histopathological detection. Written informed consent was acquired from all participants. The use of clinical specimens was reviewed and approved by the Nanfang Hospital Institutional Review Board (Guangzhou, China).

Cells and cell culture

The human immortalized normal esophageal epithelial cell line, Het-1A, was acquired from the American Type Culture Collection (ATCC, cat. no. CRL-2692). The ESCC cell lines, TE-1 (cat. no. TCHu 89), KYSE150 (cat. no. TCHu 236) and Eca-109 (cat. no. TCHu 69), were acquired from the Institute of Biochemistry and Cell Biology of the Chinese Academy of Sciences (Shanghai, China). Het-1A cells were maintained in BEBM (Lonza Group, Ltd.). ESCC cell lines were maintained in RPMI-1640 medium (Invitrogen; Thermo Fisher Scientific, Inc.). All cells were cultured in medium supplemented with 10% fetal bovine serum (Gibco; Thermo Fisher Scientific, Inc.) at a humidified incubator with 5% CO2 at 37°C.

Plasmid construction

LincIN full-length sequences were PCR-amplified using the Phusion Flash High-Fidelity PCR Master Mix (Thermo Fisher Scientific, Inc.) and the primers, 5′-CCC AAG CTT ACT CTG TAG TCA CCC AGG CT-3′ (sense) and 5′-GCT CTA GAT TCT GTA AAT TAA GTT TAA TGC TG-3′ (antisense). Subsequently, the PCR products were sub-cloned into the HindIII and XbaI sites of the pcDNA™3.1(+) plasmid (Invitrogen; Thermo Fisher Scientific, Inc.). In total, 2 independent oligonucleotides inhibiting LincIN expression were synthesized and inserted into the SuperSilencing shRNA (short hairpin RNA) expression plasmid pGPU6/Neo (GenePharma, Inc.), named LincIN-shRNA1 and LincIN-shRNA2. The targeted sequences of the LincIN shRNAs were as follows: 5′-GAC ATT ATG CAA GGA GAT GGC A-3′ (LincIN-shRNA1) and 5′-CAG TTG GTC ACT CTA CTC AGT-3′ (LincIN-shRNA2) (40). The HOXB13 3′UTR containing miR-7 target sites were PCR-amplified using the Phusion Flash High-Fidelity PCR Master Mix (Thermo Fisher Scientific, Inc.) and the primers, 5′-CGA GCT CCC CTT CCA TTA CAC CTC TCA C-3′ (sense) and 5′-GCT CTA GAT CCT CCT CCT CGT CCT CTT-3′ (antisense). The PCR products were then sub-cloned into the SacI and XbaI sites of pmirGLO plasmid (Promega Corporation).

Transfection of plasmids and miRNAs

miR-7 mimics and miRNAs negative control (NC) were purchased from GenePharma, Inc. HOXB13 specific shRNA (HOXB13-shRNA) and scrambled shRNA which was used as a negative control were purchased from GeneCopoeia, Inc. The transfections of plasmids (2.5 µg) and miRNAs (10 pmol) were carried out using Lipofectamine 3000 (Invitrogen; Thermo Fisher Scientific, Inc.) according to the protocol. At 48 h following transfection, subsequent experimentation was performed.

Construction of stable cell lines

The LincIN overexpression plasmid and control pcDNA3.1 plasmid were transfected into TE-1 cells. TE-1 cells were established from a well-differentiated human squamous cell carcinoma of the esophagus and had a male karyotype (41). TE-1 has a mutation at codon 272 of p53, leading to p53 activity impairment (42). At 48 h following transfection, the cells were selected with neomycin for 4 weeks to obtain LincIN stably overexpressed and control TE-1 cells. LincIN-specific shRNAs (2.5 µg) and control shRNA (2.5 µg) were transfected into the Eca-109 cells. The Eca-109 cells were derived from a human oesophageal carcinoma in 1973 (43). Eca-109 harbors wild-type p53 (44). At 72 h following transfection, the cells were selected with neomycin for 4 weeks to obtain LincIN stably-depleted and control Eca-109 cells. 2×106 transducing units of miR-7 overexpression lentivirus (GeneCopoeia, Inc.) was infected into LincIN stably overexpressing TE-1 cells. At 96 h following infection, the cells were selected with neomycin and puromycin for 4 weeks to obtain LincIN and miR-7 simultaneously stably overexpressed TE-1 cells. HOXB13-specific shRNAs (2.5 µg) were transfected into LincIN stably overexpressing TE-1 cells. At 72 h following transfection, the cells were selected with neomycin and puromycin for 4 weeks to obtain HOXB13 stably depleted and simultaneously LincIN stably overexpressing TE-1 cells.

Total RNA extraction and reverse transcription-quantitative polymerase chain reaction (RT-qPCR)

Total RNA was extracted using TRIzol reagent (Invitrogen; Thermo Fisher Scientific, Inc.) following the manufacturer's instructions. The RNA was then used to carry out reverse transcription using the M-MLV Reverse Transcriptase (Invitrogen; Thermo Fisher Scientific, Inc.). The first-stand cDNA generated from reverse transcription was used to carry out quantitative polymerase chain reaction (qPCR) with SYBR Premix Ex Taq™ II (Takara Biotechnology Co., Ltd.) on an ABI7500 System (Applied Biosystems; Thermo Fisher Scientific, Inc.) following the standard SYBR-Green PCR protocols. The qPCR cycling conditions were as follows: 95°C for 30 sec, then 40 cycles at the conditions of 95°C for 5 sec, 60°C for 20 sec, and 72°C for 20 sec. The RT-qPCR primers were as follows: LincIN, 5′-AAG AGG GGA GTG CGG AAC A-3′ (forward) and 5′-TCA CCA GGG ACA CGA GAT G-3′ (reverse); HOXB13, 5′-TGA CTC CCT GTT GCC TGT G-3′ (forward) and 5′-GAA CTT GTT AGC CGC ATA CTC-3′ (reverse); β-actin, 5′-GGG AAA TCG TGC GTG ACA TTA AG-3′ (forward) and 5′-TGT GTT GGC GTA CAG GTC TTT G-3′ (reverse). For the quantification of the expression of miR-7 and pri-miR-7, RT-qPCR was carried out using the TaqMan microRNA assays (Applied Biosystems; Thermo Fisher Scientific, Inc.) and TaqMan pri-miRNA assays (Applied Biosystems; Thermo Fisher Scientific, Inc.) respectively following the manufacturer's protocols. β-actin was used as a reference gene for the quantification of LincIN, HOXB13 and pri-miR-7. U6 was used as a reference gene for the quantification of miR-7. The quantification of RNA expression was calculated using the 2−ΔΔCq method (45).

Western blot analysis

Total proteins were extracted from the TE-1 and Eca-109 cells using RIPA buffer (Beyotime Institute of Biotechnology) supplemented with protease inhibitors (Beyotime Institute of Biotechnology). Protein concentrations were determined by bicinchoninic acid (BCA) assay with the BCA Protein Assay kit (Beyotime Institute of Biotechnology) following the manufacturer's protocol. Equal amounts of proteins (10 µg) were separated by 15% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and subsequently transferred onto PVDF membranes (EMD Millipore). After blocking with 5% non-fat milk in Tris-buffered saline with 0.1% Tween-20 (TBST) at room temperature for 2 h, the membranes were incubated with primary antibodies against HOXB13 (1:1,000, cat. no. 90944, Cell Signaling Technology, Inc.) or β-actin (1:10,000, cat T0022, Affinity Biosciences) overnight at 4°C. After washing with TBST buffer, the membranes were incubated with IRDye 680RD goat anti-mouse IgG secondary antibody (1:10,000, cat. no. 925-68070, LI-COR Biosciences) or IRDye 800CW goat anti-rabbit IgG secondary antibody (1:10,000, cat. no. 926-32211, LI-COR Biosciences) at room temperature for 1 h, followed by detection on an Odyssey infrared scanner (LI-COR Biosciences). β-actin was used as a loading control.

Cell growth assays

Cell growth was assessed using a Cell Counting kit-8 (CCK-8) and 5-ethynyl-2-deoxyuridine (EdU) incorporation assays (27). For CCK-8 assay, the TE-1 and Eca-109 cells were seeded at 2×103 cells/well into 96-well plates. After being incubated for 24, 48, or 72 at 37°C, 10 µl CCK-8 solution (Dojindo Laboratories, Inc.) were added to each well. The optical density values at 450 nm were measured with an enzyme-linked immune detector to indicate cell viability. EdU incorporation assay was performed using the Cell-Light™ EdU Apollo® 643 In Vitro Imaging kit (Guangzhou RiboBio Co., Ltd.) following the manufacturer's instructions. The results were detected using a Zeiss photomicroscope (Carl Zeiss AG) and counting was performed based on at least 10 random visual fields.

Transwell cell migration and invasion assays

Cell migration and invasion were assessed using Transwell assays (46). The TE-1 and Eca-109 cells suspended in serum-free medium were plated in the upper chamber of Transwell inserts (8 µm-pore size, EMD Millipore) with Matrigel (for invasion assay) or without Matrigel (for migration assay). Medium containing 20% fetal bovine serum was plated into the lower chamber. Following incubation for 48 h at 37°C, cells remaining on the upper chamber were wiped away using a cotton swab. Cells migrating into the lower surfaces of the Transwell inserts were fixed in methanol and stained with 0.5% crystal violet (Beyotime Institute of Biotechnology) at room temperature for 30 min. The results were detected using a Zeiss photomicroscope and counted based on at least 10 random visual fields.

RNA immunoprecipitation (RIP) assay

RIP assays were carried out in the TE-1 and Eca-109 cells using the EZ-Magna RIP™ RNA Binding Protein Immunoprecipitation kit (EMD Millipore) and primary antibody against NF90 (5 µg, ab131004, Abcam). The enriched RNA (LincIN and pri-miR-7) was quantified by RT-qPCR as described above.

Dual-luciferase reporter assay

Luciferase reporter plasmid pmirGLO containing HOXB13 3′UTR (50 ng) was co-transfected with LincIN overexpression plasmid (50 ng) and miR-7 mimics (1 pmol) into TE-1 cells. Luciferase reporter plasmid pmirGLO containing HOXB13 3′UTR (50 ng) was co-transfected with LincIN specific shRNAs (50 ng) into Eca-109 cells. Following incubation for 48 h at 37°C, the Firefly luciferase and Renilla luciferase activities were measured using the Dual-Luciferase Reporter Assay System (Promega Corporation).

Statistical analysis

GraphPad Prism software (version 5.0) was employed to carry out statistical analyses. For comparisons, the Wilcoxon signed-rank test, Pearson's Chi-squared test, Log-rank test, one-way ANOVA followed by Dunnett's multiple comparison test, the Student's t-test, and Pearson's correlation analysis were performed as indicated. A value of P<0.05 was considered to indicate a statistically significant difference.

Results

LincIN expression is increased in ESCC and is associated with a poor prognosis of patients with ESCC

To confirm the expression pattern of LincIN in ESCC, the expression of LincIN was examined in 56 pairs of ESCC tissues and adjacent non-cancerous tissues by RT-qPCR. The results revealed that the expression of LincIN was significantly increased in ESCC tissues compared with adjacent non-cancerous tissues (Fig. 1A). The analyses of the clinicopathological characteristics of these 56 ESCC cases revealed that the increased expression of LincIN was positively associated with an advanced tumor invasion depth, lymph node metastasis and TNM stage (Table I). Patients with ESCC with a higher expression of LincIN had a shorter survival time than those with a lower expression of LincIN (Fig. 1B). In addition, the expression of LincIN was significantly increased in the ESCC cell lines, TE-1, KYSE150 and Eca-109, compared with the normal esophageal epithelial cell line, Het-1A (Fig. 1C). Collectively, these results suggested that LincIN exprewssion was increased in ESCC tissues and cell lines, and was associated with an advanced clinical stage and a poor prognosis of patients with ESCC.

Figure 1

LincIN expression is increased in ESCC and is associated with a poor prognosis of patients with ESCC. (A) The expression of LincIN in 56 pairs of ESCC tissues and adjacent noncancerous tissues was detected by RT-qPCR. P<0.0001, determined by Wilcoxon signed-rank test. (B) Kaplan-Meier survival analysis of the association between LincIN expression level and the overall survival of 56 patients with ESCC. P=0.0160, determined by the log-rank test. (C) The expression of LincIN in immortalized normal esophageal epithelial cell line, Het-1A, and the ESCC cell lines, TE-1, KYSE150 and Eca-109, was detected by RT-qPCR. Results are shown as the means ± SD from 3 separate experiments. ***P<0.001, vs. control Het-1A cells, determined by one-way ANOVA followed by Dunnett's multiple comparison test vs. Het-1A group. ESCC, esophageal squamous cell carcinoma.

Table I

Association between LincIN expression and the clinicopathological characteristics of patients with ESCC (n=56).

Characteristics	Total	LincIN expression		P-value
		Low	High
Sex				1.000
Male	30	15	15
Female	26	13	13
Age (years)				0.788
≤60	25	13	12
>60	31	15	16
Differentiation				0.158
Well or moderate	37	21	16
Poor	19	7	12
Tumor invasion depth (T)				0.032
T1/T2	30	19	11
T3/T4	26	9	17
Lymph node metastasis (N)				0.029
N0	34	21	13
N1-N3	22	7	15
TNM stage				0.043
I	26	17	9
II	17	8	9
III	13	3	10

P-values were calculated using Pearson's Chi-squared test. Values in bold font indicate statistically significant differences (P<0.05). ESCC, esophageal squamous cell carcinoma.

Overexpression of LincIN promotes ESCC cell growth, migration and invasion

To explore the biological roles of LincIN in ESCC, LincIN was stably overexpressed in the TE-1 cells, which exhibited the lowest LincIN expression out of the 3 ESCC lines (Fig. 1C), by transfection with LincIN overexpression and control plasmids (Fig. 2A). CCK-8 assays revealed that the overexpression of LincIN promoted ESCC cell growth (Fig. 2B). EdU incorporation assays further verified the roles of LincIN overexpression in promoting ESCC cell growth (Fig. 2C). Transwell migration and invasion assays also revealed that the overexpression of LincIN promoted ESCC cell migration and invasion (Fig. 2D and E). Collectively, these results suggested that the overexpression of LincIN promoted ESCC cell growth, migration and invasion.

Figure 2

Overexpression of LincIN promotes ESCC cell growth, migration, and invasion. (A) Expression of LincIN in LincIN stably overexpressing and control TE-1 cells detected by RT-qPCR. Empty plasmid overexpressing nothing was used as a control. (B) Cell growth of LincIN stably overexpressing and control TE-1 cells detected by CCK-8 assays. (C) Cell growth of LincIN stably overexpressing and control TE-1 cells detected by EdU incorporation assays. Red color indicates EdU-positive cells. Scale bars, 100 µm. (D) Cell migration of LincIN stably overexpressing and control TE-1 cells was detected by Transwell migration assays. Scale bars, 100 µm. (E) Cell invasion of LincIN stably overexpressed and control TE-1 cells detected by Transwell invasion assays. Scale bars, 100 µm. Results are shown as the means ± SD from 3 separate experiments. **P<0.01, ***P<0.001, vs. control group determined by the Student's t-test. ESCC, esophageal squamous cell carcinoma.

Knockdown of LincIN inhibited ESCC cell growth, migration, and invasion

To further explore the biological roles of LincIN knockdown in ESCC, LincIN expression was stably depleted in the Eca-109 cells, which exhibited the highest expression of LincIN out of the 3 ESCC cell lines, by transfection with 2 independent LincIN-specific shRNAs (Fig. 3A). CCK-8 and EdU incorporation assays revealed that the knockdown of LincIN inhibited ESCC cell growth (Fig. 3B and C). Transwell migration and invasion assays also revealed that the knockdown of LincIN inhibited ESCC cell migration and invasion (Fig. 3D and E). Collectively, these results suggested that the knockdown of LincIN inhibited ESCC cell growth, migration and invasion.

Figure 3

Knockdown of LincIN inhibits ESCC cell growth, migration and invasion. (A) Expression of LincIN in LincIN stably depleted and control Eca-109 cells detected by RT-qPCR. Scrambled shRNA targeting nothing was used as a control. (B) Cell growth of LincIN stably depleted and control Eca-109 cells detected by CCK-8 assays. (C) Cell growth of LincIN stably depleted and control Eca-109 cells detected by EdU incorporation assays. Red color indicates EdU-positive cells. Scale bars, 100 µm. (D) Cell migration of LincIN stably depleted and control Eca-109 cells detected by Transwell migration assays. Scale bars, 100 µm. (E) Cell invasion of LincIN stably depleted and control Eca-109 cells detected by Transwell invasion assays. Scale bars, 100 µm. Results are shown as the means ± SD from 3 separate experiments. **P<0.01, ***P<0.001, vs. control-shRNA group determined by one-way ANOVA followed by Dunnett's multiple comparisons test. ESCC, esophageal squamous cell carcinoma.

LincIN enhances the suppressive effects of NF90 on miR-7 biogenesis

LincIN was previously identified as a NF90-binding lncRNA in breast cancer (40). Therefore, the present study further investigated whether LincIN also bound NF90 in ESCC. RIP assays revealed that LincIN was specifically enriched in the NF90 antibody group, suggesting the binding between LincIN and NF90 (Fig. 4A). NF90 has been reported to suppress miR-7 biogenesis by binding primary miR-7 (pri-miR-7) (47). miR-7 is a well-known tumor suppressor in a number of types of cancer, including ESCC (46,48,49). Therefore, the present study further explored the effects of LincIN on NF90/miR-7. RIP assays revealed that the overexpression of LincIN promoted, while the knockdown of LincIN suppressed the binding between NF90 and pri-miR-7 (Fig. 4B and C), suggesting that LincIN enhanced the binding between NF90 and pri-miR-7. Subsequently, the expression levels of pri-miR-7 and mature miR-7 in ESCC cells in which LincIN was stably overexpressed or depleted were measured by RT-qPCR. The results displayed that the overexpression of LincIN induced the accumulation of pri-miR-7 and the downregulation of mature miR-7 in ESCC cells (Fig. 4D and E). The knockdown of LincIN downregulated the accumulation of pri-miR-7 and upregulated mature miR-7 in ESCC cells (Fig. 4F and G). Collectively, these results suggested that LincIN enhanced the suppressive roles of NF90 on miR-7 biogenesis and therefore, downregulated the mature miR-7 level.

Figure 4

LincIN enhances the suppressive effects of NF90 on miR-7 biogenesis. (A) RIP assays were performed using the TE-1 cells, followed by RT-qPCR to detect LincIN expression associated with NF90. (B) RIP assays were performed in TE-1 cells transiently overexpressing LincIN or the control, followed by RT-qPCR to detect pri-miR-7 associated with NF90. (C) RIP assays were performed in Eca-109 cells transiently depleting LincIN or control, followed by RT-qPCR to detect pri-miR-7 expression associated with NF90. (D) Expression of pri-miR-7 in LincIN stably overexpressing and control TE-1 cells detected by RT-qPCR. (E) Expression of miR-7 in LincIN stably overexpressing and control TE-1 cells detected by RT-qPCR. (F) Expression of pri-miR-7 in LincIN stably depleted and control Eca-109 cells detected by RT-qPCR. (G) Expression of miR-7 in LincIN stably depleted and control Eca-109 cells detected by RT-qPCR. Results are shown as the means ± SD from 3 separated experiments. **P<0.01, ***P<0.001, vs. vs. IgG, control, or control-shRNA group; ns, not significant, determined by the Student's t-test (A, B, D and E) or one-way ANOVA followed by Dunnett's multiple comparisons test (C, F and G). ESCC, esophageal squamous cell carcinoma.

miR-7 expression is decreased and inversely correlates with LincIN expression in ESCC tissues

To determine whether the downregulation of miR-7 by LincIN exists in vivo, the expression of miR-7 was measured in the same 56 pairs of ESCC tissues and adjacent non-cancerous tissues used in Fig. 1. Conversely to the expression pattern of LincIN, miR-7 expression was significantly decreased in ESCC tissues compared with that in adjacent non-cancerous tissues (Fig. 5A). Furthermore, a statistically significant inverse correlation between miR-7 and LincIN expression levels was found in these 56 ESCC tissues (r=−0.4024, P=0.0021, Fig. 5B).

Figure 5

miR-7 expression is decreased and inversely correlates with LincIN expression in ESCC tissues. (A) Expression of miR-7 in 56 pairs of ESCC tissues and adjacent noncancerous tissues detected by RT-qPCR. P<0.0001, determined by Wilcoxon signed-rank test. (B) Correlation between miR-7 and LincIN expression level in 56 ESCC tissues was analysed; r=-0.4024, P=0.0021, determined by Pearson's correlation analysis. ESCC, esophageal squamous cell carcinoma.

LincIN upregulates HOXB13 expression by suppressing miR-7

HOXB13 was previously identified as a critical target of miR-7 and it was found to mediate the tumor suppressive roles of miR-7 (46,50). Therefore, the present study further examined the effects of LincIN on HOXB13. Dual-luciferase reporter assays revealed that the overexpression of LincIN significantly upregulated the 3′UTR activity of HOXB13, which was reversed by the concurrent overexpression of miR-7 (Fig. 6A and B). Reciprocally, the knockdown of LincIN significantly reduced the 3′UTR activity of HOXB13 (Fig. 6C). RT-qPCR assays revealed that the overexpression of LincIN significantly upregulated the HOXB13 mRNA level, which was reversed by the concurrent overexpression of miR-7 (Fig. 6D). Reciprocally, the knockdown of LincIN significantly reduced the HOXB13 mRNA level (Fig. 6E). Western blot assays revealed that the overexpression of LincIN markedly upregulated the HOXB13 protein level, which was reversed by the concurrent overexpression of miR-7 (Fig. 6F). Reciprocally, the knockdown of LincIN markedly reduced the HOXB13 protein level (Fig. 6G). Collectively, these results suggested that LincIN upregulated HOXB13 by suppressing miR-7.

Figure 6

LincIN upregulates HOXB13 expression by suppressing miR-7. (A) Expression of miR-7 in TE-1 cells transfected with miR-7 mimics or miRNAs negative control (NC) detected by RT-qPCR. (B) Luciferase activity in TE-1 cells co-transfected with LincIN overexpression plasmid, miR-7 mimics and luciferase reporter containing HOXB13 3′UTR. Results are shown as the relative ratio of Firefly luciferase activity to Renilla luciferase activity. (C) Luciferase activity in Eca-109 cells co-transfected with LincIN specific shRNAs and luciferase reporter containing HOXB13 3′UTR. Results are shown as the relative ratio of Firefly luciferase activity to Renilla luciferase activity. (D) Expression of HOXB13 mRNA in TE-1 cells co-transfected with LincIN overexpression plasmid and miR-7 mimics detected by RT-qPCR. (E) Expression of HOXB13 mRNA in LincIN stably depleted and control Eca-109 cells detected by RT-qPCR. (F) Expression of HOXB13 protein in TE-1 cells co-transfected with LincIN overexpression plasmid and miR-7 mimics detected by western blot analysis. (G) Expression of HOXB13 protein in LincIN stably depleted and control Eca-109 cells detected by western blot analysis. Results are shown as the means ± SD from 3 separate experiments. **P<0.01, ***P<0.001, vs. control or control-shRNA group; ns, not significant, determined by one-way ANOVA followed by Dunnett's multiple comparisons test. ESCC, esophageal squamous cell carcinoma.

Overexpression of miR-7 or depletion of HOXB13 attenuates the promoting effects of LincIN in ESCC cell growth, migration and invasion

To examine whether miR-7/HOXB13 mediate the roles of LincIN in promoting ESCC cell growth, migration and invasion, miR-7 was stably overexpressed or HOXB13 was depleted in TE-1 cells stably overexpressing LincIN (Fig. 7A). CCK-8 and EdU incorporation assays revealed that both miR-7 overexpression and HOXB13 depletion attenuated the promoting effects of LincIN on ESCC cell growth (Fig. 7B and C). Transwell migration and invasion assays also revealed that both miR-7 overexpression and HOXB13 depletion attenuated the promoting effects of LincIN on ESCC cell migration and invasion (Fig. 7D and E). Collectively, these results suggested that the overexpression of miR-7 or depletion of HOXB13 both attenuated the promoting effects of LincIN on ESCC cell growth, migration and invasion.

Figure 7

Overexpression of miR-7 or depletion of HOXB13 attenuates the roles of LincIN. (A) Expression of HOXB13 protein in HOXB13 stably depleted and control TE-1 cells detected by western blot analysis. (B) Cell growth of miR-7 and LincIN simultaneously stably overexpressed or HOXB13 stably depleted and simultaneously LincIN stably overexpressed TE-1 cells detected by CCK-8 assays. (C) Cell growth of miR-7 and LincIN simultaneously stably overexpressed or HOXB13 stably depleted and simultaneously LincIN stably overexpressed TE-1 cells detected by EdU incorporation assays. Red color indicates EdU-positive cells. Scale bars, 100 µm. (D) Cell migration of miR-7 and LincIN simultaneously stably overexpressed or HOXB13 stably depleted and simultaneously LincIN stably overexpressed TE-1 cells detected by Transwell migration assays. Scale bars, 100 µm. (E) Cell invasion of miR-7 and LincIN simultaneously stably overexpressed or HOXB13 stably depleted and simultaneously LincIN stably overexpressed TE-1 cells was detected by Transwell invasion assays. Scale bars, 100 µm. Results are shown as the means ± SD from 3 separate experiments. ***P<0.001, vs. control group; ns, not significant, determined by one-way ANOVA followed by Dunnett's multiple comparisons test. ESCC, esophageal squamous cell carcinoma.

Discussion

As one of the most lethal types of cancer, ESCC involves complex molecular mechanisms which remain unclear. Whole-exome sequencings have identified a number of mutations in ESCC, including TP53, CCND1, CDKN2A, NFE2L2, RB1, KMT2D, KMT2C, EP300, FAT1, NOTCH1 and others (51). Furthermore, transcriptome sequencing analyses have identified a number of differentially expressed mRNAs and lncRNAs in ESCC (52). The expression and roles of several lncRNAs in ESCC have also been investigated (50). lncRNA CCAT1 has been shown to be upregulated in ESCC tissues and to be associated with a poor survival (50). CCAT1 promotes ESCC cell proliferation and migration by downregulating SPRY4 and upregulating HOXB13 (50). CCAT1 binds and recruits EZH2 and SUV39H1 to epigenetically silence SPRY4 (45). Furthermore, CCAT1 upregulates HOXB13 by competitively binding miR-7 (45). lncRNA TTN-AS1 is highly expressed in ESCC, and promotes ESCC cell proliferation and metastasis via sponging miR-133b and inducing Snail1 and FSCN1 expression (53). lncRNA NMR promotes ESCC progression via NSUN2 and BPTF (54). Other lncRNAs have also been reported to have oncogenic or tumor suppressive roles in ESCC, such as PART1, CASC9, HOTTIP, LUCAT1, HNF1A-AS1, AFAP1-AS1, POU3F3 and others (55-61).

Transcriptome sequencing has identified >58,000 lncRNAs in human cells (62). Due to the large number of lncRNAs, other lncRNAs may also play biological roles in human ESCC. lncRNA LincIN, whose encoding gene is located on chromosome 10p11-12, was first identified in breast tumors (40). LincIN is overexpressed in human breast tumors, and is associated with a poor prognosis of patients with breast cancer, and is involved in the metastasis of breast cancer (40). However, the expression, clinical significance and role of LincIN in other diseases remain unknown. In the present study, it was found that LincIN expression was also significantly increased in ESCC tissues and cell lines, and was associated with an advanced clinical stage and a poor prognosis of patients with ESCC. Functional experiments revealed that the overexpression of LincIN promoted ESCC cell growth, migration and invasion. The knockdown of LincIN inhibited ESCC cell growth, migration and invasion. Thus, these data suggested LincIN as a promising prognostic biomarker and therapeutic target for ESCC.

The mechanisms of action of lncRNAs are complex. One of the main mechanisms is to bind proteins (13). In breast tumors, LincIN was revealed to bind NF90 (40). NF90, also known as ILF3, is an RNA binding protein (63). By binding double-stranded RNA, NF90 regulates gene expression and/or mRNA stabilities (64). In addition, by binding pri-miRNAs, NF90 modulates miRNAs biogenesis (47). The biogenesis of miR-7, which plays tumor suppressive roles in a number of types of cancer, has been reported to be suppressed by NF90 (46,47). Therefore, the present study further investigated the effects of LincIN on NF90/miR-7. It was verified that LincIN directly bound NF90 in ESCC cells. Moreover, LincIN promoted the binding between NF90 and pri-miR-7, thereby enhancing the suppressive roles of NF90 on miR-7 biogenesis. It was then verified that LincIN downregulated mature miR-7. Conversely to LincIN, miR-7 expression was significantly decreased in ESCC tissues. The expression of mature miR-7 was significantly and inversely associated with that of LincIN in ESCC tissues, supporting the negative regulation of miR-7 by LincIN. HOXB13 is a direct target of miR-7 and plays oncogenic roles in several types of cancer, including ESCC (46,50). By repressing miR-7, LincIN was verified to upregulate HOXB13 expression in ESCC cells. Functional rescue assays demonstrated that the overexpression of miR-7 or depletion of HOXB13 attenuated the promoting effects of LincIN on ESCC cell growth, migration and invasion, which further suggested that miR-7/HOXB13 are critical mediators of the oncogenic roles of LincIN in ESCC. Similar roles of lncRNAs in miRNAs biogenesis have also been reported in other types of cancer, such as the promotion of miR-145 biogenesis by lncRNA-ATB in bladder cancer and the disruption of miR-125b biogenesis by LINC01578 in breast cancer (63,65). The detailed mechanisms underlying the inverse effects of lncRNAs on the biogenesis of different miRNAs warrant further investigation.

Apart from miR-7, NF90 has also been reported to suppress p21 translation, upregulate hypoxia inducible factor-1α (HIF-1α) and VEGF-A, stabilize PARP1 mRNA, and facilitate DICER expression (40,64,66,67). Thus, whether LincIN modulates p21, HIF-1α, VEGF-A, PARP1 and DICER via interacting with NF90 in ESCC cells needs to be further investigated. Nevertheless, the present study identified a novel mechanism of lncRNA in ESCC and provides a novel target against ESCC.

In conclusion, the present study demonstrates that lncRNA LincIN is highly expressed in ESCC, and is associated with an advanced clinical stage and a poor prognosis of patients with ESCC. LincIN promotes ESCC cell growth, migration and invasion by binding NF90, suppressing the biogenesis of miR-7 and upregulating HOXB13. The data presented herein suggest that LincIN may be a promising prognostic biomarker and therapeutic target for ESCC.