Stratification of Estrogen Receptor-Negative Breast Cancer Patients by Integrating the Somatic Mutations and Transcriptomic Data

Patients with estrogen receptor-negative breast cancer generally have a worse prognosis than estrogen receptor-positive patients. Nevertheless, a significant proportion of the estrogen receptor-negative cases have favorable outcomes. Identifying patients with a good prognosis, however, remains difficult, as recent studies are quite limited. The identification of molecular biomarkers is needed to better stratify patients. The significantly mutated genes may be potentially used as biomarkers to identify the subtype and to predict outcomes. To identify the biomarkers of receptor-negative breast cancer among the significantly mutated genes, we developed a workflow to screen significantly mutated genes associated with the estrogen receptor in breast cancer by a gene coexpression module. The similarity matrix was calculated with distance correlation to obtain gene modules through a weighted gene coexpression network analysis. The modules highly associated with the estrogen receptor, called important modules, were enriched for breast cancer-related pathways or disease. To screen significantly mutated genes, a new gene list was obtained through the overlap of the important module genes and the significantly mutated genes. The genes on this list can be used as biomarkers to predict survival of estrogen receptor-negative breast cancer patients. Furthermore, we selected six hub significantly mutated genes in the gene list which were also able to separate these patients. Our method provides a new and alternative method for integrating somatic gene mutations and expression data for patient stratification of estrogen receptor-negative breast cancers.

Interrelation of androgen receptor and miR-30a and miR-30a function in ER-, PR-, AR+ MDA-MB-453 breast cancer cells

The association between androgen-induced androgen receptor (AR) activating signal and microRNA (miR)-30a was investigated, as well as the function of miR-30a in estrogen receptor-negative (ER^-), progesterone receptor-negative (PR^-), and AR-positive (AR⁺) MDA-MB-453 breast cancer cells. Androgen-induced AR activating signal upregulated the expression of AR, and downregulated the expression of miR-30a, b and c. Bioinformatics analysis indicated a putative miR-30a, b and c binding site in the 3'-untranslated region of AR mRNA. It was confirmed that the AR gene is a direct target of miR-30a, whereas AR does not target the miR-30a promoter, and AR activating signal may indirectly downregulate miR-30a through other cell signaling pathways. In this positive feedback mechanism AR is then upregulated through miR-30a. Overexpression of miR-30a inhibited cell proliferation, whereas inhibition of miR-30a expression by specific antisense oligonucleotides, increased cell growth. Previously, androgen-induced AR activating signal was demonstrated to inhibit cell proliferation in ER^-, PR^- and AR⁺ MDA-MB-453 breast cancer cells, but AR activating signal downregulated the expression of miR-30a, relieving the inhibition of MDA-MB-453 cell growth. Therefore, in MDA-MB-453 breast cancer cells, miR-30a has two different functions regarding cell growth: Inhibition of cell proliferation through a positive feedback signaling pathway; and the relative promotion of cell proliferation through downregulation of miR-30a. Thus, the association between AR activating signal and microRNAs is complex, and microRNAs may possess different functions due to different signaling pathways. Although the results of the present study were obtained in one cell line, they contribute to subsequent studies on ER^-, PR^- and AR⁺ breast cancer.