Landscape of RNAs in human lumbar disc degeneration

Accumulating evidence indicates noncoding RNAs (ncRNAs) fine-tune gene expression with mysterious machinery. We conducted a combination of mRNA, miRNA, circRNA, LncRNA microarray analyses on 10 adults' lumbar discs. Moreover, we performed additional global exploration on RNA interacting machinery in terms of in silico computational pipeline. Here we show the landscape of RNAs in human lumbar discs. In general, the RNA-abundant landscape comprises 14,635 mRNAs (37.93%), 2,059 miRNAs (5.34%), 18,995 LncRNAs (49.23%) and 2,894 (7.5%) circRNAs. Chromosome 1 contributes for RNA transcription at most (10%). Bi-directional transcription contributes evenly for RNA biogenesis, in terms of 5′ to 3′ and 3′ to 5′. Despite the majority of circRNAs are exonic, antisense (1.49%), intergenic (0.035%), intragenic (1.69%), and intronic (6.29%) circRNAs should not be ignored. A single miRNA could interact with a multitude of circRNAs. Notably, CDR1as or ciRS-7 harbors 66 consecutive binding sites for miR-7-5p (previous miR-7), evidencing our pipeline. The majority of binding sites are perfect-matched (78.95%). Collectively, global landscape of RNAs sheds novel insights on RNA interacting mechanisms in human intervertebral disc degeneration.