Abstract
The interest in chemical RNA modifications is rapidly growing in the field of molecular biology. Dynamic and reversible alterations of N6-methyladenosine (m6A) RNA modification are responsible for a platter of structural and functional changes in healthy and cancerous cell states. Although many studies reported the link between tumor initiation/progression and m6A modulators, there are few studies exploring transcriptome-wide m6A profile of non-coding RNAs. The aim of current study was to identify glioma stem cell (GSC) specific m6A landscape of long non-coding RNAs (lncRNAs) applying MeRIP-seq approach. MeRIP-seq analysis assigned 77.9% of m6A peaks to mRNAs and 8.16% to lncRNAs. GSCs and differentiated cells showed 76.4% conservation of m6A peaks, while 19.4% were unique to GSCs. Seven novel GSC-specific m6A modified lncRNAs were identified: HRAT92, SLCO4A1-AS1, CEROX1, PVT1, AGAP2-AS1, MIAT, and novel lncRNA gene ENSG00000262223. Analysis disclosed a strong negative correlation between lncRNAs m6A modification rate and expression. MeRIP-seq analysis revealed m6A modifications on previously reported glioma-associated lncRNAs: LINC000461, HOTTIP, CRNDE, TUG1, and XIST. Moreover, current study disclosed that most highly m6A modified lncRNAs primarily contain m6A modifications close to 3' and 5' ends. Our results provide basis and insight for further studies of m6A modifications in non-coding transcriptome of GSCs.