Transcriptome data Analysis
Transcriptome data were processed using Fastp (v0.23.2) for quality control, Hisat2 (v2.1.0) for accurate transcript alignment, and featureCounts (v2.0.6) for efficient quantification of transcript abundance. These raw counts were subsequently converted into Transcripts Per Million (TPM) values for normalization and comparative analysis. A weighted gene coexpression network analysis (WGCNA) was performed using the R (v4.3.3) package WGCNA (v1.72.5), following published procedure (Langfelder & Horvath, 2008).