Breast cancer ranks first in terms of mortality and incidence rates worldwide among women. The HER2+ molecular subtype is one of the most aggressive subtypes; its treatment includes neoadjuvant chemotherapy and the use of a HER2 antibody. Some patients develop resistance despite positive results obtained using this therapeutic strategy. Objective. To identify prognostic markers for treatment and survival in HER2+ patients. Methods. Patients treated with neoadjuvant chemotherapy were assigned to sensitive and resistant groups based on their treatment response. Differentially expressed genes (DEGs) were identified using RNA-seq analysis. KEGG pathway, gene ontology, and interactome analyses were performed for all DEGs. An enrichment analysis Gene set enrichment analysis was performed. All DEGs were analyzed for overall (OS) and disease-free survival (DFS). Results. A total of 94 DEGs were related to treatment resistance. Survival analysis showed that 12 genes (ATF6B, DHRS13, DIRAS1, ERAL1, GRIN2B, L1CAM, IRX3, PRTFDC1, PBX2, S100B, SLC9A3R2, and TNXB) were good predictors of disease-free survival, and eight genes (GNG4, IL22RA2, MICA, S100B, SERPINF2, HLA-A, DIRAS1, and TNXB) were good predictors of overall survival (OS). Conclusion: We highlighted a molecular expression signature that can differentiate the treatment response, overall survival, and DFS of patients with HER2+ breast cancer.