A whole genome association analysis identifies loci associated with Mycobacterium avium subsp. paratuberculosis infection status in US holstein cattle

The purpose of this study was to identify loci associated with Mycobacterium avium subspecies paratuberculosis ( Map ) infection status in US Holsteins using the Illumina BovineSNP50 BeadChip whole genome single nucleotide polymorphism (SNP) assay. Two hundred forty-five cows from dairies in New York, Pennsylvania and Vermont enrolled in longitudinal herd studies between January 1999 and November 2007 were assessed for the presence of Map in both faecal and tissue samples. An animal was considered tissue infected if any sample contained at least one colony forming unit of Map per gram of tissue (CFU/g) and the same definition was employed for faecal samples. Each animal was genotyped with the Illumina BovineSNP50 BeadChip and after quality assurance filtering, 218 animals and 45 683 SNPs remained. We sought to identify loci associated with four different case/control classifications: presence of Map in the tissue, presence of Map in faeces, presence of Map in both tissue and faeces and presence of Map in tissue but not faeces. A case–control genome wide association study was conducted to test the four different classifications of Map infection status (cases) when compared with a Map -negative control group (control). Regions on chromosomes 1, 5, 7, 8, 16, 21 and 23 were identified with moderate significance ( P  < 5 × 10⁻⁵). Two regions, one on chromosome 3 (near EDN2 ) and another on chromosome 9 (no positional gene candidates), were identified with a high level of association to the presence of Map in tissue and both tissue and faeces respectively ( P  <   5 × 10⁻⁷, genome-wide Bonferonni P  <   0.05).