The other portion of the samples was washed with nuclease-free water, surface moisture was absorbed with blotting paper, and the samples were wrapped in aluminum foil. These samples were rapidly frozen in liquid nitrogen and stored in a -80°C freezer. Using roots, stems, and leaves from two treatments, with three biological replicates as materials, each tissue was placed in a mortar and ground in liquid nitrogen. Total RNA was extracted using TRIzol reagent. Agarose gel electrophoresis (gel concentration: 1.2%; 0.5×TBE electrophoresis buffer; 150v for 15 minutes) was performed to assess the integrity of the RNA. The purity of the RNA was detected using a spectrophotometer (ThermoScientific, USA) under OD260/OD280 conditions. The VAHTS Universal V6 RNA-seq Library Prep Kit was used to construct the transcriptome libraries. Sequencing was performed on the Illumina Novaseq6000 sequencing platform, generating 150bp paired-end reads. The raw data (raw reads) in fastq format were processed using Trimmomatic[1] to remove reads containing poly-N and low-quality reads, resulting in clean reads. The clean reads were then assembled into expressed sequence tags (contigs), and the transcripts were de novo assembled using Trinity[2] software. Based on sequence similarity and length, the longest Unigene was selected, and the obtained Unigenes were annotated separately in functional databases such as the Non-Redundant Protein Database (NR), Clusters of Orthologous Groups (KOG), Gene Ontology (GO), Swiss-Prot Protein Database (Swiss-Prot), evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Protein Family Database (Pfam).