Non-coding circular RNAs repertoire and expression profile during Brassica rapa pollen development

Circular RNAs (circRNAs) are covalently closed long non-coding RNA (lncRNA) molecules generated by the back-splicing of exons from linear precursor mRNAs. Though linear lncRNAs have been shown to play important regulatory roles in diverse biological and developmental processes, little is known about the role of their circular counterparts. In this study, we have performed high-throughput RNA sequencing to delineate the expression profile and potential function of circRNAs during the five stages of pollen development in Brassica rapa. A total of 1180 circRNAs were detected in pollen development, of which 367 showed stage-specific expression patterns. Functional enrichment and metabolic pathway analysis showed that the parent genes of circRNAs were mainly involved in pollen-related molecular and biological processes such as mitotic and meiosis cell division, DNA processes, protein synthesis, protein modification, and polysaccharide biosynthesis. Moreover, by predicting the circRNA-miRNA network from our differentially expressed circRNAs, we found 88 circRNAs with potential miRNA binding sites suggesting their role in post-transcriptional regulation of the genes. Finally, we confirmed the back-splicing sites of 9 randomly selected circRNAS using divergent primers and Sanger sequencing. Our study presents the first systematic analysis of circular RNAs during pollen development and forms the basis of future studies for unlocking complex gene regulatory networks underpinning reproduction in flowering plants.