design: Total RNA was extracted using Trizol and quantified using Qubit and Qsep400. Libraries were constructed and sequenced on an Illumina platform by Metware Biotechnology Co., Ltd. Raw reads were quality-controlled using fastp, aligned to the human reference genome using HISAT2, and gene expression levels were quantified using featureCounts (FPKM method). Differential expression analysis was performed using DESeq2 with Benjamini-Hochberg correction. Functional enrichment analyses were conducted using Reactome and KEGG databases. Weighted gene co-expression network analysis (WGCNA) was applied to identify co-expressed modules.