Breast cancer has become the most common cancer that leads to women’s death. Breast cancer is a complex, highly heterogeneous disease classified into various subtypes based on histological features, which determines the therapeutic options. System identification of effective drugs for each subtype remains challenging. In this work, we present a computational network biology approach to screen precision drugs for different breast cancer subtypes by considering the impact intensity of candidate drugs on the pathway crosstalk mediated by miRNAs. Firstly, we constructed and analyzed the subtype-specific risk pathway crosstalk networks mediated by miRNAs. Then, we evaluated 36 Food and Drug Administration (FDA)-approved anticancer drugs by quantifying their effects on these subtype-specific pathway crosstalk networks and combining with survival analysis. Finally, some first-line treatments of breast cancer, such as Paclitaxel and Vincristine, were optimized for each subtype. In particular, we performed precision screening of subtype-specific therapeutic drugs and also confirmed some novel drugs suitable for breast cancer treatment. For example, Sorafenib was applicable for the basal subtype treatment, Irinotecan was optimum for Her2 subtype treatment, Vemurafenib was suitable for the LumA subtype treatment, and Vorinostat could apply to LumB subtype treatment. In addition, the mechanism of these optimal therapeutic drugs in each subtype of breast cancer was further dissected. In summary, our study offers an effective way to screen precision drugs for various breast cancer subtype treatments. We also dissected the mechanism of optimal therapeutic drugs, which may provide novel insight into the precise treatment of cancer and promote researches on the mechanisms of action of drugs.

Systematic investigation of tumor-infiltrating immune (TII) cells is important to the development of immunotherapies, and the clinical response prediction in cancers. There exists complex transcriptional regulation within TII cells, and different immune cell types display specific regulation patterns. To dissect transcriptional regulation in TII cells, we first integrated the gene expression profiles from single-cell datasets, and proposed a computational pipeline to identify TII cell type-specific transcription factor (TF) mediated activity immune modules (TF-AIMs). Our analysis revealed key TFs, such as BACH2 and NFKB1 play important roles in B and NK cells, respectively. We also found some of these TF-AIMs may contribute to tumor pathogenesis. Based on TII cell type-specific TF-AIMs, we identified eight CD8+ T cell subtypes. In particular, we found the PD1 + CD8+ T cell subset and its specific TF-AIMs associated with immunotherapy response. Furthermore, the TII cell type-specific TF-AIMs displayed the potential to be used as predictive markers for immunotherapy response of cancer patients. At the pan-cancer level, we also identified and characterized six molecular subtypes across 9680 samples based on the activation status of TII cell type-specific TF-AIMs. Finally, we constructed a user-friendly web interface CellTF-AIMs (http://bio-bigdata.hrbmu.edu.cn/CellTF-AIMs/) for exploring transcriptional regulatory pattern in various TII cell types. Our study provides valuable implications and a rich resource for understanding the mechanisms involved in cancer microenvironment and immunotherapy.