SRX20276743 GSM7316144_r1 SRP436835 PRJNA971103 SRS17603999 GSM7316144 GSM7316144 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276744 GSM7316161_r1 SRP436835 PRJNA971103 SRS17604000 GSM7316161 GSM7316161 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276745 GSM7316162_r1 SRP436835 PRJNA971103 SRS17603997 GSM7316162 GSM7316162 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276746 GSM7316163_r1 SRP436835 PRJNA971103 SRS17603998 GSM7316163 GSM7316163 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276747 GSM7316164_r1 SRP436835 PRJNA971103 SRS17604001 GSM7316164 GSM7316164 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276748 GSM7316165_r1 SRP436835 PRJNA971103 SRS17604004 GSM7316165 GSM7316165 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276749 GSM7316166_r1 SRP436835 PRJNA971103 SRS17604007 GSM7316166 GSM7316166 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276750 GSM7316167_r1 SRP436835 PRJNA971103 SRS17604002 GSM7316167 GSM7316167 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276751 GSM7316168_r1 SRP436835 PRJNA971103 SRS17604003 GSM7316168 GSM7316168 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276752 GSM7316185_r1 SRP436835 PRJNA971103 SRS17604005 GSM7316185 GSM7316185 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276753 GSM7316186_r1 SRP436835 PRJNA971103 SRS17604008 GSM7316186 GSM7316186 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276754 GSM7316187_r1 SRP436835 PRJNA971103 SRS17604006 GSM7316187 GSM7316187 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276755 GSM7316188_r1 SRP436835 PRJNA971103 SRS17604009 GSM7316188 GSM7316188 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276756 GSM7316189_r1 SRP436835 PRJNA971103 SRS17604011 GSM7316189 GSM7316189 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276757 GSM7316190_r1 SRP436835 PRJNA971103 SRS17604012 GSM7316190 GSM7316190 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276758 GSM7316191_r1 SRP436835 PRJNA971103 SRS17604010 GSM7316191 GSM7316191 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276759 GSM7316192_r1 SRP436835 PRJNA971103 SRS17604013 GSM7316192 GSM7316192 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276760 GSM7316209_r1 SRP436835 PRJNA971103 SRS17604015 GSM7316209 GSM7316209 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276761 GSM7316210_r1 SRP436835 PRJNA971103 SRS17604016 GSM7316210 GSM7316210 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276762 GSM7316211_r1 SRP436835 PRJNA971103 SRS17604014 GSM7316211 GSM7316211 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276763 GSM7316212_r1 SRP436835 PRJNA971103 SRS17604017 GSM7316212 GSM7316212 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276764 GSM7316213_r1 SRP436835 PRJNA971103 SRS17604019 GSM7316213 GSM7316213 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276765 GSM7316214_r1 SRP436835 PRJNA971103 SRS17604018 GSM7316214 GSM7316214 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276766 GSM7316215_r1 SRP436835 PRJNA971103 SRS17604020 GSM7316215 GSM7316215 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276767 GSM7316216_r1 SRP436835 PRJNA971103 SRS17604021 GSM7316216 GSM7316216 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276768 GSM7316145_r1 SRP436835 PRJNA971103 SRS17604023 GSM7316145 GSM7316145 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276769 GSM7316146_r1 SRP436835 PRJNA971103 SRS17604022 GSM7316146 GSM7316146 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276770 GSM7316147_r1 SRP436835 PRJNA971103 SRS17604024 GSM7316147 GSM7316147 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276771 GSM7316148_r1 SRP436835 PRJNA971103 SRS17604025 GSM7316148 GSM7316148 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276772 GSM7316149_r1 SRP436835 PRJNA971103 SRS17604028 GSM7316149 GSM7316149 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276773 GSM7316150_r1 SRP436835 PRJNA971103 SRS17604026 GSM7316150 GSM7316150 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276774 GSM7316151_r1 SRP436835 PRJNA971103 SRS17604029 GSM7316151 GSM7316151 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276775 GSM7316152_r1 SRP436835 PRJNA971103 SRS17604027 GSM7316152 GSM7316152 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276776 GSM7316193_r1 SRP436835 PRJNA971103 SRS17604030 GSM7316193 GSM7316193 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276777 GSM7316194_r1 SRP436835 PRJNA971103 SRS17604033 GSM7316194 GSM7316194 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276778 GSM7316195_r1 SRP436835 PRJNA971103 SRS17604034 GSM7316195 GSM7316195 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276779 GSM7316196_r1 SRP436835 PRJNA971103 SRS17604031 GSM7316196 GSM7316196 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276780 GSM7316197_r1 SRP436835 PRJNA971103 SRS17604032 GSM7316197 GSM7316197 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276781 GSM7316198_r1 SRP436835 PRJNA971103 SRS17604035 GSM7316198 GSM7316198 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276782 GSM7316199_r1 SRP436835 PRJNA971103 SRS17604036 GSM7316199 GSM7316199 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276783 GSM7316200_r1 SRP436835 PRJNA971103 SRS17604040 GSM7316200 GSM7316200 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276784 GSM7316225_r1 SRP436835 PRJNA971103 SRS17604041 GSM7316225 GSM7316225 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276785 GSM7316226_r1 SRP436835 PRJNA971103 SRS17604037 GSM7316226 GSM7316226 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276786 GSM7316153_r1 SRP436835 PRJNA971103 SRS17604038 GSM7316153 GSM7316153 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276787 GSM7316154_r1 SRP436835 PRJNA971103 SRS17604039 GSM7316154 GSM7316154 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276788 GSM7316155_r1 SRP436835 PRJNA971103 SRS17604042 GSM7316155 GSM7316155 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276789 GSM7316156_r1 SRP436835 PRJNA971103 SRS17604043 GSM7316156 GSM7316156 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276790 GSM7316157_r1 SRP436835 PRJNA971103 SRS17604047 GSM7316157 GSM7316157 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276791 GSM7316158_r1 SRP436835 PRJNA971103 SRS17604046 GSM7316158 GSM7316158 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276792 GSM7316159_r1 SRP436835 PRJNA971103 SRS17604044 GSM7316159 GSM7316159 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276793 GSM7316160_r1 SRP436835 PRJNA971103 SRS17604045 GSM7316160 GSM7316160 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276794 GSM7316177_r1 SRP436835 PRJNA971103 SRS17604049 GSM7316177 GSM7316177 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276795 GSM7316178_r1 SRP436835 PRJNA971103 SRS17604048 GSM7316178 GSM7316178 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276796 GSM7316179_r1 SRP436835 PRJNA971103 SRS17604052 GSM7316179 GSM7316179 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276797 GSM7316180_r1 SRP436835 PRJNA971103 SRS17604053 GSM7316180 GSM7316180 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276798 GSM7316181_r1 SRP436835 PRJNA971103 SRS17604050 GSM7316181 GSM7316181 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276799 GSM7316182_r1 SRP436835 PRJNA971103 SRS17604051 GSM7316182 GSM7316182 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276800 GSM7316183_r1 SRP436835 PRJNA971103 SRS17604054 GSM7316183 GSM7316183 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276801 GSM7316184_r1 SRP436835 PRJNA971103 SRS17604056 GSM7316184 GSM7316184 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276802 GSM7316201_r1 SRP436835 PRJNA971103 SRS17604064 GSM7316201 GSM7316201 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276803 GSM7316202_r1 SRP436835 PRJNA971103 SRS17604055 GSM7316202 GSM7316202 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276804 GSM7316203_r1 SRP436835 PRJNA971103 SRS17604065 GSM7316203 GSM7316203 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276805 GSM7316204_r1 SRP436835 PRJNA971103 SRS17604057 GSM7316204 GSM7316204 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276806 GSM7316205_r1 SRP436835 PRJNA971103 SRS17604058 GSM7316205 GSM7316205 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276807 GSM7316206_r1 SRP436835 PRJNA971103 SRS17604059 GSM7316206 GSM7316206 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276808 GSM7316207_r1 SRP436835 PRJNA971103 SRS17604061 GSM7316207 GSM7316207 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276809 GSM7316208_r1 SRP436835 PRJNA971103 SRS17604060 GSM7316208 GSM7316208 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276810 GSM7316169_r1 SRP436835 PRJNA971103 SRS17604062 GSM7316169 GSM7316169 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276811 GSM7316170_r1 SRP436835 PRJNA971103 SRS17604063 GSM7316170 GSM7316170 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276812 GSM7316171_r1 SRP436835 PRJNA971103 SRS17604066 GSM7316171 GSM7316171 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276813 GSM7316172_r1 SRP436835 PRJNA971103 SRS17604067 GSM7316172 GSM7316172 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276814 GSM7316173_r1 SRP436835 PRJNA971103 SRS17604072 GSM7316173 GSM7316173 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276815 GSM7316174_r1 SRP436835 PRJNA971103 SRS17604071 GSM7316174 GSM7316174 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276816 GSM7316175_r1 SRP436835 PRJNA971103 SRS17604068 GSM7316175 GSM7316175 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276817 GSM7316176_r1 SRP436835 PRJNA971103 SRS17604069 GSM7316176 GSM7316176 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276818 GSM7316217_r1 SRP436835 PRJNA971103 SRS17604070 GSM7316217 GSM7316217 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276819 GSM7316218_r1 SRP436835 PRJNA971103 SRS17604073 GSM7316218 GSM7316218 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276820 GSM7316219_r1 SRP436835 PRJNA971103 SRS17604074 GSM7316219 GSM7316219 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276821 GSM7316220_r1 SRP436835 PRJNA971103 SRS17604078 GSM7316220 GSM7316220 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276822 GSM7316221_r1 SRP436835 PRJNA971103 SRS17604077 GSM7316221 GSM7316221 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276823 GSM7316222_r1 SRP436835 PRJNA971103 SRS17604076 GSM7316222 GSM7316222 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276824 GSM7316223_r1 SRP436835 PRJNA971103 SRS17604075 GSM7316223 GSM7316223 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000 SRX20276825 GSM7316224_r1 SRP436835 PRJNA971103 SRS17604079 GSM7316224 GSM7316224 RNA-Seq TRANSCRIPTOMIC cDNA For next generation sequencing the RNA-Bee total-RNA isolation kit (Campro Scientific, Veenendaal, the Netherlands) was used to isolate total RNA from individual liver samples and RNA integrity was determined with the RNA 6000 Nano Lab-on-a-Chip kit and bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands). Total RNA's were processed into tagged random sequence libraries. Next Generation Sequencing was performed at Genome Scan (Leiden, the Netherlands). For this, strand-specific cDNA libraries were generated from total-RNA using oligo-dT magnetic beads, mRNA fragmentation, NEBNext Ultra Directional RNA Library Prep Kit from Illumina, NEB #E7760S/L resulting in 300–500 bp amplified libraries/sample. Libraries were multiplexed, clustered, and sequenced on a NovaSeq6000 v1.5 system (Illumina, New England Biolabs, Inc. Ipswich, MA, USA), with a paired-end, 2x150 bp sequencing protocol, 20 million (min-max=11-38 million) reads per sample and indexing. Clustering and DNA sequencing using the NovaSeq6000 v1.5 was performed according to manufacturer's protocols. A concentration of 1.1 nM of DNA was used yielding paired end reads (2x 150bp). NovaSeq control software NCS v1.7 was used. Image analysis, base calling, and quality check was performed with the Illumina data analysis pipeline RTA3.4.4 and Bcl2fastq v2.20. Illumina NovaSeq 6000