Genotype by environment interaction and heteroscedasticity influence the expression of parasite resistance in Katahdin sheep

Authors

Brian C. Arisman, University of Nebraska-Lincoln
Joan M. Burke, USDA, Agricultural Research Service
James L. M. Morgan, Round Mountain Consulting
Ronald M. Lewis, University of Nebraska - LincolnFollow

Document Type

Article

Date of this Version

7-5-2023

Citation

Journal of Animal Science, 2023, 101, 1–10 https://doi.org/10.1093/jas/skad228

Comments

Open access.

Abstract

Increasingly, sheep producers are choosing breeds that express resistance to gastrointestinal parasites due to reduced efficacy of anthelminthic drugs. One such breed is Katahdin. Katahdins are raised in various climates and management systems in the United States, which can be combined into eco-management clusters to describe production environments more holistically. The objectives of this study were to determine if genotype by environment interaction (G × E) and heteroscedasticity existed across these eco-management clusters for traits indicative of parasite resistance. Body weights (BW), FAMACHA scores (FAM), and fecal egg counts (FEC) were collected at around 90 d in 3,527 Katahdin lambs delineated into nine eco-management clusters. A tri-variate animal model including birth-rearing type, sex, and dam age (as a quadratic covariate) as fixed effects, and eco-management cluster, direct additive, uncorrelated maternal environmental (for BW), and residual as random effects, was fitted with ASReml. Heritability estimates for BW, FEC, and FAM were 0.36 ± 0.07, 0.31 ± 0.07, and 0.26 ± 0.05, respectively. The genetic (additive) correlation between BW with FEC was −0.26 ± 0.08 and with FAM was −0.16 ± 0.08, and thereby favorable. Heritabilities were also estimated univariately within eco-management clusters and ranged from 0.30 ± 0.05 to 0.37 ± 0.05 for BW, 0.18 ± 0.12 to 0.50 ± 0.13 for FEC, and 0.07 ± 0.06 to 0.40 ± 0.19 for FAM. Significant genetic and phenotypic heteroscedasticity among eco-management clusters was detected in FEC and FAM. A sire by eco-management cluster interaction term was added to the initial model fitted to evaluate G × E. This interaction defined substantial variation (P < 0.01) in all traits and explained 12% (FEC) to 20% (BW) of the phenotypic variation. Accounting for G × E and heteroscedasticity in the design and implementation of breeding programs may introduce operational challenges. Still, doing so would improve the efficacy of selection programs to improve parasite resistance.