GSEA was performed using the fgsea package (Korotkevich and Sukhov) in R version 3.6.3. Unbiased gene ontology analyze was performed with the C5 collection adapted for mouse that contains gene sets annotated by Gene Ontology terms that has been downloaded on http://bioinf.wehi.edu.au/software/MSigDB/. The function “fgseaMultilevel” has been used with ranked fold change values corresponding to Krt8-YFP+ cells over Krt8-YFP- cells from cervical esophagus [abs(LFC) > 2; FDR < 0.05].