Abstract Background Endocrine therapy reduces breast cancer mortality by 40%, but resistance remains a major clinical problem. In this study, we sought to investigate the impact of aromatase inhibitor (AI) therapy on gene expression and identify gene modules representing key biological pathways that relate to early AI therapy resistance. Methods Global gene expression was measured on pairs of core-cut biopsies taken at baseline and at surgery from 254 patients with ER-positive primary breast cancer randomised to receive 2-week presurgical AI (n = 198) or no presurgical treatment (control n = 56) from the POETIC trial. Data from the AI group was adjusted to eliminate artefactual process-related changes identified in the control group. The response was assessed by changes in the proliferation marker, Ki67. Results High baseline ESR1 expression associated with better AI response in HER2+ tumours but not HER2− tumours. In HER2− tumours, baseline expression of 48 genes associated with poor antiproliferative response (p