Next, we sought to understand if phenotypic changes that were measured as scores by TCGA23 show difference between our heterocellular subtypes in luminal-A tumors (Fig. 3e–l; scores in Fig. 3f–l are from reverse-phase protein microarray; RPPA as published by TCGA.23) Our analysis showed that tumor purity, hormone_a (represents signatures associated hormone receptors,30) proliferation and DNA damage response scores were significantly high in goblet-like and TA subtypes compared with the other subtypes (Fig. 3e–g and i). This evidence concerns the gene NR4A1 and neoplasm.