Abstract:Aberrant biological information occurs naturally at exposure to bisphenol A or its alternatives, which was associated with the occurrence and development of breast cancer. However, the potential molecular variation in gene expression during the breast tumor development is still unclear. Herein, high throughput RNA sequencing(RNA-Seq) and bioinformatics analysis were used to investigate the variation of tumor-mRNA profile exposed with BPS5 (5 μg/kg bw/day) or BPS50(50 μg/kg bw/day) in tumor development-associated MMTV-PyMT transgenic mouse model. Meanwhile, we analyzed the dose-effects of bisphenol S(BPS) and BPS-induced tumor development on the gene level exhaustively. In dose-effect aspects of BPS, the increased concentration of BPS significantly changed the numbers and enrichment pathway of differentially expressed genes(DEGs), especially the enrichment pathways involved in up-regulated genes including ribosome, peroxisome proliferators-activated receptor(PPAR) signaling pathway and progesterone-mediated oocyte maturation pathway. In effects of BPS exposure to tumor development, expression of IgκC, Zfp385b, Cldn10, Pgr and Snord14d has changed significantly throughout the tumor development. Gene ontology(GO) and Kyoto encyclopedia of genes and genomes(KEGG) results obtained from BPS-induced tumor development showed that the functional classifications were intensively altered with an extension of time in high-dose BPS groups. The acquired DEGs and pathway information could help with the accurate exploration of molecular mechanisms of tumor development, screening of molecular targets of breast cancer, and toxicological evaluation of environmental pollutants.

RNA-sequencing-based Gene Expression Profile Revealing Breast Tumor Development Induced by Exposure of Bisphenol S