Using biological information to analyze potential miRNA-mRNA regulatory networks in the plasma of patients with non-small cell lung cancer

Abstract

Background: Lung cancer is the most common malignant tumor, and it has a high mortality rate. However, the study of miRNA-mRNA regulatory networks in the plasma of patients with non-small cell lung cancer (NSCLC) is insufficient. Therefore, this study explored the differential expression of mRNA and miRNA in the plasma of NSCLC patients. Methods: The Gene Expression Omnibus (GEO) database was used to download microarray datasets, and the differentially expressed miRNAs (DEMs) were analyzed. We predicted transcription factors and target genes of the DEMs by using FunRich software and the TargetScanHuman database, respectively. The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was used for GO annotation and KEGG enrichment analysis of downstream target genes. We constructed protein-protein interaction (PPI) and DEM-hub gene networks using the STRING database and Cytoscape software. The GSE20189 dataset was used to screen out the key hub gene. Using The Cancer Genome Atlas (TCGA) and UALCAN databases to analyze the expression and prognosis of the key hub gene and DEMs. Then, GSE17681 and GSE137140 datasets were used to validate DEMs expression. Finally, the receiver operating characteristic (ROC) curve was used to verify the ability of the DEMs to distinguish lung cancer patients from healthy patients. Results: Four upregulated candidate DEMs (hsa-miR199a-5p, hsa-miR-186-5p, hsa-miR-328-3p, and hsa-let-7d-3p) were screened from 3 databases, and 6 upstream transcription factors and 2253 downstream target genes were predicted. These genes were mainly enriched in cancer pathways and PI3k-Akt pathways. Among the top 30 hub genes, the expression of KLHL3 was consistent with the GSE20189 dataset. Except for let-7d-3p, the expression of other DEMs and KLHL3 in tissues were consistent with those in plasma. LUSC patients with high let-7d-3p expression had poor overall survival rates (OS). External validation demonstrated that the expression of hsa-miR-199a-5p and hsa-miR-186-5p in peripheral blood of NSCLC patients was higher than the healthy controls. The ROC curve confirmed that the DEMs could better distinguish lung cancer patients from healthy people. Conclusion: The results showed that miR-199a-5p and miR-186-5p may be noninvasive diagnostic biomarkers for NSCLC patients. MiR-199a-5p-KLHL3 may be involved in the occurrence and development of NSCLC.