Even following Bonferroni correction for testing multiple cytokines (p < 0.05/48 = 1.1 × 10−3), three cytokines, including MCP1, MIP1b and IL13, still showed statistically significant associations with risk of total BC or subtype BC. The gene discussed is IL13; the disease is breast cancer.