RNA-seq libraries were prepared from 100 FACS-sorted cells/sample obtained from the pancreases of reporter
Isl1CKO-Ai14 mutant (n = 5 samples) and reporter control-
Ai14 (n = 6) from E14.5 embryos; and
Isl1CKO-Ai14 mutant (n = 6) and reporter control-
Ai14 (n = 5) from P9 mice. Each sample contained 100 tdTomato
+ endocrine cells. Following the manufacturer's instructions, the NEB Next single-cell low input RNA library prep kit for Illumina was used for cDNA synthesis, amplification, and library generation [67 (
link)] at the Gene Core Facility (Institute of Biotechnology CAS, Czechia). Fragment Analyzer assessed the quality of cDNA libraries. The libraries were sequenced on an Illumina NextSeq 500 next-generation sequencer. NextSeq 500/550 High Output kit 75 cycles (Illumina #200,024,906) were processed at the Genomics and Bioinformatics Core Facility (Institute of Molecular Genetics CAS, Czechia). RNA-Seq reads in FASTQ files were mapped to the mouse genome using STAR [version 2.7.0c [68 (
link)]] GRCm38 primary assembly and annotation version M8. The raw data of RNA sequencing were processed with a standard pipeline. Using cutadapt v1.18 [69 (
link)], the number of reads (minimum, 32 million; maximum, 73 million) was trimmed by Illumina sequencing adaptor and of bases with reading quality lower than 20, subsequently reads shorter than 20 bp were filtered out TrimmomaticPE version 0.36 [70 (
link)].
Ribosomal RNA and reads mapping to UniVec database were filtered out using bowtie v1.2.2. with parameters -S -n 1 and SortMeRNA [71 (
link)]. A count table was generated by Rsubread v2.0.1 package using default parameters without counting multi mapping reads. The raw RNA-seq data were deposited at GEO: (
https://www.ncbi.nlm.nih.gov/geo/).
DESeq2 [v1.26.0 [72 (
link)]] default parameters were used to normalize data and compare the different groups. Genes were then filtered using the criteria of an adjusted P-value P
adj < 0.05, and a base mean ≥ 50, and Fold change > 1.5 for upregulated genes and < 0.5 for downregulated genes for both E14.5 and P9 data to identify differentially expressed genes between
Isl1CKO and control endocrine cells. The enrichment of the functional categories and functional annotation clustering of the differentially expressed genes was performed using g: Profiler [73 (
link)] using version e104_eg51_p15_3922dba with g: SCS multiple testing correction methods applying a significance threshold of 0.05. Transcription factor (TF) enrichment analysis (TFEA) [41 (
link)] was used to identify the enrichment of TF target genes in our set of differentially expressed genes. The top seven enriched TFs are listed (Additional file
6: Dataset S1c).
Bohuslavova R., Fabriciova V., Lebrón-Mora L., Malfatti J., Smolik O., Valihrach L., Benesova S., Zucha D., Berkova Z., Saudek F., Evans S.M, & Pavlinkova G. (2023). ISL1 controls pancreatic alpha cell fate and beta cell maturation. Cell & Bioscience, 13, 53.