For RNA-seq analysis, total RNA was extracted from calli collected from CIM and SIM at various times (CIM 5 d, SIM 12 h, SIM 24 h, and SIM 48 h) using the TRIzol method. Sequencing was performed by LC Sciences (Hangzhou, China) on the Illumina Hiseq platform (Illumina, San Diego, CA, USA). Raw data were processed to obtain TPM data (Dobin et al., 2013). Clustering of gene expression patterns across different sampling time points was conducted using the mfuzz method in ClusterGVis (https://github.com/junjunlab/ClusterGVis). The raw RNA-seq data have been deposited in the Genome Sequence Archive (Genomics, Proteomics & Bioinformatics 2021) at the National Genomics Data Center, China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences (https://ngdc.cncb.ac.cn/gsa) under the accession number GSA: CRA023698. The RNA-seq data can also be accessed using the online tool (http://xulinlab.cemps.ac.cn/).