Comprehensive identification and characterization of lncRNAs and circRNAs reveal potential brown planthopper-responsive ceRNA networks in rice

Brown planthopper (Nilaparvata lugens Stål, BPH) is one of the most destructive pests of rice. Non-coding RNA plays an important regulatory role in various biological processes. However, comprehensive identification and characterization of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in BPH-infested rice have not been performed. Here, we performed a genome-wide analysis of lncRNAs and circRNAs in BPH6-transgenic (resistant, BPH6G) and Nipponbare (susceptible, NIP) rice plants before and after BPH feeding (early and late stage) via deep RNA-sequencing. A total of 310 lncRNAs and 129 circRNAs were found to be differentially expressed. To reveal the different responses of resistant and susceptible rice to BPH herbivory, the potential functions of these lncRNAs and circRNAs as competitive endogenous RNAs (ceRNAs) were predicted and investigated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Dual-luciferase reporter assays revealed that miR1846c and miR530 were targeted by the lncRNAs XLOC_042442 and XLOC_028297, respectively. In responsive to BPH infestation, 39 lncRNAs and 21 circRNAs were predicted to combine with 133 common miRNAs and compete for miRNA binding sites with 834 mRNAs. These mRNAs predictably participated in cell wall organization or biogenesis, developmental growth, single-organism cellular process, and the response to stress. This study comprehensively identified and characterized lncRNAs and circRNAs, and integrated their potential ceRNA functions, to reveal the rice BPH-resistance network. These results lay a foundation for further study on the functions of lncRNAs and circRNAs in the rice-BPH interaction, and enriched our understanding of the BPH-resistance response in rice.