3.1 Aberrant expression of linc00342 in various tumors
Significant expression differences between the tumor group and the normal group were observed in the adrenocortical carcinoma (ACC), breast invasive carcinoma (BRCA), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), lymphoid neoplasm diffuse large B-cell lymphoma (DLBC), glioblastoma multiforme (GBM), kidney renal clear cell carcinoma (KIRC), acute myeloid leukemia (LAML), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), ovarian serous cystadenocarcinoma (OV), pancreatic adenocarcinoma (PAAD), pheochromocytoma and paraganglioma (PCPG), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), skin cutaneous melanoma (SKCM), stomach adenocarcinoma (STAD), testicular germ cell tumors (TGCT), thyroid carcinoma (THCA), thymoma (THYM), uterine corpus endometrial carcinoma (UCEC), uterine carcinosarcoma (UCS) in GEPIA dataset. High expression in tumor relative to normal tissue was found in CHOL, DLBC, KIRC, LAML, PAAD, PCPG and THYM; low expression in tumor relative to normal tissue was found in ACC, BRCA, COAD, GBM, LUAD, LUSC, OV, PRAD, READ, SKCM, STAD, TGCT, THCA, UCEG and UCS. (Fig. 1A, 1B, 1C, P < 0.05). Meanwhile, we downloaded the gene expression and clinical information (including survival time and survival state) of linc00342 in KIRC from the TCGA database. The expression of linc00342 in KIRC was significantly higher than in normal tissues (Fig. 1D, P < 0.001), and in the receiver operating characteristic (ROC) curve, the area under the curve (AUC) was greater than 0.5, indicating that linc00342 has significant diagnostic value in KIRC (Fig. 1E, AUC = 0.889, 95%CI, 0.860–0.919). The results of survival analysis suggested that the 5-year survival rate of patients with high expression of linc00342 in KIRC was significantly lower than that of patients with low linc00342 expression (Fig. 1F, P < 0.05).
3.2 The prognostic role of linc00342 in KIRC
In order to better study the importance of linc00342 in the prognosis of KIRC, we conducted univariate and multivariate Cox analysis on the clinical data from KIRC patients (including linc00342 gene expression, age, gender, clinical stage), and constructed a nomogram prognosis prediction model based on the multivariate Cox analysis. Univariate analysis results suggested that linc00342 gene expression and clinical stage (TMN) might be independent factors for KIRC prognosis (Fig. 2A, p < 0.01), and multivariate analysis illustrated that linc00342 expression level was an independent factor for KIRC prognosis (Fig. 2B, P < 0.01). The nomogram prediction model results suggested that KIRC patients with the higher linc00342 expression had poorer 3-year and 5-year OS (Fig. 2C). We then plotted validation curves to evaluate the prediction model. The blue line was not far from the diagonal, suggesting the predicted probabilities were consistent with the observed proportions (Fig. 2D).
3.3 GO annotation and KEGG pathway enrichment analysis
Linc00342-related genes obtained from the MEM dataset (https://biit.cs.ut.ee/mem/) were used to conduct functional analysis through the David database (https://david-d.ncifcrf.gov/). Enriched GO terms and their target genes and P values are shown in Fig. 3A, 3B, including cell cycle, cell cycle arrest, M phase, cell division, mRNA metabolic process, RNA splicing, RNA processing, DNA repair, negative regulation of gene expression, and nuclear division, all of which may be associated with the tumorigenesis of KIRC (P < 0.05). KEGG pathway analysis demonstrated that linc00342 related genes were involved in a variety of significant cancer-associated pathways, such as MAPK signaling pathway, mTOR signaling pathway, insulin signaling pathway, NOD-like receptor signaling pathway, Wnt signaling pathway and GnRH signaling pathway (P < 0.05). The detailed data of KEGG analysis are listed in Table. A network was plotted by Cytoscape to present the relationship between the GO terms, KEGG terms and their correlated genes. (Fig. 3C, 3D)
3.4 Protein-protein interaction (PPI) network and Pearson’s correlation analysis
In the PPI network (Fig. 4), nodes represent proteins, edges represent protein-protein associations, and line color indicates the type of interaction. Through correlation analysis based on the Starbase database, CCAR1 (r = 0.578, P = 6.20e-49), HNRNPA1 (r = 0.089, P = 4.07e-02), HNRNPA2B1 (r = 0.649, P = 3.31e-65), HNRNPD (r = 0.410, P = 4.62e-23), HNRNPL (r = 0.422, P = 1.62e-24), RBM5 (r = 0.805, P = 5.46e-123), SRSF5 (r = 0.698, P = 2.57e-79), SRSF7 (r = 0.542, P = 3.92e-42), SRSF11 (r = 0.795, P = 8.14e-118), and WDR33 (r = 0.411, P = 4.e-02) (Fig. 5A-5J) all have significant relationships with linc00342.
3.5 High expression of linc00342 in KIRC cells
RT-qPCR analysis showed linc00342 was remarkably upregulated in OSRC2 and ACHN KIRC cells compared with normal kidney cells. (Fig. 6A, P < 0.05).
3.6 shRNA with the highest knockdown efficiency was selected
To investigate the potential mechanism of linc00342 in tumorigenesis of KIRC, linc00342 shRNA#1, #2, #3 were transfected into OSRC2 and ACHN KIRC cells to knock down the expression of linc00342. ShRNA#2 had the highest knockdown efficiency (Fig. 6B) and was used for further studies.
3.7 Depletion of linc00342 retards proliferation in KIRC cells.
Colony-formation assays revealed decreased cell proliferation in OSRC2 and ACHN KIRC cells transfected with the linc00342 shRNA#2, compared with control group with scrambled shRNA (Fig. 7, P < 0.01). The results indicated that downregulation of linc00342 could inhibit the proliferative properties of KIRC cells.
3.8 Silencing linc00342 accelerates cell apoptosis and arrests cell cycle in G0/G1
Base on the results of flow cytometry analysis, we found that along with the knockdown of linc00342 in OSRC2 and ACHN KIRC cells, the cell apoptotic rate in linc00342-shRNA-KIRC cells was elevated compared with linc00342-scramble-KIRC cells. (Fig. 8A, P < 0.05). Furthermore, the results showed linc00342-silencing KIRC cells were arrested in G0/G1. (Fig. 8B, P < 0.05)
3.9 Knocking down linc00342 exhibited significantly reduced migration and invasion ability
By comparing the results of transfection with linc00342-scramble and linc00342-shRNA, it was revealed that after knocking down linc00342, the post-healing width of the scratches as a percentage of the original scratch width was significantly lower than in the control group, indicating linc00342 can promote the migration ability of OSRC2 and ACHN KIRC cell lines. (Fig. 9A, P < 0.05) Furthermore, similar findings were obtained via a Transwell assay, in which the density of KIRC cells transfected with the linc00342-shRNA was remarkably decreased compared with controls. (Fig. 9B, P < 0.001) Accordingly, we drew the conclusion that linc00342 had the properties that promoted migration and invasion in KIRC.