Mgo+Pb
Leaf samples from MgO-primed and unprimed plants were collected after 10 days of Pb stress. The total RNA was extracted using the RNAprep Pure Plant Kit (Tiangen, Beijing, China) according the instructions provided by the manufacturer. RNA concentration and purity was measured using NanoDrop 2000(Thermo Fisher Scientific, Wilmington, DE). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA). A total amount of 1 g RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using NEBNext UltraTM RNA Library Prep Kit for Illumina (NEB, USA) following manufacturers recommendations. Briefly, mRNA was isolated by polyA selection, fragmented, and used for cDNA synthesis. Adapters were ligated, and libraries were size-selected (200300 bp) and amplified by PCR (15 cycles). Paired-end sequencing (2  125 bp) was performed using the Illumina NovaSeq6000 platform in the PE150 mode. Raw sequencing reads were processed to remove adapters and low-quality bases (Phred score < 20) using SeqPrep and Sickle. Non-coding RNAs were filtered by alignment to the Rfam database (v9.1).  Clean reads were aligned to the Vigna radiata reference genome using HISAT2