SRX12776278 GSM5655849 SRP343162 SRS10721434 GSM5655849 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655849 GSM5655849 GEO Accession GSM5655849 SRX12776279 GSM5655850 SRP343162 SRS10721435 GSM5655850 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655850 GSM5655850 GEO Accession GSM5655850 SRX12776280 GSM5655851 SRP343162 SRS10721436 GSM5655851 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655851 GSM5655851 GEO Accession GSM5655851 SRX12776281 GSM5655852 SRP343162 SRS10721437 GSM5655852 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655852 GSM5655852 GEO Accession GSM5655852 SRX12776282 GSM5655853 SRP343162 SRS10721438 GSM5655853 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655853 GSM5655853 GEO Accession GSM5655853 SRX12776283 GSM5655854 SRP343162 SRS10721439 GSM5655854 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655854 GSM5655854 GEO Accession GSM5655854 SRX12776284 GSM5655855 SRP343162 SRS10721440 GSM5655855 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655855 GSM5655855 GEO Accession GSM5655855 SRX12776285 GSM5655856 SRP343162 SRS10721441 GSM5655856 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655856 GSM5655856 GEO Accession GSM5655856 SRX12776286 GSM5655857 SRP343162 SRS10721442 GSM5655857 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655857 GSM5655857 GEO Accession GSM5655857 SRX12776287 GSM5655858 SRP343162 SRS10721443 GSM5655858 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655858 GSM5655858 GEO Accession GSM5655858 SRX12776288 GSM5655859 SRP343162 SRS10721444 GSM5655859 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655859 GSM5655859 GEO Accession GSM5655859 SRX12776289 GSM5655860 SRP343162 SRS10721445 GSM5655860 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655860 GSM5655860 GEO Accession GSM5655860 SRX12776290 GSM5655861 SRP343162 SRS10721446 GSM5655861 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655861 GSM5655861 GEO Accession GSM5655861 SRX12776291 GSM5655862 SRP343162 SRS10721448 GSM5655862 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655862 GSM5655862 GEO Accession GSM5655862 SRX12776292 GSM5655863 SRP343162 SRS10721447 GSM5655863 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655863 GSM5655863 GEO Accession GSM5655863 SRX12776293 GSM5655864 SRP343162 SRS10721449 GSM5655864 RNA-Seq TRANSCRIPTOMIC cDNA The total RNA was harvested using Trizol-LS reagent. Total RNA was extracted from liver and brain using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's directions. RNA samples were digested using DNase I to remove potential genomic DNA and stored at -80℃. RNA was quantified and quality was assessed by NanoDrop Spectrophotometer 2000c (NanoDrop Technologies, Wilmington, DE, USA) and 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and all samples were standardised to 500 ng/μl. Three samples from per group at each time point were used for RNA library construction which OD ratios of 260/280 and 260/230 greater than 1.8 and RNA Integrity Number (RIN) greater than 8. RNA-seq library preparation and sequencing were carried out by Annoroad Gene Technology Co., Ltd (Beijing, China), forty five sequencing libraries were constructed using a TruSeq™ RNA Sample Preparation Kit (Illumina) according to the manufacturer's directions. The next step was KAPA quantifcation and dilution, after which the library was sequenced on an Illumina HiSeq.2500 with 125-bp paired-end reads. After removing adaptor sequences, ambiguous N nucleotides (those an N ratio <5%) and low quality sequences (reads with <50% bases of quality value), the remaining clean reads were assembled using Trinity sofware Illumina HiSeq 2500 gds 305655864 GSM5655864 GEO Accession GSM5655864