Natural killer (NK) cell phenotype is partially mediated through binding of killer immunoglobulin-like receptors (KIR) with HLA class I ligands. The KIR gene family is highly polymorphic and not well captured by standard GWAS approaches. Here we tested the hypothesis that variations in KIR gene content combined with HLA class I ligand status is associated with keratinocyte skin cancers using a population-based study of basal cell carcinoma (BCC) and squamous cell carcinomas (SCC). We conduced an interaction analysis of KIR gene content variation and HLA-B (Bw4 vs Bw6) and HLA-C (C1 vs C2). KIR centromeric B haplotype was associated with significant risk of multiple BCC tumors (OR=2.39, 95% CI: 1.10–5.21), and there was a significant interaction between HLA-C and the activating gene KIR2DS3 for BCC (pinteraction = 0.005). Furthermore, there was significant interaction between HLA-B and telomeric KIR B haplotype (containing the activating genes KIR3DS1 and KIR2DS1) as well as HLA-B and the activating KIR gene 2DS5 (pinteraction 0.001 and 0.012, respectively). Similar but greatly attenuated associations were observed for SCC. Moreover, previous in vitro models demonstrated that p53 is required for upregulation of NK ligands, and accordingly, we observed there was a strong association between the KIR B haplotype and p53 alteration in BCC tumors, with a higher likelihood that KIR B carriers harbor abnormal p53 (p<0.004). Taken together, our data suggest functional interactions between KIR and HLA modify risks of BCC and SCC, and that KIR encoded by the B genes provide selective pressure for altered p53 in BCC tumors.