Genome-wide association study for Marek's disease mortality in layer chickens

A genome-wide association study (GWAS) using Bayesian variable selection was performed to determine genomic regions associated with mortality due to Marek's disease virus (MDV) infection in layers. Mortality (%) under experimental disease challenge (500 plaque-forming units of a very virulent plus MDV strain) was recorded for progeny groups (average 15.5 birds; range 3 to 30) of 253 genotyped sires from four generations of a brown-egg layer line. An additional generation of 43 sires with progeny data was used to validate results. Sires were genotyped with a 42K Illumina single-nucleotide polymorphism (SNP) chip. Methods BayesB (pi = 0.995) and BayesCpi, with or without weighting residuals by the size of progeny groups were applied. The proportion of genetic variance contributed by SNPs within each 1-megabase (Mb) genomic region was quantified. Average mortality was 33% but differed significantly between generations. Genetic markers explained about 11% of phenotypic variation in mortality. Correlations between genomic estimated breeding values and percentage of progeny mortality for the validation generation (sons of individuals in training) were 0.12, 0.17, 0.02, and 0.16 for BayesB, weighted BayesB, BayesCpi, and weighted BayesCpi, respectively, when using the whole genome, and 0.03, 0.20, -0.06, and 0.14, when using only SNP from the 10, 1-Mb regions, explaining the largest proportion of genetic variance according to each method. Results suggest that regions on chromosomes 2, 3, 4, 9, 15, 18, and 21 are associated with Marek's disease resistance and can be used for selection and that accounting for the size of progeny groups has a large impact on correct localization of such genomic regions.