CHGA and neoplasm: To compare lineage pathways in the tumor samples analyzed, we calculated and plotted z-scores to identify outlier genes for each lineage (color and boxes) including basal (Itga6, Cd44, and Col17a1), squamous (Tgm1), luminal (Gpx2 and Erbb2), EMT-claudin (Snail1, Snail2, Cldn7), Neuronal differentiation (Sox2 and Chga), Immune (Saa3), and luminal (Gpx2 and Upk2) shown together (Supplementary Fig. 11B) or as independent lineages (Supplementary Fig. 11C).