Total RNA was extracted from plant tissues using the Monarch Total RNA Miniprep Kit (New England Biolabs, Massachusetts, USA; Cat. No. T2010S). The tissue was disrupted in 400 µL of 1× DNA/RNA Protection Reagent using a TacoPrep Bead Beater for 13 s, and the homogenized lysate was centrifuged at 16,000 × g for 2 min to remove plant debris. The supernatant was mixed with 400 µL of RNA lysis buffer, and genomic DNA was removed using a gDNA removal column by centrifugation at 16,000 × g for 30 s. An equal volume of 95% ethanol was added to the flow-through, and the mixture was loaded onto an RNA purification column. The column was washed with RNA wash buffer and centrifuged, and total RNA was eluted with 50 µL of nuclease-free water. The isolated RNA was stored at −80 °C until further analysis. RNA quantity and quality were assessed using a Nanodrop spectrophotometer, an Agilent 5400 Fragment Analyzer System, agarose gel electrophoresis, Qubit, real-time PCR, and a bioanalyzer. RNA sequencing libraries were constructed by purifying mRNA from total RNA using poly-T oligo-attached magnetic beads, followed by mRNA fragmentation, first- and second-strand cDNA synthesis, end repair, A-tailing, adapter ligation, size selection, amplification, and purification. Qualified libraries were sequenced on an Illumina NovaSeq X Plus platform using 150-bp paired-end sequencing. RNA-Seq