RNA-Seq Analysis
Total RNA was extracted from the tissue using TRIzol(R) Reagent according to the manufacturers instructions. Then RNA quality was determined by 5300 Bioanalyser (Agilent) and quantified using the ND-2000 (NanoDrop Technologies). The sequencing library was performed on NovaSeq X Plus platform (PE150) using NovaSeq Reagent Kit. Or the sequencing library was performed on DNBSEQ-T7 platform (PE150) using DNBSEQ-T7RS Reagent Kit (FCLPE150) version 3.0. Functional-enrichment analysis including GO and KEGG were performed to identify which DEGs were significantly enriched in GO terms and metabolic pathways at Bonferroni-corrected P-value<0.05 compared with the whole-transcriptome background. GO functional enrichment and KEGG pathway analysis were carried out by Goatools and Python scipy software, respectively.