BackgroundAcquired resistance to 5-fluorouracil (5-FU) frequently results in chemotherapy failure and disease recurrence in advanced colorectal cancer (CRC) patients. Research has demonstrated that dysregulation of long non-coding RNAs (lncRNAs) mediates the development of chemotherapy resistance in cancerous cells. The present study aims to identify key lncRNAs associated with 5-FU resistance in CRC using bioinformatic and experimental validation approaches.
MethodsThe Gene Expression Omnibus (GEO) dataset GSE119481, which contains miRNA expression profiles of the parental CRC HCT116 cell line (HCT116/P) and its in-vitro established 5-FU-resistant sub-cell line (HCT116/FUR), was downloaded. Firstly, differentially expressed microRNAs (DEmiRNAs) between the parental and 5-FU resistance cells were identified. LncRNAs and mRNAs were then predicted using online databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to uncover relevant biological mechanisms and pathways. Networks integrating lncRNAs, miRNAs, and mRNAs interactions were constructed, and topological analyses were used to identify key lncRNAs associated with 5-FU resistance. An in-vitro model of the HCT116/FUR sub-cell line was developed by exposing the HCT116/P cell line to increasing concentrations of 5-FU. Finally, real-time quantitative PCR (RT-qPCR) was performed on total RNA extracted from the HCT116/P cell line and the HCT116/FUR sub-cell line to validate the in-silico predictions of key lncRNAs.
ResultsA total of 32 DEmiRNAs were identified. Enrichment analysis demonstrated that these DEmiRNAs were mainly enriched in several cancer hallmark pathways that regulate cell growth, cell cycle, cell survival, inflammation, immune response, and apoptosis. The predictive analysis identified 237 unique lncRNAs and 123 mRNAs interacting with these DEmiRNAs. The pathway analysis indicated that most of these predicted genes were enriched in the cellular response to starvation, protein polyubiquitination, chromatin remodeling, and negative regulation of gene expression. Topological analyses of the lncRNA-miRNA-mRNA network highlighted the nuclear enriched abundant transcript 1 (NEAT1), metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), and Opa interacting protein 5 antisense RNA 1 (OIP5-AS1) as central lncRNAs. Experimental analysis by RT-qPCR confirmed that the expression levels of NEAT1 and MALAT1 were significantly increased in HCT116/FUR cells compared to HCT116/P cells. However, no significant difference was observed in the OIP5-AS1 expression level between the two cells.
ConclusionOur findings specifically highlight MALAT1 and NEAT1 as significant contributors to 5-FU resistance in CRC. These lncRNAs are promising biomarkers for diagnosing and predicting outcomes in CRC.
