Differences in individual drug responses are obstacles in breast cancer (BRCA) treatment, so predicting responses would help to plan treatment strategies. The accumulation of cancer molecular profiling and drug response data provide opportunities and challenges to identify novel molecular signatures and mechanisms of tumor responsiveness to drugs in BRCA. This study evaluated drug responses with a multi-omics integrated system that depended on long non-coding RNAs (lncRNAs). We identified drug response-related lncRNAs (DRlncs) by combining expression data of lncRNA, microRNA, messenger RNA, methylation levels, somatic mutations, and the survival data of cancer patients treated with drugs. We constructed an integrated and computational multi-omics approach to identify DRlncs for diverse chemotherapeutic drugs in BRCA. Some DRlncs were identified with Adriamycin, Cytoxan, Tamoxifen, and all samples for BRCA patients. These DRlncs showed specific features regarding both expression and computational accuracies. The DRlnc-gene co-expression networks were constructed and analyzed. Key DRlncs, such as HOXA-AS2 (Ensembl: ENSG00000253552), in the drug Adriamycin were characterized. The experimental analysis also suggested that HOXA-AS2 (Ensembl: ENSG00000253552) was a key DRlnc in Adriamycin drug resistance in BRCA patients. Some DRlncs were associated with survival and some specific functions. A possible mechanism of DRlnc HOXA-AS2 (Ensembl: ENSG00000253552) in the Adriamycin drug response for BRCA resistance was inferred. In summary, this study provides a framework for lncRNA-based evaluation of clinical drug responses in BRCA. Understanding the underlying molecular mechanisms of drug responses will facilitate improved responses to chemotherapy and outcomes of BRCA treatment.Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.