Recent studies suggest that tumor cells overexpressing aldoketoreductases (AKRs) exhibit increased resistance to DNA damaging agents such as anthracyclines. AKRs may induce resistance to the anthracycline doxorubicin by catalyzing its conversion to the less toxic 13-hydroxy metabolite doxorubicinol. However, it has not been established whether during selection for anthracycline resistance, AKR overexpression in tumor cells can be correlated with the onset or magnitude of drug resistance and with appreciable conversion of anthracyclines to 13-hydroxy metabolites. Through microarray and quantitative polymerase chain reaction studies involving rigid selection criteria and both correlative discriminate statistics and time-course models, we have identified several genes whose expression can be correlated with the onset and/or magnitude of anthracycline resistance, including AKR1C2 and AKR1C3. Also associated with the onset or magnitude of anthracycline resistance were genes involved in drug transport (ABCB1, ABCC1), cell signaling and transcription (RDC1, CXCR4), cell proliferation or apoptosis (BMP7, CAV1), protection from reactive oxygen species (AKR1C2, AKR1C3, FTL, FTH, TXNRD1, MT2A), and structural or immune system proteins (IFI30, STMN1). As expected, doxorubicin-resistant and epirubicin-resistant cells exhibited higher levels of doxorubicinol than wild-type cells, although at insufficient levels to account for significant drug resistance. Nevertheless, an inhibitor of Akr1c2 (5beta-cholanic acid) almost completely restored sensitivity to doxorubicin in ABCB1-deficient doxorubicin-resistant cells, while having no effect on ABCB1-expressing epirubicin-resistant cells. Taken together, we show for the first time that a variety of genes (particularly redox genes such as AKR1C2 and AKR1C3) can be temporally and causally correlated with the acquisition of anthracycline resistance in breast tumor cells.