Abstract
RNA-binding proteins (RBPs) play essential roles in regulating gene expression. However, the RNA ligands of RBPs are poorly understood in plants, not least due to the lack of efficient tools for genome-wide identification of RBP-bound RNAs. An RBP-fused adenosine deaminase acting on RNA (ADAR) can edit RBP-bound RNAs, which allows efficient identification of RNA ligands of RBPs in vivo. Here, we report the RNA editing activities of the ADAR deaminase domain (ADARdd) in plants. Protoplast experiments indicated that RBP-ADARdd fusions efficiently edited adenosines within 41 nucleotides (nt) of their binding sites. We then engineered ADARdd to profile the RNA ligands of rice (Oryza sativa) Double-stranded RNA-Binding Protein 1 (OsDRB1). Overexpressing the OsDRB1-ADARdd fusion protein in rice introduced thousands of A-to-G and T-to-C RNA‒DNA variants (RDVs). We developed a stringent bioinformatic approach to identify A-to-I RNA edits from RDVs, which removed 99.7% to 100% of background single-nucleotide variants in RNA-seq data. This pipeline identified a total of 1,798 high-confidence RNA editing (HiCE) sites, which marked 799 transcripts as OsDRB1-binding RNAs, from the leaf and root samples of OsDRB1-ADARdd-overexpressing plants. These HiCE sites were predominantly located in repetitive elements, 3'-UTRs, and introns. Small RNA sequencing also identified 191 A-to-I RNA edits in miRNAs and other sRNAs, confirming that OsDRB1 is involved in sRNA biogenesis or function. Our study presents a valuable tool for genome-wide profiling of RNA ligands of RBPs in plants and provides a global view of OsDRB1-binding RNAs.