Comprehensive mRNA expression analysis uncovers HLA polymorphism-dependent responses to abacavir in keratinocytes. a) Scatter plot showing 3D-Gene DNA microarray hybridization intensities for abacavir-exposed keratinocytes of B*57:01-Tg compared to keratinocytes of B*57:03-Tg. b) Positively enriched GO gene sets in abacavir-exposed keratinocytes of B*57:01-Tg relative to B*57:03-Tg. In this panel, gene sets with significant differences (nominal [NOM] P < 0.01). c) Positively enriched KEGG gene sets in abacavir-exposed keratinocytes of B*57:01-Tg relative to B*57:03-Tg. KEGG pathway gene sets are listed according to their normalized enrichment score. Significant differences (* NOM P < 0.05) were observed between abacavir-exposed keratinocytes of B*57:01-Tg and B*57:03-Tg. d–i) mRNA expression levels of Hspa1a d), Ube2j1 e), Tmem129 f), Hspa1 l g), Pdgfra h), and Itpr1 i) relative to GAPDH in keratinocytes from B*57:01-Tg, their LM, and B*57:03-Tg. Keratinocytes were incubated with (+) or without (−) 100 μM abacavir for 12 h. mRNA levels are presented relative to the mean value in keratinocytes from B*57:01-Tg without abacavir exposure. Data are expressed as the mean ± SEM (n = 4/group). Significant differences (* P < 0.05, ** P < 0.01, *** P < 0.001) were observed compared with another group (one-way ANOVA, followed by Bonferroni’s multiple comparisons tests; d–i).