Three separate trifoliate leaves from 8-weeks old receiver plants were harvested in liquid N2, homogenized with pestle and mortar for RNA extraction using ReliaPrepTM RNA Miniprep System (Promega, cat. Z6012). RNA quantification was performed with Qubit™ RNA BR Assay Kit and Agilent Bioanalyzer RNA 6000 Pico Kit. Samples from 4 different plants were processed as biological replicates for each treatment. Total 400ng of each RNA sample (RIN≥7, OD260/280 between 1.9-2.1, OD260/230 between 1.5-2.0) was processed for library preparation using MGIEasy RNA Directional Library Prep Set following manufacturer's instructions. QIAseq FastSelect –rRNA Plant Kit (QIAGEN) was used to deplete rRNA prior to the library fragmentation step. Qubit® dsDNA HS Assay and Agilent HS DNA Kit was used to determine library concentration and length (average 326 nt) prior to multiplexing and circularization. The final pooled libraries (1.94 ng/µL) were sequenced on DNB G400RS high throughput sequencing machine (MGI) using FLC PE150 flowcell.