The present invention relates to methods of predicting therapeutic efficacy of thiopurines in an individual by determining the presence of one or more risk variants. In one embodiment, the effective therapeutic efficacy of thiopurines is determined by the presence of risk variants at the genetic loci of HLA-DRB1, CREM, TAGAP, PLCL1, GPX4, SBNO2, MEF2A and / or LYSMD4. In another embodiment, the risk variants are located at the genetic loci of ARL4C, IL1R2, JAK2, 19q13, CARD9, SNAPC4, and / or 8q24. In another embodiment, the individual is has been diagnosed with inflammatory bowel disease.