We present a new model of ESR1 network regulation based on analysis of Doxorubicin, Estradiol, and TNFα combination treatment in MCF-7. We used Doxorubicin as a therapeutic agent, TNFα as marker and mediator of an inflammatory microenvironment and 17β-Estradiol (E2) as an agonist of Estrogen Receptors, known predisposing factor for hormone-driven breast cancer, whose pharmacological inhibition reduces the risk of breast cancer recurrence. Based on the results of transcriptomics analysis, we found 71 differentially expressed genes that are specific for the combination treatment with Doxorubicin + Estradiol + TNFα in comparison with single or double treatments. The responsiveness to the triple treatment was examined for seven genes by qPCR, of which six were validated, and then extended to four additional cell lines differing for p53 and/or ER status. The results of differential regulation enrichment analysis highlight the role of the ESR1 network that included 36 of 71 specific differentially expressed genes. We propose that the combined activation of p53 and NF-kB transcription factors significantly influences ligand-dependent, ER-driven transcriptional responses, also of the ESR1 gene itself. These results provide a model of coordinated interaction of TFs to explain the Doxorubicin, E2 and TNFα induced repression mechanisms.

Unlike other tumours, TP53 is rarely mutated in melanoma; however, it fails to function as a tumour suppressor. We assume that its functions might be altered through interactions with several families of proteins, including p53/p73, NME and GLI. To elucidate the potential interplay among these families we analysed the expression profiles of aforementioned genes and proteins in a panel of melanoma cell lines, metastatic melanoma specimens and healthy corresponding tissue. Using qPCR a higher level of NME1 gene expression and lower levels of Δ40p53β, ΔNp73, GLI1, GLI2 and PTCH1 were observed in tumour samples compared to healthy tissue. Protein expression of Δ133p53α, Δ160p53α and ΔNp73α isoforms, NME1 and NME2, and N'ΔGLI1, GLI1FL, GLI2ΔN isoforms was elevated in tumour tissue, whereas ∆Np73β was downregulated. The results in melanoma cell lines, in general, support these findings. In addition, we correlated expression profiles with clinical features and outcome. Higher Δ133p53β and p53α mRNA and both GLI1 mRNA and GLI3R protein expression had a negative impact on the overall survival. Shorter overall survival was also connected with lower p53β and NME1 gene expression levels. In conclusion, all examined genes may have implications in melanoma development and functional inactivity of TP53.