Search

Background: Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (N e ), runs of homozygosity (ROH) and Wright's fixation index (F ST ) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs)., Results: Katahdin sheep had the largest current N e at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the F ST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the F ST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all F ST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms., Conclusions: In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity., (© 2024. The Author(s).)

The Katahdin hair breed gained popularity in the United States as low input and prolific, with a propensity to exhibit parasite resistance. With the introduction of genomically enhanced estimated breeding values (GEBV) to the Katahdin genetic evaluation, defining the diversity present in the breed is pertinent. Utilizing pedigree records (n = 92,030) from 1984 to 2019 from the National Sheep Improvement Program, our objectives were to (i) estimate the completeness and quality of the pedigree, (ii) calculate diversity statistics for the whole pedigree and relevant reference subpopulations and (iii) assess the impact of current diversity on genomic selection. Reference 1 was Katahdins born from 2017 to 2019 (n = 23,494), while reference 2 was a subset with at least three generations of Katahdin ancestry (n = 9327). The completeness of the whole pedigree, and the pedigrees of reference 1 and reference 2, were above 50% through the fourth, fifth and seventh generation of ancestors, respectively. Effective population size (N e ) averaged 111 animals with a range from 42.2 to 451.0. The average generation interval was 2.9 years for the whole pedigree and reference 1, and 2.8 years for reference 2. The mean individual inbreeding and average relatedness coefficients were 1.62% and 0.91%, 1.74% and 0.90% and 2.94% and 1.46% for the whole pedigree, reference 1, and reference 2, respectively. There were over 300 effective founders in the whole pedigree and reference 1, with 169 in reference 2. Effective number of ancestors were over 150 for the whole pedigree and reference 1, while there were 67 for reference 2. Prediction accuracies increased as the reference population grew from 1k to 7.5k and plateaued at 15k animals. Given the large number of founders and ancestors contributing to the base genetic variation in the breed, the N e is sufficient to maintain diversity while achieving progress with selection. Stable low rates of inbreeding and relatedness suggest that incorporating genetic conservation in breeding decisions is currently not of high priority. Current N e suggests that with limited genotyping, high levels of accuracy for genomic prediction can be achieved. However, intense selection on GEBV may cause loss of genetic diversity long term., (© 2023 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

Gastrointestinal nematodes (GIN) negatively affect the performance and well-being of sheep. Due to anthelmintic resistance, GIN are difficult to control leading producers to choose breeds that can exhibit resistance to parasitism. An example is Katahdin sheep. Katahdins are raised in various climates and management systems in the United States. These environmental factors can be combined to form eco-management groupings or clusters. We hypothesized that GIN challenge varies predictably based on the characteristics of these environmental clusters. Forty Katahdin producers from across the United States were surveyed for management information, with body weights (BW), fecal egg counts (FEC), and FAMACHA scores (FAM) available from 17 of the 40 flocks. The performance data included 3,426 lambs evaluated around 90 d of age. Management and climate data were combined into clusters using multiple correspondence and principal component (PC) analysis. Performance data were aligned with their corresponding cluster. Depending on the trait, eco-management cluster, birth-rearing type, sex, and, as a covariate, dam age, were fitted as systematic effects with ANOVA. Clusters also were formed based on climate or management data alone. When compared with fitting the eco-management clusters, they defined less variation in each of the traits based on Akaike and Bayesian information criterion, and adjusted r2 values. To further examine variation defined by eco-management clusters, residuals from an ANOVA model excluding eco-management cluster were retained, and their correlation with PC loadings calculated. All PC loadings were included as potential independent variables and tested for significance using backward stepwise regression. The PC loadings with a correlation |≥0.49| explained significant variation in each trait and were included in the final models chosen; adjusted r2 values for BW, FEC, and FAM were 0.90, 0.81, and 0.97, respectively. When analyzing GIN challenge, eco-management clusters corresponding with hotter temperatures and greater rainfall, and with pasture-born lambs, suffered greater parasitism. Conversely, the eco-management clusters with lambs turned out to pasture at older ages benefited from reduced parasitism. Through the formation of eco-management clusters, an environmental variable can be defined to study interactions of genotypes to their environment, providing a potentially useful tool for identifying parasite-resistant sheep., (© The Author(s) 2023. Published by Oxford University Press on behalf of the American Society of Animal Science.)

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of the American Society of Animal Science.)

Estimated breeding values (EBV) for fecal egg counts (FEC) at 42-90 days of age (WFEC) and 91-150 days of age (PFEC) for 84 progeny-tested Katahdin sires were used to identify associations of deregressed EBV with single-nucleotide polymorphisms (SNP) using 388,000 SNP with minor-allele frequencies ≥0.10 on an Illumina high-density ovine array. Associations between markers and FEC EBV were initially quantified by single-SNP linear regression. Effects of linkage disequilibrium (LD) were minimized by assigning SNP to 2,535 consecutive 1-Mb bins and focusing on the effect of the most significant SNP in each bin. Bonferroni correction was used to define bin-based (BB) genome- and chromosome-wide significance. Six bins on chromosome 5 achieved BB genome-wide significance for PFEC EBV, and three of those SNP achieved chromosome-wide significance after Bonferroni correction based on the 14,530 total SNP on chromosome 5. These bins were nested within 12 consecutive bins between 59 and 71 Mb on chromosome 5 that reached BB chromosome-wide significance. The largest SNP effects were at 63, 67, and 70 Mb, with LD among these SNP of r 2 ≤ 0.2. Regional heritability mapping (RHM) was then used to evaluate the ability of different genomic regions to account for additive variance in FEC EBV. Chromosome-level RHM indicated that one 500-SNP window between 65.9 and 69.9 Mb accounted for significant variation in PFEC EBV. Five additional 500-SNP windows between 59.3 and 71.6 Mb reached suggestive ( p < 0.10) significance for PFEC EBV. Although previous studies rarely identified markers for parasite resistance on chromosome 5, the IL12B gene at 68.5 Mb codes for the p40 subunit of both interleukins 12 and 23. Other immunoregulatory genes are also located in this region of chromosome 5, providing opportunity for additive or associative effects., Competing Interests: Author JLM is the co-owner of Round Mountain Consulting. Author TS is employed by Acceligen. Author MM is the owner of Delta G. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Notter, Heidaritabar, Burke, Shirali, Murdoch, Morgan, Morota, Sonstegard, Becker, Spangler, MacNeil and Miller.)

Gastrointestinal nematodes (GIN) pose a severe threat to sheep production worldwide. Anthelmintic drug resistance coupled with growing concern regarding potential environmental effects of drug use have demonstrated the necessity of implementing other methods of GIN control. The aim of this study was to test for genetic variants associated with resistance or susceptibility to GIN in Katahdin sheep to improve the current understanding of the genetic mechanisms responsible for host response to GIN. Linear regression and case-control genome-wide association studies were conducted with high-density genotype data and cube-root transformed weaning fecal egg counts (tFEC) of 583 Katahdin sheep. The case-control GWAS identified two significant SNPs ( P -values 1.49e-08 to 1.01e-08) within introns of the gene adhesion G protein-coupled receptor B3 ( ADGRB3 ) associated with lower fecal egg counts. With linear regression, four significant SNPs ( P -values 7.82e-08 to 3.34e-08) were identified within the first intron of the gene EGF-like repeats and discoidin domains 3 ( EDIL3 ). These identified SNPs were in very high linkage disequilibrium ( r 2 of 0.996-1), and animals with alternate homozygous genotypes had significantly higher median weaning tFEC phenotypes compared to all other genotypes. Significant SNPs were queried through public databases to identify putative transcription factor binding site (TFBS) and potential lncRNA differences between reference and alternate alleles. Changes in TFBS were predicted at two SNPs, and one significant SNP was found to be within a predicted lncRNA sequence with greater than 90% similarity to a known lncRNA in the bovine genome. The gene EDIL3 has been described in other species for its roles in the inhibition and resolution of inflammation. Potential changes of EDIL3 expression mediated through lncRNA expression and/or transcription factor binding may impact the overall immune response and reduce the ability of Katahdin sheep to control GIN infection. This study lays the foundation for further research of EDIL3 and ADGRB3 towards understanding genetic mechanisms of susceptibility to GIN, and suggests these SNPs may contribute to genetic strategies for improving parasite resistance traits in sheep., Competing Interests: Author JLM was the co-owner of Round Mountain Consulting. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Becker, Burke, Lewis, Miller, Morgan, Rosen, Van Tassell, Notter and Murdoch.)

Resistance to gastrointestinal nematodes has previously been shown to be a moderately heritable trait in some breeds of sheep, but the mechanisms of resistance are not well understood. Selection for resistance currently relies upon faecal egg counts (FEC), blood packed cell volumes and FAMACHA visual indicator scores of anaemia. Identifying genomic markers associated with disease resistance would potentially improve the selection process and provide a more reliable means of classifying and understanding the biology behind resistant and susceptible sheep. A GWAS was conducted to identify possible genetic loci associated with resistance to Haemonchus contortus in Katahdin sheep. Forty animals were selected from the top and bottom 10% of estimated breeding values for FEC from a total pool of 641 sires and ram lambs. Samples were genotyped using Applied Biosystems™ Axiom™ Ovine Genotyping Array (50K) consisting of 51 572 SNPs. Following quality control, 46 268 SNPs were included in subsequent analyses. Analyses were conducted using a linear regression model in plink v1.90 and a single-locus mixed model in snp and variation suite. Genome-wide significance was determined by a Bonferroni correction for multiple testing. Using linear regression, loci on chromosomes 2, 3, 16, 23 and 24 were significantly associated at the genome level with FEC estimated breeding values, and we identified a region on chromosome 2 that was significant using both statistical analyses. We suggest a potential role for the gene DIS3L2 for gastrointestinal nematode resistance in Katahdin sheep, although further research is needed to validate these findings., (© 2020 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

This study estimated genetic parameters for ewe reproductive traits [number of lambs born (NLB) and weaned (NLW) per ewe lambing] and fecal egg counts (FEC) during the peri-parturient rise (PPR) for use in genetic evaluation of Katahdin sheep. Data included NLB and NLW for 23,060 lambings by 9,295 Katahdin ewes, 1,230 PPR at lambing (PPR0) for 750 ewes, 1,070 PPR at approximately 30 d postpartum (PPR30) for 611 ewes, BW at birth, weaning, and (or) post-weaning for 12,869 lambs, and FEC at weaning and (or) post-weaning for 4,676 lambs. Direct additive, permanent environmental, and residual (co)variances were estimated in univariate and bivariate animal models. Fixed effects included effects of ewe management group and ewe age for all traits, and, for PPR, a continuous effect of days between lambing and measurement. Effects of litter size on PPR0 and number of lambs suckled on PPR30 were included in univariate models but excluded from bivariate models for PPR and NLB or NLW. Heritability estimates in univariate models for NLB, NLW, PPR0, and PPR30 were 0.09 ± 0.01, 0.06 ± 0.01, 0.35 ± 0.06, and 0.24 ± 0.07, respectively. Estimates of permanent environmental variance as a proportion of total phenotypic variance were 0.02 ± 0.01 for NLB, 0.03 ± 0.01 for NLW, 0.05 ± 0.06 for PPR0, and 0.13 ± 0.07 for PPR30. Direct additive, phenotypic, permanent environmental, and residual correlations between NLB and NLW were 0.88 ± 0.03, 0.74 ± 0.004, 0.54 ± 0.15, 0.74 ± 0.003, respectively; corresponding correlations between PPR0 and PPR30 were 0.96 ± 0.07, 0.46 ± 0.03, 0.98 ± 0.50, 0.18 ± 0.05, respectively. The additive genetic correlation (rd) between ewe reproductive traits and PPR ranged from 0.12 to 0.18. Estimates of rd between lamb BW and subsequent ewe NLB and NLW ranged from 0.07 to 0.20, and those between PPR and lamb BW ranged from -0.03 to 0.29. The rd between ewe reproductive traits and lamb FEC ranged from 0.27 to 0.40, and those between PPR and lamb FEC ranged from 0.56 to 0.77. Correlations between maternal additive effects on BW and direct additive effects on PPR were low (-0.08 to 0.10), and those between maternal additive effects on BW and direct additive effects on ewe reproductive traits were variable (-0.36 to 0.11). We conclude that FEC in growing lambs and peri-parturient ewes are controlled by similar genes and that modest, but manageable, genetic antagonisms may exist between FEC and ewe productivity.

Reliance on anthelmintic drugs to control internal parasites in sheep is no longer sustainable because of the development of resistance to these drugs in parasite populations. Genetic selection may offer an alternative long-term solution, as differences in parasite resistance exist both within and among sheep breeds. However, selection for parasite resistance may have correlated effects on other production traits. The objectives of this study were to estimate genetic parameters for weaning (WFEC) and postweaning (PWFEC) fecal egg counts (FEC) and assess their relationship with birth (BWT), weaning (WWT), and postweaning (PWWT) BW in Katahdin lambs. The study used WFEC (n = 2,537), PWFEC (n = 3.421), BWT (n = 12,869), WWT (n = 10,961), and PWWT (n = 7,812) from 12,869 lambs measured between 2003 and 2015 in 13 flocks enrolled in the U.S. National Sheep Improvement Program. Animal and sire models were fitted to the data using the ASReml statistical package. Records were corrected for fixed effects of dam age, joint effect of type of birth and rearing, and management group (defined by joint effects of flock, sex, and birth year and season); lamb age in days at each measurement time was fitted as a covariate. Maternal additive and maternal permanent environmental effects were not significant (P > 0.05), but litter effects influenced (P < 0.01) both WFEC and PWFEC. Heritability estimates ranged from 0.18 to 0.26 for WFEC and 0.23 to 0.46 for PWFEC, depending on the model used. Heritability estimates from sire models were higher than estimates from animal models. Direct additive, litter, residual, and phenotypic correlations between WFEC and PWFEC were 0.82, 0.25, 0.15, and 0.29, respectively. Bivariate analyses revealed low to moderate correlations between BW and FEC. Moderate heritabilities for FEC in this study indicated that genetic progress for this trait can be achieved in Katahdin lambs and that selection for low FEC should have little or no effect on BW.

The aim of this study was to estimate genetic parameters for BW in Katahdin lambs. Six animal models were used to study direct and maternal effects on birth weight (BWT), weaning weight (WWT), and postweaning weight (PWWT) using 41,066 BWT, 33,980 WWT, and 22,793 PWWT records collected over 17 yr in 100 flocks. Models included fixed effects of management group, dam age, type of birth (for BWT) or birth and rearing (for WWT and PWWT), and lamb age at weighing (fitted as a covariate for WWT and PWWT; all < 0.05). Variance components for random effects were estimated in sequentially more complex models and tested for significance with likelihood-ratio tests. A model that fitted only an additive animal effect overestimated additive variance for all BW, resulting in larger estimates of direct heritability than models that included maternal effects. Maternal effects explained variation ( < 0.05) in all BW. Heritability estimates for optimal models were 0.15 ± 0.01 for BWT, 0.18 ± 0.02 for WWT, and 0.20 ± 0.02 for PWWT. Estimates of maternal heritabilities were 0.14 for BWT, 0.10 for WWT, and 0.06 for PWWT, with SE = 0.01. Permanent environmental maternal effects explained 4 to 6% (±1%) of total phenotypic variances for these BW. Litter effects included temporary environmental effects common to littermates and a proportion of the dominance genetic variance and accounted for 16 to 19% (±1%) of phenotypic variance. Correlations between additive direct and maternal genetic effects were -0.14 for BWT, -0.23 for WWT, and -0.04 for PWWT but differed from 0 ( < 0.05) only for WWT. The total heritability predicted the total response in direct and maternal genetic effects from mass selection and was 0.23 for BWT, 0.20 for WWT, and 0.23 for PWWT. Direct and maternal additive, maternal permanent environmental, residual, and phenotypic correlations between BWT and WWT were 0.53 ± 0.05, 0.58 ± 0.06, 0.51 ± 0.06, 0.39 ± 0.01, and 0.44 ± 0.01, respectively; those between BWT and PWWT were 0.45 ± 0.06, 0.58 ± 0.08, 0.36 ± 0.08, 0.33 ± 0.01, and 0.37 ± 0.01 respectively; and those between WWT and PWWT were 0.85, 0.99, 0.92, 0.77, and 0.81, respectively, with SE ≤ 0.02. Therefore, both direct and maternal effects had an important impact on BW in Katahdin lambs. Models that included both additive and permanent environmental maternal effects as well as a temporary environmental litter effect should result in more accurate estimates of breeding values and better selection decisions.

Hair sheep of Caribbean origin have become an important part of the U.S. sheep industry. Their lack of wool eliminates a number of health concerns and drastically reduces the cost of production. More importantly, Caribbean hair sheep demonstrate robust production performance even in the presence of drug-resistant gastrointestinal nematodes, a rising concern to the industry. Despite the growing importance of hair sheep in the Americas their genetic origins have remained speculative. Prior to this report no genetic studies were able to identify a unique geographical origin of hair sheep in the New World. Our study clarifies the African and European ancestry of Caribbean hair sheep. Whole-genome structural analysis was conducted on four established breeds of hair sheep from the Caribbean region. Using breeds representing Africa and Europe we establish an objective measure indicating Caribbean hair sheep are derived from Iberian and West African origins. Caribbean hair sheep result from West African introgression into established ecotypes of Iberian descent. Genotypes from 47,750 autosomal single nucleotide polymorphism markers scored in 290 animals were used to characterize the population structure of the St. Croix, Barbados Blackbelly, Morada Nova, and Santa Ines. Principal components, admixture, and phylogenetic analyses results correlate with historical patterns of colonization and trade. These patterns support co-migration of these sheep with humans.