Molecular characteristics analysis
According to the lnCAR database, AL365181.3 is an intronless transcript of 3222 nucleotides (nt)[20]. Figure S1A shows the genomic information for AL365181.3. AL365181.3 is located at chr1:156641666.156644887. According to the online database, AL365181.3 did not exhibit coding potential (Fig. S1B-C). Furthermore, we found that LUAD cells contain AL365181.3, which is located primarily in the cytoplasm (Fig. S1D).
The expression of AL365181.3 varies in human tumors
To validate the involvement of AL365181.3 in regulating human tumor development, we verified the RNA expression pattern of AL365181.3 in a variety of cancers and showed that the expression level of AL365181.3 varied greatly among different tumors, with AL365181.3 RNA being significantly highly expressed in 15 cancers in the TCGA database (Fig. 1A). We also confirmed that AL365181.3 was highly expressed in 13 paired adjacent normal tissues (Fig. 1B). These findings suggest that AL365181.3 can play a procancer role in different types of tumors.
Prognostic value of AL365181.3 in human tumors
AL365181.3 expression differs among different types of cancer, so its prognostic value in human cancer was investigated. We found that the expression level of AL365181.3 was related to the overall survival (OS) time of patients with four kinds of tumors (Fig. 2A), disease-specific survival (DSS) of patients with 4 kinds of tumors (Fig. 2B) and progression-free interval (PFI) of patients with 8 kinds of tumors (Fig. 2C).
ROC curve analysis of AL365181.3 in human tumors
We investigated whether AL365181.3 can serve as a biomarker for human tumors. According to the ROC curve analysis, AL365181.3 could be used to diagnose 20 types of tumors with high sensitivity and specificity (AUC > 0.75) (Fig. S2A-E).
AL365181.3 is highly expressed in LUAD tissues and is correlated with adverse clinical parameters in LUAD patients
Using AL365181.3 expression data from TCGA, we examined AL365181.3 expression in LUAD, and we found that AL365181.3 was more highly expressed in LUAD (Fig. 3A-B). Moreover, Gene Expression Omnibus (GEO) data provided similar results (Fig. 3C-D). We also examined the clinical relevance of AL365181.3 in LUAD using the TCGA LUAD dataset. The pathological stage and TNM stage were significantly correlated with AL365181.3 expression (Fig. 3E-H, Table 1). Importantly, in LUAD patients with many clinical features, including TNM stage, pathological stage, residual tumor, cancer location, anatomic neoplasm subdivision, age, smoking status, and number of peak years smoked, those with high AL365181.3 expression had significantly worse OS (Fig. 3I-R).
Table 1
Correlation of AL365181.3 expression with clinicopathologic features in LUAD patients.
Characteristics | Low expression of AL365181.3 | High expression of AL365181.3 | p |
n | 269 | 270 | | |
Pathologic T stage, n (%) | | | 0.009 | |
T1&T2 | 244 (45.5%) | 224 (41.8%) | | |
T3&T4 | 24 (4.5%) | 44 (8.2%) | | |
Pathologic N stage, n (%) | | | 0.007 | |
N0 | 187 (35.8%) | 163 (31.2%) | | |
N1&N2&N3 | 71 (13.6%) | 102 (19.5%) | | |
Pathologic M stage, n (%) | | | 0.066 | |
M0 | 186 (47.7%) | 179 (45.9%) | | |
M1 | 8 (2.1%) | 17 (4.4%) | | |
Pathologic stage, n (%) | | | < 0.001 | |
Stage I | 172 (32.4%) | 124 (23.4%) | | |
Stage II | 47 (8.9%) | 78 (14.7%) | | |
Stage III&Stage IV | 47 (8.9%) | 63 (11.9%) | | |
A nomogram was also constructed using the TCGA-LUAD cohort to predict OS, DSS, and the PFI. The prognostic indicators in the nomogram included AL365181.3 expression and pathological stage (Fig. 4A–C). LUAD OS, DSS, and the PFI were reliably predicted by the nomogram based on calibration curves (Fig. 4D–F). In conclusion, AL365181.3 can be used as a sensitive diagnostic indicator and can be used as a promising biomarker for LUAD.
Analysis of the AL365181.3-related signaling pathways in LUAD
Using the “clusterProfiler” R package, we performed functional annotation of the AL365181.3-associated differentially expressed genes (DEGs) in LUAD patients to elucidate the mechanism of the effect of AL365181.3, and 418 DEGs (mRNAs and lncRNAs) were identified with threshold values of |log2-fold change (FC)|>2 and adjusted p value < 0.05; these genes included 29 upregulated and 282 downregulated lncRNAs and 71 upregulated and 36 downregulated mRNAs (Fig. S3A–B). The KEGG enrichment results showed that the DEGs were involved mainly in pentose and glucuronate interconversions, steroid hormone biosynthesis, ascorbate and aldarate metabolism, drug metabolism-cytochrome P450, and metabolism of xenobiotics by cytochrome P450 (Fig. S3C). Next, we performed Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the GEO dataset. The results revealed that AL365181.3 is involved in the regulation of metabolic metabolism, MAPK signaling and other tumor regulatory signaling pathways (Fig. S3D-E).
GSEA indicated that AL365181.3 was involved mainly in arginine and proline metabolism, phenylalanine metabolism and other pathways, other enzymes involved in drug metabolism, fructose and mannose metabolism, glycosaminoglycan degradation, the pentose phosphate pathway, the metabolism of cytochromes to xenobiotics, linoleic acid metabolism, the pentose phosphate pathway, nicotinic acid and nicotinamide metabolism, phenylalanine metabolism, retinol metabolism, ascorbic and aldehyde metabolism, and steroid hormone biosynthesis and tyrosine metabolism (Fig. 5).
The enrichment of hallmark genes indicated the following: early estrogen response, late estrogen response, glycolysis, heme metabolism, reactive oxygen species pathway, xenobiotic metabolism, and pancreatic beta cells (Fig. 6). These findings strongly demonstrated that AL365181.3 is involved mainly in the regulation of LUAD metabolism-related pathways.
AL365181.3 knockdown inhibits the proliferation of LUAD cells in vitro
The above results suggest that AL365181.3 expression is significantly upregulated in LUAD tissues and that AL365181.3 may affect the progression of LUAD. To investigate the biological function of AL365181.3 in LUAD, two siRNA sequences targeting AL365181.3, named siRNA-AL365181.3#1 and siRNA-AL365181.3#2, were designed and synthesized, and random nonsense siRNA sequences were selected as negative controls. The target cells were transfected separately by RNAi, and the knockdown efficiency of the transfected cell lines was examined via qRT‒PCR analysis (Fig. 7A). CCK8 and cell colony formation assays confirmed that in both A549 and H1299 cells, AL365181.3 knockdown reduced the proliferative capacity (Fig. 7B-C). In addition, knocking down AL365181.3 enhanced LUAD cell apoptosis (Fig. 7D).
AL365181.3 knockdown inhibits the migration and invasion of LUAD cells
The effect of AL365181.3 on the metastasis of LUAD cells was subsequently confirmed. Transwell and wound healing assays demonstrated that cell migration was markedly inhibited by downregulation of AL365181.3 expression (Fig. 7E-F). In conclusion, these results demonstrated that AL365181.3 has an essential function in regulating the migratory capacity of LUAD cells.
AL365181.3 knockdown inhibits LUAD cell xenograft tumor growth
Using a xenograft tumor model, we examined the effect of AL365181.3 on LUAF in vivo. Beginning on day 15, the mice were injected with NC or AL365181.3 siRNA every three days, and on day 30, they were euthanized. Compared to those of the NC cells, the tumors of the A549 cells grown with lower AL365181.3 concentrations grew at a significantly slower rate (Fig. 8A-C). Compared with those in the NC group, the weights of the si-AL365181.3#2-treated xenograft tumors were lower (Fig. 8D). qPCR revealed that AL365181.3 levels were significantly lower in si-AL365181.3-injected xenograft tumors than in control tumors (Fig. 8E). Compared with those in the NC group, the injection of si-AL365181.3#2 into xenograft tumors resulted in a decrease in the level of Ki67 expression, as observed through IHC staining (Fig. 8F). These results suggest that AL365181.3 could contribute to the oncogenic impact of LUAD.