Genome-wide identification and characterization of circular RNAs by high throughput sequencing in soybean

Circular RNAs (circRNAs) arise during pre-mRNA splicing, in which the 3' and 5' ends are linked to each other by a covalent bond. Soybean is an ancient tetraploid, which underwent two whole genome duplications. Most of soybean genes are paralogous genes with multiple copies. Although many circRNAs have been identified in animals and plants, little is known about soybean circRNAs, especially about circRNAs derived from paralogous genes. Here, we used deep sequencing technology coupled with RNase R enrichment strategy and bioinformatic approach to uncover circRNAs in soybean. A total of 5,372 circRNAs were identified, approximately 80% of which were paralogous circRNAs generated from paralogous genes. Despite high sequence homology, the paralogous genes could produce different paralogous circRNAs with different expression patterns. Two thousand and one hundred thirty four circRNAs were predicted to be 92 miRNAs target mimicry. CircRNAs and circRNA isoforms exhibited tissue-specific expression patterns in soybean. Based on the function of circRNA-host genes, the soybean circRNAs may participate in many biological processes such as developmental process, multiorganism process, and metabolic process. Our study not only provided a basis for research into the function of circRNAs in soybean but also new insights into the plant circRNA kingdom.