Clustering analyses of DEGs in the four temperature treatments. Clusters were obtained using the hierarchical-means method based on the 5592 DEGs (the intersection set of the DEGs in the four treatments). (A) Heatmap of hierarchical clustering. Each column represents a treatment, and each row represents a unigene. Differences in expression were shown in different colors. Negative numbers indicate down-regulated and positive number means up-regulated. (B) Hierarchical clustering figure of differentially expressed genes. The “x” axis represents the treatments. The “y” axis represents the value of the relative expression level (log2 (HT, CS, FS/ NT)). In the subcluster, the gray line represents the value of the relative expression levels, the blue line represents the mean value of the relative expression level, and the red line is the reference. The subclusters in Figure B represent the models in Figure A from the top down.