J. Anim Sci.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Pollott, G. E.
Right arrow Articles by Greeff, J. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Pollott, G. E.
Right arrow Articles by Greeff, J. C.
J. Anim. Sci. 2004. 82:2840-2851
© 2004 American Society of Animal Science


ANIMAL GENETICS

Genotype x environment interactions and genetic parameters for fecal egg count and production traits of Merino sheep1

G. E. Pollott2 and J. C. Greeff

Great Southern Agricultural Research Institute, Katanning, Western Australia 6317, Australia


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Breeding for host resistance to parasites has become an imperative in many sheep industries. Because of the widespread use of AI in sheep breeding schemes, it is important to understand how the performance of offspring from rams varies in different flock environments, both for resistance to parasites and key production traits. This study used both variance component and reaction norm models to investigate the level of genotype x environment interaction for fecal egg count (FEC) and important Merino production traits in a range of flock environments in Australia. These flocks were linked by the use of common rams in a sire-referencing scheme. Both linear and quadratic polynomial reaction norm models were used. The heritability of these traits and the genetic correlation between them and FEC also was investigated using the reaction norm model. A contemporary group (CG) was defined by a flock, year, age class, sex, and paddock combination. Each CG environment was characterized by the mean value of any given trait for that CG. The recorded data used in the study were analyzed in a standardized form. Standardization for each trait was achieved within a CG by subtracting the CG mean from each observation and dividing by the CG SD. The genotype x environment effect accounted for < 0.05 of the phenotypic variance for all traits. In most traits the heritability varied little across environments. The exceptions were FEC, BW, and both greasy and clean fleece weights, which had a higher heritability at the lower end of the environmental range. Fecal egg count also had a higher heritability in high-FEC environments. Genetic correlations between FEC and several key production traits were similar in the flock environments studied. Quadratic polynomial models and models with a variable residual fitted the data better than linear models. The genotype x environment effect for FEC and the genetic correlations between FEC and production traits were effectively zero; thus, sheep breeding programs for increased parasite resistance can be run effectively by ignoring these factors. Some account should be taken of the high heritabilities of FEC and fleece and BW in different flock environments.

Key Words: Fecal Egg Count • Gastrointestinal Parasites • Genetic Parameters • Genotype x Environment Interaction • Reaction Norm • Sheep


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Genotype x environment interaction is a term used to describe the phenomenon that occurs when a set of genotypes change their relative performance in different environments (Falconer, 1981Go). The most commonly described situation is where two breeds are kept in two environments and change their ranking for a given trait in those two environments. However, several studies have extended this concept to consider how a group of genotypes varies in several environments. The estimation of a variance component to describe the additional variance in the population due to the genotype x environment effect is one approach (Hickman and Henderson, 1955Go; Meyer, 1987Go). Recent developments in the use of random regression models have led to the suggestion that the genotype x environment effect can be systematically investigated by regressing sire breeding values on some measure of the environments in which their offspring are found (Kolmodin et al., 2002Go). This paper describes an investigation of the genotype x environment effect for fecal egg count (FEC) and production traits in Australian Merino sheep using a range of methods. The use of the random regression approach to estimate genotype x environment effects also facilitated the calculation of both heritabilities and genetic correlations as they varied across a range of environments.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Flock Data

The data were collected between 1988 and 2002 from 55 commercial Merino flocks located in the four main sheep producing states of Australia: Western Australia, New South Wales, Victoria, and South Australia. They were derived by combining databases from four major sources of records: Merino Genetic Services (Meat and Livestock Australia); Merino Benchmark (New South Wales Agriculture); and Sire Evaluation Test data from New England and Yardstick (Western Australia). Flocks were selected from all four sources on the basis of having a high proportion of their animals recorded for FEC over recent years. All available animal records were used from these "FEC flocks" regardless of whether the individual animal had a FEC record or not. The final amalgamated database comprised records from over 250,000 animals, of which 127,723 had recorded sire information.

The distribution of gastrointestinal worm species is commonly described with reference to the annual rainfall pattern in different parts of Australia. In the winter rainfall regions (Western Australia, South Australia, and parts of Victoria), the dominant species are Trichostrongylus spp. and Ostertagia spp., with Nematodirus spp., Oesophagostomum venulosum, and Chabertia ovina of lesser importance. In the summer rainfall regions (Queensland and New South Wales), Haemonchus contortus (in summer) and Trichostrongylus spp. (in winter) are the main species found, with Oesophagostomum and Nematodirus of secondary importance (Cole, 1986Go). In uniform rainfall areas, species prevalence varies on an annual basis (parts of Victoria, Tasmania). The exact parasite status of the contemporary groups used in this study was largely unknown, apart from the mean FEC level calculated from the data.

The traits analyzed in this study are described in Table 1Go. Fecal egg count was adjusted using the cube-root transformation (Woolaston and Piper, 1996Go). Standardized values for all traits were used in these analyses. There were two reasons for transforming the data in this way. First, heterogeneous variances in contemporary groups have resulted in spurious genotype x environment interactions (Maniatis and Pollott, 2002Go). This possibility was removed by standardizing the data in such a way as to make within-contemporary group variances have a value of one. Second, the use of random regression methods to model the way sire variance changed across different environments was found to give unrealistic values for heritability with untransformed data. Standardized values were calculated for each trait by adjusting the phenotypic value using the relevant contemporary group (CG) mean and standard deviation for that trait. For example, FECstd = (FEC –MeanFECCG)/SDFEC CG, where MeanFECCG and SDFECCG refer to the mean and standard deviation of the CG of the individual animal whose data is being standardized for FEC (Hill, 1984Go). This was repeated for each trait using its own CG means and standard deviations to standardize that trait.


View this table:
[in this window]
[in a new window]
 
Table 1. The traits analyzed in this study and the categorical descriptors used to characterize the traits
 
For any particular analysis, records with complete trait, covariate, and categorical information (fixed effects) were selected from the database. The number of animals recorded for each trait and categorical effect was varied. Consequently, the number of animals used in the various analyses differed considerably. Age class of the recorded animals was incorporated in the contemporary group fixed effect and records from all ages were treated as if they were from the same trait.

It was possible to explore the data for genotype xenvironment effects because the use of reference sires linked flocks by common sires, and so the offspring of certain sires were born in a wide range of flocks. Table 2Go shows the distribution of sires across flocks. Of the 2,186 sires in the dataset, 2,004 were only used in one flock and the remaining 182 were used in several flocks. Two rams were used very widely and 35 rams were used in more than three flocks. However, CG (flock, year, sex, age class, paddock) was the main grouping of animals used in these analyses. Linkage between CG varied according to the trait being analyzed. Of the 291 CG in the FEC dataset, all used at least two rams to mate their ewes and only 28 sires were used in a single CG. Thus, there would be a reasonable data structure with which to estimate any genotype xenvironment effect if all traits were recorded on all animals. However, this was clearly not the case (Table 3Go), and so the data structure for each trait varied. The number of CG for each trait is shown in Table 3Go and ranged from 43 (staple strength) to 593 (fiber diameter). Mean CG size ranged between 78 animals for loin muscle depth to 185 for staple strength.


View this table:
[in this window]
[in a new window]
 
Table 2. The distribution of sires between flocks and fecal egg count contemporary groups (FEC CG)
 

View this table:
[in this window]
[in a new window]
 
Table 3. Number of records, means, and fixed effect models used to analyze the 12 standardized traits
 
Basic Statistical Model

A sire model (Model 1) was fitted to all standardized traits using ASREML (VSN Int., Hemel Hempstead, U.K.), and variances were considered to have been estimated when the log likelihood in successive iterations differed by less than 0.002x the current iteration number.


[Model 1]

In Model 1, Yijkl was a record from one of the 12 standarized traits analyzed for the lth animal, in the ith birth type, the jth rearing type, and from the kth sire. Age and Age2 were fitted as covariates. The sire and residual effects were fitted as random variables. The number of animals, CG, and sires varied for each trait, depending on the level of recording of the trait, fixed effects, and covariates.

Genotype xEnvironment Effect

Random Effects for Sire xContemporary Group. The first approach to estimating the genotype x environment effect for FEC was to extend Model 1 by fitting sire x CG as an additional random effect. Mixed models were run for all standardized traits using the final fixed effect structure derived from fitting Model 1 (Table 3Go), with sire, residual and sire x contemporary group as random effects. The likelihood ratio test was used to test for significant additions/subtractions to the model against a {chi}2 value, the number of degrees of freedom being the difference in the number of sets of effects in the two models being compared. A genotype x environment variance ratio (ge2) was calculated from the sire x CG variance as a proportion of the total phenotypic variance.

Reaction Norm Approach. The second approach to estimating genotype x environment effects used the so-called "reaction norm" approach (Kolmodin et al., 2002Go) by fitting random regression terms for each sire (Model 2). The CG average for the trait was used as the independent variable in all random regression models. This created an environmental gradient across the range of environments in the dataset, in this case each contemporary group being defined as a different environment. The environmental gradient was redefined for each trait according to the CG mean for the trait being analyzed. An example random regression model, for FEC, is shown in Model 2.


[Model 2]

In Model 2 Yij = the jth animal’s recorded standardized FEC value, Sai = the additive effect of the ith sire and Sbi the regression coefficient of the ith sire’s breeding value on Xk, the kth CG mean FEC value. The variance of the sire "intercept" (Sai), slope of a linear regression of sire CG values on CG mean (Sbi), and the covariance between them (covSab) were estimated as random effects. For other traits a different set of fixed effects were used, depending on which effects were found to significantly affect the trait from previous analyses, and the Xk variable was the CG mean for the trait being analyzed from the kth CG.

Model 2 estimated the sire variance as the linear regression of sire breeding value in a specific environment on the mean CG value for that environment. Further analyses were carried out using a quadratic polynomial in place of the linear random regression for FEC (Model 3) to allow some flexibility in the sire terms. This was achieved by adding the term Sci(2) to the model, plus its covariances with the other two random regression terms. The sire x CG term was added to the random part of Model 3 to allow an estimate of how much of the Sx CG variation was accounted for by the systematic effects described by the random regressions (Model 4).


[Model 3]


[Model 4]

Analyses using these models were repeated for each trait using standardized data with the appropriate fixed effect model for the trait.

Heritability of a Trait in Different Environments

The use of reaction norm models, using random regression across a range of environments, allows the estimation of the way in which the heritability of a trait varies with the environment. As described above, for each trait, the CG mean for that trait was used as the environmental variable. The method described by Kolmodin et al. (2002)Go was used to estimate the sire variance ({sigma}2s) of a trait at different values of the environmental variable. In a general matrix form this was:


where V was a 1 x n vector of the environmental variable value (X) at which the sire variance was being estimated (1, X, X2, ..., Xn), n was the order of polynomial being fitted, and G was the n x n variance/covariance matrix of the random regression terms being fitted in the model. The additive genetic variance, at any given X value, was calculated as four times the sire variance and the phenotypic variance ({sigma}2p) calculated as the residual plus the sire variance. Standard errors of the heritability were calculated as described by Gilmour et al. (2002)Go for each value of X.

Genetic Correlations across Contemporary Groups with Different Mean Fecal Egg Count Values

Genetic correlations between standardized FEC and six standardized production traits were calculated across the range of FEC environments using the reaction norm approach described by Kolmodin et al. (2002)Go. The six traits were chosen for a number of reasons. Fiber diameter and greasy fleece weight are widely recorded and important traits in Merino flocks but had a zero genetic correlation with FEC. Body weight was a widely recorded trait, whereas fat depth, loin muscle depth, and staple strength all had point estimate genetic correlations different from zero. The sire covariance between FEC and each production trait was calculated using VGV', as described above, but in this case the G matrix was the variance/covariance matrix of all possible combinations of the polynomial terms fitted to both traits. Standard errors were calculated as described by Gilmour et al. (2002)Go. In all cases standardized data were used for both traits and a quadratic polynomial was fitted for each trait.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Genotype x Environment Interactions

A summary of the analyses adding sire x CG to the basic model (Table 3Go) for each trait is shown in Table 4Go. The likelihood ratio test showed that in all traits, except staple length and fiber curvature, the addition of sire x CG resulted in a better fit of the model to the data (P < 0.01). The proportion of phenotypic variance accounted for by this additional term ranged from 0.6 (fiber curvature) to 4% (grease fleece weight and loin muscle depth). The genotype x environment effect for FEC was 3%.


View this table:
[in this window]
[in a new window]
 
Table 4. Comparison of models fitted with and without the sire x contemporary group effect to standardized data (unitless values)
 
The results of fitting linear and quadratic polynomial regressions of FEC on the environmental variable, CG average FEC, are summarized as likelihood ratio test comparisons between various models in Table 5Go. The linear random regression model reduced the log likelihood (P < 0.001) compared with the sire model (Model 1) and the polynomial random regression improved the fit of the model further (P < 0.001). There was little difference between the fit of the models containing either the sire x CG or the linear random regression effects.


View this table:
[in this window]
[in a new window]
 
Table 5. Likelihood ratio test (LRT) values comparing a range of models used to investigate genotype x environment effects in standardized fecal egg count
 
The results of fitting the linear random regression model to the 12 standardized traits are summarized in Table 6Go. Staple strength was the only trait where the use of the random regression model did not significantly decrease the log likelihood from the use of the sire model alone. Greasy fleece weight, clean fleece weight, and fiber diameter coefficient of variation showed a very high correlation between the slope of the sire random regression lines and the level of the sire breeding values. The slope of the individual sire regression lines accounted for more than 1% of the phenotypic variation in clean fleece weight, fiber diameter, fiber curvature, BW, fat depth, and loin muscle depth.


View this table:
[in this window]
[in a new window]
 
Table 6. Variances, covariances, correlations, and proportion of phenotypic variance attributed to slope variance from fitting a linear random regression model to each of the 12 standardized traits (unitless values)
 
The most accurately estimated individual sire regressions are shown in Figure 1Go from the quadratic polynomial random regression model. The rams demonstrating the greatest range of breeding values over the different environments were Ram 386 (–1.03 epg0.33) and Ram 925 (+1.08 epg0.33). Another indication of a genotype x environment effect would be a change in relative ranking of sires in different environments. This is difficult to summarize in a single statistic, but inspection of Figure 1Go clearly demonstrates that the 10 sires ranked differently over the range of the environmental variable shown. Rams 925 and 536 are good examples of this; they had identical values of FEC in CG with a mean FEC value of 5 epg0.33, but were two units of FEC different in environments with a mean FEC level of 13 epg0.33.



View larger version (25K):
[in this window]
[in a new window]
 
Figure 1. The change in breeding value of 10 rams over the range of fecal egg count (FEC) environments analyzed using a polynomial random regression model. Fecal egg count environments were characterized by mean FEC in contemporary groups (CG).

 
Heritability of Fecal Egg Count in Different Environments

The calculation of the heritability of FEC from the use of quadratic random regression model to model the way the sires’ breeding values changed across the range of FEC environments is shown in Figure 2Go. In low- and high-FEC environments heritability was higher than in the middle of the range of environments. Because of the use of standardized data, this variation in heritability was almost entirely due to the change in sire variance across the environmental range.



View larger version (13K):
[in this window]
[in a new window]
 
Figure 2. The heritability of standardized fecal egg count (FEC), standard error of heritability, and number of contemporary groups across a range of environments. Sire model heritability was 0.24 ± 0.016. Model used was the quadratic regression of sire’s breeding value on mean contemporary group (CG) fecal egg count value.

 
The results of fitting the linear random regression model (not shown) to the standardized FEC data produced a similar estimate of heritability across the range of environments.

Heritability of Production Traits in Different Environments

The heritability results for the production traits are presented as a series of graphs (Figures 3Go and 4Go) derived from using the polynomial random regression. In Figures 3Go and 4Go, the environmental variable was the CG mean for the original data of each trait but expressed as a standard deviation of the CG means. This allowed a much easier comparison between traits on a common scale. The range ± 2 SD was used since this encompasses 95% of environments and represents the range of environments over which the heritability was most accurately estimated. Figures 2GoGo to 4Go also show the heritability standard errors and the distribution of CG means across the range of environments.



View larger version (38K):
[in this window]
[in a new window]
 
Figure 3. The heritability of standardized greasy fleece weight (GFW), clean fleece weight (CFW), fiber diameter (FD), coefficient of variation of fiber diameter (FDCV), comfort factor (FDPK), and curvature (CRV), standard error of heritability, and number of contemporary groups (CG) across a range of environments. Model used was the quadratic regression of sire’s breeding value on mean CG value. Results are shown using standard deviation units of the different environmental variables to aid comparison between traits.

 


View larger version (29K):
[in this window]
[in a new window]
 
Figure 4. The heritability of standardized staple strength (SS), staple length (SL), BW, fat depth (FAT), and loin muscle depth (EMD), standard error of heritability, and the number of contemporary groups (CG) across a range of environments. Model used was the quadratic regression of sire’s breeding value on mean CG value. Results shown using standard deviation units of the different environmental variables to aid comparison between traits.

 
A range of patterns of heritability is apparent from Figures 2GoGo to 4Go. Greasy and clean fleece weight (Figure 3Go) had a similar pattern to body weight (Figure 4Go) and fiber curvature (Figure 3Go). In environments characterized by a low environmental variable, those four traits had a high heritability, which fell rapidly as the environment improved and leveled off by ± 1 SD on either side of the mean. In the best environments, their heritability declined further. The fleece weight traits and BW were estimated from a large number of CG across the range of environments, as shown by the small standard errors at all but the lowest environmental values (–2 SD). The CG means were skewed slightly toward the lower end of the scale and there were no groups with a mean less than –2 SD from the overall mean, but there were groups with up to 3 SD above the mean.

Some other traits showed very little variation in heritability in different environments (fiber diameter, fiber diameter CV, staple length). Because the phenotypic variance is about unity across the range of environments with standardized data, these results reflect a stable additive variance across the range as well. All three traits showed a normal distribution of CG means.

The remaining traits all had different patterns of heritability depending on a range of factors. The heritability of comfort factor (Figure 3Go) decreased from 0.62 (–2 SD) to 0.24 (0.5 SD). Because of the way comfort factor is calculated, it has a truncated normal distribution; no value can be higher than 100%. This is reflected in the distribution of CG means shown in Figure 3Go.

Staple strength (Figure 4Go) was recorded on the smallest number of CG, and its heritability was low (approximately 0.2) at both extremes of the environmental scale (+2 and –2 SD) and highest in the middle (0.43). The standard errors were large and indicate that this trait may have a constant heritability across the range of environments.

The distribution of fat depth CG means (Figure 4Go) was also skewed toward the lower end of the scale, and fat depth was not widely recorded. Its heritability was effectively the same across the range of values shown.

The distribution of loin muscle depth CG means (Figure 4Go) was skewed toward the lower end of the scale and ranged between –1.7 and 2.5 SD. It was not widely recorded, and the distribution of CG means above 0.25 SD was very low. This may account for the irregular pattern of loin muscle depth heritability shown in Figure 4Go, with low values at either end of the scale and a rise on the positive side of the scale.

Genetic Correlations

The pattern of genetic correlations between FEC and six production traits are summarized in Figure 5Go. In all cases, the environmental variable used was the mean CG FEC value. The two widely recorded wool traits, fiber diameter and FEC, showed very little correlation with FEC across the range of FEC environments studied; the highest correlation was < 0.1 at the highest level of FEC (fiber diameter). The two body measurements (fat depth and loin muscle depth) were found to vary with changing FEC environments; as the environment declined, the genetic correlation became more negative. This implies that under high levels of parasite challenge, genetically resistant animals will have a greater breeding value for fat and loin muscle depth than less resistant animals. Staple strength appeared to have a zero correlation with FEC across the range of environments, but this may have been due to it being recorded on a relatively small number of animals. Body weight also seemed to have a low and constant genetic correlation with FEC.



View larger version (18K):
[in this window]
[in a new window]
 
Figure 5. The genetic correlation between standardized fecal egg count and six production traits over a range of fecal egg count environments. The model used was the quadratic regression of sire’s breeding value on mean contemporary group (CG) fecal egg count value. Point estimate genetic correlations using standardized data: loin muscle depth (EMD) = –0.19; fat depth (FAT) = –0.25; fiber diameter (FD) = –0.04; greasy fleece weight (GFW) = 0.05; BW = –0.09; and staple strength (SS) = –0.20.

 
In general, traits tended to have very little correlation with FEC in low-FEC environments, but this increased in high-FEC environments, fat depth being the only exception. Figure 5Go also shows the point estimates of the genetic correlations estimated from these standardized data. All these point estimates could be found on their respective random regression lines, with the exception of staple strength, but there was no consistency about where on the environmental scale they were located. Most were below or at the mean FEC value of 8.85 epg0.33


    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
Genotype xEnvironment Interaction for Fecal Egg Count

There was clear evidence of a genotype xenvironment interaction for FEC in this Australian Merino data. This was shown in a number of ways: by the ge2 value being different from zero, and by the slope variance of the random regression model being different from zero—the ram breeding values changing in different environments and the better fit of the models, including a genotype xenvironment effect (or random regression term), compared with those without.

The general level of the effect was rather low. On the basis of the results from the current dataset, the level of the genotype xenvironment effect for FEC would seem to be of little practical value. However, there were clearly some rams that did change their breeding values in different environments, and some rams appeared to have a negative random regression slope (i.e., they had a higher breeding value in a low-FEC environment than in a high-FEC environment).

The low level of genotype xenvironment effect may have been due to a range of factors. Firstly, parasite resistance, as measured by FEC, may be due to a general immunity to gastrointestinal parasites rather than a resistance to specific parasites. Thus, rams with a high level of genetic resistance to gastrointestinal parasites, used in a range of FEC environments with a different mix of parasites or environmental factors, may confer the same level of a general immunity to their offspring in all environments. This would result in no genotype xenvironment effect being found.

Secondly, FEC measurements do not indicate which parasites are present, and so the offspring of a particular ram may have immunity to different parasites in different environments. In Australia, internal parasite species vary between regions and seasons (Cole, 1986Go). This variation was not measured in this study, but would be useful information to collect to help the interpretation of this type of data.

There is very little evidence from other studies to help interpret these results. Morris et al. (1993)Go and our own unpublished results failed to detect a genotype xenvironment effect in trials designed to investigate the use of two contrasting FEC genotypes in two distinct FEC environments. More recent studies in New Zealand, looking at the genetic correlation of FEC for Trichostrongylus and Nematodirus spp. (Morris et al., 2004Go), found the genetic correlation to be approximately 0.43. This indicates that different genes are controlling the resistance to the two types of parasites and that a genotype xenvironment effect may exist.

Genotype xEnvironment Interactions in Production Traits

The results for the wool traits indicate that fleece weight is subjected to the highest level of genotype xenvironment interaction, with fiber diameter traits having a lower level of effect and the other fiber traits generally unaffected by genotype xenvironment influences. Very few reports of genotype xenvironment effects on wool traits are available. Dominik et al. (1999)Go quoted genetic correlations between related animals kept in two specific environments and concluded that there was no genotype xenvironment effect present for a range of production traits, despite some correlations being 0.62 ± 0.20 (clean fleece weight), 0.77 ±0.11 (fiber diameter), and 0.80 ±0.13 (staple strength).

Genotype xenvironment interactions of body traits in sheep have been reported by Hagger (1998)Go and Maniatis and Pollott (2002)Go. Hagger (1998)Go analyzed pre-weaning growth rate in Swiss Black-Brown Mountain sheep and found ge2 values ranging from 4 to 6%. The results of Maniatis and Pollott (2002)Go were slightly lower at 2 to 3% for pre- and postweaning weights of Suffolk sheep in a U.K. sire-referencing scheme. The value of 2.6% reported in the current study for BW was very similar. Maniatis and Pollott (2002)Go also reported values of 2% for ultrasonically measured loin muscle depth and 2 to 4% for fat depth, which were similar to the values of 4 and 3.2% found in this study. From this limited amount of evidence, there appears to be a real but low level of genotype xenvironment effect for fleece and body traits in sheep, with less evidence for an effect on fiber characteristics.

Heritability of Traits in Different Environments

One by-product from the use of random regression models to investigate genotype xenvironment effects is that they allow the heritability of a trait to be estimated in a wide range of environments. The traditional way of estimating heritability, a point estimate for any given dataset, has become the norm, and authors tend to look for similarities between their estimates as a means of "verifying" their results. However, heritability was originally defined to relate only to the population under study in the environment in which it was found (Falconer 1981Go), and a range of values is to be expected. The use of reaction norm models makes it possible to revisit the idea of the heritability of a trait being related to a particular flock environment in a more systematic and comparable way than has been possible previously.

The heritability of FEC for this dataset was 0.23 ± 0.016. From Figure 2Go, it is apparent that this value is not that calculated by random regression at the mean FEC value (8.85 epg0.33), but may represent heritability calculated at an average weighted mean environmental value. Inspection of Figures 3Go and 4Go confirms this. In some cases, the point estimate heritability was at the mean environmental value (clean fleece weight, fiber diameter CV, staple length, fat depth), or within the standard error shown (comfort factor, staple strength, fiber curvature, BW, loin muscle depth). In other traits, this is clearly not the case (FEC, fiber diameter), but heritability tends toward the middle of the range of values found.

Figures 2GoGo to 4Go show how the heritabilities of the different traits vary in the range of environments analyzed in this study. These fall into several distinct patterns. Fleece and BW characteristics have a high heritability in environments characterized by a low value of the environmental variable, a stable heritability in the middle range of environments, and then a decreasing heritability as the environment improves at the upper extreme. The standard errors of the estimates at the extremes of the environmental range suggest that these are real differences. It is not surprising that these traits have a similar pattern because they are all related to size of the animal in some way. At low values of the environmental variable, the genetic variability of the traits was high, implying that some animals have the genetic ability to cope with poor conditions better than others. Under a range of relatively benign environments, the variability of breeding variables is stable and less than in the poor environments. Under good regimes, the genetic variability of growth and fleece traits is low and most genotypes respond similarly to the environmental conditions. It is impossible to tell from the current dataset how many, and which, genes are influencing these traits. However, it does seem logical to assume that different genes contribute to body and fleece weight under different environmental conditions.

A number of traits showed no difference in heritability across the range of environments in which they were recorded (fiber diameter CV, staple length, fat depth), or their standard errors were large at the environmental extremes and so no difference could be detected across the range, despite a nonlinear pattern (staple strength, fiber curvature, loin muscle depth). In these traits, the genetic variability remained stable in all environments, implying that the breeding values for the animals were invariant to the environmental conditions.

The remaining traits (FEC, fiber diameter, and comfort factor) demonstrated a variety of patterns for the way heritability varied across environments. Fiber diameter heritability was well estimated and only varied a small amount over the environmental range, and so it can be considered to be almost similar across the range. Comfort factor was a trait with a very skewed distribution and, because of the way it was calculated, higher variance is expected at lower trait values. Fecal egg count (Figure 2Go) was found to have a high heritability at the extremes of the environmental range, but a moderate level in the middle. This implies that in low-FEC environments, the genetic variation for FEC was high and therefore some rams have the genetic predisposition to have high egg counts even when the challenge is relatively low (i.e., they have no resistance to the parasites at all). Under a moderate challenge regimen, animals are more similar in their genetic control of FEC; presumably the animals with lower resistance in good environments increase their FEC values toward the animals with no resistance. At high challenges, some animals have the ability to resist the parasites more than others. Inspection of the small number of rams shown in Figure 1Go adds weight to this interpretation. Some rams appear to have a reducing breeding value for FEC as the environment declines.

Several authors have reported the heritability of Merino production traits in two environments using common rams as the link between environments. Lewer and Ritchie (1992)Go compared two groups of ewes kept in simulated commercial (harsh environment) and stud (less harsh environment) environments and linked by common rams. The heritability of clean fleece weight, staple strength, and fiber diameter were lower in the commercial environment, although the standard errors were large. Dominik et al. (1999)Go analyzed a dataset that included that of Lewer and Ritchie (1992)Go and found no significant difference between heritabilities for FEC, clean fleece weight, fiber diameter, staple strength, fiber diameter CV, staple strength, and BW. Similar conclusions were reached by Erasmus (1988)Go, who also compared the heritabilities of Merino production traits in two different environments. These results represent differences in heritabilities between the two environments over a very narrow range of environments compared with the ranges shown in Figures 2GoGo to 4Go. Thus, the results reported in this article and those of Dominik et al. (1999)Go and Erasmus (1988)Go are in agreement.

The variety of patterns for the way heritability varies across the range of trait environments may have some implications for selection programs. For example, selection schemes could be designed to select against the animals that have high FEC in poor environments, as well as selecting those that have a low FEC in all environments.

Genetic Correlations

The pattern of the way that genetic correlations between FEC and some production traits change as the FEC environment changes suggest that in low-FEC environments, there will be no effect on most other traits of selection for low-FEC animals, fat depth being the only exception. However, in medium- to high-FEC environments, some increase in fat depth and muscle depth can be expected.

Heritability and the Genotype x Evironment Effect

Using the reaction norm approach to describe the systematic way that a genotype may vary in different environments has also allowed for a description of the way that the heritability of a trait varies across those same environments. The results were explored to investigate the relationship between the pattern of heritability across environments and any genotype xenvironment effect.

Using Figure 1Go as an example, if no genotype xenvironment effect existed, then the regression lines of sire EBV on environmental value would be parallel. This would result in a heritability of the same value in all environments (i.e., the heritability line would be parallel to the x-axis on graphs such as Figure 2Go). Any deviation from the situation of parallel regression lines implies that a genotype xenvironment effect is present and would be reflected in a changing additive variance across the environmental range. Thus, a genotype xenvironment effect is expected for any standardized trait that does not have a constant heritability across the range of environmental variables.

Looking at Figures 2GoGo to 4Go, some traits have a constant heritability across the environmental range (fiber diameter CV, staple length) and some traits may be considered to have constant heritability due to large standard errors (fiber diameter, fiber curvature, staple strength, fat depth, loin muscle depth). The other traits have a variable heritability (fecal egg count, clean fleece weight, comfort factor, BW). If we ignore comfort factor, because it is a trait with a truncated normal distribution, then the four traits with the most variable heritability patterns have the highest ge2 values; FEC (0.030), FEC (0.039), clean fleece weight (0.031), and body weight (0.024). The traits with either a flat heritability or large standard errors, indicating the possibility of a flat heritability line, all had ge2 values ≤0.01 and were mostly not different from 0. The only exception to this was fiber diameter (ge2 = 0.024) which was well estimated and shows a small rise in heritability in the middle of the environmental range.

Relationship between Different Measures of the Genotype xEnvironment Effect

In this article, a range of measures of the genotype xenvironment effect has been used to investigate its presence in a number of Merino traits. These different measures were compared to see if any relationship between them existed. The use of ge2 estimates the variance in the trait due to deviations from the sire plus CG main effects. Thus, if a particular sire-CG combination has a value that differs from the sum of the two main effects, it will contribute to the s2 value. The variance due to the slope of the sire’s breeding value, calculated from the random regression model, shows how similarly all sires respond in the range of environments. If all sires responded in the same way to each environment, then the variance of slopes would be zero (i.e., they would all have the same slope). Thus, ge2 gives an insight into the existence of a genotype xenvironment effect, and the slope variance tells us how consistent that effect is across the environmental range; a low slope variance indicating a consistent systematic effect across all sires.


    Implications
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 
The significant but small genotype xenvironment interactions for fecal egg count and fleece and body traits should allow rams to be used widely in Merino breeding programs without a need to consider the flock environment in which they are bred. The heritability of several of these traits is invariant to the flock environments studied and can be used widely in selection programs. However, the heritability of body and fleece weights and fecal egg count varies with the flock environment, and selection programs may need to be adjusted accordingly. Breeding for increased resistance to parasites should not lead to any adverse correlated effects in other traits in any of the flock environments studied.


    Footnotes
 
1 This research was carried out while the first author was on sabbatical leave from Imperial College London, U.K. It was funded by the Sheep Cooperative Research Centre’s program while the first author was based at the Dept. of Agric. of Western Australia. All three institutions are acknowledged for their various contributions to this work. Data were supplied by Merino Genetic Services, Merino Benchmark, New England Sire Test Evaluation, and Yardstick (Western Australia), and we wish to thank them all, plus D. J. Brown, A. Casey, and A. Swan, for their assistance in getting the dataset together. Back

2 Correspondence: Dept. of Anim. Sci., Imperial College London, Wye Campus, Ashford, Kent, TN25 5AH, U.K. (phone: +44 (0)20 759 42707; fax: +44 (0)20 759 42919; e-mail: g.pollott{at}imperial.ac.uk).

Received for publication December 23, 2003. Accepted for publication June 1, 2004.


    Literature Cited
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 Literature Cited
 


Cole, V. G. 1986. Animal Health in Australia. Vol. 8. Helminth Parasites of Sheep and Cattle. Australian Agric. Health and Quarantine Serv., Dept. of Primary Industry, Canberra, Australia.

Dominik, S., B. J. Crook, and B. P. Kinghorn. 1999. Genotype xmanagement interaction on wool production traits and body weight in Western Australian merino sheep. Proc. Aust Assoc. Anim. Breed. Genet. 13:98–101.

Erasmus, G. J. 1988. A mixed model analysis of a selection experiment with merino sheep in an arid environment. Ph.D. Diss., Univ. of the Orange Free State, Bloemfontein, South Africa.

Falconer, D. S. 1981. Introduction to Quantitative Genetics. 2nd ed. Longman, London, U.K.

Gilmour, A. R., B. J. Gogel, B. R. Cullis, S. J. Welham, and R. Thompson. 2002. ASReml Users Guide Release 1.0 VSN Int. Ltd., Hemel Hempstead, U.K.

Hagger, C. 1998. Litter, permanent environmental, ram-flock and genetic effects on early weight gain of lambs. J. Anim. Sci. 76:452–457.[Abstract/Free Full Text]

Hickman, C. G., and C. R. Henderson. 1955. Components of the relationship between level of production and rate of mortality in dairy cattle. J. Dairy Sci. 38:883–890.[Abstract/Free Full Text]

Hill, W. G. 1984. On selection among groups with heterogeneous variance. Anim. Prod. 39:473–478.

Kolmodin, R., E. Strandberg, P. Madsen, J. Jensen, and H. Jorjani. 2002. Genotype by environment interaction in Nordic dairy cattle studied using reaction norms. Acta Agric. Scand. Sect. A. 52:11–24.

Lewer, R. P., and A. J. M. Ritchie. 1992. Genetics for staple strength. Pages 121–126 in Proc. Natl. Workshop Manag. Wool Quality in Mediterranean Environ. Western Australia Dept. of Agric., Perth, Australia.

Maniatis, N., and G. E. Pollott. 2002. Genotype xenvironment interactions in lamb weight and composition traits. Anim. Sci. 75:3–14.

Meyer, K. 1987. Estimates of variances due to sire xherd interactions and environmental covariances between paternal half-sibs for first lactation dairy production. Livest. Prod. Sci. 17:95–115.

Morris, C. A., S. A. Bisset, R. L. Baker, T. G. Watson, D. L. Johnson, and M. Wheeler. 1993. An investigation of sire by location interactions for fecal nematode egg counts in lambs. Proc. N. Z. Soc. Anim. Prod. 53:231–233.

Morris, C. A., S. A. Bisset, A. Vlassoff, C. J. West, and M. Wheeler. 2004. Genetic parameters of Nematodirus spp. egg counts in Romney lambs in New Zealand. Anim. Sci. 79:33–39.

Woolaston, R. R., and L. R. Piper. 1996. Selection of merino sheep for resistance to Haemonchus contortus: Genetic variation. Anim. Sci. 62:451–460.


This article has been cited by other articles:


Home page
J ANIM SCIHome page
B. Zumbach, I. Misztal, S. Tsuruta, J. P. Sanchez, M. Azain, W. Herring, J. Holl, T. Long, and M. Culbertson
Genetic components of heat stress in finishing pigs: Parameter estimation
J Anim Sci, September 1, 2008; 86(9): 2076 - 2081.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
G. D. Snowder, L. D. Van Vleck, L. V. Cundiff, G. L. Bennett, M. Koohmaraie, and M. E. Dikeman
Bovine respiratory disease in feedlot cattle: Phenotypic, environmental, and genetic correlations with growth, carcass, and longissimus muscle palatability traits
J Anim Sci, August 1, 2007; 85(8): 1885 - 1892.
[Abstract] [Full Text] [PDF]


Home page
J DAIRY SCIHome page
J. J. Windig, M. P. L. Calus, B. Beerda, and R. F. Veerkamp
Genetic Correlations Between Milk Production and Health and Fertility Depending on Herd Environment
J Dairy Sci, May 1, 2006; 89(5): 1765 - 1775.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
N. R. Adams and J. R. Briegel
Multiple effects of an additional growth hormone gene in adult sheep
J Anim Sci, August 1, 2005; 83(8): 1868 - 1874.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
G. E. Pollott, L. J. E. Karlsson, S. Eady, and J. C. Greeff
Genetic parameters for indicators of host resistance to parasites from weaning to hogget age in Merino sheep
J Anim Sci, October 1, 2004; 82(10): 2852 - 2864.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Pollott, G. E.
Right arrow Articles by Greeff, J. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Pollott, G. E.
Right arrow Articles by Greeff, J. C.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS