PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of gseBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleGenetics, Selection, Evolution : GSEJournal Front Page
 
Genet Sel Evol. 2010; 42(1): 36.
Published online 2010 September 16. doi:  10.1186/1297-9686-42-36
PMCID: PMC2949606

Mapping Quantitative Trait Loci (QTL) in sheep. III. QTL for carcass composition traits derived from CT scans and aligned with a meta-assembly for sheep and cattle carcass QTL

Abstract

An (Awassi × Merino) × Merino single-sire backcross family with 165 male offspring was used to map quantitative trait loci (QTL) for body composition traits on a framework map of 189 microsatellite loci across all autosomes. Two cohorts were created from the experimental progeny to represent alternative maturity classes for body composition assessment. Animals were raised under paddock conditions prior to entering the feedlot for a 90-day fattening phase. Body composition traits were derived in vivo at the end of the experiment prior to slaughter at 2 (cohort 1) and 3.5 (cohort 2) years of age, using computed tomography. Image analysis was used to gain accurate predictions for 13 traits describing major fat depots, lean muscle, bone, body proportions and body weight which were used for single- and two-QTL mapping analysis. Using a maximum-likelihood approach, three highly significant (LOD ≥ 3), 15 significant (LOD ≥ 2), and 11 suggestive QTL (1.7 ≤ LOD < 2) were detected on eleven chromosomes. Regression analysis confirmed 28 of these QTL and an additional 17 suggestive (P < 0.1) and two significant (P < 0.05) QTL were identified using this method. QTL with pleiotropic effects for two or more tissues were identified on chromosomes 1, 6, 10, 14, 16 and 23. No tissue-specific QTL were identified.

A meta-assembly of ovine QTL for carcass traits from this study and public domain sources was performed and compared with a corresponding bovine meta-assembly. The assembly demonstrated QTL with effects on carcass composition in homologous regions on OAR1, 2, 6 and 21.

Background

Sheep production is a major contributor to global food production and sheep are one of the few sources of meat with little cultural and religious restriction in consumption. Body composition traits in sheep, primarily muscle mass and fatness, are economically important to the sheep meat industry. There are numerous methods to predict body composition in sheep. Much of the variation that exists in sheep body composition is expressed as between- and within-breed differences. In order to understand the genetic architecture of these economically important traits it is essential to accurately define the phenotypes which describe carcass composition [1].

Live-weight is considered as a standard measurement of body mass, but is a poor indicator of body composition due to the inability to distinguish between different stages of physiological maturity. Body weight may be used as indicator of body composition in animals of similar genetic backgrounds and at the same physiological maturity, however, at different maturity stages the accuracy is greatly reduced [2,3]. Improved predictions of carcass composition can be determined by using ultrasound. Such scans provide a basis to estimate breeding values for eye muscle area and subcutaneous fat depth [3-5]. Increased accuracy and prediction of full body carcass characteristics can be achieved using computed tomography (CT) [6,7] but this is not routinely implemented due to cost constraints.

In addition to the difficulties in obtaining accurate carcass measurements, generation intervals are large, time to assessment is long and therefore the response to selection is slow. Therefore, the use of marker assisted selection or MAS is seen as an attractive aid to increase the efficiency of selection for these traits expensive to measure.

Linkage studies indicate the presence of one or a few major genes for increased muscling and fatness in different sheep populations [8-10]. Two full and 12 partial genome scans have reported QTL for carcass composition including bone density on chromosomes 1-6, 8, 18, 20, 21, and 24 in populations of Coopworth, Scottish Blackface, British Texel, Charollais, Suffolk, Texel and different cross-breed sheep populations [8,11-18]. At present two DNA tests (LoinMax and MyoMax; http://www.pfizeranimalgenetics.com.au/sites/PAG/aus/Pages/sheep.aspx[19]) are commercially available, which test for genetic variants in the Carwell and Myostatin genes [8,10,16,17,20-25].

This study uses CT imaging to accurately determine body composition in vivo in relation to body weight at two different stages of maturity. For the first time, a full genome scan was conducted to identify genomic regions associated with CT-derived parameters in an ovine backcross resource population.

Methods

Resource population

A resource population from crosses between fat-tail Awassi (A) and small-framed Merino (M) sheep was established. Further details of the development of the resource population can be found in Raadsma et al. [26,27]. In the QTL study reported here, only phenotypic and genotypic information from the second generation male backcross (AMM) progeny from one of four F1 sires was analysed in full.

Carcass traits

The backcross progeny were weighed approximately bi-monthly until 83 weeks of age. Weights were recorded as non-fasted body weights immediately off pasture on the same day. At 83 weeks of age, male animals were randomly allocated to two management cohorts. Cohort 1 (n = 86) was lot fed for 90 days after which time all animals were CT scanned prior to slaughter at two years of age. Cohort 2 (n = 79) were grazed under paddock conditions for a further 18 months and then lot fed for 90 days followed by CT scanning and slaughter at 3.5 years of age. Both cohorts were fed ad libitum on a grain and lucerne pelleted ratio with a metabolisable energy content of 12.1 MJ/kg during the feedlot period. The two cohorts were created to capture the differences in fat deposition due to changes in maturity.

At the end of the ad libitum phase and three days prior slaughter, CT scanning was used to estimate lean, fat and bone quantities for individual sheep. Animals were fasted overnight, body weights were recorded and animals were scanned using a Hitachi CT-W400 scanner located in the Meat Science Group at the University of New England, Armidale. Animals were restrained in the supine position using three adjustable belts over the abdomen, chest and neck during the scans at 120 kV tube voltages and 150 mA current. Cross-section images were collected every 40 mm starting proximal to the articulatio genus (rear knee joint) and finishing at the first cervical vertebra. Between 24 and 28 images were collected from each animal depending on their length. The carcass weight was estimated from the CT images. Three sets of data (images) were derived from each image by cropping restraining equipment, internal organs and hooves, distal portion of leg, internal fat and kidney, using AUTOCAT [28]. These images provided an estimate of total body composition including hooves, internal organs and abdominal fat (first set), internal fat - comprising kidney, pelvic, mesenteric and heart fat (second set minus third set) and typical carcass components including total lean, carcass lean and total amount of bone (third set). Furthermore AUTOCAT was used to calculate the area, mean pixel value and variance of each tissue group for each animal from the three sets of images. Subcutaneous fat depth was measured over the eye muscle at the first lumbar two thirds ventral to the vertebrae. The area of fat surrounding the eye muscle (M. longissimus dorsi) was termed the subcutaneous fat area. The eye muscle area was estimated by averaging the area of muscle at the closest image to the first lumbar and the next caudal image. Percentages of lean, fat and bone were calculated as a percentage of the carcass weight estimated by CT (i.e. the sum of individual components estimated by CT). A list of all traits used in this study is provided in Table Table11.

Table 1
Summary statistics of traits used in this in this study

A linear model was fitted using SAS (version 9.2) to adjust the scanning results for final body weight and cohort. For some of the traits, a scatter plot of the trait versus final body weight revealed a linear association for the first cohort but a nonlinear association for the second cohort. To allow for this nonlinearity, a quadratic term was included for the second cohort only. The full model allowing for this takes the form

Trait=β0+β1Cohort2+β2FBW+β3Cohort2×FBW+β4Cohort2×FBW2+ε

where Trait is the measurement to be adjusted for, Cohort2 is a 0-1 indicator variable taking the value 1 for the second cohort, FBW is the final body weight of the sheep, and ε is the random error. Non-significant terms from the above model were dropped, with quadratic terms retained for all traits except dressing percentage, carcass bone, percentage fat in carcass, percentage lean in carcass. Carcass weight and final body weight were adjusted only for cohort effects (Additional file 1). Residuals from the fitted models were obtained, and these were treated as the adjusted traits for subsequent QTL mapping.

Marker analysis QTL mapping procedure

A genome scan using 189 polymorphic microsatellite markers covering all 26 sheep autosomes was conducted in 510 backcross animals. For the linkage analysis, genotypic and phenotypic information from the CT scan of 165 animals was used. The procedure of DNA extraction, genotyping, allele calling and map positions has been outlined previously [26].

QTL analyses were performed for all traits using two methods. Based on a type I error of 0.05, the design (n = 160 animals) had a predicted power of 0.88 to detect QTL with 0.5 SD effect [29]. Solutions were obtained using the QTL-MLE procedure for normally distributed traits in 'R' [26]. As described in previous papers [26,27], when using QTL-MLE, a QTL with LOD ≥ 3.0 was deemed highly significant, significant if LOD ≥ 2.0, and suggestive for QTL with 1.75 ≤ LOD < 2.0.

The second method involved regression analysis for a half-sib design implemented using the web-based program QTL Express [30]. QTL with chromosome-wide significance (P < 0.05) were described as suggestive QTL, whereas QTL exceeding the P < 0.01 chromosome-wide levels and P < 0.05 experiment-wide levels were labelled as significant and highly significant QTL, respectively. A two-QTL model was also fitted to the data using a full two-dimensional scan of each chromosome in QTL Express [30].

Meta-assembly

A meta-assembly of QTL identified in this study was conducted by collating all known QTL from public sources for matched traits based on individual QTL locations and meta-scores as described previously [27]. The positions and confidence intervals of ovine and bovine QTL and blocks of conserved synteny across both species were identified and aligned to the genomes of both species. The individual QTL locations and their scores, and meta-score profiles can be browsed at http://crcidp.vetsci.usyd.edu.au/cgi-bin/gbrowse/oaries_genome/. In addition to the lactation traits, QTL profiles for growth, body weight and carcass composition can now be browsed on this website. Growth and body weight meta-scores from the first paper of this series [26] were also loaded into the website. The carcass composition traits were summarised into four trait classes: bone (percentage bone, bone weight, bone yield), fat (fat yield, back fat, fat depth, marbling, fat thickness, subcutaneous fat thickness), muscle (longissimus muscle area, rib eye area, carcass yield, retail product yield, shear force, lean meat yield) and weight (hot and cold carcass weight, yearling, weaning and slaughter weight). Single and aggregated bars, heat maps and plots can be selected for sheep and cattle as well as meta-scores for both species. Hyperlinks to the original manuscript reference are given.

Results

Analysis of carcass data

The summary statistics for each phenotype are shown in Table Table1.1. For the second cohort, carcass weight and predicted carcass weight from the scan were highly correlated (r = 0.90, P < 0.01) and both traits were also highly correlated with final body weight (r = 0.92 and 0.89, for both cohorts respectively, P < 0.01) (Additional file 2). Across both cohorts, the average body weight at scanning was 51 kg, with an average carcass weight of 28 kg (dressing percentage 55%). Animals from cohort 2 were significantly (P < 0.01) heavier, with a higher mass of total bone, fat and lean compared to cohort 1. However, they had a significantly (P < 0.01) lower percentage bone in the carcass (Additional file 3). Within tissue groups, lean, fat (except internal fat and subcutaneous fat depth) and bone parameters were significantly correlated (r = 0.27 to 0.81, all P < 0.01) (Additional file 4). Significant correlations (P < 0.05) were also detected between many traits among fat and lean tissue groups, with the highest correlation between percentage lean and fat (r = -0.97, P < 0.01). No significant correlations were detected between carcass bone, total bone and eye muscle area and most of the other traits (Additional file 4).

Putative QTL identified

In total, three highly significant (LOD ≥ 3.0), 15 significant (LOD ≥ 2.0) and 12 suggestive (1.7 < LOD < 2.0) QTL were detected on chromosomes 1 to 3, 6, 7, 9-11, 14, 16 and 23 across the 13 traits using QTL-MLE. A summary of the suggestive and significant QTL positions, effect sizes, and 1-LOD support intervals is shown in Table Table2.2. The genome-wide LOD score profiles for all traits are shown in Figures Figures1,1, ,2,2, ,33 and and4.4. With the exception of one suggestive QTL on chromosome 6, all QTL detected by QTL-MLE were confirmed by the QTL regression analysis of QTL Express. A total of five highly significant (experiment-wide P < 0.05), six significant (chromosome-wide P < 0.01) and 34 suggestive (chromosome-wide P < 0.05) QTL were identified on chromosomes 1-3, 6, 7, 9, 10, 11, 14, 16, 19, 23 and 26 using QTL Express (Additional file 5). Among these, two significant (chromosome-wide P < 0.01) and 16 suggestive (chromosome-wide P < 0.05) QTL on chromosomes 6, 8-14, 16, 23 and 26 were not detected using QTL-MLE. Confidence intervals and 1-LOD support intervals for QTL locations extended across a large proportion of each of the chromosomes (Table (Table2,2, additional file 5).

Table 2
Summary of QTL for carcass traits using QTL-MLE
Figure 1
QTL map of the entire genome for body and carcass weight and dressing percentage.
Figure 2
QTL map of the entire genome for carcass lean, total lean, eye muscle area and lean percentage.
Figure 3
QTL map of the entire genome for carcass fat, total fat, internal fat, subcutaneous fat depth, subcutaneous fat area and percentage fat.
Figure 4
QTL map of the entire genome for total bone, carcass bone and bone percentage.

Common QTL for body and carcass weight were identified on chromosomes 2, 6 and 11 using both QTL analysis methods, in addition to the QTL for body weight on chromosome 16 and for dressing percentage on chromosome 14. For muscle traits, eight QTL were detected on seven chromosomes, for fat traits ten QTL on seven chromosomes and for bone traits only two QTL. There were no QTL which solely contributed to traits related to a single tissue i.e. QTL just for muscle, fat or bone. For chromosomes 1, 6, 10, 14, 16 and 23, the QTL for different tissue groups acted pleiotropically, with the same QTL describing traits for different tissue groups. Among the six QTL identified on chromosome 6, two were for weight and three for fat parameters, although the peak positions of the QTL for these two traits groups differed. Similarly, the QTL regions for final body weight, percent lean and subcutaneous fat area were all on chromosome 16, but the peak positions varied. The effect sizes of the QTL ranged from 0.73 to 0.99 SD (Table (Table2)2) and accounted for 3.8 to 9.4% of the phenotypic variance (Additional file 5). Three of the QTL identified here were deemed cryptic QTL, with an effect opposite to what was expected based on breed of origin.

The two-QTL model implemented in QTL Express showed four pairs of QTL which were separated by at least one marker; carcass lean (OAR1), percent bone (OAR1), percent fat (OAR18) and internal fat (OAR19). QTL for carcass lean on chromosome 1 were in coupling phase, whereas all other QTL pairs were in repulsion phase. The QTL in repulsion phase were not identified using the single QTL model since the opposite sign of the QTL effects may have prevented detection under the single QTL model. Details describing QTL positions and effect sizes, and comparisons with single and no QTL models are in Table Table33.

Table 3
Summary of significant QTL for carcass traits using QTL Express under a two-QTL model

Meta-assembly

Published QTL reports for carcass traits in sheep, comprising four genome-wide linkage studies [26,31-33] and 13 partial genome scans [8,11,13-18,34-36] were used for the meta-assembly. QTL for a wide range of carcass traits, including traits not measured in our study (muscle growth, muscle depth, and meat colour), were reported on chromosomes 1-6, 8, 11, 18, 20, 21, 23, 24 and 26 in various sheep populations [8,13,15-18,31-33,35,36]. For two of the studies, the locations of the QTL were not given [11,34]. No QTL were reported on chromosomes 7, 9, 10, 12-17, 19, 22, and 25, but these results might be biased due to partial genome scans, favouring chromosomes with known QTL or candidate genes. The meta-scores showed consistency on six regions of interest across multiple studies for fat, muscle and weight traits, specifically for fat on OAR2 (BTA2) and OAR6 (BTA6), for muscle QTL on OAR2 (BTA2) and for weight on OAR1 (BTA1), 6 (BTA6) and 21 (BTA29) (Figure (Figure55).

Figure 5
Comparative genome map of aggregated meta-scores for carcass-related QTL derived from sheep and cattle studies.

The results of the ovine and bovine meta-assembly are shown as a comparative meta-score plot against the ovine genome in Figure Figure55 and are visualised on the ovine genome browser http://crcidp.vetsci.usyd.edu.au/cgi-bin/gbrowse/oaries_genome/. The very broad range of traits describing carcass and body composition in cattle resulted in QTL being reported on all chromosomes. Furthermore, in contrast to studies in sheep, the majority of studies in cattle reviewed here refer to genome-wide genome scans (n = 14) [37-39]. In addition, eight partial genome scans or candidate gene analyses in cattle were included here [40-47].

Discussion

This study is interesting in that it is the fourth full genome scan for mapping QTL in sheep with respect to carcass traits, and the first where carcass traits were determined from data derived by CT scan which can provide highly accurate profiles of tissue distribution.

Analysis of carcass data

CT scanning was first developed for medical applications and has been extended to animal applications since the 1980s, firstly in pigs and subsequently in sheep [48]. Experiments in sheep and lambs showed that the correlation between CT measures of carcass composition and those derived from manual dissection is very high, but CT or virtual dissection is more precise and reliable [48]. Our study confirmed the high correlation between carcass weight and estimates of carcass weight from scanning [49]. Compared to ultrasound, the standard errors of the predicted values are lower [48,50]. Vester-Christensen et al. [51] and Young et al. [48] have proposed that CT scanning should be an essential reference tool for body and carcass composition. The use of the more precise phenotypes derived from CT measures will also lead to better phenotypes for genetic analysis. Heritabilities for CT-derived traits have been found to be moderate to high [48,52,53]. Theoretical predictions of the genetic progress by incorporating CT traits into selection indices suggest increases in response by 50% or even 100% when combining different measurement methods [6].

The sheep in our study were managed as two cohorts. These cohorts differed significantly in carcass weight and stage of maturity and were considerably heavier than animals in studies published previously [49]. Animals investigated here were taken to a greater stage of maturity to measure specific effects on fat and fat distribution. Sheep from cohort 1 had similar muscle/carcass lean weights compared to meat sheep [54] and Norwegian lambs [49]. However, for both these studies, the proportion of muscle was higher than in our study, largely due to differences in fatness and stage of development (age, maturity). For the same reasons, the proportion of bone in the carcass was lower in our study than in studies presented by Young et al. [54] and Kongsro et al. [49].

The main focus of our project was the study of fat characteristics in the carcass. Therefore, older and consequently more mature sheep were used. Adjusting body composition traits for body weight at the time of scanning was considered the best method to accurately measure tissue groups independently of their body mass. Animals from the second cohort had higher fat content and total percent fat compared to animals from cohort 1. There were significant correlations between the major tissue groups (lean, fat and bone). Fat traits tended to be significantly and negatively correlated with lean traits, as reported by Lambe et al. [55]. Without adjusting for body weight, the correlations would have been strongly positive [55,56], as was also the case here (results not shown). The importance of adjustment for body weight is that properties of body tissue can be investigated free from the effects of body mass. The differences in stage of maturity resulted in different adjustments for body weight, namely a linear effect for cohort 1 and a curvilinear effect for cohort 2, suggesting a plateau of growth had been reached and animals were in the mature fattening phase of development.

QTL analysis

Genome-wise error rates were controlled by adjustment of P-values through the use of a chromosome- and experiment-wide permutation test in the case of QTL Express, therefore the number of false positive QTL was assumed to be minimal. For the maximum-likelihood analysis we chose thresholds for a LOD statistic which was deemed to be conservative at LOD of 2 (P ≈ 0.01) and LOD of 3 (P ≈ 0.001). The close agreement between the number of QTL detected in each method suggests that the likelihood of random false positives is expected to be small.

For body and carcass weight, QTL were identified on chromosomes 2, 6, 11 and 16. The QTL on chromosomes 6 and 11 were consistent with those reported in the same study population at earlier time points [26]. The QTL for final and carcass weight on chromosome 2 was the only one that corresponded to a QTL for live weight in Scottish Blackface and Suffolk, Texel sheep [13,17]. A total of eight QTL across seven chromosomal regions were identified for muscle. QTL on chromosomes 1 and 6 were consistent with other studies in Suffolk and Texel populations [11,16,17], whereas QTL on chromosomes 7, 9, 10, 16 and 23 can be considered novel.

QTL for fat have previously been reported on OAR 1-4, 18 and 20 in different sheep populations [14,16,17,31,33,34]. Within the confidence interval of our QTL, we confirmed QTL on chromosome 1 and 3, and novel QTL were identified on OAR 6, 10, 14, 16 and 23. QTL for fatness have consistently been reported on chromosomes 2, 3 and 18 [14,16,17], but the QTL on OAR18 was only identified using the two-QTL model and no QTL on OAR2 was detected in the current study despite the emphasis on fat traits.

Few reports are available for bone-related traits in sheep, and no QTL study on bone yield in the carcass has been reported. Previous QTL studies have analysed bone density and cross sectional area in Scottish Blackface and Coopworth sheep [13,31,32]. The two QTL detected here for bone yield suggest that the QTL landscape is rather featureless for this trait.

In summary, the first interesting discovery of this paper was the identification of novel QTL with small to moderate effects on body composition and body weight on chromosomes 1, 6, 7, 9, 10, 14, 16 and 23. This may in part be due to an increase in accuracy of phenotyping using CT image analysis.

A notable finding of this study was that there were no QTL which exclusively affected multiple measures of the same tissue group, i.e. fat, lean or bone. The effect of measuring fat at individual or a limited number of sites was discussed by Thompson [57], who proposed that individual depots may not reflect total body fat in the animal. This is seen in the correlations of non-unity between traits indicative of fatness at different sites (Table (Table1).1). This suggests that different measures of the same tissue reflect different traits with different QTL. This may have implications for QTL detection and application. For instance, QTL used to reduce subcutaneous fat content may not necessarily result in a reduction of total carcass fat.

Despite many QTL reports and a significant association on chromosome 18, we could not confirm the effects of the important loci such as the Carwell and Callipyge genes with known effects on muscle lean in sheep [20,58,59]. These and other genes on the same chromosomal region are known to be imprinted with paternally expressed protein-coding genes, as well as several maternally expressed non-coding RNA genes [20], which may have prevented their detection in our study, which used only one paternal half-sib family. Single-marker association analyses revealed significant associations for markers close these genes (results not shown) but this requires more detailed analysis. In future studies, the use of multi-sire families and linkage disequilibrium among maternal alleles should focus on the identification of these imprinting effects.

We considered an interdependency of QTL for body weight on OAR11 and 16. These chromosomes contain the growth hormone (GH) and growth hormone receptor (GHR) genes, with known effects on body weight and carcass composition across species [60-66]. Even though we identified QTL for final body weight on both chromosomes, we could not detect an interaction between these two QTL and assume that the genes underlying these two QTL act in a simple additive fashion.

Although we examined the likely importance of two QTL for all traits and report on five cases with significant support for QTL pairs, the interpretation of the results warrant caution, especially where the QTL are closely located or no corresponding QTL were detected under the single QTL model.

Meta-assembly and comparative analysis

A meta-assembly of QTL identified for carcass traits was conducted by collating all known ovine QTL from public sources for matched traits, as previously described [27]. Additionally, studies in cattle were summarised using the same methodology. A summary of the carcass meta-scores from cattle and sheep that is incorporated into the ovine genome browser http://crcidp.vetsci.usyd.edu.au/cgi-bin/gbrowse/oaries_genome/ is shown in Figure Figure5.5. Furthermore, due to the lack of studies based on CT-derived phenotypes, and the different methods, models, and population types used across studies, we considered them to be the same carcass trait if they described the same carcass characteristic.

Despite the large number of QTL detected in cattle, relatively few QTL in sheep were found in comparative locations. However, for some traits, especially for muscle and fat, loci were mapped to homologous regions; these QTL may possibly describe the same gene.

Among the 11 studies summarised using the meta-assembly in cattle, four reported QTL for body weight on chromosome 14 [39,46,67,68], but no QTL was found in sheep in the homologous position. The highest ovine meta-scores for body weight were derived for chromosomes 1 [13,14,17], 6 [26] and 21 [13,26]. QTL were reported for body weight in cattle on the homologous chromosomes 1 [37,39,67], 6 [40,67,69] and 29 [39,70,71].

Amongst the six ovine and 13 bovine QTL studies reporting linkage regions for muscle related traits, we found one region in common between sheep and cattle, namely OAR2 and BTA2 for carcass weight, eye muscle area and retail product yield [8,16,36-38,67,70,72] (Figure (Figure55).

High meta-scores for fat QTL were derived for sheep and cattle on homologous chromosomes OAR2/BTA2 and OAR6/BTA6 [37,38,45,67,68,70,71,73]. However, no QTL for fat traits were identified in sheep, which align to bovine chromosome 14, where the highest meta-score was calculated for corresponding traits in cattle [46,67,69,70,74].

No homologous regions were found between sheep and cattle for bone traits. This is likely due to the limited number of studies conducted and QTL reported for these traits in either species.

Conclusion

This is the first study using CT-derived carcass measures for a full-genome scan in sheep. To our knowledge this is the only study with a focus on carcass fat characteristics in mature sheep. We present evidence for a significant number of new QTL for muscle, fat and bone traits in sheep. We also confirm and support the presence of previously published QTL in breeds other than those studied here. Finally, homology in QTL regions between sheep and cattle for muscle, bone, fat and body weight was demonstrated.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

CRC analysed the CT scans of the second cohort and ran the QTL analyses, was responsible for the data assembly, phenotypic analysis, contributed to genotyping and preparing the manuscript. EJ ran the final models and QTL analyses, participated in the manuscript preparation. MH calculated the QTL scores and designed the online data source. PCT developed the QTL methodology, implemented the QTL-MLE program, and contributed to manuscript preparation and the overall design. IT was involved in the project management and participated in the preparation of the manuscript, HWR was responsible for the overall design, project management, and was involved in analysis of the data and writing the manuscript. All authors read and approved the final manuscript.

Supplementary Material

Additional file 1:

Summary of linear models for trait pre-correction used in this study. R-square, overall P-value (Model P-value) and P-values for the regression coefficients in the fitted models. Models were adjusted for final body weight (FBW) and cohort allowing for nonlinearity if observed in a scatter plot of the trait versus final body weight

Additional file 2:

Phenotypic correlation between body weight and carcass weight measures. Phenotypic correlation between body weight and carcass weight measures. Correlations exceeded the P < 0.01 threshold using with n = 160 animals and 1n = 72 animals

Additional file 3:

Summary statistics of traits used in this study. Summary statistics of traits used in this in this study within the two cohorts

Additional file 4:

Phenotypic correlation between the carcass traits. Phenotypic correlation between carcass lean (CL), total lean (TL), eye muscle area (EMA), carcass fat (CF), total fat (TF), internal fat (IF), subcutaneous fat depth (SFD), carcass bone (B), total bone (TB), percent lean in carcass (PL), percent fat in carcass (PF), percent bone in carcass (PB); r > 0.2 corresponds to P < 0.05, and r > 0.3 corresponds to P < 0.01 with n = 160

Additional file 5:

QTL for body weight and carcass traits using QTL Express. QTL for body weight and carcass traits using QTL Express; *chromosome-wide P < 0.05; **chromosome-wide P < 0.01; ***experiment-wide P < 0.05; ****experiment-wide P < 0.01; variance or QTL heritability (Heritab) as a proportion of the phenotypic variance accounted for by the QTL [variance explained by the QTL effect = 1-(MS of full model/MS of reduced model)]

Acknowledgements

The authors are particularly gratefully to Mrs Marilyn Jones and Mrs Gina Attard for their assistance in the genome scan, Mr Dave Palmer, Mr Joe Davis, Ms Renee Deever and the many casual staff who have contributed to the rearing of the lambs, the husbandry of the sheep, and the collation of data. Dr Mary Lam is thanked for her support in the first part of the analysis and the establishment of the database. The research was approved by the University of Sydney Animal Ethics Committee. The study was in part funded with an ARC-Linkage Grant. The resource flock and part of the sheep dairy was established with contributions from Mr Tom and Phillip Grant from Awassi Australia. The animals were scanned at The University of New England, the authors thank the staff of the facility and the visiting students for the analysis of the images for the first cohort in 1999.

References

  • Lawrie RA. Lawrie's meat science. Woodhead Publishing; 1998.
  • Fortin A, Shrestha JNB. In vivo estimation of carcass meat by ultrasound in ram lambs slaughtered at an average live weight of 37 Kg. Anim Prod. 1986;43:469–475. doi: 10.1017/S0003356100002683. [Cross Ref]
  • Stanford K, Clark I, Jones SDM. Use of uUltrasound in prediction of carcass characteristics in lambs. Can J Anim Sci. 1995;75:185–189. doi: 10.4141/cjas95-027. [Cross Ref]
  • Miles KA, Hayball M, Dixon AK. Color perfusion imaging - a new application of computed-tomography. Lancet. 1991;337:643–645. doi: 10.1016/0140-6736(91)92455-B. [PubMed] [Cross Ref]
  • Gilmour AR, Luff AF, Fogarty NM, Banks R. Genetic-parameters for ultrasound fat depth and eye muscle measurements in live poll Dorset sheep. Aust J Agric Res. 1994;45:1281–1291. doi: 10.1071/AR9941281. [Cross Ref]
  • Simm G, Dingwall WS. Selection indexes for lean meat production in sheep. Livest Prod Sci. 1989;21:223–233. doi: 10.1016/0301-6226(89)90052-3. [Cross Ref]
  • Jopson NB, McEwan JC, Fennessy PF, Dodds KG, Nicoll GB, Wade CM. Economic benefits of including computed tomography measurements in a large terminal sire breeding programme. Proceedings of the Twelfth Conference of the Association for the Advancement of Animal Breeding and Genetics: 6-10 April 1997; Dubbo. 1997. pp. 72–76.
  • Marcq F, Larzul C, Marot V, Bouix J, Eychenne F, Laville E, Bibe B, Leroy PL, Georges M, Elsen JM. Preliminary results of a whole-genome scan targeting QTL for carcass traits in a Texel * Romanov intercross. Proceedings of the 7th World Congress on Genetics Applied to Livestock Production: 19-23 August 2002; Montpellier. 2002. pp. 0–4.
  • Marshall K, Henshall J, Banks RG, Werf JHJvd. Finding major gene effects in Australian meat sheep - feasibility study for a Texel dataset. Proceedings of the Thirteenth Conference Association for the Advancement of Animal Breeding and Genetics: 4-7 July 1999; Mandurah. 1999. pp. 86–89.
  • Marcq F, Elsen JM, El Barkouki S, Bouix J, Eychenne F, Grobet L, Karim L, Laville E, Nezer C, Royo L, Sayd T, Bibe B, Le Roy PL, Georges M. Investigating the role of myostatin in the determinism of double muscling characterizing Belgian Texel sheep. Anim Genet. 1998;29(Suppl 1):52–53.
  • Karamichou E, Nute GR, Richardson RI, McLean K, Bishop SC. Identifying QTL for meat quality and carcass composition traits in Blackface sheep. Proceedings of the British Society of Animal Science Annual Conference; 4-6 April 2005; York. 2005. p. 45.
  • Karamichou E, Richardson RI, Nute GR, Bishop SC. Identification of quantitative trait loci affecting fatty acid composition in sheep. Proceedings of the 8th World Congress on Genetics Applied to Livestock Production; 13-18 August 2006; Belo Horizonte. 2006. pp. 13–12.
  • Karamichou E, Richardson RI, Nute GR, McLean KA, Bishop SC. A partial genome scan to map quantitative trait loci for carcass composition, as assessed by X-ray computer tomography, and meat quality traits in Scottish Blackface Sheep. Anim Sci. 2006;82:301–309.
  • McRae AF, Bishop SC, Walling GA, Wilson AD, Visscher PM. Mapping of multiple quantitative trait loci for growth and carcass traits in a complex commercial sheep pedigree. Anim Sci. 2005;80:135–141. doi: 10.1079/ASC41040135. [Cross Ref]
  • Johnson PL, McEwan JC, Dodds KG, Purchas RW, Blair HT. A directed search in the region of GDF8 for quantitative trait loci affecting carcass traits in Texel sheep. J Anim Sci. 2005;83:1988–2000. [PubMed]
  • Walling GA, Visscher PM, Simm G, Bishop SC. Confirmed linkage for qtls affecting muscling in texel sheep on chromosomes 2 and 18. Proceedings of the 52nd Annual Meeting of the European Association for Animal Production; 26-29 August 2001; Budapest. 2001. p. G5.6.
  • Walling GA, Visscher PM, Wilson AD, McTeir BL, Simm G, Bishop SC. Mapping of quantitative trait loci for growth and carcass traits in commercial sheep populations. J Anim Sci. 2004;82:2234–2245. [PubMed]
  • McEwan JC, Gerard EM, Jopson NB, Nicoll GB, Greer GJ, Dodds KG, Bain WE, Burkin HR, Lord EA, Broad TE. Localization of a QTL for rib-eye muscling on OAR18. Anim Genet. 1998;29(Suppl 1):66. [PubMed]
  • Pfizer. http://www.pfizeranimalgenetics.com.au/sites/PAG/aus/Pages/sheep.aspx
  • Cockett NE, Smit MA, Bidwell CA, Segers K, Hadfield TL, Snowder GD, Georges M, Charlier C. The callipyge mutation and other genes that affect muscle hypertrophy in sheep. Genet Sel Evol. 2004;36(Suppl 1):S65–S81. [PMC free article] [PubMed]
  • Smit M, Segers K, Carrascosa LG, Shay T, Baraldi F, Gyapay G, Snowder G, Georges M, Cockett N, Charlier C. Mosaicism of solid gold supports the causality of a noncoding A-to-G transition in the determinism of the callipyge phenotype. Genetics. 2003;163:453–456. [PubMed]
  • Forrest RH, Hickford JGH, Hogan A, Frampton C. Polymorphism at the ovine beta(3)-adrenergic receptor locus: associations with birth weight, growth rate, carcass composition and cold survival. Anim Genet. 2003;34:19–25. doi: 10.1046/j.1365-2052.2003.00936.x. [PubMed] [Cross Ref]
  • Masri AY, Lambe NR, Macfarlane JM, Brotherstone S, Haresign W, Rius-Vilarrasa E, Bunger L. The effects of a loin muscling quantitative trait locus (LoinMAX (TM)) on carcass and VIA-based traits in crossbred lambs. Animal. pp. 407–416. [PubMed] [Cross Ref]
  • Campbell AW, McLaren RJ. LoinMAXTM and MyoMAXTM: taking DNA marker tests from the research environment to commercial reality. Proc N Z Soc Anim Prod. 2007;67:160–162.
  • Bain WE, Johnson PL, Greer GJ, Dodds KG, McLean NJ, McLaren RJ, Galloway SM, Stijn TCv, McEwan JC. Effect of MyoMAX on carcass lean and fat. Proc N Z Soc Anim Prod. 2008;68:43–44.
  • Raadsma HW, Thomson PC, Zenger KR, Cavanagh C, Lam MK, Jonas E, Jones M, Attard G, Palmer D, Nicholas FW. Mapping quantitative trait loci (QTL) in sheep. I. A new male framework linkage map and QTL for growth rate and body weight. Genet Sel Evol. 2009;41:34. doi: 10.1186/1297-9686-41-34. [PMC free article] [PubMed] [Cross Ref]
  • Raadsma HW, Jonas E, McGill D, Hobbs M, Lam MK, Thomson PC. Mapping quantitative trait loci (QTL) in sheep. II. Meta-assembly and identification of novel QTL for milk production traits in sheep. Genet Sel Evol. 2009;41:45. doi: 10.1186/1297-9686-41-45. [PMC free article] [PubMed] [Cross Ref]
  • Jopson NB, Kolstad K, Sehested E, Vangen O. Computer tomography as an accurate and cost effective alternative to carcass dissection. Proc Aust Assoc Anim Breed Genet. 1995. pp. 635–638.
  • Van der Werf J. Power calculation to determine sample size for given 2-way contrast and power/type 1 error. http://www-personal.une.edu.au/~jvanderw/software.htm
  • Seaton GGR, Haley CS, Knott SA, Visscher PM. QTL express: rapid and user-friendly mapping of quantitative trait loci in livestock. Proceedings of the 7th World Congress on Genetics Applied to Livestock Production: 19-23 August 2002; Montpellier. 2002. pp. 1–2.
  • Campbell AW, Bain WE, McRae AF, Broad TE, Johnstone PD, Veenvliet BA, Greer GJ, Glass BC, Beattie AE, Jopson NB, Detection of quantitative trait loci for bone mineral density in Coopworth sheep. Proceedings of the 7th World Congress on Genetics Applied to Livestock Production: 19-23 August 2002; Montpellier. 2002. pp. 0–4.
  • Campbell AW, Bain WE, McRae AF, Broad TE, Johnstone PD, Dodds KG, Veenvliet BA, Greer GJ, Glass BC, Beattie AE. et al. Bone density in sheep: genetic variation and quantitative trait loci localisation. Bone. 2003;33:540–548. doi: 10.1016/S8756-3282(03)00228-X. [PubMed] [Cross Ref]
  • Walling GA, Wilson AD, McTeir BL, Visscher PM, Simm G, Bishop SC. A candidate region approach allows efficient QTL detection in UK Suffolk and Texel populations. Proceedings of the 7th World Congress on Genetics Applied to Livestock Production: 19-23 August 2002; Montpellier. 2002. pp. 0–4.
  • Broad TE, Veenvliet BA, Glass BC, Johnstone PD, Jopson NB, Bain WE, Greer GJ, McEwan JC, Clarke JN. A partial genome scan reveals three QTL affecting fat deposition in sheep. HUGO - Human genome meeting. Edinburgh. 2001.
  • Karamichou E, Richardson RI, Nute GR, Gibson KP, Bishop SC. Genetic analyses and quantitative trait loci detection, using a partial genome scan, for intramuscular fatty acid composition in Scottish Blackface sheep. J Anim Sci. 2006;84:3228–3238. doi: 10.2527/jas.2006-204. [PubMed] [Cross Ref]
  • Laville E, Bouix J, Sayd T, Bibe B, Elsen JM, Larzul C, Eychenne F, Marcq F, Georges M. Effects of a quantitative trait locus for muscle hypertrophy from Belgian Texel sheep on carcass conformation and muscularity. J Anim Sci. 2004;82:3128–3137. [PubMed]
  • Stone RT, Keele JW, Shackelford SD, Kappes SM, Koohmaraie M. A primary screen of the bovine genome for quantitative trait loci affecting carcass and growth traits. J Anim Sci. 1999;77:1379–1384. [PubMed]
  • Abe T, Saburi J, Hasebe H, Nakagawa T, Kawamura T, Saito K, Nade T, Misumi S, Okumura T, Kuchida K. et al. Bovine quantitative trait loci analysis for growth, carcass, and meat quality traits in an F-2 population from a cross between Japanese Black and Limousin. J Anim Sci. 2008;86:2821–2832. doi: 10.2527/jas.2007-0676. [PubMed] [Cross Ref]
  • Kim JJ, Farnir F, Savell J, Taylor JF. Detection of quantitative trait loci for growth and beef carcass fatness traits in a cross between Bos taurus (Angus) and Bos indicus (Brahman) cattle. J Anim Sci. 2003;81:1933–1942. [PubMed]
  • Setoguchi K, Furuta M, Hirano T, Nagao T, Watanabe T, Sugimoto Y, Takasuga A. Cross-breed comparisons identified a critical 591-kb region for bovine carcass weight QTL (CW-2) on chromosome 6 and the Ile-442-Met substitution in NCAPG as a positional candidate. BMC Genetics. 2009;10:43. doi: 10.1186/1471-2156-10-43. [PMC free article] [PubMed] [Cross Ref]
  • Taylor JF, Coutinho LL, Herring KL, Gallagher DS, Brenneman RA, Burney N, Sanders JO, Turner RV, Smith SB, Miller RK. et al. Candidate gene analysis of GH1 for effects on growth and carcass composition of cattle. Anim Genet. 1998;29:194–201. doi: 10.1111/j.1365-2052.1998.00317.x. [PubMed] [Cross Ref]
  • Moore SS, Li C, Basarab J, Snelling WM, Kneeland J, Murdoch B, Hansen C, Benkel B. Fine mapping of quantitative trait loci and assessment of positional candidate genes for backfat on bovine chromosome 14 in a commercial line of Bos taurus. J Anim Sci. 2003;81:1919–1925. [PubMed]
  • Imai K, Matsughige T, Watanabe T, Sugimoto Y, Ihara N. Mapping of a quantitative trait locus for beef marbling on bovine chromosome 9 in purebred Japanese Black cattle. Anim Biotechnol. 2007;18:75–80. doi: 10.1080/10495390601090992. [PubMed] [Cross Ref]
  • Marques E, Nkrumah JD, Sherman EL, Moore SS. Polymorphisms in positional candidate genes on BTA14 and BTA26 affect carcass quality in beef cattle. J Anim Sci. 2009;87:2475–2484. doi: 10.2527/jas.2008-1456. [PubMed] [Cross Ref]
  • Li C, Basarab J, Snelling WM, Benkel B, Kneeland J, Murdoch B, Hansen C, Moore SS. Identification and fine mapping of quantitative trait loci for backfat on bovine chromosomes 2, 5, 6, 19, 21, and 23 in a commercial line of Bos taurus. J Anim Sci. 2004;82:967–972. [PubMed]
  • Mizoshita K, Watanabe T, Hayashi H, Kubota C, Yamakuchi H, Todoroki J, Sugimoto Y. Quantitative trait loci analysis for growth and carcass traits in a half-sib family of purebred Japanese Black (Wagyu) cattle. J Anim Sci. 2004;82:3415–3420. [PubMed]
  • Casas E, Stone RT, Keele JW, Shackelford SD, Kappes SM, Koohmaraie M. A comprehensive search for quantitative trait loci affecting growth and carcass composition of cattle segregating alternative forms of the myostatin gene. J Anim Sci. 2001;79:854–860. [PubMed]
  • Young MJ, Nsoso SJ, Logan CM, Beatson PR. Prediction of carcass tissue weight in vivo using live weight, ultrasound or X-ray CT measurements. New Zealand Society of Animal Production: 1996. 1996. pp. 205–211.
  • Kongsro J, Roe M, Aastveit AH, Kvaal K, Egelandsdal B. Virtual dissection of lamb carcasses using computer tomography (CT) and its correlation to manual dissection. J Food Eng. 2008;88:86–93. doi: 10.1016/j.jfoodeng.2008.01.021. [Cross Ref]
  • Sehested E. Evaluation of carcass composition of live lambs based on computed tomography. Proceedings of the 35th Annual Meeting of the EAAP: 6-9 August 1984; The Hague. 1984. G5.21.
  • Vester-Christensen M, Erbou SGH, Hansen MF, Olsen EV, Christensen LB, Hviid M, Ersboll BK, Larsen R. Virtual dissection of pig carcasses. Meat Sci. 2009;81:699–704. doi: 10.1016/j.meatsci.2008.11.015. [PubMed] [Cross Ref]
  • Kvame T, Vangen O. In-vivo composition of carcass regions in lambs of two genetic lines, and selection of CT positions for estimation of each region. Small Rumin Res. 2006;66:201–208. doi: 10.1016/j.smallrumres.2005.09.014. [Cross Ref]
  • Jones HE, Lewis RM, Young MJ, Wolf BT, Warkup CC. Changes in muscularity with growth and its relationship with other carcass traits in three terminal sire breeds of sheep. Anim Sci. 2002;74:265–275.
  • Young MJ, Simm G, Glasbey CA. Computerised tomography for carcass analysis. British Society of Animal Science: 2001. 2001. pp. 250–254.
  • Lambe NR, Conington J, McLean KA, Navajas EA, Fisher AV, Bünger L. In vivo prediction of internal fat weight in Scottish Blackface lambs, using computer tomography. J Anim Breed Genet. 2006;123:105–113. doi: 10.1111/j.1439-0388.2006.00570.x. [PubMed] [Cross Ref]
  • Theriault M, Pomar C, Castonguay FW. Accuracy of real-time ultrasound measurements of total tissue, fat, and muscle depths at different measuring sites in lamb. J Anim Sci. 2009;87:1801–1813. doi: 10.2527/jas.2008-1002. [PubMed] [Cross Ref]
  • Thompson WR, Theuninck DH, Meiske JC, Goodrich RD, Rust JR, Byers FM. Linear Measurements and Visual Appraisal as Estimators of Percentage Empty Body-Fat of Beef-Cows. J Anim Sci. 1983;56:755–760.
  • Cockett NE, Jackson SP, Snowder GD, Shay TL, Berghmans S, Beever JE, Carpenter C, Georges M. The callipyge phenomenon: Evidence for unusual genetic inheritance. J Anim Sci. 1999;77:221–227. [PubMed]
  • Cockett NE, Snowder GD. Genetic influences on carcass merit of sheep. J Dairy Sci. 2001;84(Suppl 1):99.
  • Sherman EL, Nkrumah JD, Murdoch BM, Li C, Wang Z, Fu A, Moore SS. Polymorphisms and haplotypes in the bovine neuropeptide Y, growth hormone receptor, ghrelin, insulin-like growth factor 2, and uncoupling proteins 2 and 3 genes and their associations with measures of growth, performance, feed efficiency, and carcass merit in beef cattle. J Anim Sci. 2008;86:1–16. doi: 10.2527/jas.2006-799. [PubMed] [Cross Ref]
  • Fang MX, Nie QH, Luo CL, Zhang DX, Zhang XQ. Associations of GHSR gene polymorphisms with chicken growth and carcass traits. Mol Biol Rep. pp. 423–428. [PubMed] [Cross Ref]
  • Garrett AJ, Rincon G, Medrano JF, Elzo MA, Silver GA, Thomas MG. Promoter region of the bovine growth hormone receptor gene: Single nucleotide polymorphism discovery in cattle and association with performance in Brangus bulls. J Anim Sci. 2008;86:3315–3323. doi: 10.2527/jas.2008-0990. [PubMed] [Cross Ref]
  • Tatsuda K, Oka A, Iwamoto E, Kuroda Y, Takeshita H, Kataoka H, Kouno S. Relationship of the bovine growth hormone gene to carcass traits in Japanese black cattle. J Anim Breed Genet. 2008;125:45–49. doi: 10.1111/j.1439-0388.2007.00688.x. [PubMed] [Cross Ref]
  • de Faria DA, Guimaraes SEF, Lopes PS, Pires AV, Paiva SR, Sollero BP, Wenceslau AA. Association between G316A growth hormone polymorphism and economic traits in pigs. Genet Mol Biol. 2006;29:634–640.
  • Curi RA, Palmieri DA, Suguisawa L, Ferraz ALJ, de Oliveira HN, Furlan LR, Silveira AC, Lopes CR. Effects of GHR gene polymorphisms on growth and carcass traits in Zebu and crossbred beef cattle. Livest Sci. 2006;101:94–100. doi: 10.1016/j.livprodsci.2005.09.015. [Cross Ref]
  • Maj A, Oprzadek J, Dymnicki E, Zwierzchowski L. Association of the polymorphism in the 5 '-noncoding region of the bovine growth hormone receptor gene with meat production traits in Polish Black-and-White cattle. Meat Sci. 2006;72:539–544. doi: 10.1016/j.meatsci.2005.09.001. [PubMed] [Cross Ref]
  • Takasuga A, Watanabe T, Mizoguchi Y, Hirano T, Ihara N, Takano A, Yokouchi K, Fujikawa A, Chiba K, Kobayashi N. et al. Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping. Mamm Genome. 2007;18:125–136. doi: 10.1007/s00335-006-0096-5. [PubMed] [Cross Ref]
  • Mizoshita K, Takano A, Watanabe T, Takasuga A, Sugimoto Y. Identification of a 1.1-Mb region for a carcass weight QTL on bovine Chromosome 14. Mamm Genome. 2005;16:532–537. doi: 10.1007/s00335-005-0024-0. [PubMed] [Cross Ref]
  • Casas E, Shackleford SD, Keele JW, Stone RT, Kappes SM, Koohmaraie M. Quantitative trait loci affecting growth and carcass composition of cattle segregating alternate forms of myostatin. J Anim Sci. 2000;78:560–569. [PubMed]
  • Casas E, Shackelford SD, Keele JW, Koohmaraie M, Smith TPL, Stone RT. Detection of quantitative trait loci for growth and carcass composition in cattle. J Anim Sci. 2003;81:2976–2983. [PubMed]
  • MacNeil MD, Grosz MD. Genome-wide scans for QTL affecting carcass traits in Hereford × composite double backcross populations. J Anim Sci. 2002;80:2316–2324. [PubMed]
  • Morris CA, Pitchford WS, Cullen NG, Esmailizadeh AK, Hickey SM, Hyndman D, Dodds KG, Afolayan RA, Crawford AM, Bottema CDK. Quantitative trait loci for live animal and carcass composition traits in Jersey and Limousin back-cross cattle finished on pasture or feedlot. Anim Genet. 2009;40:648–654. doi: 10.1111/j.1365-2052.2009.01895.x. [PubMed] [Cross Ref]
  • Gutierrez-Gil B, Williams JL, Homer D, Burton D, Haley CS, Wiener P. Search for quantitative trait loci affecting growth and carcass traits in a cross population of beef and dairy cattle. J Anim Sci. 2009;87:24–36. doi: 10.2527/jas.2008-0922. [PubMed] [Cross Ref]
  • Marquez GC, Enns RM, Grosz MD, Alexander LJ, MacNeil MD. Quantitative trait loci with effects on feed efficiency traits in Hereford × composite double backcross populations. Anim Genet. 2009;40:986–988. doi: 10.1111/j.1365-2052.2009.01946.x. [PubMed] [Cross Ref]

Articles from Genetics, Selection, Evolution : GSE are provided here courtesy of BioMed Central