Heritability of susceptibility to Johne's disease in cattle has been shown to vary from 0.041 to 0.159. Although the presence of genetic variation involved in susceptibility to Johne's disease has been demonstrated, the understanding of genes contributing to the genetic variance is far from complete. The objective of this study was to contribute to further understanding of genetic variation involved in susceptibility to Johne's disease by identifying associated chromosomal regions using a genome-wide association approach. Log-transformed ELISA test results of 265,290 individual Holstein-Friesian cows from 3,927 herds from the Netherlands were analyzed to obtain sire estimated breeding values for Mycobacterium avium subspecies paratuberculosis (MAP)-specific antibody response in milk using a sire-maternal grandsire model with fixed effects for parity, year of birth, lactation stage, and herd; a covariate for milk yield on test day; and random effects for sire, maternal grandsire, and error. For 192 sires with estimated breeding values with a minimum reliability of 70%, single nucleotide polymorphism (SNP) typing was conducted by a multiple SNP analysis with a random polygenic effect fitting 37,869 SNP simultaneously. Five SNP associated with MAP-specific antibody response in milk were identified distributed over 4 chromosomal regions (chromosome 4, 15, 18, and 28). Thirteen putative SNP associated with MAP-specific antibody response in milk were identified distributed over 10 chromosomes (chromosome 4, 14, 16, 18, 19, 20, 21, 26, 27, and 29). This knowledge contributes to the current understanding of genetic variation involved in Johne's disease susceptibility and facilitates control of Johne's disease and improvement of health status by breeding.