library_construction_protocol: Total RNA extraction method:①RNA degradation and contamination was monitored on 1%agarose gels. ② RNA purity was checked using the NanoPhotometer®spectrophotometer (IMPLEN, CA, USA) . ③RNA integrity was assessed using the RNA Nano 6000 Assay Kit of theBioanalyzer 2100 system (Agilent Technologies, CA, USA). 1. Library preparation for Transcriptome sequencing: A total amount of 1 µg RNA per sample was used as input material for the RNAsample preparations. Sequencing libraries were generated using NEBNext®UltraTMRNA Library Prep Kit for Illumina® (NEB, USA) following manufacturer'srecommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attachedmagnetic beads. Fragmentation was carried out using divalent cations under elevatedtemperature in NEBNext First Strand Synthesis Reaction Buffer(5X). First strandcDNA was synthesized using random hexamer primer and M-MuLVReverseTranscriptase(RNase H - ). Second strand cDNA synthesis was subsequentlyperformed using DNA Polymerase I and RNase H. Remaining overhangs wereconverted into blunt ends via exonuclease/polymerase activities. After adenylationof 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure wereligated to prepare for hybridization. In order to select cDNAfragments of preferentially 250~300 bp in length, the library fragments were purified with AMPureXP system (Beckman Coulter, Beverly, USA). Then 3 µl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37°C for 15 min followedby5min at 95 °Cbefore PCR. Then PCR was performed with Phusion High -Fidelity DNApolymerase, Universal PCR primers and Index (X) Primer. At last, PCRproductswere purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. 2. Clustering and sequencing (GenechemExperimental Department): The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia) accordingtothe manufacturer's instructions. After cluster generation, the library preparations weresequenced on an Illumina Novaseq platform and 150 bp paired-end reads weregenerated.