Differentiated organoids were dissociated to single cell, and sorted into two populations (Venus-positive and -negative) using an Influx Cell Sorter (BD Biosciences) at the Cambridge Institute of Medical Research (CIMR) Flow Cytometry Core Facility. DAPI and DRAQ5 staining, side scatter, forward scatter and pulse width gates were used to remove debris and cell clusters. RNA extraction and isolation were performed using RNAeasy Micro Plus kit (Qiagen). RNA concentration and quality (RIN 7-9) were assessed using RNA6000 Pico Kit and Bioanalyser 2000 (Agilent). The library was constructed with using the SMARTer Stranded Total RNA-Seq v2 Pico Input Mammalian kit (Takara Bio).