design: The RNA-seq library was prepared using the Illumina NovaSeq 6000 platform (Novogene Co. Ltd.) or GENEWIZ. Illumina adapter sequences and low-quality sequences (phred score < 20) at the 3 end were trimmed, and paired reads were removed if either of the two reads did not meet the minimum length (30 bp). Cleaned reads after trimming were aligned to a human genome (GRCH38) using HISAT2 (v.2.2.1). For each sample, we count the reads corresponding to each gene using htseq-count. Differential analysis among cell groups was performed using limma, a count-based method implemented in R with voom for normalization. BenjaminiHochberg adjusted p < 0.06 was the criteria used to identify significant differentially expressed genes.