SRX18487385 GSM6784420_r1 SRP411259 PRJNA908254 SRS15960300 GSM6784420 GSM6784420 RNA-Seq TRANSCRIPTOMIC cDNA Once organoids reached confluence (~3 days after first passage), the organoids were extracted with Cell Recovery Solution (Corning, Tewksbury, MA). RNA was then extracted utilizing the RNeasy micro kit (Qiagen, Redwood, CA). Total RNA was isolated by Trizol reagent (Invitrogen, Carlsbad, CA) according to manufacturer's instructions. RNA purity was assessed by the 260/280 and 260/230 nm absorbance ratio of 2 or greater. 1 µg of total RNA was used for cDNA synthesis using the Bio-Rad iScript cDNA synthesis kit (Bio-Rad Life Science, Hercules, CA). Analysis of gene expression via Reverse-Transcription Polymerase Chain Reaction (RTPCR) was performed on Bio-Rad C1000 Thermal cycler using Taq PCR Master Mix Kit (Qiagen, Valencia, CA) and intron spanning gene specific primers (Table S2). Quantitative Real-Time PCR (qPCR) was performed using cDNA using Light-Cycler Taqman Master (Roche Applied Science, Indianapolis, IN) and probes from the Universal Probe Library (Roche Applied Science) using intron spanning, gene specific primers. Relative expression levels were calculated by Δ–Δ Ct method. Actin was used to normalize the gene expression. Illumina HiSeq 3000 SRX18487386 GSM6784421_r1 SRP411259 PRJNA908254 SRS15960301 GSM6784421 GSM6784421 RNA-Seq TRANSCRIPTOMIC cDNA Once organoids reached confluence (~3 days after first passage), the organoids were extracted with Cell Recovery Solution (Corning, Tewksbury, MA). RNA was then extracted utilizing the RNeasy micro kit (Qiagen, Redwood, CA). Total RNA was isolated by Trizol reagent (Invitrogen, Carlsbad, CA) according to manufacturer's instructions. RNA purity was assessed by the 260/280 and 260/230 nm absorbance ratio of 2 or greater. 1 µg of total RNA was used for cDNA synthesis using the Bio-Rad iScript cDNA synthesis kit (Bio-Rad Life Science, Hercules, CA). Analysis of gene expression via Reverse-Transcription Polymerase Chain Reaction (RTPCR) was performed on Bio-Rad C1000 Thermal cycler using Taq PCR Master Mix Kit (Qiagen, Valencia, CA) and intron spanning gene specific primers (Table S2). Quantitative Real-Time PCR (qPCR) was performed using cDNA using Light-Cycler Taqman Master (Roche Applied Science, Indianapolis, IN) and probes from the Universal Probe Library (Roche Applied Science) using intron spanning, gene specific primers. Relative expression levels were calculated by Δ–Δ Ct method. Actin was used to normalize the gene expression. Illumina HiSeq 3000 SRX18487387 GSM6784422_r1 SRP411259 PRJNA908254 SRS15960302 GSM6784422 GSM6784422 RNA-Seq TRANSCRIPTOMIC cDNA Once organoids reached confluence (~3 days after first passage), the organoids were extracted with Cell Recovery Solution (Corning, Tewksbury, MA). RNA was then extracted utilizing the RNeasy micro kit (Qiagen, Redwood, CA). Total RNA was isolated by Trizol reagent (Invitrogen, Carlsbad, CA) according to manufacturer's instructions. RNA purity was assessed by the 260/280 and 260/230 nm absorbance ratio of 2 or greater. 1 µg of total RNA was used for cDNA synthesis using the Bio-Rad iScript cDNA synthesis kit (Bio-Rad Life Science, Hercules, CA). Analysis of gene expression via Reverse-Transcription Polymerase Chain Reaction (RTPCR) was performed on Bio-Rad C1000 Thermal cycler using Taq PCR Master Mix Kit (Qiagen, Valencia, CA) and intron spanning gene specific primers (Table S2). Quantitative Real-Time PCR (qPCR) was performed using cDNA using Light-Cycler Taqman Master (Roche Applied Science, Indianapolis, IN) and probes from the Universal Probe Library (Roche Applied Science) using intron spanning, gene specific primers. Relative expression levels were calculated by Δ–Δ Ct method. Actin was used to normalize the gene expression. Illumina HiSeq 3000 SRX18487388 GSM6784423_r1 SRP411259 PRJNA908254 SRS15960303 GSM6784423 GSM6784423 RNA-Seq TRANSCRIPTOMIC cDNA Once organoids reached confluence (~3 days after first passage), the organoids were extracted with Cell Recovery Solution (Corning, Tewksbury, MA). RNA was then extracted utilizing the RNeasy micro kit (Qiagen, Redwood, CA). Total RNA was isolated by Trizol reagent (Invitrogen, Carlsbad, CA) according to manufacturer's instructions. RNA purity was assessed by the 260/280 and 260/230 nm absorbance ratio of 2 or greater. 1 µg of total RNA was used for cDNA synthesis using the Bio-Rad iScript cDNA synthesis kit (Bio-Rad Life Science, Hercules, CA). Analysis of gene expression via Reverse-Transcription Polymerase Chain Reaction (RTPCR) was performed on Bio-Rad C1000 Thermal cycler using Taq PCR Master Mix Kit (Qiagen, Valencia, CA) and intron spanning gene specific primers (Table S2). Quantitative Real-Time PCR (qPCR) was performed using cDNA using Light-Cycler Taqman Master (Roche Applied Science, Indianapolis, IN) and probes from the Universal Probe Library (Roche Applied Science) using intron spanning, gene specific primers. Relative expression levels were calculated by Δ–Δ Ct method. Actin was used to normalize the gene expression. Illumina HiSeq 3000 SRX18487389 GSM6784424_r1 SRP411259 PRJNA908254 SRS15960304 GSM6784424 GSM6784424 RNA-Seq TRANSCRIPTOMIC cDNA Once organoids reached confluence (~3 days after first passage), the organoids were extracted with Cell Recovery Solution (Corning, Tewksbury, MA). RNA was then extracted utilizing the RNeasy micro kit (Qiagen, Redwood, CA). Total RNA was isolated by Trizol reagent (Invitrogen, Carlsbad, CA) according to manufacturer's instructions. RNA purity was assessed by the 260/280 and 260/230 nm absorbance ratio of 2 or greater. 1 µg of total RNA was used for cDNA synthesis using the Bio-Rad iScript cDNA synthesis kit (Bio-Rad Life Science, Hercules, CA). Analysis of gene expression via Reverse-Transcription Polymerase Chain Reaction (RTPCR) was performed on Bio-Rad C1000 Thermal cycler using Taq PCR Master Mix Kit (Qiagen, Valencia, CA) and intron spanning gene specific primers (Table S2). Quantitative Real-Time PCR (qPCR) was performed using cDNA using Light-Cycler Taqman Master (Roche Applied Science, Indianapolis, IN) and probes from the Universal Probe Library (Roche Applied Science) using intron spanning, gene specific primers. Relative expression levels were calculated by Δ–Δ Ct method. Actin was used to normalize the gene expression. Illumina HiSeq 3000 SRX18487390 GSM6784425_r1 SRP411259 PRJNA908254 SRS15960305 GSM6784425 GSM6784425 RNA-Seq TRANSCRIPTOMIC cDNA Once organoids reached confluence (~3 days after first passage), the organoids were extracted with Cell Recovery Solution (Corning, Tewksbury, MA). RNA was then extracted utilizing the RNeasy micro kit (Qiagen, Redwood, CA). Total RNA was isolated by Trizol reagent (Invitrogen, Carlsbad, CA) according to manufacturer's instructions. RNA purity was assessed by the 260/280 and 260/230 nm absorbance ratio of 2 or greater. 1 µg of total RNA was used for cDNA synthesis using the Bio-Rad iScript cDNA synthesis kit (Bio-Rad Life Science, Hercules, CA). Analysis of gene expression via Reverse-Transcription Polymerase Chain Reaction (RTPCR) was performed on Bio-Rad C1000 Thermal cycler using Taq PCR Master Mix Kit (Qiagen, Valencia, CA) and intron spanning gene specific primers (Table S2). Quantitative Real-Time PCR (qPCR) was performed using cDNA using Light-Cycler Taqman Master (Roche Applied Science, Indianapolis, IN) and probes from the Universal Probe Library (Roche Applied Science) using intron spanning, gene specific primers. Relative expression levels were calculated by Δ–Δ Ct method. Actin was used to normalize the gene expression. Illumina HiSeq 3000