After batch correction (Supplementary Fig. 5a), we merged the K5:BRAFV600E and K14:BRAFV600E tumour samples as they clustered independently of the model and our analysis returned no differentially expressed genes (Supplementary Fig. 5b, c). This evidence concerns the gene KRT5 and neoplasm.