SRX5388166 GSM3611362 SRP186173 SRS4376609 GSM3611362 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611362 GSM3611362 GEO Accession GSM3611362 SRX5388167 GSM3611363 SRP186173 SRS4376610 GSM3611363 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611363 GSM3611363 GEO Accession GSM3611363 SRX5388168 GSM3611364 SRP186173 SRS4376611 GSM3611364 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611364 GSM3611364 GEO Accession GSM3611364 SRX5388169 GSM3611365 SRP186173 SRS4376612 GSM3611365 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611365 GSM3611365 GEO Accession GSM3611365 SRX5388170 GSM3611366 SRP186173 SRS4376613 GSM3611366 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611366 GSM3611366 GEO Accession GSM3611366 SRX5388171 GSM3611367 SRP186173 SRS4376614 GSM3611367 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611367 GSM3611367 GEO Accession GSM3611367 SRX5388172 GSM3611368 SRP186173 SRS4376615 GSM3611368 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611368 GSM3611368 GEO Accession GSM3611368 SRX5388173 GSM3611369 SRP186173 SRS4376616 GSM3611369 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611369 GSM3611369 GEO Accession GSM3611369 SRX5388174 GSM3611370 SRP186173 SRS4376617 GSM3611370 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611370 GSM3611370 GEO Accession GSM3611370 SRX5388175 GSM3611371 SRP186173 SRS4376618 GSM3611371 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611371 GSM3611371 GEO Accession GSM3611371 SRX5388176 GSM3611372 SRP186173 SRS4376619 GSM3611372 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611372 GSM3611372 GEO Accession GSM3611372 SRX5388177 GSM3611373 SRP186173 SRS4376620 GSM3611373 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611373 GSM3611373 GEO Accession GSM3611373 SRX5388178 GSM3611374 SRP186173 SRS4376621 GSM3611374 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611374 GSM3611374 GEO Accession GSM3611374 SRX5388179 GSM3611375 SRP186173 SRS4376622 GSM3611375 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611375 GSM3611375 GEO Accession GSM3611375 SRX5388180 GSM3611376 SRP186173 SRS4376623 GSM3611376 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611376 GSM3611376 GEO Accession GSM3611376 SRX5388181 GSM3611377 SRP186173 SRS4376624 GSM3611377 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611377 GSM3611377 GEO Accession GSM3611377 SRX5388182 GSM3611378 SRP186173 SRS4376625 GSM3611378 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611378 GSM3611378 GEO Accession GSM3611378 SRX5388183 GSM3611379 SRP186173 SRS4376626 GSM3611379 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611379 GSM3611379 GEO Accession GSM3611379 SRX5388184 GSM3611380 SRP186173 SRS4376627 GSM3611380 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611380 GSM3611380 GEO Accession GSM3611380 SRX5388185 GSM3611381 SRP186173 SRS4376628 GSM3611381 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611381 GSM3611381 GEO Accession GSM3611381 SRX5388186 GSM3611382 SRP186173 SRS4376629 GSM3611382 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611382 GSM3611382 GEO Accession GSM3611382 SRX5388187 GSM3611383 SRP186173 SRS4376630 GSM3611383 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611383 GSM3611383 GEO Accession GSM3611383 SRX5388188 GSM3611384 SRP186173 SRS4376631 GSM3611384 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611384 GSM3611384 GEO Accession GSM3611384 SRX5388189 GSM3611385 SRP186173 SRS4376632 GSM3611385 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611385 GSM3611385 GEO Accession GSM3611385 SRX5388190 GSM3611386 SRP186173 SRS4376633 GSM3611386 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611386 GSM3611386 GEO Accession GSM3611386 SRX5388191 GSM3611387 SRP186173 SRS4376634 GSM3611387 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611387 GSM3611387 GEO Accession GSM3611387 SRX5388192 GSM3611388 SRP186173 SRS4376635 GSM3611388 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611388 GSM3611388 GEO Accession GSM3611388 SRX5388193 GSM3611389 SRP186173 SRS4376636 GSM3611389 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611389 GSM3611389 GEO Accession GSM3611389 SRX5388194 GSM3611390 SRP186173 SRS4376637 GSM3611390 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611390 GSM3611390 GEO Accession GSM3611390 SRX5388195 GSM3611391 SRP186173 SRS4376638 GSM3611391 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611391 GSM3611391 GEO Accession GSM3611391 SRX5388196 GSM3611392 SRP186173 SRS4376639 GSM3611392 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611392 GSM3611392 GEO Accession GSM3611392 SRX5388197 GSM3611393 SRP186173 SRS4376640 GSM3611393 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611393 GSM3611393 GEO Accession GSM3611393 SRX5388198 GSM3611394 SRP186173 SRS4376641 GSM3611394 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611394 GSM3611394 GEO Accession GSM3611394 SRX5388199 GSM3611395 SRP186173 SRS4376642 GSM3611395 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611395 GSM3611395 GEO Accession GSM3611395 SRX5388200 GSM3611396 SRP186173 SRS4376643 GSM3611396 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611396 GSM3611396 GEO Accession GSM3611396 SRX5388201 GSM3611397 SRP186173 SRS4376644 GSM3611397 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611397 GSM3611397 GEO Accession GSM3611397 SRX5388202 GSM3611398 SRP186173 SRS4376645 GSM3611398 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611398 GSM3611398 GEO Accession GSM3611398 SRX5388203 GSM3611399 SRP186173 SRS4376646 GSM3611399 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611399 GSM3611399 GEO Accession GSM3611399 SRX5388204 GSM3611400 SRP186173 SRS4376647 GSM3611400 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611400 GSM3611400 GEO Accession GSM3611400 SRX5388205 GSM3611401 SRP186173 SRS4376648 GSM3611401 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611401 GSM3611401 GEO Accession GSM3611401 SRX5388206 GSM3611402 SRP186173 SRS4376649 GSM3611402 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611402 GSM3611402 GEO Accession GSM3611402 SRX5388207 GSM3611403 SRP186173 SRS4376650 GSM3611403 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611403 GSM3611403 GEO Accession GSM3611403 SRX5388208 GSM3611404 SRP186173 SRS4376651 GSM3611404 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611404 GSM3611404 GEO Accession GSM3611404 SRX5388209 GSM3611405 SRP186173 SRS4376652 GSM3611405 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611405 GSM3611405 GEO Accession GSM3611405 SRX5388210 GSM3611406 SRP186173 SRS4376653 GSM3611406 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611406 GSM3611406 GEO Accession GSM3611406 SRX5388211 GSM3611407 SRP186173 SRS4376654 GSM3611407 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611407 GSM3611407 GEO Accession GSM3611407 SRX5388212 GSM3611408 SRP186173 SRS4376655 GSM3611408 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611408 GSM3611408 GEO Accession GSM3611408 SRX5388213 GSM3611409 SRP186173 SRS4376656 GSM3611409 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611409 GSM3611409 GEO Accession GSM3611409 SRX5388214 GSM3611410 SRP186173 SRS4376657 GSM3611410 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611410 GSM3611410 GEO Accession GSM3611410 SRX5388215 GSM3611411 SRP186173 SRS4376658 GSM3611411 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611411 GSM3611411 GEO Accession GSM3611411 SRX5388216 GSM3611412 SRP186173 SRS4376659 GSM3611412 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611412 GSM3611412 GEO Accession GSM3611412 SRX5388217 GSM3611413 SRP186173 SRS4376660 GSM3611413 ATAC-seq GENOMIC other Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611413 GSM3611413 GEO Accession GSM3611413 SRX5388218 GSM3611414 SRP186173 SRS4376661 GSM3611414 RNA-Seq TRANSCRIPTOMIC cDNA Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611414 GSM3611414 GEO Accession GSM3611414 SRX5388219 GSM3611415 SRP186173 SRS4376662 GSM3611415 RNA-Seq TRANSCRIPTOMIC cDNA Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611415 GSM3611415 GEO Accession GSM3611415 SRX5388220 GSM3611416 SRP186173 SRS4376663 GSM3611416 RNA-Seq TRANSCRIPTOMIC cDNA Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611416 GSM3611416 GEO Accession GSM3611416 SRX5388221 GSM3611417 SRP186173 SRS4376664 GSM3611417 RNA-Seq TRANSCRIPTOMIC cDNA Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611417 GSM3611417 GEO Accession GSM3611417 SRX5388222 GSM3611418 SRP186173 SRS4376665 GSM3611418 RNA-Seq TRANSCRIPTOMIC cDNA Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611418 GSM3611418 GEO Accession GSM3611418 SRX5388223 GSM3611419 SRP186173 SRS4376666 GSM3611419 RNA-Seq TRANSCRIPTOMIC cDNA Nulcei were isolated from FACSorted cells and incubated in the tagmentation reaction. Tagmented-chromatin fragmented were PCR amplified and purified. The quality of ATAC-seq libraries was assessed by using 2100 Bioanalyzer (Agilent Technologies). ATAC-seq was performed with minor modifications to the published protocol (Buenrostro et al., 2013). 3000-5000 cells obtained by FACS were centrifuged at 800 x g for 4 min at 4 °C for each sample. Cells were resuspended in 5 µL of cold lysis buffer (10 mM Tris-HCl [pH 7.4], 10 nM NaCL, 3 nM MgCl2, 0,1% v/v Igepal CA-360 [Sigma-Aldrich] and incubated on ice for 5 min. After centrifugation of the cells at 500 x g for 10 min at 4 °C, the supernatant was discarded.  The tagmentation reaction was performed at 37 °C for 30 min using Nextera Sample Preparation kit (Illumina). After tagmentation, transposed DNA fragments were amplified using the following PCR condition: 1 cycle of 72°C for 5 min and 98°C for 30 s, followed by 12 to 14 cycles of 98°C for 10 s, 63°C for 30 s and 72°C for 1 min. Amplified libraries were purified using PCR purification MinElute kit (Qiagen), and the quality of of the purified library was assessed by using 2100 Bioanalyzer (Agilent Technologies) and a High Sensitivity DNA analysis Kit (Agilent Technologies) to confirm a period pattern in the size of the amplified DNA. After size selection by using Agencourt Ampure XP beads (Beckman Coulter, Cat#A633880) with a bead-to-sample ratio of 1.8:1, the libraries were sequenced as paired-end 50 bp by using the HiSeq 2500 platform (Illumina) according to the manufacturer's instructions. An input control was also generated by using 10 ng of C. robusta genomic DNA. Illumina HiSeq 2500 gds 303611419 GSM3611419 GEO Accession GSM3611419