Dataset Information

Identification of a novel prognosis-associated ceRNA network in lung adenocarcinoma via bioinformatics analysis.

ABSTRACT:

Background

Lung adenocarcinoma (LUAD) is the most common subtype of nonsmall-cell lung cancer (NSCLC) and has a high incidence rate and mortality. The survival of LUAD patients has increased with the development of targeted therapeutics, but the prognosis of these patients is still poor. Long noncoding RNAs (lncRNAs) play an important role in the occurrence and development of LUAD. The purpose of this study was to identify novel abnormally regulated lncRNA-microRNA (miRNA)-messenger RNA (mRNA) competing endogenous RNA (ceRNA) networks that may suggest new therapeutic targets for LUAD or relate to LUAD prognosis.

Methods

We used the SBC human ceRNA array V1.0 to screen for differentially expressed (DE) lncRNAs and mRNAs in four paired LUAD samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to annotate the DE lncRNAs and mRNAs. R bioinformatics packages, The Cancer Genome Atlas (TCGA) LUAD database, and Kaplan-Meier (KM) survival analysis tools were used to validate the microarray data and construct the lncRNA-miRNA-mRNA ceRNA regulatory network. Then, quantitative real-time PCR (qRT-PCR) was used to validate the DE lncRNAs in 7 LUAD cell lines.

Results

A total of 2819 DE lncRNAs and 2396 DE mRNAs (P < 0.05 and fold change ≥ 2 or ≤ 0.5) were identified in four paired LUAD tissue samples. In total, 255 of the DE lncRNAs were also identified in TCGA. The GO and KEGG analysis results suggested that the DE genes were most enriched in angiogenesis and cell proliferation, and were closely related to human cancers. Moreover, the differential expression of ENST00000609697, ENST00000602992, and NR_024321 was consistent with the microarray data, as determined by qRT-PCR validation in 7 LUAD cell lines; however, only ENST00000609697 was associated with the overall survival of LUAD patients (log-rank P = 0.029). Finally, through analysis of ENST00000609697 target genes, we identified the ENST00000609697-hsa-miR-6791-5p-RASL12 ceRNA network, which may play a tumor-suppressive role in LUAD.

Conclusion

ENST00000609697 was abnormally expressed in LUAD. Furthermore, downregulation of ENST00000609697 and its target gene RASL12 was associated with poor prognosis in LUAD. The ENST00000609697-hsa-miR-6791-5p-RASL12 axis may play a tumor-suppressive role. These results suggest new potential prognostic and therapeutic biomarkers for LUAD.

SUBMITTER: Li Y

PROVIDER: S-EPMC8611860 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundThe tumor microenvironment plays an important role in the progression and malignancy of lung adenocarcinoma and affects the immunotherapy response. There is increasing evidence that long non-coding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) have significant functions in the development and treatment response of various kinds of cancer. This study aimed to explore the association between immune-related lncRNA-microRNA (miRNA)-messenger RNA (mRNA)-ceRNA networks, and the prognosis of and immunotherapy response in lung adenocarcinoma.MethodsRNA-sequencing (RNA-seq) and miRNA-seq data from The Cancer Genome Atlas (TCGA) were used to evaluate the infiltration of immune cells in lung adenocarcinoma samples by undertaking a single-sample gene set enrichment analysis (ssGSEA) to divide the cells into high and low immune cell infiltration groups. The differentially expressed mRNA (DEmRNA) was further analyzed by a weighted gene co-expression network analysis (WGCNA), search tool for recurring instances of neighboring genes (STRING), and Cytoscape to select hub genes. The ceRNA network was constructed using Cytoscape. Additionally, survival analyses were conducted to screen out prognostic candidate genes.ResultsSeven thousand five hundred and thirty-eight mRNAs, 540 lncRNAs, and 138 miRNAs were found to be differentially expressed between the high and low immune cell infiltration groups. The two DEmRNA modules most significantly associated with immune cell infiltration were further analyzed, and four clusters, including 179 DEmRNAs, were selected based on Molecular Complex Detection (MCODE) scores. The selected DEmRNAs in the four clusters were mainly enriched in pathways involved in regulating the immune response. Ultimately, a ceRNA network of SNHG6-hsa-miR-30e-5p-CYSLTR1 was identified as being associated with the prognosis of and immunotherapy response in lung adenocarcinoma.ConclusionsThe present study extends understandings of immune-related lncRNA-miRNA-mRNA-ceRNA networks and identifies novel targets and a regulatory pathway for anti-tumor immunotherapy.

Project description:Pancreatic ductal adenocarcinoma (PDAC) is a class of the commonest malignant carcinomas. The present study aimed to elucidate the potential biomarker and prognostic targets in PDAC. The array data of GSE41368, GSE43795, GSE55643, and GSE41369 were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) and differentially expressed microRNAs (DEmiRNAs) in PDAC were obtained by using GEO2R, and overlapped DEGs were acquired with Venn Diagrams. Functional enrichment analysis of overlapped DEGs and DEmiRNAs was conducted with Metascape and FunRich, respectively. The protein-protein interaction (PPI) network of overlapped DEGs was constructed by STRING and visualized with Cytoscape. Overall survival (OS) of DEmiRNAs and hub genes were investigated by Kaplan-Meier (KM) plotter (KM plotter). Transcriptional data and correlation analyses among hub genes were verified through GEPIA and Human Protein Atlas (HPA). Additionally, miRNA targets were searched using miRTarBase, then miRNA-DEG regulatory network was visualized with Cytoscape. A total of 32 DEmiRNAs and 150 overlapped DEGs were identified, and Metascape showed that DEGs were significantly enriched in cellular chemical homeostasis and pathways in cancer, while DEmiRNAs were mainly enriched in signal transduction and Glypican pathway. Moreover, seven hub genes with a high degree, namely, V-myc avian myelocytomatosis viral oncogene homolog (MYC), solute carrier family 2 member 1 (SLC2A1), PKM, plasminogen activator, urokinase (PLAU), peroxisome proliferator activated receptor γ (PPARG), MET proto-oncogene, receptor tyrosine kinase (MET), and integrin subunit α 3 (ITGA3), were identified and found to be up-regulated between PDAC and normal tissues. miR-135b, miR-221, miR-21, miR-27a, miR-199b-5p, miR-143, miR-196a, miR-655, miR-455-3p, miR-744 and hub genes predicted poor OS of PDAC. An integrative bioinformatics analysis identified several hub genes that may serve as potential biomarkers or targets for early diagnosis and precision target treatment of PDAC.

Project description:IntroductionThe regulatory mechanisms of super enhancers (SEs) and ceRNA networks in LUAD progression are not well understood. We aimed to discover the prognostic-related ceRNA network regulated by SEs in metastatic LUAD.MethodsRNA-seq data were extracted from The Cancer Genome Atlas (TCGA) database. Differentially expressed (DE) RNAs were identified by edgeR. CeRNA network was predicted and visualized using starBase and Cytoscape. H3K27ac ChIP-seq data were derived from the Gene Expression Omnibus (GEO) database, and used for SE identification. Kaplan-Meier curve and multivariate Cox model were applied for prognostic analysis. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein-protein interaction (PPI) network were performed for functional analysis. SEs of AC074117.1 were verified by ChIP-qPCR in A549 and H1299 cells. MTT assay was performed to analyze cell proliferation. Luciferase activity assay was carried out to validate the target targeting relationships of ceRNA network.ResultsA total of 2355 DEmRNA, 483 DElncRNA and 155 DEmiRNA were identified between metastatic LUAD and adjacent normal tissues. CeRNA network consisting of 7 DElncRNAs, 18 DEmiRNAs and 15 DEmRNAs was constructed. Among the seven DElncRNAs in ceRNA network, only AC074117.1 was regulated by SEs. SE-regulated prognostic ceRNA sub-network consisting of FKBP3, E2F2, AC074117.1 and hsa-let-7c-5p was screened and verified. The overlapping co-expressed mRNAs of FKBP3, E2F2, AC074117.1 and hsa-let-7c-5p were mainly related to cell division and Fanconi anemia pathway. Genes in the ceRNA sub-network were correlated with DNA mismatch repair markers. Functional experiments proved that AC074117.1 was highly expressed in LUAD cells. AC074117.1 silencing notably inhibited proliferation of A549 and H1299 cells. Luciferase activity assay confirmed the direct relationship in AC074117.1-hsa-let-7c-5p-FKBP3/E2F2 network.ConclusionA novel prognostic ceRNA sub-network regulated by SEs was identified in metastatic LUAD. This study provided potential therapeutic targets and prognostic markers for further study of metastatic LUAD.