Abstract
Circular RNAs (circRNAs) are covalently closed RNA molecules generated by the back-splicing of exons from linear precursor mRNAs. Though various linear RNAs have been shown to play important regulatory roles in many biological and developmental processes, little is known about the role of their circular counterparts. In this study, we 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 meiotic 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 nine selected circRNAs using divergent primers and Sanger sequencing. Our study presents the 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.