9 (i.e. total KLH-binding NAb titers) will result in similar responses of the other NAb (iso)types. The GWAS were performed with 57,636 single nucleotide polymorphisms (SNP) for total KLH-binding NAb titers, and for the isotypes IgM, IgA, and IgG on 1,628 white purebred laying chickens of the same line as the base population. One genomic region was significantly associated to KLH-binding IgM NAb titers, and to a lesser extent to total KLH-binding NAb titers. The region showed full dominance, and had a major effect on IgM. This region is located on chromosome 4, and contained two Toll-like receptors (TLR). To further characterize the found association, total antibody concentration, and antibody concentrations of the isotypes IgM, IgA, and IgG were measured as well, and GWAS were performed. One genomic region, the same as identified before, was significantly associated to IgM antibody concentration. Full sequence data of key ancestors of the study population allowed imputation to full genome sequence for further association: 16 candidate genes were identified. SNP located in coding regions of these candidate genes were checked for predicted changes in protein functioning. One SNP, a C/G polymorphism at 69,965,939 base pairs (Gallus_gallus-5.0), received the maximum impact score from two independent prediction tools, which makes this SNP the most likely causal variant. The C-variant had an allele frequency of 0.45, showed a dominant mode of gene action, and was associated with high IgM levels. The G-variant had an allele frequency of 0.55, and was associated with low IgM levels. This SNP is located in TLR1A, which suggests a fundamental role of TLR1A on regulation of IgM levels, or B cells, or both.