Genetic parameters for resistance to the Salmonella abortusovis vaccinal strain Rv6 in sheep

An experimental population (1216 lambs from 30 sires) of the Inra401 sheep was created in an Inra flock to allow QTL detection for susceptibility to Salmonella infection, wool and carcass traits. The Inra401 is a sheep composite line developed from two breeds: Berrichon du Cher and Romanov. At 113 days of age on average, the lambs were inoculated intravenously with 10⁸ Salmonella abortusovis Rv6 (vaccinal strain). They were slaughtered 10 days after the inoculation. Several traits were measured at inoculation and/or slaughtering to estimate the genetic resistance of the lambs to Salmonella infection: specific IgM and IgG1 antibody titres, body weight loss, spleen and pre-scapular node weights and counts of viable Salmonella persisting in these organs. This paper presents a quantitative analysis of the genetic variability of the traits related to salmonellosis susceptibility. The heritabilities of the traits varied between 0.10 and 0.64 (significantly different from zero). Thus, in sheep as well as in other species, the determinism of resistance to Salmonella infection is under genetic control. Moreover, the correlations between the traits are in agreement with the known immune mechanisms. The genetic variability observed should help QTL detection.     

Genetic analysis for mastitis resistance and milk somatic cell score in French Lacaune dairy sheep

Genetic analysis for mastitis resistance was studied from two data sets. Firstly, risk factors for different mastitis traits, i.e. culling due to clinical or chronic mastitis and subclinical mastitis predicted from somatic cell count (SCC), were explored using data from 957 first lactation Lacaune ewes of an experimental INRA flock composed of two divergent lines for milk yield. Secondly, genetic parameters for SCC were estimated from 5 272 first lactation Lacaune ewes recorded among 38 flocks, using an animal model. In the experimental flock, the frequency of culling due to clinical mastitis (5%) was lower than that of subclinical mastitis (10%) predicted from SCC. Predicted subclinical mastitis was unfavourably associated with the milk yield level. Such an antagonism was not detected for clinical mastitis, which could result, to some extent, from its low frequency or from the limited amount of data. In practice, however, selection for mastitis resistance could be limited in a first approach to selection against subclinical mastitis using SCC. The heritability estimate of SCC was 0.15 for the lactation mean trait and varied from 0.04 to 0.12 from the first to the fifth test-day. The genetic correlation between lactation SCC and milk yield was slightly positive (0.15) but showed a strong evolution during lactation, i.e. from favourable (-0.48) to antagonistic (0.27). On a lactation basis, our results suggest that selection for mastitis resistance based on SCC is feasible. Patterns for genetic parameters within first lactation, however, require further confirmation and investigation.     

Genetic parameters for litter size in sheep: natural versus hormone-induced oestrus

The litter size in Suffolk and Texel-sheep was analysed using REML and Bayesian methods. Litters born after hormonal induced oestrus and after natural oestrus were treated as different traits in order to estimate the genetic correlation between the traits. Explanatory variables were the age of the ewe at lambing, period of lambing, a year*flock-effect, a permanent environmental effect associated with the ewe, and the additive genetic effect. The heritability estimates for litter size ranged from 0.06 to 0.13 using REML in bi-variate linear models. Transformation of the estimates to the underlying scale resulted in heritability estimates from 0.12 to 0.17. Posterior means of the heritability of litter size in the Bayesian approach with bi-variate threshold models varied from 0.05 to 0.18. REML estimates of the genetic correlations between the two types of litter size ranged from 0.57 to 0.64 in the Suffolk and from 0.75 to 0.81 in the Texel. The posterior means of the genetic correlation (Bayesian analysis) were 0.40 and 0.44 for the Suffolk and 0.56 and 0.75 for the Texel in the sire and animal model respectively. A bivariate threshold model seems appropriate for the genetic evaluation of prolificacy in the breeds concerned.     

Quantitative trait loci for resistance to infection in sheep using a live Salmonella Abortusovis vaccine

Quantitative trait loci (QTL) mapping for susceptibility to a Salmonella Abortusovis vaccinal strain was performed using an experimental design involving 30 Romane sheep sire families (1216 progenies). Nine QTL corresponding to bacterial load, weight variations and antibody response criteria were mapped on eight chromosomes, including the major histocompatibility complex area on chromosome 20. Surprisingly, none was found to be significant in the SLC11A1 region (formerly NRAMP1) that has been shown to influence Salmonella susceptibility in other species.     

Quantitative trait loci for resistance to trichostrongylid infection in Spanish Churra sheep

Background

For ruminants reared on grazing systems, gastrointestinal nematode (GIN) parasite infections represent the class of diseases with the greatest impact on animal health and productivity. Among the many possible strategies for controlling GIN infection, the enhancement of host resistance through the selection of resistant animals has been suggested by many authors. Because of the difficulty of routinely collecting phenotypic indicators of parasite resistance, information derived from molecular markers may be used to improve the efficiency of classical genetic breeding.

Methods

A total of 181 microsatellite markers evenly distributed along the 26 sheep autosomes were used in a genome scan analysis performed in a commercial population of Spanish Churra sheep to detect chromosomal regions associated with parasite resistance. Following a daughter design, we analysed 322 ewes distributed in eight half-sib families. The phenotypes studied included two faecal egg counts (LFEC0 and LFEC1), anti-Teladorsagia circumcincta LIV IgA levels (IgA) and serum pepsinogen levels (Peps).

Results

The regression analysis revealed one QTL at the 5% genome-wise significance level on chromosome 6 for LFEC1 within the marker interval BM4621-CSN3. This QTL was found to be segregating in three out of the eight families analysed. Four other QTL were identified at the 5% chromosome-wise level on chromosomes 1, 10 and 14. Three of these QTL influenced faecal egg count, and the other one had an effect on IgA levels.

Conclusion

This study has successfully identified segregating QTL for parasite resistance traits in a commercial population. For some of the QTL detected, we have identified interesting coincidences with QTL previously reported in sheep, although most of those studies have been focused on young animals. Some of these coincidences might indicate that some common underlying loci affect parasite resistance traits in different sheep breeds. The identification of new QTL may suggest the existence of complex host-parasite relationships that have unique features depending on the host-parasite combination, perhaps due to the different mechanisms underlying resistance in adult sheep (hypersensitivity reactions) and lambs (immunity). The most significant QTL identified on chromosome 6 for LFEC₁ may be the target for future fine-mapping research efforts.     

Genetic variations in blood proteins within and between Spanish dairy sheep breeds

Eleven flocks of three Spanish sheep breeds have been studied: four Churra, four Lacha and three Manchega. Genetic variations between flocks of the same and different breeds have been analysed, using eight genetic blood systems and a wide range of statistical methods. The results show that the Churra, Lacha and Manchega breeds present significant differences between each other. Significant intraracial variations have been observed in Churras and Lachas, although they were almost non-existent in Manchegas. In some cases the genetic distances between flocks of the same breed were greater than the distances between flocks of different breeds. The data indicate that the closest breeds genetically are Lacha and Manchega.     

Environmental and genetic variation factors of artificial insemination success in French dairy sheep

Artificial inseminations (n = 678 168) recorded during 5 years in five French artificial insemination (AI) centres (2 'Lacaune', 1 'Manech tête rousse', 1 'Manech tête noire' and 1 'Basco béarnaise') were analysed to determine environmental and genetic factors affecting the insemination results. Analyses within centre-breed were performed using a linear model, which jointly estimates male and female fertility. This model combined four categories of data: the environmental effects related to the female, those related to the male, the non-sex-specific effects and finally the pedigree data of these males and females. After selection, the environmental female effects considered were age, synchronisation (0/1) on the previous year, total number of synchronisations during the female reproductive life, time interval between previous lambing and insemination, already dry or still lactating (0/1) when inseminated, and milk quantity produced during the previous year expressed as quartiles intra herd * year. The environmental male effects were motility and concentration of the semen. The non-sex-specific effects were the inseminator, the interaction herd * year nested within the inseminator, considered as random effects and the interaction year * season considered as a fixed effect. The main variation factors of AI success were relative to non-sex-specific effects and to female effects. Heritability estimates varied from 0.001 to 0.005 for male fertility and from 0.040 to 0.078 for female fertility. Repeatability estimates varied from 0.007 to 0.015 for male fertility and from 0.104 to 0.136 for female fertility. These parameters indicate that genetic improvement of AI results through a classical polygenic selection would be difficult. Moreover, in spite of the large quantity of variation factors fitted by the joint model, a very large residual variance remained unexplained.     

Response to helminth infection of sheep selected for resistance to Haemonchus contortus

Lines of Merino sheep selected for increased (IRH) and decreased (DRH) resistance to Haemonchus contortus were compared with an unselected (CH) line, after approximately four generations of selection. Measurements were recorded on 69 IRH, 47 DRH and 84 CH animals. Following artificial challenge with H. contortus, the IRH line had significantly (P < 0.001) lower faecal egg counts than the CH and DRH lines (2730,12,720 and 17,400 epg, respectively). Significant differences (P<0.05) were found between all lines in the minimum packed cell volumes during artificial infection (25.7, 22.0 and 20.3%) and in faecal egg counts after natural infection (140, 3590 and 8750 epg). Differences were also recorded (P<0.05) following artificial challenge with Trichostrongylus colubriformis (490, 840 and 1340 epg). On a percentage basis, faecal egg counts in the IRH line deviated less from the CH line following artificial infection with T. colubriformis (42%) than with H. contortus (79%). The reverse was true for the DRH line (60 and 37%, respectively). Differences in egg output of this magnitude should have a marked effect on requirements for anthelmintic treatment, rate of development of drug resistance and level of pasture contamination when the lines are grazed separately.     

Genetic parameters of milk fatty acid profile in sheep: comparison between gas chromatographic measurements and Fourier-transform IR spectroscopy predictions

Fatty acid (FA) composition is a key component of sheep milk nutritional quality. However, breeding for FA composition in dairy sheep is hampered by the logistic and phenotyping costs. This study was aimed at estimating genetic parameters for sheep milk FA and to test the feasibility of their routine measurement by using Fourier-transform IR (FTIR) spectroscopy. Milk FA composition of 989 Sarda ewes farmed in 48 flocks was measured by gas chromatography (FA_GC). Moreover, FTIR spectrum was collected for each sample, and it was used to predict FA composition (FA_FTIR) by partial least squares regression: 700 ewes were used for estimating model parameters, whereas the remaining 289 animals were used to validate the model. One hundred replicates were performed by randomly assigning animals to estimation and validation data set, respectively. Variance components for both measured and predicted traits were estimated with an animal model that included the fixed effects of parity, days in milking interval, lambing month, province, altitude of flock location, the random effects of flock-test-date and animal genetic additive. Genetic correlations among FA_GC, and between corresponding FA_GC and FA_FTIR were estimated by bivariate analysis. Coefficients of determination between FA_GC and FA_FTIR ranged from moderate (about 0.50 for odd- and branched-chain FA) to high (about 0.90 for de novo FA) values. Low-to-moderate heritabilities were observed for individual FA (ranging from 0.01 to 0.47). The largest value was observed for GC measured C16:0. Low–to-moderate heritabilities were estimated for FA groups. In most of cases, heritabilites were slightly larger for FA_GC than FA_FTIR. Estimates of genetic correlations among FA_GC showed a large variability in magnitude and sign. The genetic correlation between FA_FTIR and FA_GC was higher than 60% for all investigated traits. Results of the present study confirm the existence of genetic variability of the FA composition in sheep and suggest the feasibility of using FA_FTIR as proxies for these traits in large scale breeding programs.     

Identification of quantitative trait loci affecting resistance to gastrointestinal parasites in a double backcross population of Red Maasai and Dorper sheep

A genome‐wide scan for quantitative trait loci (QTL) affecting gastrointestinal nematode resistance in sheep was completed using a double backcross population derived from Red Maasai and Dorper ewes bred to F₁ rams. This design provided an opportunity to map potentially unique genetic variation associated with a parasite‐tolerant breed like Red Maasai, a breed developed to survive East African grazing conditions. Parasite indicator phenotypes (blood packed cell volume – PCV and faecal egg count – FEC) were collected on a weekly basis from 1064 lambs during a single 3‐month post‐weaning grazing challenge on infected pastures. The averages of last measurements for FEC (AVFEC) and PCV (AVPCV), along with decline in PCV from challenge start to end (PCVD), were used to select lambs (N = 371) for genotyping that represented the tails (10% threshold) of the phenotypic distributions. Marker genotypes for 172 microsatellite loci covering 25 of 26 autosomes (1560.7 cm ) were scored and corrected by Genoprob prior to qxpak analysis that included Box–Cox transformed AVFEC and arcsine transformed PCV statistics. Significant QTL for AVFEC and AVPCV were detected on four chromosomes, and this included a novel AVFEC QTL on chromosome 6 that would have remained undetected without Box–Cox transformation methods. The most significant P‐values for AVFEC, AVPCV and PCVD overlapped the same marker interval on chromosome 22, suggesting the potential for a single causative mutation, which remains unknown. In all cases, the favourable QTL allele was always contributed from Red Maasai, providing support for the idea that future marker‐assisted selection for genetic improvement of production in East Africa will rely on markers in linkage disequilibrium with these QTL.     

Phenotypical characteristics of Shiga-like toxin Escherichia coli isolated from sheep dairy products

AIMS: To analyse phenotypical characteristics of Shiga toxin-producing Escherichia coli (STEC) strains from ovine origin. METHODS AND RESULTS: A total of 13 STEC strains (eight O157 and five non-O157) isolated from sheep dairy products were used in this study. Biochemical traits, motility, haemolytic activity, resistance to tellurite-cefixime, maximum growth temperature and antibiotic resistance were determined. The STEC strains were grouped into nine biochemical and physiological biotypes (five for the O157 and four for the non-O157 strains). All STEC strains showed resistance to bacitracin, cloxacilin, penicillin and tylosin. CONCLUSIONS: Different biotypes and antibiotic resistance patterns of STEC isolated from sheep dairy products were observed. SIGNIFICANCE AND IMPACT OF THE STUDY: This work will be a contribution to the better characterization of STEC isolated from sheep dairy products, which have, to date, been scarcely studied, and to the better understanding of the risks associated with its consumption.     

Comparative physiological responses of stress in sheep selected for resistance or susceptibility to Haemonchus contortus in semi-arid tropical region during different seasons

The present study was conducted to assess the comparative physiological responses of stress in sheep selected for resistance or susceptibility to Haemonchus contortus in the semi-arid region during different seasons. A total of 20 adult Malpura sheep from resistant and susceptibility line were selected for study. The intensity of strongyle infection, blood biochemical and cortisol levels were estimated monthly basis while meteorological data were collected everyday throughout the year. The temperature-humidity index was the highest during dry-hot season. The faecal egg count increased significantly (p < 0.05) during June–August (rainfall season) in susceptible sheep as compared to resistant sheep. The plasma cortisol level increased significantly (p < 0.05) in both the groups during dry hot (March–May) season and it remained higher in the susceptible line till June–August. The total plasma protein reduced significantly during hot-humid with scantly rainfall (September–November) season as compared to other season in susceptible line. It is pertinent to conclude that though sheep remained in stress during dry hot (March–May) season, but effects of stress on the intensity of strongyle infection and blood biochemical were better tolerated by animals of resistant line during hot-humid season compared to susceptible line under semi-arid region.     

Genetic parameters and response to selection for growth in tambaqui

Although the tambaqui (Colossoma macropomum) is the most cultivated native fish species in Brazil, estimated breeding values for growth traits are rarely used for selection of superior individuals in commercial fingerling production. This study aimed to estimate the (co)variance components of growth traits. Body weight, length and width of 2500 tambaqui were determined at tagging and at 6 and 12 months after tagging in a commercial breeding programme in Brazil. Heritability estimates were low for traits measured at tagging (0.10 to 0.19) and moderate to high for traits measured at 6 and 12 months (0.23 to 0.81). Common full-sib effects were high at tagging (>73%), low at 6 months and negligible at 12 months. Positive genetic correlations were found among growth traits at 12 months (0.84 to 0.99) and between growth traits at 6 and 12 months (0.80 to 0.92). These results show that animal selection can be performed at 6 months after tagging. Expected genetic gains for growth traits ranged from 8% to 31%. A simulation of the sex ratio was performed, as individuals did not reach sexual maturity during the experimental period. Because of the sexual dimorphism, more accurate heritability estimates were obtained when considering the female proportion to be 90% in the high-weight group. The findings indicate that it is possible to obtain considerable genetic gains in growth by selecting for growth traits. The development of a tool to determine the sex of animals at early stages can improve the response to selection in tambaqui.     

Genetics of milking characteristics in dairy cows

Genetic selection for milking speed is feasible. The existence of a correlation structure between milking speed and milk yield, however, necessitates a selection strategy to increase milking speed with no repercussion on genetic merit for milk yield. Residual milking duration (RMD) and residual milking duration including somatic cell score (RMDS), defined as the residuals from a regression model of milking duration on milk yield or milk yield plus somatic cell score (SCS) have been advocated. The objective of this study was to undertake a first ever genetic analysis of these novel traits. Data on electronically recorded milking duration and other milking characteristics from 235 005 test-day records on 74 608 cows in 1075 Irish dairy herds were available. Variance components for the milking characteristic traits were estimated using animal linear mixed models and covariances with other performance traits, including udder-related type traits, were estimated using sire models. The heritability of milking duration, RMD and RMDS was 0.20, 0.22 and 0.18, respectively. There were little differences in the heritability of RMD or RMDS when defined using genetic regression. The genetic standard deviation of RMDS defined on the phenotypic or genetic level was 36.8 s and 37.6 s, respectively, clearly indicating considerable exploitable genetic variation in milking duration independent of both milk yield and SCS. The genetic correlation between phenotypically derived RMDS and milk yield was favourable (−0.43), but RMDS was unfavourably genetically correlated with SCS (−0.30); the genetic correlations with both traits when RMDS was defined at a genetic level were zero. RMDS defined at the phenotypic level was negatively (i.e. unfavourable) genetically correlated (−0.35; s.e. = 0.15) with mastitis; however, when defined using genetic regression, shorter RMDS was not associated with greater expected incidence of mastitis. RMDS, defined at the genetic level, is a useful heritable trait with ample genetic variation for inclusion in a national breeding strategy without influencing genetic gain in either milk yield or udder health.     

湖羊、同羊12个同工酶座位的检测

Gene frequencies of Hu sheep and Tong sheep were obtained with “Random sampling in typical colonies of a central area“.of the 12 loci tested in Hu shpee.11 loci were polymorphic.Reliability of the estimated frequencies of 27 alleles reached 0.95 except for Po^F,Tf^A,Tfd,Hb-β^A and CAF which had reliabilities of 0.5222,0.7478,0.5222,0.6212 and 0.899,respectively,Of the 12 loci tested in Tong sheep,11 loci were polymorphic.Reliability of the estimated frequency of 25 alleles reached 0.95 except for Tf^A,Tf^E and CAF which had reliabilities of 0.931,0.6922 and 0.7924,respectively.The average heterozygosity(H) and average homozygosity(J) was computed and the J of the two sheep colonies was 0.6619 and 0.6448,respectively.Consistent with our conclusions based on genetic data,previous research divided the native sheep populations of East and South Central Asia into three group:the “mongolian group“ “south-Asian group“and “European group“ .Consequently,the degree of genetic similarity between populations and known groups would seem to provide a reliable means of determining the genetic relationships between populations and may reflect the true genetic origin of Hu sheep and Tong sheep in China.     

Exploring the assumptions underlying genetic variation in host nematode resistance (Open Access publication)

The wide range of genetic parameter estimates for production traits and nematode resistance in sheep obtained from field studies gives rise to much speculation. Using a mathematical model describing host – parasite interactions in a genetically heterogeneous lamb population, we investigated the consequence of: (i) genetic relationships between underlying growth and immunological traits on estimated genetic parameters for performance and nematode resistance, and (ii) alterations in resource allocation on these parameter estimates. Altering genetic correlations between underlying growth and immunological traits had large impacts on estimated genetic parameters for production and resistance traits. Extreme parameter values observed from field studies could only be reproduced by assuming genetic relationships between the underlying input traits. Altering preferences in the resource allocation had less pronounced effects on the genetic parameters for the same traits. Effects were stronger when allocation shifted towards growth, in which case worm burden and faecal egg counts increased and genetic correlations between these resistance traits and body weight became stronger. Our study has implications for the biological interpretation of field data, and for the prediction of selection response from breeding for nematode resistance. It demonstrates the profound impact that moderate levels of pleiotropy and linkage may have on observed genetic parameters, and hence on outcomes of selection for nematode resistance.     

Using the choice experiment method in the design of breeding goals in dairy sheep

Market failures are the main cause of poor acknowledgement of the true impact of functional sheep traits on the management and economic performance of farms, which results in their omission from the breeding goal or the estimation of non-representative economic weights in the breeding goal. Consequently, stated-preference non-market valuation techniques, which recently emerged to mitigate these problems, are necessary to estimate economic weights for functional traits. The purpose of this paper is to present an example of the use of a choice experiment (CE) in the estimation of economic weights for sheep traits for the design of breeding goals. Through a questionnaire survey the preferences of sheep farmers are recorded and their marginal willingness to pay (MWTP) for 10 production and functional traits is estimated. Data are analysed using random parameter logit models. The results reveal unobserved preference heterogeneity for fertility, adaptability to grazing and resistance to disease, thus highlighting that these traits are appreciated differently by farmers, because their needs are diverse. Positive MWTP is found for Greek breeds, high milk production and lambs with low fat deposition, for which there is high demand in Greek markets. On the other hand, MWTP for the cheese-making ability of milk is negative, stemming from the fact that sheep milk prices in Greece are not formulated according to milk composition. In addition, farmers seem to understand differences between udder shapes and attribute different values to various types. This application of the CE method indicates that communication channels among farmers and breeders should be established in order to enhance market performance and to provide orientation to the design of breeding programmes. Non-market valuation can be used complementarily to market valuation techniques, in order to provide accurate estimates for production and functional traits.     

Genetic parameters for growth, reproductive and maternal traits in a multibreed meat sheep population

The genetic parameters for growth, reproductive and maternal traits in a multibreed meat sheep population were estimated by applying the Average Information Restricted Maximum Likelihood method to an animal model. Data from a flock supported by the Programa de Melhoramento Genético de Caprinos e Ovinos de Corte (GENECOC) were used. The traits studied included birth weight (BW), weaning weight (WW), slaughter weight (SW), yearling weight (YW), weight gain from birth to weaning (GBW), weight gain from weaning to slaughter (GWS), weight gain from weaning to yearling (GWY), age at first lambing (AFL), lambing interval (LI), gestation length (GL), lambing date (LD - number of days between the start of breeding season and lambing), litter weight at birth (LWB) and litter weight at weaning (LWW). The direct heritabilities were 0.35, 0.81, 0.65, 0.49, 0.20, 0.15 and 0.39 for BW, WW, SW, YW, GBW, GWS and GWY, respectively, and 0.04, 0.06, 0.10, 0.05, 0.15 and 0.11 for AFL, LI, GL, LD, LWB and LWW, respectively. Positive genetic correlations were observed among body weights. In contrast, there was a negative genetic correlation between GBW and GWS (-0.49) and GBW and GWY (-0.56). Positive genetic correlations were observed between AFL and LI, LI and GL, and LWB and LWW. These results indicate a strong maternal influence in this herd and the presence of sufficient genetic variation to allow mass selection for growth traits. Additive effects were of little importance for reproductive traits, and other strategies are necessary to improve the performance of these animals.     

Genetic mapping of quantitative trait loci for resistance to Haemonchus contortus in sheep

This paper presents results from a mapping experiment to detect quantitative trait loci (QTL) for resistance to Haemonchus contortus infestation in merino sheep. The primary trait analysed was faecal worm egg count in response to artificial challenge at 6 months of age. In the first stage of the experiment, whole genome linkage analysis was used for broad-scale mapping. The animal resource used was a designed flock comprising 571 individuals from four half-sib families. The average marker spacing was about 20 cM. For the primary trait, 11 QTL (as chromosomal/family combinations) were significant at the 5% chromosome-wide level, with allelic substitution effects of between 0.19 and 0.38 phenotypic standard deviation units. In general, these QTL did not have a significant effect on faecal worm egg count recorded at 13 months of age. In the second stage of the experiment, three promising regions (located on chromosomes 1, 3 and 4) were fine-mapped. This involved typing more closely spaced markers on individuals from the designed flock as well as an additional 495 individuals selected from a related population with a deeper pedigree. Analysis was performed using a linkage disequilibrium–linkage approach, under additive, dominant and multiple QTL models. Of these, the multiple QTL model resulted in the most refined QTL positions, with resolutions of <10 cM achieved for two regions. Because of the moderate size of effect of the QTL, and the apparent age and/or immune status specificity of the QTL, it is suggested that a panel of QTL will be required for significant genetic gains to be achieved within industry via marker-assisted selection.     

Genetic structure and admixture in sheep from terminal breeds in the United States

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright’s F_ST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright’s F_ST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture , and based on cross-validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.