BACKGROUND: Breast cancer is one of the most common malignancies in women. In a previous study, we found that for two patients who had a high risk of lymphatic metastasis, lymphatic metastasis did not occur; whereas, for two patients who had a low risk of lymphatic metastasis, lymphatic metastasis did occur.

MATERIAL AND METHODS: We analyzed the differential gene expressions of these four patients by RNA-sequence. The data (HRNM_T versus HRNM_N, LRYM_T versus LRYM_N, and HRNM_T versus LRYM_T) was then processed using differentially expressed genes (DEGs) analysis, functional analysis for DEGs, and PPI network construct.

RESULTS: For HRNM_T versus HRNM_N, there were 224 DEGs. There were 504 DEGs for LRYM_T versus LRYM_N, and 88 DEGs for LRYM_T versus LRYM_N. For HRNM_T versus HRNM_N, DEGs were up-regulated mainly in the cell cycle, the IL-17 signaling pathway, and the progesterone-mediated oocyte maturation; DEGs were down-regulated mainly in the IL-17 signaling pathway. For LRYM_T versus LRYM_N, DEGs were up-regulated mainly in protein digestion and absorption, and cytokine-cytokine receptor interaction; DEGs were down-regulated mainly in ECM-receptor interaction. For HRNM_T versus LRYM_T, DEGs were up-regulated mainly in the PPAR signaling pathway; DEGs were downregulated mainly in the adipocytokine signaling pathway. The DEGs were screened to construct PPI networks.

CONCLUSIONS: The GO and KEGG functional enrichments of HRNM_T versus HRNM_N, and LRYM_T versus LRYM_N were consistent with earlier studies. For HRNM_T versus LRYM_T, DEGs were up-regulated mainly in PPAR signaling; DEGs were down-regulated mainly in the adipocytokine pathway.

Keywords: Gene Expression Profiling, Protein Interaction Maps, RNA 3' Polyadenylation Signals