Abstract
Long non-coding RNAs (lncRNAs) are emerging in molecular biology as crucial regulators of cancer. The efficacy of doxorubicin--based chemotherapy in osteosarcoma (OS) is usually limited by acquired drug resistance. To explore the mechanism of chemoresistance of OS in terms of lncRNA, using a human lncRNA-mRNA combined microarray, we identified 3,465 lncRNAs (1,761 up and 1,704 down) and 3,278 mRNAs (1,607 up and 1,671 down) aberrantly expressed in all three sets of doxorubicin-resistant MG63/DXR and their paired parental MG63 cells (fold-change >2.0, P<0.05 and FDR <0.05). Fifteen randomly selected lncRNAs were dysregulated in MG63/DXR cells relative to MG63 cells by qRT-PCR detection, which were consistent with our microarray data. Bioinformatics analysis identified that classical genes and pathways involved in cell proliferation, apoptosis, and drug metabolism were differently expressed in these cell lines. A lncRNA-mRNA co-expression network identified lncRNAs, including ENST00000563280 and NR-036444, may play a critical role in doxorubicin-resistance of OS by interacting with important genes such as ABCB1, HIF1A and FOXC2. Besides, we found that lncRNA ENST00000563280 was distinctly increased in specimens of OS patients with a poor chemoresponse compared to those with a good chemoresponse and the patients of lower expression of it may survive longer than those of higher expression, which suggest that it may serve as a biomarker to predict the chemoresponse and prognosis of osteosarcoma patients. These results provide important insights about the lncRNAs involved in osteosarcoma chemoresistance and lay a solid foundation for uncovering the mechanism ultimately.