library_construction_protocol: RNA was harvested with the Aurum™ Total RNA Mini Kit according to manufacturers' instructions. RNA sequencing was performed by Brooks Life Sciences Genewiz® with PolyA selection for RNA removal, 2x150 bp sequencing configuration and 20-30 million reads per sample. RNA Library Preparation and NovaSeq Sequencing by AzentaLifeSciences Genewiz: RNA samples were quantified using Qubit 4.0 Fluorometer (Life Technologies) and RNA integrity was checked with an RNA Kit on Agilent 5300 Fragment Analyzer (Agilent Technologies). RNA sequencing libraries were prepared using the NEBNext Ultra RNA Library Prep Kit for Illumina following the manufacturer's instructions (NEB, Ipswich). Briefly, mRNAs were first enriched with Oligo(dT) beads. Enriched mRNAs were fragmented for 15 minutes at 94 °C. First-strand and second strand cDNAs were subsequently synthesized. cDNA fragments were end repaired and adenylated at 3'ends, and universal adapters were ligated to cDNA fragments, followed by index addition and library enrichment by limited-cycle PCR. Sequencing libraries were validated using NGS Kit on the Agilent 5300 Fragment Analyzer (Agilent Technologies), and quantified by using Qubit 4.0 Fluorometer (Invitrogen). The sequencing libraries were multiplexed and loaded on the flowcell on the Illumina NovaSeq 6000 instrument according to manufacturer's instructions. The samples were sequenced using a 2x150 Pair-End (PE) configuration v1.5. Image analysis and basecalling were conducted by the NovaSeq Control Software v1.7 on the NovaSeq instrument. Raw sequence data (.bcl files) generated from Illumina NovaSeq was converted into fastq files and de-multiplexed using Illumina bcl2fastq program version 2.20. One mismatch was tolerated for index sequence identification.