- Open Access
Mapping quantitative trait loci (QTL) in sheep. I. A new male framework linkage map and QTL for growth rate and body weight
Genetics Selection Evolution volume 41, Article number: 34 (2009)
A male sheep linkage map comprising 191 microsatellites was generated from a single family of 510 Awassi-Merino backcross progeny. Except for ovine chromosomes 1, 2, 10 and 17, all other chromosomes yielded a LOD score difference greater than 3.0 between the best and second-best map order. The map is on average 11% longer than the Sheep Linkage Map v4.7 male-specific map. This map was employed in quantitative trait loci (QTL) analyses on body-weight and growth-rate traits between birth and 98 weeks of age. A custom maximum likelihood program was developed to map QTL in half-sib families for non-inbred strains (QTL-MLE) and is freely available on request. The new analysis package offers the advantage of enabling QTL × fixed effect interactions to be included in the model. Fifty-four putative QTL were identified on nine chromosomes. Significant QTL with sex-specific effects (i.e. QTL × sex interaction) in the range of 0.4 to 0.7 SD were found on ovine chromosomes 1, 3, 6, 11, 21, 23, 24 and 26.
Over the past few decades, a number of quantitative trait loci (QTL) analyses have been conducted on many livestock breeds. These studies have provided very useful genetic information and enriched our knowledge on the underlying biology and genetic architecture of complex traits. A general review of QTL mapping can be found in Weller .
An important input to be considered in QTL studies is the availability of a robust framework map of the genome. The initial work by Crawford et al.  has resulted in the first extensive ovine genetic linkage map covering 2,070 cM of the sheep genome and comprising 246 polymorphic markers . It has been followed by second  and third generation updates . The latest update of the ovine linkage map has been recently published and is available on the Australian Sheep Gene Mapping website http://rubens.its.unimelb.edu.au/~jillm/jill.htm. Several QTL studies have established independent linkage maps to position QTL, e.g. Beh et al. , Crawford et al. , Beraldi et al. , Murphey et al.  and Gutierrez-Gil et al. , using independent populations of Merino, Coopworth, Soay, Suffolk, and Churra sheep, respectively.
In sheep, growth rate and body mass represent economically important traits, which are under moderate genetic control and respond to directional selection . Despite extensive background information, relatively few QTL studies have been reported for growth in sheep and furthermore they have been mostly restricted to partial genome scans, limiting the discovery of and reports on new QTL. QTL studies contribute to the understanding of the genetic basis of a biologically complex trait such as growth because they can identify positional candidate genes. Walling et al.  have reported QTL affecting muscle depth and live weight at eight weeks of age in Texel sheep from partial genome scans in candidate gene regions on Ovis aries chromosome 2 (OAR2) and OAR18. Using candidate regions on OAR1, 2, 3, 5, 5, 6, 11, 18 and 20 in Suffolk and Texel commercial sheep populations, Wallinget al. [13, 14] have revealed suggestive QTL for body weight. Based on previous studies in sheep and other livestock species, McRae et al.  have analysed results of partial scans on selected autosomes (OAR1, 2, 3, 18 and 20) and identified QTL for body weight at eight and 20 weeks of age on OAR1. A whole genome linkage study, conducted in an Indonesian Thin Tail × Merino sheep population, has revealed QTL for birth weight on OAR5 and for body weight at yearling on OAR18 .
Combining results from QTL analyses in different livestock species and functional and positional candidate gene studies have shown that the myostatin gene on OAR2, the insulin-like growth factor-1 gene on OAR3, the callipyge gene and the Carwell rib eye muscling locus on OAR18 and the MHC locus on OAR20 are linked to growth or muscularity QTL in sheep and/or cattle [13, 17–29]. However, incomplete genome scans and positional candidate gene studies give an incomplete picture of the whole genome and of the location of growth and body weight QTL.
In this paper, we report the development of a framework map for male sheep, derived from a paternal half-sib design within an Awassi × Merino resource population. We use this map to search for putative QTL for growth rate and body weight in this resource population. In subsequent papers, we will report other putative QTL for economically important production traits such as milk yield and milk persistency, fleece/wool production, carcass characteristics, reproduction, behaviour, feed intake, and type traits. The range of phenotypes collected during this study is listed in the additional file 1.
As described by Raadsma et al. , a resource population from crosses between Awassi and Merino sheep was established to exploit the extreme differences between these two types of sheep in a range of production characteristics. Awassi sheep is a large-frame fat-tailed breed, which has its origins in the Middle East as a multi-purpose breed for milk, carpet wool and meat production and where it is dominant. From this source, the modern milking Awassi sheep was developed in Israel , which is the breed used in the present resource. Merino sheep is known for high-quality apparel wool but poor maternal characteristics . The Australian Merino breed, which is dominant in Australia, was derived from Spanish and Saxon Merinos crossed with meat breeds imported from Capetown and Bengal . Both super-fine and medium-wool Merinos were used in the present resource: they have a much smaller frame size than the milking Awassi breed and a very different fat distribution.
This resource population was developed in three phases, coinciding with different stages of research. A diagrammatic representation of the mating structure is shown in Figure 1 for one of the sire families and the other families have similar mating structures. In Phase 1, four sires from an imported strain of improved dairy Awassi , were crossed with 30 super-fine and medium-wool Merino ewes. Four resulting F1 sires (AM) were backcrossed to 1650 fine and medium-wool Merino ewes, resulting in approximately 1000 generation-2 (G2) backcrosses (AMM). In Phase 2, 280 AMM G2 ewes were mated to the four AM F1 sires so that matings were both within family (F1 sire mated with his daughters) and across families (F1 sire mated with daughters of other F1 sires) to produce approximately 900 G3 animals (AM_AMM). In Phase 3, 280 of the available G3 ewes were mated to three of the AM F1 sires (both within and across sire families) to produce G4a animals (AM_AM_AMM). In addition, four G3 males (each replacing one of the F1 sires) were mated to G3 ewes, resulting in 490 G4b animals (AM_AMM_AM_AMM). A total of 2,700 progeny were produced over 10 years, representing four generations. A broad range of phenotypes was collected from the progeny, as well as a DNA and tissue (blood, milk, fat, muscle, wool) repository for each available animal. In the initial QTL study reported here, only phenotypic and genotypic information from the G2 backcross progeny of the first F1 sire were analysed in detail, as this was the only family where a genome-wide scan was performed. The additional families will be used for confirmation of QTL effects and, when combined with high-density marker analysis, for fine mapping of confirmed QTL.
Progeny were reared in typical Australian paddock conditions for a NSW Southern Tablelands environment. Supplementary feeding occurred at times when feed availability from pasture was limited and corresponded to periods of negative growth (approximately 12 months of age). From 83 to 98 weeks (at which time the growth study was terminated), only the males were maintained on pasture as a single cohort till separate feed intake and carcass studies were undertaken. Ewes were relocated to a separate farm for lambing and milk recording.
DNA was extracted from blood using a modification of the protocol described by Montgomery and Sise . Purity of all extracted DNA was assessed by calculating the 260/280 nm ratios determined with an Eppendorf BioPhotometer. All DNA samples were dispensed to 96-well plates using a robotic workstation (Beckman Biomek 2000 with integrated MJ research DNA Engine PCR cycler).
Two hundred previously published polymorphic microsatellite markers covering all 26 autosomes were used in the construction of the map. They comprised 112 cattle (Bos taurus) markers, 73 sheep (Ovis aries) markers, and 15 other bovidae markers sourced from Prof. Yoshikazu Sugimoto (pers. comm.). All markers were screened for phase-known heterozygosity for the sire genotype. Markers were chosen on their Polymorphic Information Content  (PIC; > 0.6 if possible), and ease of scoring. Five hundred and ten animals were genotyped, comprising the Awassi grandsire, the Merino grand dam, and 510 AMM backcross G2 progeny (246 ewes and 264 wethers).
PCR was performed in 10 μL reactions containing 50 ng DNA, 1 × PCR buffer, 1 × 2.5 mM MgCl2, 200 μM of each dNTP, 0.8 pmol of each forward primer (with M13-29 tail) and reverse primer, 0.2 pmol of M13-29 primer labelled with either IRD 700 or IRD800 dye, and 0.5 units of Taq polymerase. PCR amplifications were carried out using one of the following three MJ Research (Watertown, Massachusetts, USA) 96 well PCR machines, namely, PTC-100, PTC-200, and PTC-200 Gradient Cycler.
The touchdown program (Licor-50) was used for the majority of the PCR, and a second program (Cav-low) was used for markers with a lower annealing temperature if amplification was unsuccessful using the Licor-50 program. The Licor-50 thermocycler touchdown cycles were as follows: initial denaturation for 5 min at 95°C, 5 cycles of 95°C for 45 s, 68°C for 1.5 min (-2°C per cycle), 72°C for 1 min, followed by 4 cycles of 95°C for 45 s, 58°C for 1 min (-2°C per cycle), 72°C for 1 min, followed by 25 cycles of 95°C for 45 s, 50°C for 1 min, 72°C for 1 min and a final 5 min extension at 72°C. The Cav-low cycles were as follows: initial denaturation for 5 min at 95°C, 5 cycles of 95°C for 30 s, 55°C for 1.5 min, 72°C for 45 s, followed by 5 cycles of 95°C for 30 s, 50°C for 30 s, 72°C for 45 s and a final 5 min extensions at 72°C.
Microsatellite PCR products were separated by polyacrylamide electrophoresis (PAGE) and detected using a Licor 4200 semi-automated sequencer.
Scoring of genotypes
The following description applies to the genotype scoring of the AMM backcross only as mentioned previously. All genotypes were scored by at least two independent scorers. To facilitate linkage analysis, only the F1 allele source was scored (Awassi or Merino origin), rather than the actual allele size. The Awassi allele was scored as '1', while the Merino allele was scored as '2', giving a genotype for the F1 sires. Only the identities of the alleles that were in the F1 sire were scored in the G2 AMM backcrosses, their genotypes identified as '1', '2' or '12'. A score of 1 can be homozygous '11' or 1x, where x is not equal to 2. Similarly a score of 2 can be homozygous '22' or 2x, where x is not equal to 1. Since information of the maternal allele was not available, heterozygous '12' in the backcross progeny was only semi-informative, as one cannot determine which allele originated from the F1 sire or from the Merino dam. The QTL mapping methodology used here exploited the semi-informative marker information (additional file 2).
Using the genotype information from our Awassi-Merino resource population, we generated an independent sheep linkage framework map comprising the 200 microsatellites genotyped in this resource. Carthagene version 4.0 [36, 37] and Multipoint http://www.multiqtl.com/ were used for the construction and validation of the map. These two programs use a multipoint maximum likelihood estimation method. Carthagene was used for the initial map construction, and Multipoint was used to test and validate marker orders. Only markers showing consistent results from both programs were included in the final framework map.
We used information from the Sheep Linkage Map v4.7 http://rubens.its.unimelb.edu.au/~jillm/jill.htm to group markers according to their chromosomal location as a prior to the construction of the framework map. Marker ordering and validation were performed for each linkage (chromosome) group separately. A minimum LOD score of 3.0 and a maximum recombination fraction of 0.4 were used as thresholds for linkage and sub-linkage grouping within the same chromosome. The Kosambi map function  was used to convert recombination fractions to distances. A framework map was considered satisfactory for the marker positions within a linkage group if the LOD score difference between the best and next-best map order was greater than or equal to 3.0.
Analysis of growth data
Non-fasted body-weight measurements were taken at weeks 2, 15, 25, 32, 37, 43, 48, 50, 56, 60, 67, 74, 79, and 83 for 510 G2 AMM backcrosses (246 ewes and 264 wethers). Birth weight was recorded for some animals, and body weights at weeks 90 and 98 were recorded for males only. The analysis of these data indicated distinct changes in growth rate at weeks 43, 56, and 86, presumably as a result of seasonal influences. Thus, growth rates were divided into four growth phases: week 0 to week 43, week 43 to week 56, week 56 to week 83, and week 83 to week 98. To accommodate these distinct changes, a piecewise-linear mixed model was used to model growth of each animal. Linear mixed models were fitted with separate slopes in each phase, but constrained to connect at each breakpoint (spline knot). While, arguably, a non-linear growth model may have been more applicable, the major purpose of the modelling was to capture the main features of the growth data. A full description of the piecewise-linear mixed model can be found in the additional file 2.
QTL mapping procedure
A maximum likelihood procedure, named QTL-MLE, suitable for the backcross design of the present resource (in which only the paternal allele was identified in G2 animals) was developed and programmed using R  by one of us (PCT). The software allows easy modification for the identification of QTL for most types of traits, including binary (e.g. disease presence-absence), ordinal (e.g. 5-point disease severity scale), or survival-time traits. Details of the algorithm are provided below, in terms of the models used to analyze body weight and growth data.
For a normally distributed trait, a linear model may be appropriate, i.e. y i = β'x i + γq i + ε i , where y i = observed trait value of animal i, i = 1, ... n; x i = set of covariates and fixed effects for animal i; β = corresponding set of regression parameters; γ = sire family allelic QTL effect (Q relative to q); q i = unobserved QTL allele of animal i, = 1 if Q, 0 if q; and ε i = random error, assumed N(0,σ2). Note the Merino dam effects will be absorbed into this last term. The genotype of the F1 sire is assumed to be Qq, with Q originating from the Awassi line and q from the Merino line.
Since there are only two types of QTL alleles in backcross animals, the phenotype distribution is a mixture of two distributions. We calculate the QTL transmission probability (π i ) as the probability of the sire transmitting QTL allele Q = π i = p(q i = 1 | m i ), while the probability of transmitting the other allele q is 1 - π i = p(q i = 0 | m i ), where m i is the "flanking" marker genotype information. Probabilities depend on the distance from the putative QTL to the marker(s) calculated via Haldane's mapping function. If the immediate flanking markers are "informative" (genotyped as '1' or '2'), they provide all possible information. Wherever a "semi-informative" marker ('12') is encountered adjacent to a putative QTL, the minimal set of markers that contains all the information for that QTL comprises the smallest set of contiguous markers flanked by "informative" markers.
At regular distances (typically 1 cM) along the length of the chromosome, the log-likelihood is constructed assuming a QTL at that position (d), i.e.
where f(·) is the probability density function (PDF) for a normal distribution (assuming that is the appropriate model for the data type). The log-likelihood is maximized using the E-M algorithm, which allows standard linear model software to be used, in an iterative manner. This requires computation at each iteration of the posterior probabilities (τ i ) that the sire transmits allele Q, conditional on its phenotype,
At the peak log-likelihood position (i.e. estimated QTL location), these τ i values can be used to classify backcross animals with high probability of having received the Q (or q) allele. Also at the peak, a 1-LOD support interval for estimated QTL position was determined by determining the range of map positions that are within one LOD of the peak.
Implementation of the program in R has the advantage that the QTL mapping procedure can be extended within other modelling and graphical capabilities of this package. For normally distributed traits, the linear model function lm() is used, and this easily allows model extension to include interactions between the QTL and other fixed effects, such as sex-specific QTL effects: most other QTL analysis programs do not allow such extensions. Another advantage of the R system is the relative ease to model traits of different types. This is achieved by changing only a few lines of code, primarily (1) replacing the lm() call by another function call, and (2) replacing the normal PDF in the τ i calculation (dnorm()) by the appropriate PDF (or discrete probability function) for the required distribution.
Using QTL-MLE, separate genome scans were conducted for single QTL on the bodyweights at the start and end of the four growth phases. For these traits, the model-based predictions from the piecewise-linear mixed model output were analysed rather than the raw data. The stages analysed were at weeks 2, 43, 56, 83, and 98. Note that week-2 bodyweights were selected in preference to week-0 (start of Phase I) due to the relatively few birth weights available. The model fitted to these values was as follows:
Weight i = model-based bodyweight at week i (2, 43, 56, 83, and 98);
Sex = 1 if ram/wether; 0 if ewe;
QTL = 1 if Awassi allele, Q; 0 if Merino allele, q (allele type is unobserved); and
ε = residual random error term.
Note that the unobserved QTL term is taken into account using the E-M algorithm of the interval mapping procedure. The interaction term was added to allow for sex-specific QTL effects.
Similarly, the average growth rates during each growth phase were analysed as separate traits. Again, model-based growth rates were used, as obtained from the piecewise-linear mixed model, and the model-based bodyweight at the start of each growth phase was used as a covariate. (As in the growth rate QTL model, the week-2 predicted bodyweights were used in preference to week-0 predicted ones). The model fitted for this QTL analysis took the following form:
GR i = model-based average growth rate in growth phase i and
Weight i = model-based bodyweight at start of growth phase i.
Since data for only wethers were available for the last growth phase (83–98 weeks), a term for sex was not included in either the week-98 body weight analysis, or the growth rate analysis. An additional series of analyses was performed without inclusion of the initial weight as a covariate.
Because of the large number of analyses, we adopted the false discovery rate (FDR) method of Benjamini and Hochberg  to adjust P-values for all traits to control for genome-wise error rates. Results were concluded to be significant when the adjusted P-values were less than 0.05. In all of these cases, LOD scores generated by QTL-MLE were larger than 2; QTL are described as suggestive where the F-value exceeds chromosome wide P < 0.05 threshold but not the 0.01 threshold. Based on a type I error of 0.01, the design had a power of 0.80 to detect QTL with 0.3 SD effect with 510 animals and an average marker spacing of 20 cM .
QTL mapping using QTL Express
For comparative purposes, all traits were analysed using the half-sib applet in QTL Express . With the exception of the QTL × fixed effect interaction, the same fixed effects as in the MLE analysis were fitted. Chromosome-wide significance thresholds were assessed using permutation tests , and bootstrap procedures  were used to obtain confidence intervals, both implemented in QTL Express using 1,000 re-samplings.
Methods for mapping a single QTL can be biased by the presence of other QTL [47, 48]. To address this situation, two-QTL models were also fitted for all traits using QTL Express . To control for false-positive QTL due to multiple testing, the permutation thresholds obtained in the single-QTL analyses were used to test for the significance of the two-versus one-QTL for a particular trait. Corresponding F-values for the two-versus zero-QTL test are included for comparison and additional support, although the same significance thresholds would not be applicable (given it would be a two numerator df test rather than a one df test).
Sheep framework map
From the 200 markers used, 194 markers showed significant linkage with at least one other marker at a LOD score of 3 or greater within their assigned linkage group (chromosome). The six markers that did not show significant linkage with other markers on their assigned chromosome were DIK4933 and OARFCB129 on OAR3, TGLA116 on OAR4, MCM185 on OAR7, BM6108 on OAR10 and RM024 on OAR24. All these markers were excluded from the framework map. A further three markers were excluded because their inclusion did not improve the overall LOD score of the framework map, even though they had a LOD of 3 or greater with one other marker within their linkage group. These three markers were KAP8 on OAR1, TGLA67 and OARFCB5 on OAR3. The final map contains 191 markers.
For the framework map, both Carthagene and Multipoint produced the same linkage and map order results. The additional file 3 presents the LOD score differences between the best and second-best map order for each chromosome generated by Carthagene. Except for OAR1, 2, 10 and 17, all other chromosomes yield a LOD score difference greater than 3.0 between the best and second-best map order. Thus the framework map can be considered fixed for the majority of the chromosomes. A detailed higher resolution order and length can be found in additional file 4.
In our framework map, we have also included four bovine microsatellite markers (DIK4572, DIK4527, DIK4612, and DIK2269) that are presently not included on the Sheep Linkage v4.7 Best Position Map. DIK4572 has been mapped to BTA2  and in the present study is placed on OAR2 with a two-point LOD score of 4.8 with its closest marker INRA135. DIK 4527, DIK4612 and DIK2269 all map on BTA20 , and in the present study are placed on OAR16 with respective two-point LOD scores of 28.2, 14.7 and 11.8 with their closest neighbouring markers. These bovine and ovine positions are consistent with the cattle-sheep comparative map as shown on the Sheep Linkage Map web site http://rubens.its.unimelb.edu.au/~jillm/jill.htm.
Apart from a slight difference in marker position, the marker order of the ReproGen Framework Map is the same as the Sheep Linkage Map Best Position Map v4.7. Sixteen chromosomes had a length at least a 7 cM greater than that in Sheep Linkage Map v 4.7, indicating slightly more recombination in the ReproGen map population. Six chromosomes (OAR4, 6, 12, 13, 23, 26) showed a similar length (within 3 cM) in both maps.
Overall growth performance
Table 1 presents the number of observations, the mean and the standard deviation of body weight at each of the measurement weeks. The plot of the weights (Figure 2A) indicates distinct changes at weeks 43, 56, and 86, suggesting growth phases. The fitted piecewise-linear mixed models for individual sheep are shown in Figure 2B.
All fixed effect terms in the piecewise-linear mixed model are significant (Table 2) indicating different growth profiles for both sexes, and support for the change in growth rate across the four phases. Table 2 also shows the estimated variance components, with their approximate standard errors. These represent individual animal variation in birth weights, and also in their individual growth rates, across the different phases.
Putative QTL identified for growth rate and body weight
Single QTL Analysis
Table 3 presents detailed results of the genome scan for QTL of body weight (BW) at the critical weeks separating the growth phases. Table 4 shows the corresponding information for growth rate (GR) during each of the four phases, whilst Table 5 shows the same information for growth rate traits, but after adjustment for body weight at the start of the growth phase. The 1-LOD support intervals generated by QTL-MLE are also reported. Figure 3 presents a QTL map showing the alignment of the QTL for all body weight traits along the genome, and Figures 4 and 5 show similar scans for growth rate QTL, unadjusted and adjusted for initial body weights. The additional file 5 contains all results using QTL-MLE and QTL Express showing the relative positions of the peaks along the genome for the different traits.
With the exception of BW02, QTL for body weight traits have been identified across the sheep genome (OAR1, 3, 6, 11, 21, 23, 24, and 26). Importantly, examination of the 1-LOD support intervals suggests that the same QTL are involved in various body weight traits (OAR3 for BW43, BW56, and BW83, OAR6 for BW43, BW56, and BW83, OAR11 for BW43, BW56, and BW83, OAR21 for BW43, BW56, and BW83 and OAR24 for BW43, and BW83). In addition, the QTL effects for males were almost always greater in absolute value than for females, and for males in particular, the effect of the Awassi allele led to an increase in body weight relative to the Merino allele.
Multiple QTL were also detected for the growth rate traits, and in general, these correspond to the QTL identified for the critical body weight traits, in terms of map position and also effect. All the body weight QTL also mapped to growth rate QTL, but in addition a suggestive QTL was found on OAR8 for GR00-43. While the growth rate QTL are in general the same as the body weight QTL, the analysis of growth rate QTL adjusting for the body weight at the start of the growth phase shows quite different results. Note that for the first growth phase, the body weight covariate adjusted for was BW02, since there were relatively few animals with birth weights data. After adjusting for initial body weight, QTL were identified for the first growth phase, GR00-43, corresponding to many of the regions previously identified for body weight and unadjusted growth rate traits, and an additional suggestive QTL was mapped on OAR16. However, no QTL were detected for GR43-56 after adjusting for BW43 (this period corresponding to a period of weight loss). Three QTL (on OAR3, 7 and 18) were detected for GR56-83, and only one QTL (on OAR1) for GR83-98.
Note that OAR1 is involved in body weight and growth rate QTL on three chromosomal locations, namely 32–68 cM (GR83-98 adj for BW83, positive effect of Awassi allele), 95–154 cM (BW43, GR00-43, both positive effects), 346–380 cM (BW83, GR43-56, GR56-83, GR00-43 adj for BW02, all negative effects).
Mapping results obtained by QTL Express were consistent with those obtained by QTL-MLE, particularly for those with greater effects (additional file 5). QTL Express also identified additional QTL on OAR6, 16 (GR02 in week 2) and OAR3 and 26 (GR4 in week 42) (but as noted earlier, it was not possible to fit sex-specific QTL effects in QTL Express).
Significant results for the two-QTL model are presented in Table 6. Overall, the two-QTL procedure detected far fewer QTL compared with the single-QTL methods, as QTL were detected for only three traits. For adjusted GR56-83, two QTL were detected in coupling phase on OAR3, one at 104 cM and the other at 284 cM, both with effects of ~0.35 SD (Table 6), the 3D surface plot of the variance ratio test statistic is shown in the additional file 6. In the single-QTL analyses, the second peak at 284 is clearly visible on the interval map. For the adjusted GR00-43, two QTL were mapped in very close proximity on OAR4 (108 cM and 112 cM) in repulsion phase. However, given that both these positions are flanked by the same markers (OARCP26 and OARHH35), and that both QTL have extremely large estimated effects (~ ± 2.8 SD), it is unlikely that these are real QTL. For adjusted GR83-98, a similar situation occurred on OAR4 (24 cM and 28 cM), and on OAR22 for the same trait (68 cM: -1.75 SD, and 88 cM: +1.90 SD) with the QTL being mapped to separate marker bracket intervals.
This paper reports the construction of a male distance framework map for sheep and its application in the identification of QTL for body weight and growth. There are several advantages in developing a separate framework map. First, it provides an independent verification to the Sheep Linkage Map v4.7, since it originates from a single sheep resource. It would be possible to integrate the data in this map with the data of other Sheep Linkage Maps to create an integrated framework map for sheep. The ReproGen framework map and the Sheep Linkage Map v4.7 agree well, with no changes in marker order. With the exception of OAR4, 6, 12, 13, 14, 18, 23 and 26, the ReproGen map is 11% larger than the male-averaged Sheep Linkage Map v4.7. Perhaps one explanation for this is that the current map was constructed using a single sire family with a large number of progeny, compared with the Sheep Linkage Map v4.7, which is derived from a large number of smaller families. The second advantage in constructing an experiment specific framework map is that the QTL can be unambiguously mapped to a fixed location since the markers are in a fixed order. The use of a framework map not only allows integration of markers in a consensus map, but also alignment of QTL in integrated maps for future meta-analyses such as those undertaken in dairy cattle by Khatkar et al. .
The pattern of growth in this flock is consistent for sheep maintained on semi-improved pasture in a temperate Australian tablelands climate. Such grazing systems are characterized by low pasture availability in the colder winter months, and abundant pasture growth in spring with a residual pasture carry over in summer and autumn. The availability of pasture is reflected in the growth curve over a period of 52 weeks, with a rapid growth following birth in spring till the end of autumn, a period of no growth or decline coinciding with winter, and a subsequent cycle between 52 and 98 weeks.
Consistent with findings in previous studies using European sheep, we have identified a number of significant QTL for body weight and growth rates on different chromosomes. In order to minimize the large number of possible QTL detected for single-point estimates of body weight, all data were combined through a growth curve over a period where changes in body weight were similar for all sheep as shown in the stick-point graph (Figure 2B). The points of change in growth (break points) were identified by single-point body weights, which were used in the QTL analyses as reference body weights. Since growth between each break point was strongly influenced by starting body weights at each time, true growth rate was analysed by adjusting for starting body weights. The final outcome of summarizing all body weights in relatively few growth and body weight indicators was that 17 body weight time points were collapsed to 12 core traits instead of the 49 possible correlated traits. It allows for a greater consistency in QTL reporting for traits related to growth.
Despite the economic importance of body weight and growth rates in sheep and the extensive studies reporting genetic variation in this trait [12, 51, 52], relatively few QTL studies have been reported to date. In this study, we report 54 QTL involved in body weight and growth rate. The majority of the QTL (49 out of 54) are co-located across nine chromosomal regions, suggesting QTL with a general effect on growth and body weight (OAR1, 3, 6, 11, 21, 23, 24 and 26). Only five QTL (OAR7, 8, 12, 16 and 18) were observed for a single growth trait, and notably four of these were for growth rate adjusted for starting body weight. Our study identified new QTL for body weight/growth (OAR6, 7, 8, 11, 12, 16, 21, 24 and 26) since our study covers a full autosomal genome scan as distinct from the previous partial genome scans, to the best of our knowledge. Our study confirmed QTL for growth and body weight previously reported by Walling et al. [13, 14], and McRae  and slaughter live weight  on OAR1, 3, and 18. In our study, no QTL for growth was observed on OAR2, and OAR5 where Walling et al. , Karamichou et al.  and Margawati et al.  reported QTL.
The effects of the QTL detected in this study are relatively large (in the order of 0.4 to 0.7 phenotypic standard deviations) and originated predominantly from the larger-framed Awassi grandsire, as expected with the favourable allele for growth and body size. In five out of 54 cases, the reverse was observed with the favourable allele originating from the Merino breed, suggesting the presence of cryptic QTL. These estimated QTL effects are consistent with the significant genetic variation in growth that has been observed within this breed (e.g. Merino body weight heritabilities are in the range of 0.39–0.76, based on a recent comprehensive study by Huisman et al. . The current study has the power to detect QTL effects of 0.3 SD or larger. However, surprisingly for most traits, significant sex by QTL interactions were observed, with most QTL expressed in males but not females (Tables 3, 4, 5). No immediate or obvious explanation can be given for this, since the ewes and wethers were managed together for most of the period. A notable feature of this study was that the QTL locations have relatively small 1-LOD intervals, most likely due to the relatively large family size. This will facilitate positional candidate gene analyses and make a significant contribution to possible future meta-analyses. Some of the QTL reported in this paper are consistent with possible candidate genes that have been suggested for growth and body weight. The additional file 7 shows a comparison between the QTL region detected here with linkage and association studies of other groups as well as possible candidate genes possibly included in the regions.
The maximum likelihood approaches to QTL detection employed in this analysis were in general robust and their results were in good agreement with those of QTL Express , which is based on the least-squares methodology. Overall, 33 QTL were detected by both methods, with an additional 20 QTL identified by QTL-MLE alone, and another nine detected by QTL Express alone. However of the additional 29 QTL detected by only one of these methods, the majority (15/20 for QTL-MLE and 9/9 for QTL Express) had suggestive significance levels, as indicated by P-values for each method. However, it should be pointed out that the QTL-MLE could detect sex-specific QTL, which may have increased the power of the analysis. This is considered a notable advantage of our approach in that QTL × fixed effects can be fitted, which at the time of development was not possible in QTL Express. Obviously, this is a major consideration in our study since most QTL for growth and body weight appear to have differential effects in both sexes. An additional benefit of QTL-MLE is the flexibility of analysis offered, due to working in the R environment, with the ready potential for further enhancements in the future. However, QTL Express has the benefit of being able to fit two-QTL models. While three chromosomes containing two QTL were identified, two were unlikely given the size of the effect estimates, but strong support is indicated for two QTL for growth rate between weeks 56 and 83, adjusted for week-56 body weight on OAR3.
In conclusion, we present significant evidence for an independent framework map, which is in very good agreement with the previously published Sheep Linkage Map v4.7 framework map for sheep. We also present evidence for a significant number of new QTL for body weight and growth rate in sheep, and confirm and support the presence of previously published QTL in breeds other than those studied here.
As indicated earlier, a large range of traits have been recorded from this sheep mapping study. Subsequent papers will describe QTL for these traits, using the current paper as a foundation to describe the overall approach. As SNP chips become more widely available for research in sheep genetics, this series of papers should provide a reference point to QTL studies in sheep using microsatellite markers.
Weller JI: Quantitative trait loci analysis in animals. 2001, New York: CABI Publishing
Crawford AM, Dodds KG, Ede AJ, Pearson CA, Garmonsway HG, Beattie AE, Montgomery GW, Nguyen TC, Davies K, Maddox JF: An autosomal genetic linkage map of the sheep genome. Genetics. 1995, 140 (2): 703-724.
Crawford AM, Dodds KG, Ede AJ, Pierson CA, Montgomery GW, Garmonsway HG, Beattie AE, Davies K, Maddox JF, Kappes SW: An Autosomal Genetic-Linkage Map of the Sheep Genome. Genetics. 1995, 140 (2): 703-724.
Gortari MJd, Freking BA, Cuthbertson RP, Kappes SM, Keele JW, Stone RT, Leymaster KA, Dodds KG, Crawford AM, Beattie CW: A second-generation linkage map of the sheep genome. Mamm Genome. 1998, 9 (3): 204-209. 10.1007/s003359900726.
Maddox JF, Davies KP, Crawford AM, Hulme DJ, Vaiman D, Cribiu EP, Freking BA, Beh KJ, Cockett NE, Kang N: An enhanced linkage map of the sheep genome comprising more than 1000 loci. Genome Res. 2001, 11 (7): 1275-1289. 10.1101/gr.GR-1350R.
Maddox JF, Cockett NE: An update on sheep and goat linkage maps and other genomic resources. Small Rum Res. 2007, 70 (1): 4-20. 10.1016/j.smallrumres.2007.01.008.
Beh KJ, Hulme DJ, Callaghan MJ, Leish Z, Lenane I, Windon RG, Maddox JF: A genome scan for quantitative trait loci affecting resistance to Trichostrongylus colubriformis in sheep. Anim Genet. 2002, 33 (2): 97-106. 10.1046/j.1365-2052.2002.00829.x.
Crawford AM, Paterson KA, Dodds KG, Tascon CD, Williamson PA, Thomson MR, Bisset SA, Beattie AE, Greer GJ, Green RS: Discovery of quantitative trait loci for resistance to parasitic nematode infection in sheep: I. Analysis of outcross pedigrees. BMC Genomics. 2006, (7): 178-10.1186/1471-2164-7-178.
Beraldi D, McRae AF, Gratten J, Slate J, Visscher PM, Pemberton JM: Development of a linkage map and mapping of phenotypic polymorphisms in a free-living population of Soay sheep (Ovis aries). Genetics. 2006, 173 (3): 1521-1537. 10.1534/genetics.106.057141.
Murphy AM, MacHugh DE, Park SDE, Scraggs E, Haley CS, Lynn DJ, Boland MP, Doherty ML: Linkage mapping of the locus for inherited ovine arthrogryposis (IOA) to sheep Chromosome 5. Mamm Genome. 2007, 18 (1): 43-52. 10.1007/s00335-006-0016-8.
Gutierrez-Gil B, Arranz JJ, El-Zarei MF, Alvarez L, Pedrosa S, Primitivo FS, Bayon Y: A male linkage map constructed for QTL mapping in Spanish Churra sheep. J Anim Breed Genet. 2008, 125 (3): 201-204. 10.1111/j.1439-0388.2008.00729.x.
Safari E, Fogarty NM, Gilmour AR: A review of genetic parameter estimates for wool, growth, meat and reproduction traits in sheep. Livest Prod Sci. 2005, 92 (3): 271-289. 10.1016/j.livprodsci.2004.09.003.
Walling GA, Visscher PM, Simm G, Bishop SC: Confirmed linkage for QTLs affecting muscling in Texel sheep on chromosomes 2 and 18. Proc 52nd Ann Meeting Europ Assoc Anim Prod, 26–29 August 2001, Budapest, Hungary, Paper G56: 2001. 2001
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 (8): 2234-2245.
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. 10.1079/ASC41040135.
Margawati ET, Raadsma HW, Martojo H, Muladmo S: Quantitative Trait Loci (QTL) Analysis for Production Traits of Birth Weight and Weight 360 days in Backcross Sheep. Hayati J Biosci. 2006, 13 (1): 6-
Broad TE, Glass BC, Greer GJ, Robertson TM, Bain WE, Lord EA, McEwan JC: Search for a locus near to myostatin that increases muscling in Texel sheep in New Zealand. Proc New Zeal Soc Anim Prod. 2000, 60: 110-112.
Kmiec M: Transferrin Polymorphism versus Growth Rate in Lambs, Polish Long-wool Sheep – II. Analysis of relation between transferrin polymorphism of lamb blood serum versus growth rate of lambs up to age of 5 months. Arch Tierzucht. 1999, 42 (5): 469-479.
Marcq F, Elsen JM, El Barkouki S, Bouix J, Eychenne F, Grobet L, Karim L, Laville E, Nezer C, Royo L: Investigating the role of myostatin in the determinism of double muscling characterizing Belgian Texel sheep. Anim Genet. 1998, 29 (Suppl 1): 52-
Moody DE: Phenotypic variation and identification of quantitative trait loci for traits influencing energy balance in mice. Dissertation. 1998, University of Nebraska
Moody DE, Pomp D, Newman S, MacNeil MD: Characterization of DNA polymorphisms in three populations of hereford cattle and their associations with growth and maternal EPD in line 1 herefords. J Anim Sci. 1996, 74 (8): 1784-1793.
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 (6): 1379-1384.
Bot J, Karlsson LJE, Greef J, Witt C: Association of the MHC with production traits in Merino ewes. Livest Prod Sci. 2004, 86 (1–3): 85-91. 10.1016/S0301-6226(03)00146-5.
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 (4): 854-860.
Cockett NE, Jackson SP, Shay TL, Nielsen D, Moore SS, Steele MR, Barendse W, Green RD, Georges M: Chromosomal localization of the callipyge gene in sheep (Ovis aries) using bovine DNA markers. PNAS. 1994, 91 (8): 3019-3023. 10.1073/pnas.91.8.3019.
Freking BA, Murphy SK, Wylie AA, Rhodes SJ, Keele JW, Leymaster KA, Jirtle RL, Smith TPL: Identification of the single base change causing the callipyge muscle hypertrophy phenotype, the only known example of polar overdominance in mammals. Genome Res. 2002, 12 (10): 1496-1506. 10.1101/gr.571002.
Nicoll GB, Burkin HR, Broad TE, Jopson NB, Greer GJ, Bain WE, Wright CS, Dodds KG, Fennessy PF, McEwan JC: Genetic linkage of microsatellite markers to the Carwell locus for rib-eye muscling in sheep. Proc 6th World Cong Genet Appl Livest Prod, Armidale, Australia. 1998, 26: 529-532.
Paterson AH, DeVerna JW, Lanini B, Tanksley SD: Fine Mapping of Quantitative Trait Loci Using Selected Overlapping Recombinant Chromosomes, in an Interspecies Cross of Tomato. Genetics. 1990, 124 (3): 735-742.
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. Proc 7th World Cong Genet Appl Livest Prod, Montpellier, 19–23 August 2002, Inra, Castanet-Tolosan, France, Abstract 02–16. 2002, 0-4.
Raadsma HW, Nicholas FW, Jones M, Attard G, Palmer D, Abbott K, Grant T: QTL mapping for production traits in amulti-stage AwassixMerino back-inter-cross design. Proc Assoc Advmt Anim Breed Genet. 1999, 13: 361-364.
Gootwine E, Goot H: Lamb and milk production of Awassi and East-Friesian sheep and their crosses under Mediterranean environment. Small Rum Res. 1996, 20 (3): 255-260. 10.1016/0921-4488(95)00807-1.
Bencini R, Pulina G: The quality of sheep milk: a review. Wool Tech Sheep Breed. 1997, 45 (3): 182-220.
McLeish W: The beginnings of the sheep and wool industry in Australia. 2004, [http://www2.dpi.qld.gov.au/sheep/6570.html]
Montgomery GW, Sise JA: Extraction of DNA from Sheep White Blood-Cells. NZ J Agr Res. 1990, 33 (3): 437-441.
Botstein D, White RL, Skolnick M, Davis RW: Construction of a Genetic-Linkage Map in Man Using Restriction Fragment Length Polymorphisms. Am J Hum Genet. 1980, 32 (3): 314-331.
de Givry S, Bouchez M, Chabrier P, Milan D, Schiex T: CAR(H)(T)AGene: multipopulation integrated genetic and radiation hybrid mapping. Bioinformatics. 2005, 21 (8): 1703-1704. 10.1093/bioinformatics/bti222.
Schiex T, Gaspin C: CARTHAGENE: constructing and joining maximum likelihood genetic maps. Proc Int Conf Intell Syst Mol Biol. 1997, 5: 258-267.
Mester DI, Ronin YI, Nevo E, Korol AO: Fast and high precision algorithms for optimization in large-scale genomic problems. Comput Biol Chem. 2004, 28 (4): 281-290. 10.1016/j.compbiolchem.2004.08.003.
Kosambi DD: The estimation of map distances from recombination values. Ann Eugen. 1944, 12: 172-175.
R Development Core Team: R: A language and environment for statistical computing. 2008, R Foundation for Statistical Computing V, Austria, [http://www.R-project.org/]
Jansen RC: A General Mixture Model for Mapping Quantitative Trait Loci by Using Molecular Markers. Theor Appl Genet. 1992, 85 (2–3): 252-260.
Benjamini Y, Hochberg Y: Controlling the False Discovery Rate – a Practical and Powerful Approach to Multiple Testing. J R Stat Soc B. 1995, 57 (1): 289-300.
Werf van der J: 2009, [http://www-personal.une.edu.au/~jvanderw/software.htm]
Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in of simple and complex pedigrees. Bioinformatics. 2002, 18 (2): 339-340. 10.1093/bioinformatics/18.2.339.
Churchill GA, Doerge RW: Empirical Threshold Values for Quantitative Trait Mapping. Genetics. 1994, 138 (3): 963-971.
Visscher PM, Thompson R, Haley CS: Confidence Intervals in QTL Mapping by Bootstrapping. Genetics. 1996, 143 (2): 1013-1020.
Meuwissen THE, Karlsen A, Lien S, Olsaker I, Goddard ME: Fine mapping of a quantitative trait locus for twinning rate using combined linkage and linkage disequilibrium mapping. Genetics. 2002, 161 (1): 373-379.
Meuwissen THE, Goddard ME: Mapping multiple QTL using linkage disequilibrium and linkage analysis information and multitrait data. Genet Sel Evol. 2004, 36 (3): 261-279. 10.1051/gse:2004001.
Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, Hirano T, Itoh T, Watanabe T, Reed KM: A Comprehensive Genetic Map of the Cattle Genome Based on 3802 Microsatellites. Genome Res. 2004, 14 (10a): 1987-1998. 10.1101/gr.2741704.
Khatkar MS, Thomson PC, Tammen I, Raadsma HW: Quantitative trait loci mapping in dairy cattle: review and meta-analysis. Genet Sel Evol. 2004, 36 (2): 163-190. 10.1051/gse:2003057.
Snyman MA, Olivier JJ, Erasmus GJ, vanWyk JB: Genetic parameter estimates for total weight of lamb weaned in Afrino and Merino sheep. Livest Prod Sci. 1997, 48 (2): 111-116. 10.1016/S0301-6226(96)01418-2.
Huisman AE, Brown DJ, Ball AJ, Graser HU: Genetic parameters for bodyweight, wool, and disease resistance and reproduction traits in Merino sheep. 1. Description of traits, model comparison, variance components and their ratios. Aust J Exp Agr. 2008, 48 (9): 1177-1185. 10.1071/EA08119.
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.
We would like to thank Stephen Brown for assistance with the development of components of the QTL-MLE R coding, and Sally Isberg for conducting testing of the QTL-MLE program and the initial QTL analyses. Animals were housed and maintained in accordance with University of Sydney Animal Ethics Committee specifications. The foundation of this project (1996–1999) was supported with an ARC linkage grant and support from Awassi Australia. Research during 2002 and 2007 was supported by the Cooperative Research Centre for Innovative Dairy Products, Melbourne Australia. The support given by Tom Grant, Phillip Grant and Dr Evan Hunt during the establishment phase is gratefully acknowledged.
The authors declare that they have no competing interests.
HR drafted the overall design, did the project management, and was involved in analysis of the data and writing the manuscript. PT developed the statistical methodology for the growth trait analysis and QTL methodology, implemented the QTL-MLE program, and contributed to manuscript preparation and the overall design. KZ and CC carried out the genetic marker analysis, genetic map construction, and ran the early stage analyses. ML ran the early stage QTL analyses, was responsible for the data assembly, and participated in the growth curve analyses. EJ helped run the QTL and data analyses, participated in the manuscript preparation and final response to referees. MJ and GA performed the genetic marker analyses, and did the DNA bank assembly. DP participated in the analysis, he further collected and provided the phenotypes and was responsible for the experimental animals. FN supervised the study, coordinated the design and initiation of the project and contributed to the manuscript preparation. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Traits measured in the Awassi-Merino resource flock. A list of all traits included the traits used for the analysis in the present paper, which were recorded in the sheep resource population used in the present study. (DOC 50 KB)
Additional file 2: Methodology QTL-MLE. A detailed description of the QTL-MLE methodology which was developed by PCT and was used for the analysis shown here. (DOC 69 KB)
Additional file 3: LOD score difference between the best and second best map order. This tables shows the differences of the LOD scored between the best and the second best map order, shown are the results on each autosome. (DOC 32 KB)
Additional file 4: Marker information for the Awassi × Merino map. The figures shown here summarize the markers used in the study and compare the framework map described in the Awassi × Merino resource population with the published map. Further the identified linkage regions from the present study and QTL and candidate genes from various references are shown. All figures were designed using the MapChart software described by Voorrips R.E. in 2002 (MapChart: Software for the graphical presentation of linkage maps and QTLs. The Journal of Heredity 93 (1): 77–78). (DOC 160 KB)
Additional file 5: Summary of results using QTL-MLE and QTL Express. Values in table are estimated QTL positions (cM). This table shows the average QTL positions of all QTL presented here using both programs QTL-MLE and QTL Express. It summarizes and compares all results observed in the resource population for the growth and weight traits. (DOC 100 KB)
Additional file 6: 3D surface plot of the variance ratio test statistic for a two-QTL model. This figure shows the 3D surface plot for a two-QTL model for GR56-83 adj for BW56 on OAR3. The values of the plot are generated by QTL Express. (DOC 232 KB)
Additional file 7: Comparison of the QTL analysis results in the present studies with other results and possible candidate genes. This table describes QTL identified in other studies which are located in comparable regions to the linkage regions identified using animals of the Awassi × Merino population. Further possible candidate genes are described. The list of references for these data is attached following the table. (DOC 72 KB)
About this article
Cite this article
Raadsma, H.W., Thomson, P.C., Zenger, K.R. et al. 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 41, 34 (2009). https://doi.org/10.1186/1297-9686-41-34
- Quantitative Trait Locus
- Quantitative Trait Locus Analysis
- Quantitative Trait Locus Effect
- Quantitative Trait Locus Allele
- Putative Quantitative Trait Locus