Abstract
Chordoma is a rare bone cancer originating from embryologic notochordal remnants. Clival chordomas show different dural penetration ability, with serious dural penetration exhibiting poorer prognosis. The molecular mechanism of dural penetration is not clear. We analyzed lncRNA and mRNA profiles in 12 chordoma patients with different degrees of dural penetration using expression microarrays. The differentially expressed lncRNAs and mRNAs were used to construct a lncRNA-mRNA co-expression network. LncRNAs were classified into lincRNA, enhancer-like lncRNA, or antisense lncRNA. Biological functions for lncRNAs were predicted according to the lncRNA-mRNA network and adjacent coding genes by pathway analysis. The 2760 lncRNAs and 3988 mRNAs were differentially expressed in chordomas between two groups of patients with and without dural penetration. Possible pathway involvement of the significance among the 55 lncRNAs located in the lncRNA-mRNA network, 24 lincRNAs, 7 enhancer-like lncRNAs, and 14 antisense lncRNAs include cell adhesion, metastasis, invasion, proliferation, and apoptosis. Expression of 10 lncRNAs and mRNAs, and epidermal growth factor mRNA with two identified lncRNAs were subsequently verified by qRT-PCR in chordoma tissues. Our report predicts the biological functions of many lncRNAs which may be used as diagnostic and prognostic biomarkers as well as therapeutic targets during the process of dural penetration in chordoma.