RNA was extracted using QIAsymphony PAXgene Blood RNA Kit following manufacturer's instructions. The library was constructed using NEBnext ultra II RNA library preparation kit with rRNA and globin depletion. We processed the raw RNA-seq data by nf-core rnaseq (version 3.12.0) (doi: 10.1038/s41587-020-0439-x), a community-curated best practice workflow. In brief, the sequence read quality was checked by fastqc. Then, the adapter and low-quality reads were trimmed by “Trim Galore!”. The reads were mapped to Macaca fascicularis reference genome 6.0 by HISAT2 (doi: 10.1038/s41587-019-0201-4), and annotated by Emsembl 109. The mapped counts were counted by Salmon (doi: 10.1038/s41587-019-0201-4).