Low-quality sequences and adapters were removed by the function of trim_galore implanted in TrimGalore (version 0.5.0). Clean reads were further mapped to the GRCh38 genome using hisat2 (version 2.1.0). Gene expression was quantified at the gene level using featureCounts tool. Downstream analyses were performed in R (version 3.6.1).