BCL files generated from the Illumina HiSeq2500 platform were converted to fastq format and demultiplexed according to index primer sequences. The data was then imported into Partek Genomics Suite (Flow ver. 10.0.21.0328; Copyright 2009, Partek, St. Louis, MO, USA) for further processing. Read quality for each sample was assessed using FastQC, and high-quality reads were aligned to the mouse genome assembly (mm10, Ensemble transcripts release 99). The aligned reads were normalized to the median ratio and quantified for analysis with DeSeq2. Differential analysis was performed using DeSeq2, and pathway analysis was conducted with reference to the KEGG database.