SRX1655977 GSM2096686 SRP072214 SRS1356519 GSM2096686 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096686 GEO Accession GSM2096686 SRX1655978 GSM2096687 SRP072214 SRS1356518 GSM2096687 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096687 GEO Accession GSM2096687 SRX1655979 GSM2096688 SRP072214 SRS1356517 GSM2096688 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096688 GEO Accession GSM2096688 SRX1655980 GSM2096689 SRP072214 SRS1356516 GSM2096689 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096689 GEO Accession GSM2096689 SRX1655981 GSM2096690 SRP072214 SRS1356515 GSM2096690 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096690 GEO Accession GSM2096690 SRX1655982 GSM2096691 SRP072214 SRS1356514 GSM2096691 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096691 GEO Accession GSM2096691 SRX1655983 GSM2096692 SRP072214 SRS1356513 GSM2096692 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096692 GEO Accession GSM2096692 SRX1655984 GSM2096693 SRP072214 SRS1356512 GSM2096693 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096693 GEO Accession GSM2096693 SRX1655985 GSM2096694 SRP072214 SRS1356511 GSM2096694 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096694 GEO Accession GSM2096694 SRX1655986 GSM2096695 SRP072214 SRS1356510 GSM2096695 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096695 GEO Accession GSM2096695 SRX1655987 GSM2096696 SRP072214 SRS1356509 GSM2096696 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096696 GEO Accession GSM2096696 SRX1655988 GSM2096697 SRP072214 SRS1356319 GSM2096697 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096697 GEO Accession GSM2096697 SRX1655989 GSM2096698 SRP072214 SRS1356507 GSM2096698 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096698 GEO Accession GSM2096698 SRX1655990 GSM2096699 SRP072214 SRS1356508 GSM2096699 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096699 GEO Accession GSM2096699 SRX1655991 GSM2096700 SRP072214 SRS1356506 GSM2096700 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096700 GEO Accession GSM2096700 SRX1655992 GSM2096701 SRP072214 SRS1356505 GSM2096701 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096701 GEO Accession GSM2096701 SRX1655993 GSM2096702 SRP072214 SRS1356503 GSM2096702 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096702 GEO Accession GSM2096702 SRX1655994 GSM2096703 SRP072214 SRS1356504 GSM2096703 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096703 GEO Accession GSM2096703 SRX1655995 GSM2096704 SRP072214 SRS1356501 GSM2096704 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096704 GEO Accession GSM2096704 SRX1655996 GSM2096705 SRP072214 SRS1356502 GSM2096705 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096705 GEO Accession GSM2096705 SRX1655997 GSM2096706 SRP072214 SRS1356499 GSM2096706 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096706 GEO Accession GSM2096706 SRX1655998 GSM2096707 SRP072214 SRS1356500 GSM2096707 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096707 GEO Accession GSM2096707 SRX1655999 GSM2096708 SRP072214 SRS1356498 GSM2096708 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096708 GEO Accession GSM2096708 SRX1656000 GSM2096709 SRP072214 SRS1356497 GSM2096709 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096709 GEO Accession GSM2096709 SRX1656001 GSM2096710 SRP072214 SRS1356496 GSM2096710 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096710 GEO Accession GSM2096710 SRX1656002 GSM2096711 SRP072214 SRS1356495 GSM2096711 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096711 GEO Accession GSM2096711 SRX1656003 GSM2096712 SRP072214 SRS1356493 GSM2096712 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096712 GEO Accession GSM2096712 SRX1656004 GSM2096713 SRP072214 SRS1356494 GSM2096713 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096713 GEO Accession GSM2096713 SRX1656005 GSM2096714 SRP072214 SRS1356491 GSM2096714 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096714 GEO Accession GSM2096714 SRX1656006 GSM2096715 SRP072214 SRS1356492 GSM2096715 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096715 GEO Accession GSM2096715 SRX1656007 GSM2096716 SRP072214 SRS1356490 GSM2096716 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096716 GEO Accession GSM2096716 SRX1656008 GSM2096717 SRP072214 SRS1356489 GSM2096717 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096717 GEO Accession GSM2096717 SRX1656009 GSM2096718 SRP072214 SRS1356488 GSM2096718 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096718 GEO Accession GSM2096718 SRX1656010 GSM2096719 SRP072214 SRS1356486 GSM2096719 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096719 GEO Accession GSM2096719 SRX1656011 GSM2096720 SRP072214 SRS1356487 GSM2096720 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096720 GEO Accession GSM2096720 SRX1656012 GSM2096721 SRP072214 SRS1356485 GSM2096721 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096721 GEO Accession GSM2096721 SRX1656013 GSM2096722 SRP072214 SRS1356484 GSM2096722 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096722 GEO Accession GSM2096722 SRX1656014 GSM2096723 SRP072214 SRS1356481 GSM2096723 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096723 GEO Accession GSM2096723 SRX1656015 GSM2096724 SRP072214 SRS1356482 GSM2096724 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096724 GEO Accession GSM2096724 SRX1656016 GSM2096725 SRP072214 SRS1356483 GSM2096725 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096725 GEO Accession GSM2096725 SRX1656017 GSM2096726 SRP072214 SRS1356480 GSM2096726 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096726 GEO Accession GSM2096726 SRX1656018 GSM2096727 SRP072214 SRS1356479 GSM2096727 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096727 GEO Accession GSM2096727 SRX1656019 GSM2096728 SRP072214 SRS1356477 GSM2096728 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096728 GEO Accession GSM2096728 SRX1656020 GSM2096729 SRP072214 SRS1356478 GSM2096729 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096729 GEO Accession GSM2096729 SRX1656021 GSM2096730 SRP072214 SRS1356476 GSM2096730 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096730 GEO Accession GSM2096730 SRX1656022 GSM2096731 SRP072214 SRS1356474 GSM2096731 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096731 GEO Accession GSM2096731 SRX1656023 GSM2096732 SRP072214 SRS1356475 GSM2096732 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096732 GEO Accession GSM2096732 SRX1656024 GSM2096733 SRP072214 SRS1356473 GSM2096733 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096733 GEO Accession GSM2096733 SRX1656025 GSM2096734 SRP072214 SRS1356472 GSM2096734 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096734 GEO Accession GSM2096734 SRX1656026 GSM2096735 SRP072214 SRS1356471 GSM2096735 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096735 GEO Accession GSM2096735 SRX1656027 GSM2096736 SRP072214 SRS1356470 GSM2096736 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096736 GEO Accession GSM2096736 SRX1656028 GSM2096737 SRP072214 SRS1356465 GSM2096737 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096737 GEO Accession GSM2096737 SRX1656029 GSM2096738 SRP072214 SRS1356462 GSM2096738 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096738 GEO Accession GSM2096738 SRX1656030 GSM2096739 SRP072214 SRS1356461 GSM2096739 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096739 GEO Accession GSM2096739 SRX1656031 GSM2096740 SRP072214 SRS1356457 GSM2096740 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096740 GEO Accession GSM2096740 SRX1656032 GSM2096741 SRP072214 SRS1356459 GSM2096741 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096741 GEO Accession GSM2096741 SRX1656033 GSM2096742 SRP072214 SRS1356458 GSM2096742 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096742 GEO Accession GSM2096742 SRX1656034 GSM2096743 SRP072214 SRS1356455 GSM2096743 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096743 GEO Accession GSM2096743 SRX1656035 GSM2096744 SRP072214 SRS1356454 GSM2096744 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096744 GEO Accession GSM2096744 SRX1656036 GSM2096745 SRP072214 SRS1356456 GSM2096745 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096745 GEO Accession GSM2096745 SRX1656038 GSM2096746 SRP072214 SRS1356449 GSM2096746 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096746 GEO Accession GSM2096746 SRX1656039 GSM2096747 SRP072214 SRS1356448 GSM2096747 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096747 GEO Accession GSM2096747 SRX1656040 GSM2096748 SRP072214 SRS1356447 GSM2096748 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096748 GEO Accession GSM2096748 SRX1656041 GSM2096749 SRP072214 SRS1356446 GSM2096749 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096749 GEO Accession GSM2096749 SRX1656042 GSM2096750 SRP072214 SRS1356443 GSM2096750 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096750 GEO Accession GSM2096750 SRX1656043 GSM2096751 SRP072214 SRS1356445 GSM2096751 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096751 GEO Accession GSM2096751 SRX1656044 GSM2096752 SRP072214 SRS1356444 GSM2096752 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096752 GEO Accession GSM2096752 SRX1656045 GSM2096753 SRP072214 SRS1356442 GSM2096753 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096753 GEO Accession GSM2096753 SRX1656046 GSM2096754 SRP072214 SRS1356441 GSM2096754 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096754 GEO Accession GSM2096754 SRX1656047 GSM2096755 SRP072214 SRS1356440 GSM2096755 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096755 GEO Accession GSM2096755 SRX1656048 GSM2096756 SRP072214 SRS1356439 GSM2096756 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096756 GEO Accession GSM2096756 SRX1656049 GSM2096757 SRP072214 SRS1356438 GSM2096757 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096757 GEO Accession GSM2096757 SRX1656050 GSM2096758 SRP072214 SRS1356437 GSM2096758 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096758 GEO Accession GSM2096758 SRX1656051 GSM2096759 SRP072214 SRS1356436 GSM2096759 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096759 GEO Accession GSM2096759 SRX1656052 GSM2096760 SRP072214 SRS1356435 GSM2096760 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096760 GEO Accession GSM2096760 SRX1656053 GSM2096761 SRP072214 SRS1356434 GSM2096761 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096761 GEO Accession GSM2096761 SRX1656054 GSM2096762 SRP072214 SRS1356433 GSM2096762 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096762 GEO Accession GSM2096762 SRX1656055 GSM2096763 SRP072214 SRS1356432 GSM2096763 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096763 GEO Accession GSM2096763 SRX1656056 GSM2096764 SRP072214 SRS1356431 GSM2096764 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096764 GEO Accession GSM2096764 SRX1656057 GSM2096765 SRP072214 SRS1356430 GSM2096765 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096765 GEO Accession GSM2096765 SRX1656058 GSM2096766 SRP072214 SRS1356429 GSM2096766 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096766 GEO Accession GSM2096766 SRX1656059 GSM2096767 SRP072214 SRS1356428 GSM2096767 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096767 GEO Accession GSM2096767 SRX1656060 GSM2096768 SRP072214 SRS1356427 GSM2096768 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096768 GEO Accession GSM2096768 SRX1656061 GSM2096769 SRP072214 SRS1356426 GSM2096769 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096769 GEO Accession GSM2096769 SRX1656062 GSM2096770 SRP072214 SRS1356425 GSM2096770 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096770 GEO Accession GSM2096770 SRX1656063 GSM2096771 SRP072214 SRS1356423 GSM2096771 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096771 GEO Accession GSM2096771 SRX1656064 GSM2096772 SRP072214 SRS1356424 GSM2096772 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096772 GEO Accession GSM2096772 SRX1656065 GSM2096773 SRP072214 SRS1356422 GSM2096773 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096773 GEO Accession GSM2096773 SRX1656066 GSM2096774 SRP072214 SRS1356421 GSM2096774 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096774 GEO Accession GSM2096774 SRX1656067 GSM2096775 SRP072214 SRS1356419 GSM2096775 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096775 GEO Accession GSM2096775 SRX1656068 GSM2096776 SRP072214 SRS1356420 GSM2096776 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096776 GEO Accession GSM2096776 SRX1656069 GSM2096777 SRP072214 SRS1356417 GSM2096777 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096777 GEO Accession GSM2096777 SRX1656070 GSM2096778 SRP072214 SRS1356418 GSM2096778 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096778 GEO Accession GSM2096778 SRX1656071 GSM2096779 SRP072214 SRS1356416 GSM2096779 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096779 GEO Accession GSM2096779 SRX1656072 GSM2096780 SRP072214 SRS1356415 GSM2096780 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096780 GEO Accession GSM2096780 SRX1656073 GSM2096781 SRP072214 SRS1356414 GSM2096781 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096781 GEO Accession GSM2096781 SRX1656074 GSM2096782 SRP072214 SRS1356413 GSM2096782 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096782 GEO Accession GSM2096782 SRX1656075 GSM2096783 SRP072214 SRS1356412 GSM2096783 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096783 GEO Accession GSM2096783 SRX1656076 GSM2096784 SRP072214 SRS1356411 GSM2096784 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096784 GEO Accession GSM2096784 SRX1656077 GSM2096785 SRP072214 SRS1356410 GSM2096785 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096785 GEO Accession GSM2096785 SRX1656078 GSM2096786 SRP072214 SRS1356409 GSM2096786 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096786 GEO Accession GSM2096786 SRX1656079 GSM2096787 SRP072214 SRS1356408 GSM2096787 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096787 GEO Accession GSM2096787 SRX1656080 GSM2096788 SRP072214 SRS1356407 GSM2096788 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096788 GEO Accession GSM2096788 SRX1656081 GSM2096789 SRP072214 SRS1356406 GSM2096789 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096789 GEO Accession GSM2096789 SRX1656082 GSM2096790 SRP072214 SRS1356405 GSM2096790 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096790 GEO Accession GSM2096790 SRX1656083 GSM2096791 SRP072214 SRS1356404 GSM2096791 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096791 GEO Accession GSM2096791 SRX1656084 GSM2096792 SRP072214 SRS1356403 GSM2096792 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096792 GEO Accession GSM2096792 SRX1656085 GSM2096793 SRP072214 SRS1356402 GSM2096793 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096793 GEO Accession GSM2096793 SRX1656086 GSM2096794 SRP072214 SRS1356401 GSM2096794 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096794 GEO Accession GSM2096794 SRX1656087 GSM2096795 SRP072214 SRS1356400 GSM2096795 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096795 GEO Accession GSM2096795 SRX1656088 GSM2096796 SRP072214 SRS1356399 GSM2096796 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096796 GEO Accession GSM2096796 SRX1656089 GSM2096797 SRP072214 SRS1356397 GSM2096797 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096797 GEO Accession GSM2096797 SRX1656090 GSM2096798 SRP072214 SRS1356396 GSM2096798 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096798 GEO Accession GSM2096798 SRX1656091 GSM2096799 SRP072214 SRS1356398 GSM2096799 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096799 GEO Accession GSM2096799 SRX1656092 GSM2096800 SRP072214 SRS1356395 GSM2096800 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096800 GEO Accession GSM2096800 SRX1656093 GSM2096801 SRP072214 SRS1356393 GSM2096801 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096801 GEO Accession GSM2096801 SRX1656094 GSM2096802 SRP072214 SRS1356394 GSM2096802 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096802 GEO Accession GSM2096802 SRX1656095 GSM2096803 SRP072214 SRS1356391 GSM2096803 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096803 GEO Accession GSM2096803 SRX1656096 GSM2096804 SRP072214 SRS1356392 GSM2096804 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096804 GEO Accession GSM2096804 SRX1656097 GSM2096805 SRP072214 SRS1356389 GSM2096805 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096805 GEO Accession GSM2096805 SRX1656098 GSM2096806 SRP072214 SRS1356390 GSM2096806 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096806 GEO Accession GSM2096806 SRX1656099 GSM2096807 SRP072214 SRS1356388 GSM2096807 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096807 GEO Accession GSM2096807 SRX1656100 GSM2096808 SRP072214 SRS1356387 GSM2096808 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096808 GEO Accession GSM2096808 SRX1656101 GSM2096809 SRP072214 SRS1356386 GSM2096809 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096809 GEO Accession GSM2096809 SRX1656102 GSM2096810 SRP072214 SRS1356385 GSM2096810 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096810 GEO Accession GSM2096810 SRX1656103 GSM2096811 SRP072214 SRS1356384 GSM2096811 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096811 GEO Accession GSM2096811 SRX1656104 GSM2096812 SRP072214 SRS1356383 GSM2096812 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096812 GEO Accession GSM2096812 SRX1656105 GSM2096813 SRP072214 SRS1356382 GSM2096813 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096813 GEO Accession GSM2096813 SRX1656106 GSM2096814 SRP072214 SRS1356381 GSM2096814 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096814 GEO Accession GSM2096814 SRX1656107 GSM2096815 SRP072214 SRS1356380 GSM2096815 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096815 GEO Accession GSM2096815 SRX1656108 GSM2096816 SRP072214 SRS1356379 GSM2096816 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096816 GEO Accession GSM2096816 SRX1656109 GSM2096817 SRP072214 SRS1356378 GSM2096817 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096817 GEO Accession GSM2096817 SRX1656110 GSM2096818 SRP072214 SRS1356377 GSM2096818 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096818 GEO Accession GSM2096818 SRX1656111 GSM2096819 SRP072214 SRS1356376 GSM2096819 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096819 GEO Accession GSM2096819 SRX1656112 GSM2096820 SRP072214 SRS1356375 GSM2096820 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096820 GEO Accession GSM2096820 SRX1656113 GSM2096821 SRP072214 SRS1356374 GSM2096821 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096821 GEO Accession GSM2096821 SRX1656114 GSM2096822 SRP072214 SRS1356373 GSM2096822 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096822 GEO Accession GSM2096822 SRX1656115 GSM2096823 SRP072214 SRS1356372 GSM2096823 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096823 GEO Accession GSM2096823 SRX1656116 GSM2096824 SRP072214 SRS1356371 GSM2096824 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096824 GEO Accession GSM2096824 SRX1656117 GSM2096825 SRP072214 SRS1356370 GSM2096825 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096825 GEO Accession GSM2096825 SRX1656118 GSM2096826 SRP072214 SRS1356369 GSM2096826 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096826 GEO Accession GSM2096826 SRX1656119 GSM2096827 SRP072214 SRS1356368 GSM2096827 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096827 GEO Accession GSM2096827 SRX1656120 GSM2096828 SRP072214 SRS1356367 GSM2096828 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096828 GEO Accession GSM2096828 SRX1656121 GSM2096829 SRP072214 SRS1356366 GSM2096829 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096829 GEO Accession GSM2096829 SRX1656122 GSM2096830 SRP072214 SRS1356365 GSM2096830 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096830 GEO Accession GSM2096830 SRX1656123 GSM2096831 SRP072214 SRS1356364 GSM2096831 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096831 GEO Accession GSM2096831 SRX1656124 GSM2096832 SRP072214 SRS1356363 GSM2096832 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096832 GEO Accession GSM2096832 SRX1656125 GSM2096833 SRP072214 SRS1356362 GSM2096833 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096833 GEO Accession GSM2096833 SRX1656126 GSM2096834 SRP072214 SRS1356361 GSM2096834 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096834 GEO Accession GSM2096834 SRX1656127 GSM2096835 SRP072214 SRS1356360 GSM2096835 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096835 GEO Accession GSM2096835 SRX1656128 GSM2096836 SRP072214 SRS1356359 GSM2096836 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096836 GEO Accession GSM2096836 SRX1656129 GSM2096837 SRP072214 SRS1356358 GSM2096837 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096837 GEO Accession GSM2096837 SRX1656130 GSM2096838 SRP072214 SRS1356357 GSM2096838 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096838 GEO Accession GSM2096838 SRX1656131 GSM2096839 SRP072214 SRS1356356 GSM2096839 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096839 GEO Accession GSM2096839 SRX1656132 GSM2096840 SRP072214 SRS1356355 GSM2096840 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096840 GEO Accession GSM2096840 SRX1656133 GSM2096841 SRP072214 SRS1356354 GSM2096841 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096841 GEO Accession GSM2096841 SRX1656134 GSM2096842 SRP072214 SRS1356353 GSM2096842 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096842 GEO Accession GSM2096842 SRX1656135 GSM2096843 SRP072214 SRS1356352 GSM2096843 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096843 GEO Accession GSM2096843 SRX1656136 GSM2096844 SRP072214 SRS1356351 GSM2096844 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096844 GEO Accession GSM2096844 SRX1656137 GSM2096845 SRP072214 SRS1356350 GSM2096845 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096845 GEO Accession GSM2096845 SRX1656138 GSM2096846 SRP072214 SRS1356349 GSM2096846 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096846 GEO Accession GSM2096846 SRX1656139 GSM2096847 SRP072214 SRS1356348 GSM2096847 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096847 GEO Accession GSM2096847 SRX1656140 GSM2096848 SRP072214 SRS1356347 GSM2096848 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096848 GEO Accession GSM2096848 SRX1656141 GSM2096849 SRP072214 SRS1356346 GSM2096849 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096849 GEO Accession GSM2096849 SRX1656142 GSM2096850 SRP072214 SRS1356345 GSM2096850 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096850 GEO Accession GSM2096850 SRX1656143 GSM2096851 SRP072214 SRS1356320 GSM2096851 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096851 GEO Accession GSM2096851 SRX1656144 GSM2096852 SRP072214 SRS1356344 GSM2096852 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096852 GEO Accession GSM2096852 SRX1656145 GSM2096853 SRP072214 SRS1356343 GSM2096853 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096853 GEO Accession GSM2096853 SRX1656146 GSM2096854 SRP072214 SRS1356342 GSM2096854 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096854 GEO Accession GSM2096854 SRX1656147 GSM2096855 SRP072214 SRS1356341 GSM2096855 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096855 GEO Accession GSM2096855 SRX1656148 GSM2096856 SRP072214 SRS1356340 GSM2096856 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096856 GEO Accession GSM2096856 SRX1656149 GSM2096857 SRP072214 SRS1356339 GSM2096857 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096857 GEO Accession GSM2096857 SRX1656150 GSM2096858 SRP072214 SRS1356338 GSM2096858 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096858 GEO Accession GSM2096858 SRX1656151 GSM2096859 SRP072214 SRS1356337 GSM2096859 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096859 GEO Accession GSM2096859 SRX1656152 GSM2096860 SRP072214 SRS1356336 GSM2096860 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096860 GEO Accession GSM2096860 SRX1656153 GSM2096861 SRP072214 SRS1356335 GSM2096861 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096861 GEO Accession GSM2096861 SRX1656154 GSM2096862 SRP072214 SRS1356334 GSM2096862 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096862 GEO Accession GSM2096862 SRX1656155 GSM2096863 SRP072214 SRS1356333 GSM2096863 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096863 GEO Accession GSM2096863 SRX1656156 GSM2096864 SRP072214 SRS1356332 GSM2096864 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096864 GEO Accession GSM2096864 SRX1656157 GSM2096865 SRP072214 SRS1356331 GSM2096865 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096865 GEO Accession GSM2096865 SRX1656158 GSM2096866 SRP072214 SRS1356323 GSM2096866 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096866 GEO Accession GSM2096866 SRX1656159 GSM2096867 SRP072214 SRS1356321 GSM2096867 RNA-Seq TRANSCRIPTOMIC cDNA Isolated mouse islets were dissociated into single cells by enzymatic digestion using 0.05% Trypsin (25300-120; Life Technologies). 7-aminoactinomycin D (7-AAD, eBiosciences) was used at 1:500 dilution as Live/Dead stain to exclude dead cells. Sort gates were adjusted with reference to negative controls (wild-type islets without YFP labeling). We collected cells from YFP+, lineage-traced population and also YFPNeg, non-labeled population. Cells were sorted on a special order 5-laser FACS Aria II (BD Biosciences) using a 100 µm nozzle at a flow rate of 1 following doublet removal. Sorted single-cells were collected directly into 96-well plates (Bio-Rad cat #: HSP9601) containing 4 µL of lysis buffer with dNTPs37 for downstream single-cell RNA-Seq assays. Single-cells were collected in lysis buffer in 96-well plates, followed by reverse transcription with template-switch using an LNA-modified template switch oligo to generate cDNA. After 21 cycles of pre-amplification, DNA was purified and analyzed on an automated Fragment Analyzer (Advanced Analytical). Each cell’s cDNA fragment profile was individually inspected and only wells with successful amplification products (concentration higher than 0.06 ng/ul) and with no detectable RNA degradation were selected for final library preparation. Tagmentation assays and barcoded sequencing libraries were prepared using Nextera XT kit (FC-131-1024; Illumina) according to the manufacturer’s instructions. Barcoded libraries were pooled and subjected to 75 bp paired-end sequencing on the Illumina NextSeq instrument. NextSeq 500 gds 302096867 GEO Accession GSM2096867