|
|
||||||||
ANIMAL GENETICS |



* Roslin Institute (Edinburgh), Roslin, Midlothian EH25 9PS, UK; and
and
Dept. of Clinical Veterinary Science, Division of Farm Animal Science, University of Bristol, Langford, Bristol BS40 5DU
| Abstract |
|---|
|
|
|---|
Key Words: fatty acid genetic parameter intramuscular fat quantitative trait locus sheep
| INTRODUCTION |
|---|
|
|
|---|
Sheep meat has been criticized because of high concentration of SFA, presumed to increase the risk of heart disease, and a low PUFA to SFA ratio (Enser et al., 1996
). Recently, nutritionists have focused on the types of PUFA that lower blood cholesterol concentrations, and the balance in the diet between n-3 PUFA formed from linolenic acid (18:3) and n-6 PUFA formed from linoleic acid (18:2; Williams, 2000
).
Fatty acid composition is influenced by both genetic and environmental factors. Genetic factors have not been widely investigated in sheep, and it is not known whether fatty composition can be genetically manipulated toward a more favorable profile. In this paper, we studied muscle fatty acid composition because i.m. fat cannot be easily removed and, thus, inevitably has an impact on human health.
The aims of this study were, first, to investigate the inheritance of fatty acid composition in sheep meat; second, to estimate genetic relationships between fatty acids and computer tomograph-assessed muscle density (a visual predictor of i.m. fatness); and third, to identify QTL for fatty acid composition, as a means to determine options for genetically improving the fatty acid profile of sheep meat.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Blackface lambs born at Roslin Institutes Blythbank farm, Scotland, during the period 2000 to 2003, were used as the experimental resource. These lambs were from a flock that had been divergently selected for carcass lean content from 1988 to 1996, creating LEAN and FAT lines, after which the lines were maintained as closed populations with no further selection. A description of the selection procedures and initial responses to selection is given by Bishop (1993)
.
The flock consisted of 200 ewes, split approximately equally between the LEAN and FAT lines. A small proportion of reciprocal LEAN x FAT-line crosses were made at the 1999 matings, so that a cohort of F1 lambs was born in April 2000, along with a majority of purebred LEAN- and FAT-line lambs. The male F1 lambs were then backcrossed to the purebred LEAN- and FAT-line ewes to create a population of F1 x LEAN and F1 x FAT lambs from 2001 to 2003, for the purpose of QTL detection. This design created 9 half-sib families for QTL detection. On average, families contained 33 male offspring for measurement of fatty acids traits, as described below, with a range from 12 to 46.
Standard husbandry procedures were applied in this flock; all lambs were tagged at birth, with parentage, day of birth, gender, and mortalities recorded. Each year the lambs were kept in 2 groups (i.e., on 2 separate fields) for ease of management. Parentage information was maintained for all animals born after 1986, giving a total of 4,847 known animals in the flock pedigree. Male lambs were not castrated.
Phenotypic Measurements
Phenotypic measurements of fatty acids were made on three hundred fifty 8-mo-old, male, grass-fed lambs, composed of 300 male lambs from the 9 backcross families, plus 25 LEAN- and 25 FAT-line male lambs born in 2000. Measurements were performed at the University of Bristol, on cohorts of 20 animals treated identically during their growth, transportation, and pre-slaughter periods.
Lipids were extracted from duplicate 10-g samples of LM, essentially as per Folch et al. (1957)
, separated into neutral and phospholipid, saponified and methylated, and individual fatty acids were separated by column chromatography and quantified as described by Demirel et al. (2004)
. Total fatty acids (as used in this study) was taken as the sum of all the phospholipid and neutral lipid fatty acids quantified. Total fatty acids included some unassigned fatty acids.
Fatty acid results are reported for major fatty acids and minor identifiable fatty acids relevant to the study. As a result of incomplete resolution, the trans-18:1 isomers are reported as a single value that does not include minor isomers (trans-13, trans-16-18:1) not resolved from cis-18:1n-9 and cis-18:1n-7. In addition to the minor crosscontamination of the latter 2 fatty acids, the fatty acid listed as cis-16:1 consists of n-9 and n-7 isomers. Fatty acid data are presented as weight of fatty acids, in milligrams, per 100 g of muscle tissue.
Additional phenotypic measurements of carcass composition and meat quality were also taken on this population of lambs (Karamichou et al., 2006a
). In particular, computerized tomography (CT) assessments of carcass composition were obtained on these lambs at 5 mo of age, plus 350 females from the same families. From each scan image, the areas and image densities were obtained for the fat, muscle, and bone components of the carcass. Our previous study (Karamichou et al., 2006a
) showed that muscle density was highly heritable and also was correlated with many meat quality traits influenced by i.m. fat content. Thus, we assessed phenotypic and genetic correlations of muscle density with fatty acids.
Genotyping
The DNA was extracted from blood samples for the 9 sires and their 300 progeny for a partial genome scan covering chromosomes 1, 2, 3, 5, 14, 18, 20, and 21. These 8 chromosomes were chosen because previous work had revealed QTL affecting meat and carcass traits on these chromosomes in sheep or on equivalent regions in other species (Nicoll et al., 1998
; Walling et al., 1998a
; Elo et al., 1999
; Freking et al., 1999
; Kmiec, 1999
; Stone et al., 1999
; Broad et al., 2000
; Freking et al., 2002
).
Informative marker panels were developed separately for each sire, containing an average of 16, 8, 10, 6, 7, 6, 7, and 4 informative microsatellite markers per sire on chromosomes 1, 2, 3, 5, 14, 18, 20, and 21, respectively. Initially, each sire was genotyped for all available microsatellite markers across each candidate region, and heterozygous markers were then selected at approximately 10-cM intervals, wherever possible. All offspring were subsequently genotyped for selected markers that were heterozygous in their sire. In total, 139 markers were included in this study.
Data Analysis
All fatty acid measurement traits were right-skewed and were therefore log-transformed before further analysis. Restricted maximal likelihood methods were used to estimate variance components using an animal model, fitting the complete pedigree structure (4,847 animals), using the ASReml package (Gilmour et al., 2004
). The fixed effects included in this analysis were year of birth (2000, 2001, 2002, and 2003), management group (1 or 2), litter size (1 = single or 2 = twins and triplets), and year x slaughter day (15 classes). Triplets comprised proportionately less than 0.02 of the data, and hence were combined with the twin lambs in the litter size classification. The only significant interaction was between year and group.
Heritability estimates were calculated for each log-transformed fatty acid, for total SFA (i.e., the sum of the saturated fatty acids), for total MUFA, and for total PUFA, using a bivariate analysis with CT-assessed muscle density. These analyses also yielded genetic and phenotypic correlations between each fatty acid measurement and muscle density.
Line Effects.
To estimate genetic differences between the LEAN and FAT lines for each trait, using all of the data rather than just the 50 pure-line lambs, true line effects were estimated as the generalized least squares solutions to equations describing the genetic composition of each line or cross. Details of the methodology for predicting means and variances for line categories have been provided elsewhere (Karamichou et al., 2006a
). Line effects were estimated using untransformed fatty acid data.
QTL Analysis
The QTL analysis was performed on the log-transformed data for each fatty acid, and for the 3 categories of total fatty acids (SFA, MUFA, and PUFA).
Map Construction.
Linkage analysis was carried out with the CriMap program, option build (Green et al., 1990
). Marker locations were in close agreement with previous studies (Maddox et al., 2001
). In cases in which there was a disagreement with the published linkage map, the marker order was checked using the CriMap flips option. The marker order with the greatest likelihood was chosen to create a consensus linkage map that was used in subsequent QTL analyses.
QTL Detection.
Marker information contents were obtained as in Knott et al. (1998)
. A regression method for QTL detection was applied according to Knott et al. (1996)
. The analysis used the multimarker approach for interval mapping in half-sib families. The method contained the following steps: for each offspring the probability of inheriting a particular sire haplotype was calculated at 1-cM intervals conditional on the linkage phase of the sire and marker genotypes of the individual and its sire. For a given position, the conditional probabilities of the offspring provided independent variables on which the trait score can be regressed. Therefore, the regression model, in which fixed effects were simultaneously fitted, may be represented as
![]() | [1] |
where yij denotes the trait score of the jth individual originating from sire i; ai is the average effect for half-sib family i; bik is the effect of one of the paternal haplotypes at interval k within half-sib family i; pijk is the probability for individual j of inheriting the first paternal haplotype at interval k conditional on the marker genotypes; and eij is the residual effect for individual j. Fixed effects fitted in these analyses were those found to be significant, for each trait, by the variance analysis. To remove possible body size effects, fatty acids were also corrected for carcass weight at slaughter and the QTL reestimated.
For each regression, an F-ratio of the full model including the inheritance probabilities vs. the same model without the inheritance probabilities was calculated across families. In a 1-QTL model, the location with the largest F-ratio was taken to be the best estimated position for a QTL for each trait.
Size of QTL Effects.
For cases where QTL effects were significant, the within-sire substitution effects were obtained for each sire. The average substitution effect was calculated for those sires that showed significant evidence of a segregating QTL (i.e., for which the sire-specific t-statistic was nominally significant at P < 0.05).
For single-QTL analyses, the proportion of the phenotypic variance explained by the QTL (
2QTL/
2P) was calculated as 4 x [(1 Mean Squarefull)/Mean Squarereduced)] (Knott et al., 1996
), where
2QT is the additive variance at the QTL.
Significance Threshold.
Chromosome-wide empirical threshold values of the test statistics from the regression analysis, at and
= 0.05 chromosome level, were estimated with the permutation test procedure of Churchill and Doerge (1994)
. The thresholds varied between chromosomes depending on their length and the markers they contain.
Confidence Intervals.
Bootstrap confidence intervals were estimated, as described by Visscher et al. (1996)
. Additionally, 95% support intervals for QTL location were obtained using the 1-LOD drop approximation (Lander and Botstein, 1989
).
Two-QTL Model.
Once a single QTL on a chromosome had been identified, the presence of a second QTL was investigated by performing a grid search at 1-cM intervals. The 2-QTL model was accepted if there was a significant improvement over the best possible 1-QTL model at P < 0.05 using a variance ratio (F) test with 2 df (for the additional additive effect and position estimated for the second QTL) as an approximate significance test.
| RESULTS |
|---|
|
|
|---|
|
|
QTL Results
A total of 21 QTL, significant at chromosome-wide level, were detected in 6 out of 8 chromosomal regions, for 14 out of 17 fatty acids, and for every category of fatty acid. All families produced evidence for significant QTL in 1 or more regions. A summary of the chromosome-wide significant QTL, along with the proportions of variance attributable to the QTL and the confidence intervals for QTL location are presented in Table 3
. The profile for the QTL affecting linolenic acid is shown in Figure 1
. Results for fatty acids corrected for live weight were essentially identical to those for uncorrected fatty acids. Additionally, the 2-QTL analyses never gave a significantly better fit than the single QTL analyses; hence, all results presented in Table 3
are from the single QTL analyses.
|
|
The most remarkable results were on chromosome 21 where 10 significant QTL were identified. Five of these were at position 58 cM, and a further 3 at position 57 cM. Moreover, apart from arachidonic acid, these QTL all segregated in the same families (S3 and S7). However, these individual fatty acids were not necessarily strongly correlated with each other. Residual correlations between the individual fatty acids were calculated as an approximation of phenotypic correlations. The correlations were seldom strong, ranging from 0.03 to 0.50.
| DISCUSSION |
|---|
|
|
|---|
Line Differences
Selection for increased leanness changed some aspects of fatty acid composition; the FAT line had a significantly higher i.m. fat content (total fatty acids), and a higher content of stearic (18:0), oleic (18:1 n-9), linoleic (18:2 n-6), and DPA (22:5 n-3) fatty acids than the LEAN line. Stearic acid accounts for 18% of the total fatty acid in lamb meat (Rhee, 1992
; Enser et al., 1996
). Several studies have shown that stearic acid is essentially neutral in its effects on serum total cholesterol, similar to oleic acid (Bonanome and Grundy, 1988
; Zock and Katan, 1992
; Kris-Etherton et al., 1993
; Grundy, 1994
; Judd et al., 2002
). It is not clear why dietary stearic acid does not raise human serum cholesterol level as do other SFA. Oleic acid is the primary MUFA in lamb meat and accounts for about 32% of the total fatty acids (Rhee, 1992
; Enser et al., 1996
). Available evidence indicates that, whereas the shorter chain saturated fatty acids raise human serum cholesterol concentrations, the monounsaturated oleic acid does not (Denke, 1994
). The fatty acid composition of muscle influences important components of meat quality, such as meat flavor, flavor liking, and overall acceptability. In the study of Cameron et al. (2000)
, where the effects of genotype, diet, and their interaction on fatty acid composition of i.m. fat in pigs were examined, it was found that linoleic acid was positively correlated with pork flavor (0.33), taste panel-assessed flavor liking (0.23), and overall acceptability (0.23), whereas DPA was negatively correlated with pork flavor (0.30), flavor liking (0.33) and overall acceptability (0.30). In particular, linoleic acid (18:2 n-6), which is low in ruminant tissues (Wood et al., 1999
; about 3.5% of total fatty acids), is a plant fatty acid that can be transformed to conjugate linolenic acid by bacteria in the rumen (Kepler et al., 1966
). This low level of linoleic acid causes the polyunsaturated:saturated fatty acid ratio (an important nutritional index) to be below the recommended dietary value, which is 0.45 (Department of Health, 1994
).
The LEAN line had more CLA than the FAT line, in agreement with the results of Wachira et al. (2002)
in a comparison between (lean) Soay and (relatively fat) Suffolk breeds. Interest in CLA research has increased in the past few years as a result of reports of CLA consumption providing several health benefits (Kramer et al., 1998
). Because plants do not synthesize CLA, ruminant fats in meat are the primary dietary source of CLA for humans (Herbein et al., 2000
). It has been found that CLA has positive effects of reducing cardiovascular risk, protecting against atherosclerosis, is anticarcinogenic, reduces intake, reduces body content of adipose tissue and lipid, and enhances the immune system (Cook et al., 1993
; Ip et al., 1994
; Lee et al., 1994
; Nicolosi et al., 1997
; West et al., 1998
). In summary, these results show that, as the lamb becomes leaner, the PUFA of muscle increase and, in particular, quantities of the beneficial CLA increase.
Inheritance of Fatty Acids
The heritability estimates for most of the fatty acids measured in this study indicate that there is substantial genetic variation in these traits, such that fatty acid composition can potentially be improved, probably through indirect selection. In particular, SFA were moderately heritable. The additive genetic coefficient of variation for these traits was 0.45 for myristic and 0.09 for palmitic and stearic acids (results not shown), indicating that genetic alteration of these fatty acids is feasible. Due to lack of information on heritabilities of fatty acids in sheep, we will compare our results to published studies in other species. In pigs, heritability estimates for the content of palmitic and stearic fatty acids, obtained by Fernandez et al. (2003)
for Iberian pigs, were 0.31 for palmitic and 0.41 for stearic acid. The estimate for palmitic acid is similar to the one obtained in the current study. The review by Sellier (1998)
of 3 previous studies indicated that stearic acid had a higher heritability (0.51) than found in the current study.
The heritability estimates for monounsaturated fatty acids were high for cis-vaccenic (18:1 n-7), and vaccenic (trans 18:1 n-7) acid, and moderate for palmitoleic (16:1), gadoleic (20:1), and oleic (18:1 n-9) acid. In the study of Fernandez et al. (2003)
, the heritability of oleic acid was almost the same (0.30) as in our study.
The heritability of polyunsaturated fatty acids was variable with an average of 0.23. In pigs, Fernandez et al. (2003)
and Sellier (1998)
presented heritabilities for linoleic acid of 0.29 and 0.58, respectively. Linoleic acid appears to be more heritable in pigs than in sheep, although it should be noted that linoleic acid is much lower in ruminant tissues than in pigs and in all species is obtained from the diet rather than being synthesized in the animal. The rumen biohydrogenates a high proportion of linoleic and linolenic acid, reducing concentrations in body tissues (Wood et al., 1994
; 1999
). Significant heritabilities for fatty acids not synthesized by the animal may be a result of different rates of fatty acid catabolism between animals.
Intramuscular fat (total fatty acids) was moderately heritable, although the value reported here is slightly lower than results from studies done in pigs. Malmfors and Nilsson (1979)
obtained heritability estimates of 0.58 and 0.68 for Swedish Landrace and Large White pigs, respectively. Scheper (1979)
reported a heritability estimate of 0.35 for German Landrace pigs, and a heritability of 0.55 was reported for Landrace pigs in Denmark (Just et al., 1986
) and in Switzerland (Schworer et al., 1987
). The average heritability weighted by number of sires for i.m. fat content from previous reports is 0.53 (Sellier, 1998
).
The heritabilities of trait combinations were surprisingly high in our data set. One possible reason for this result is that, by combining correlated traits, the random measurement error or imprecision associated with individual fatty acids is reduced.
Relationships between Fatty Acids and Muscle Density
The results of our previous study (Karamichou et al., 2006a
), where we estimated genetic parameters for carcass composition, assessed in the live animal, and meat quality traits in these same sheep, demonstrated that altering carcass fatness will simultaneously change muscle density and aspects of i.m. fat content. Hence, muscle density was chosen as the trait to be used for estimating phenotypic and genetic correlations with fatty acids, and we observed that muscle density was indeed strongly correlated with i.m. fat content.
There is no published information on genetic relationships between muscle density, assessed by CT, and fatty acids in any species. Genetic correlations of SFA and MUFA with muscle density were negative, with correlations involving vaccenic and cis-vaccenic acids being the strongest. Genetic correlations of polyunsaturated fatty acids with muscle density were mostly negative, except for arachidonic, adrenic, and DPA fatty acids. The CLA, docosahexaenoic acid, and linoleic fatty acids were strongly negatively correlated with muscle density. This result indicates that selection for decreased muscle density increases the concentration of CLA, docosahexaenoic acid, and linoleic fatty acids in lambs. In summary, selection to decrease muscle density, and hence increase i.m. fat, is expected to result in lamb meat with more total fatty acids, but proportionately increased MUFA and decreased PUFA. This is in accordance with the review of meat fatty acid composition by De Smet et al. (2004)
, where they reported that variation in fat content has an effect on fatty acid composition, independent of species or breed and dietary factors. Hence, the content of total SFA and total MUFA increases more rapidly with increasing fatness than does the content of total PUFA, leading to a decrease in the relative proportion of PUFA. This is shown at the phenotypic level in Figure 2
.
|
Chromosome 2 showed significant QTL for linolenic acid (18:3 n-3), DPA (22:5 n-3), EPA (20:5 n-3), and arachidonic acid (20:4 n-6). The QTL for linolenic acid mapped to the same position as slaughter live weight in these sheep (Karamichou et al., 2006b
), and it was located 23 cM distal to the myostatin locus, which is responsible for the double muscling phenotype in cattle. Chromosome 2 was chosen for the mounting evidence of 1 or several QTL for carcass composition segregating around the myostatin locus (Broad et al., 2000
; Walling et al., 2001
; Johnson et al., 2005a
,b
; Clop et al., 2006
). In the study of Clop et al. (2003)
, a significant QTL for linolenic fatty acid was identified on pig chromosome 12 (conservation of synteny with sheep chromosome 11). In the current study, the QTL on chromosome 2 affecting arachidonic fatty acid mapped to the same position as a QTL for muscle density in Blackface sheep (Karamichou et al., 2006b
).
One QTL for CLA was found on chromosome 3, at 159 cM, excluding the candidate locus this chromosome was chosen for, IGF-I at 227 cM. The QTL for muscle density was also significant and mapped to a position (172 cM) close to the CLA fatty acid content QTL. The concentration of CLA in lamb muscle is important in human nutrition because it has been linked to a multitude of metabolic effects, including inhibition of carcinogenesis, reduced fat deposition, altered immune response, and reduced serum lipids (Mulvihill, 2001
).
Another significant region was detected on chromosome 21. Remarkably, 8 significant QTL (Figure 3
) for arachidonic (20:4 n-6), cis-vaccenic (18:1 n-7), stearic (18:0), oleic (18:1 n-9), linoleic (18:2 n-6), palmitic (16:1), linolenic (18:3 n-3), and myristic (14:0) fatty acids were mapped to the same position (5758 cM), and gadoleic (20:1) acid and DPA (22:5 n-3) were mapped to 21 and 0 cM, respectively. Additionally, the families segregating for 8 out of 10 QTL identified, at the genome level, were the same. Hence, these effects likely correspond to a single QTL. Results from other studies in Iberian x Landrace pigs have identified QTL affecting linoleic fatty acid on pig chromosome 4 (Perez-Enciso et al., 2000
; Clop et al., 2003
; conservation of synteny with sheep chromosome 25). Also, Perez-Enciso et al. (2000)
have found that effects on backfat thickness, backfat weight, and LM area in pigs were also significant and mapped to the same position as the linoleic acid content QTL. Linoleic acid is an essential fatty acid for mammals because they lack desaturase capacity beyond the ninth carbon atom (Vance and Vance, 1996
). It is a key component of cellular membranes and a precursor of prostaglandins and thromboxanes. It is also stored in adipose tissue or ß-oxidized for energy production. In fact, it is highly digestible and preferentially deposited compared with other fatty acids (Lawrence and Fowler, 1997
). High linoleic acid content is also associated with toughness and low consumer acceptability in pig meat (Whittington et al., 1986
; Cameron and Enser, 1991
; Lawrence and Fowler, 1997
).
|
This study was also successful in detecting significant QTL related to fatty acid composition on sheep chromosomes 1, 2, 3, 5, 14, 18, and 21, and the same QTL were detected irrespective of whether the data were corrected for live weight or not. These QTL are clearly of potential importance to the sheep industry; however, they first need to be confirmed in independent populations, and more precise genetic markers would be required. Several genes might be selected as functional candidate genes to explain the QTL found in this study, as fatty acid metabolism can be influenced by a large number of genes involved in complex metabolic routes (Clop et al., 2003
). Association analyses between allelic variants of these genes and fatty acid content would need to be performed to find the necessary genetic markers.
| IMPLICATIONS |
|---|
|
|
|---|
| Footnotes |
|---|
2 Corresponding author: Stephen.Bishop{at}bbsrc.ac.uk
Received for publication April 4, 2006. Accepted for publication July 26, 2006.
| LITERATURE CITED |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |