Identification of microRNAs as novel biomarkers for esophageal squamous cell carcinoma: a study based on The Cancer Genome Atlas (TCGA) and bioinformatics

BACKGROUND: MicroRNAs (miRNAs) have played important roles in the regulation of gene expression in many cancers, but their roles in esophageal squamous cell carcinoma (ESCC) are still unclear. The aim of this study was to determine the potential ESCC-specific key miRNAs from a large sample dataset in The Cancer Genome Atlas (TCGA). METHODS: Integrative bioinformatics analysis was used to identify key ESCC-specific miRNAs related to the ESCC patients’ tumor histological grade and lymphatic metastasis from TCGA. Next, these key miRNA potential gene regulatory functions and relationships with ESCC patients’ clinical characteristics and overall survival were analyzed. Finally, three key miRNAs were selected randomly and quantificational real-time polymerase chain reaction (qRT-PCR) was used to validate in 51 newly diagnosed ESCC patients’ tissues samples (collected from Nov. 2017 to Feb. 2019, in Wuwei, China) whether the bioinformatics analyses results were reliable and valid. Two-tailed Student's t test, Pearson Chi-squared test and Kaplan-Meier survival analysis were used in this study. RESULTS: Thirty-five ESCC-specific miRNAs from TCGA database were investigated (fold-change > 2.0, P < 0.05), and 28 participated in the miRNAs-mRNAs co-expression network construction, while 17 were related with ESCC patients’ tumor histological grade, TNM stage, and lymphatic metastasis (P < 0.05). Meanwhile, six miRNAs (including miR-200b-3p, miR-31-5p, miR-15b-5p, miR-141-3p, miR-135b-5p, and miR-195-5p) were correlated with overall survival of ESCC patients (log-rank, P < 0.05). MiR-135b-5p, miR-15b-5p, and miR-195-5p were selected for verification of the expression levels in 51 ESCC patients’ tissue samples by using qRT-PCR. We found that the fold-changes between qRT-PCR and TCGA were completely consistent. The results also suggested that miR-135b-5p, miR-15b-5p, and miR-195-5p were significantly correlated with tumor differentiation degrees (P < 0.05), miR-195-5p was significantly correlated with tumor TNM stage (P < 0.05), and miR-135b-5p was significantly correlated with lymph-node metastasis (P < 0.05). MiR-135b-5p, miR-15b-5p, and miR-195-5p expression levels, ESCC patient clinical features association analysis results and the aforementioned TCGA bioinformatics analyses were similar. CONCLUSION: This study identified key ESCC-related miRNAs. The key miRNAs are worthy of further investigation as potential novel biomarkers for diagnosis, classification, and prognosis of ESCC.