Heat stress alters genome-wide profiles of circular RNAs in Arabidopsis

1599 novel circRNAs and 1583 heat stress-specific circRNAs were identified in Arabidopsis. Heat stress enhanced accumulation of circRNAs remarkably. Heat stress altered the sizes of circRNAs, numbers of circularized exons and alterative circularization events. A putative circRNA-mediated ceRNA networks under heat stress was established. Heat stress retards plant growth and destabilizes crop yield. The noncoding RNAs were demonstrated to be involved in plant response to heat stress. As a newly-characterized class of noncoding RNAs, circular RNAs (circRNAs) play important roles in transcriptional and post-transcriptional regulation. A few recent investigations indicated that plant circRNAs were differentially expressed under abiotic stress. However, little is known about how heat stress mediates biogenesis of circRNAs in plants. Here, we uncovered 1599 previously-unknown circRNAs and 1583 heat-specific circRNAs, by RNA-sequencing and bioinformatic analysis. Our results indicated that much more circRNAs were expressed under heat stress than in control condition. Besides, heat stress also increased the length of circRNAs, the quantity of circularized exons, and alternative circularization events. Moreover, we observed a positive correlation between expression patterns of some circRNAs and their parental genes. The prediction of ceRNA (competing endogenous RNA) networks indicated that differentially-expressed circRNAs could influence expression of many important genes, that participate in response to heat stress, hydrogen peroxide, and phytohormone signaling pathways, by interacting with the corresponding microRNAs. Together, our observations indicated that heat stress had great impacts on the biogenesis of circRNAs. Heat-induced circRNAs might participate in plant response to heat stress through the circRNA-mediated ceRNA networks.