SRX18377218 GSM6755746_r1 SRP409752 PRJNA905359 SRS15863249 GSM6755746 GSM6755746 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377219 GSM6755747_r1 SRP409752 PRJNA905359 SRS15863250 GSM6755747 GSM6755747 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377220 GSM6755748_r1 SRP409752 PRJNA905359 SRS15863251 GSM6755748 GSM6755748 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377221 GSM6755749_r1 SRP409752 PRJNA905359 SRS15863252 GSM6755749 GSM6755749 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377222 GSM6755750_r1 SRP409752 PRJNA905359 SRS15863253 GSM6755750 GSM6755750 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377223 GSM6755751_r1 SRP409752 PRJNA905359 SRS15863254 GSM6755751 GSM6755751 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377224 GSM6755752_r1 SRP409752 PRJNA905359 SRS15863255 GSM6755752 GSM6755752 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377225 GSM6755753_r1 SRP409752 PRJNA905359 SRS15863256 GSM6755753 GSM6755753 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377226 GSM6755754_r1 SRP409752 PRJNA905359 SRS15863257 GSM6755754 GSM6755754 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377227 GSM6755755_r1 SRP409752 PRJNA905359 SRS15863258 GSM6755755 GSM6755755 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377228 GSM6755756_r1 SRP409752 PRJNA905359 SRS15863259 GSM6755756 GSM6755756 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000 SRX18377229 GSM6755757_r1 SRP409752 PRJNA905359 SRS15863260 GSM6755757 GSM6755757 RNA-Seq TRANSCRIPTOMIC cDNA Total cellular RNA was extracted from cultured cells using RNAiso Plus (Takara, 9108, Japan) according to the manufacturer's protocol. The RNA concentra-tion and purity were quantified using NanoDrop 2000 ultrafine spectrophotometer (Thermo Scientific, Massachusetts, America). And then the genome DNA was elimi-nated with the RNase-free DNase H. The RNA was then reverse-transcribed with 5× FastKing-RT SuperMix kit (TIANGEN, KR118–01, Beijing, China) according to the manufacturer's protocols. The reverse-transcribed cDNA was used as template for qPCR using Genious 2X SYBR Green Fast qPCR Mix (ABclonal, RK21204) according to the manufacturer's protocol. LightCycler 480 II real-time fluorescence quantitative PCR instrument was used for PCR reaction. The primer sequences used for amplifying the target fragments and GAPDH as internal control were listed in Table1. Gene's rela-tive expression levels were calculated by the 2-ΔΔCt (threshold cycle) method. Illumina HiSeq 2000