Email updates

Keep up to date with the latest news and content from Genetics Selection Evolution and BioMed Central.

Open Access Research

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Elisabeth Jonas, Peter C Thomson, Evelyn JS Hall, David McGill, Mary K Lam and Herman W Raadsma*

Author Affiliations

ReproGen-Animal Bioscience Group, Faculty of Veterinary Science, University of Sydney, 425 Werombi Road, Camden NSW 2570, Australia

For all author emails, please log on.

Genetics Selection Evolution 2011, 43:22  doi:10.1186/1297-9686-43-22

The electronic version of this article is the complete one and can be found online at:

Received:15 December 2010
Accepted:21 June 2011
Published:21 June 2011

© 2011 Jonas et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.


An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.


Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.


This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.


Lactation performance plays an important role for the productivity and therefore economic value of dairy cattle and dairy sheep. Persistency of lactation is a trait of considerable importance as it reflects the ability of an animal to maintain milk production at a high level after the peak yield [1]. It has been found that cows with greater lactation persistency tend to incur lower feed, health, and reproductive costs [2], and cows affected with mastitis tend to have less persistent lactations [3]. The animal health benefits of increased lactation persistency are speculative at present; but it has been hypothesised that increased lactation persistency is associated with lower incidences of peri-parturient metabolic diseases [4]. With the recent interest in year-round calving as opposed to the conventional 305-day cycle in dairy cattle, identifying the genetic basis for lactation persistency and extended lactation is of interest.

Lactation curve models have had a long history in dairy science, dating back at least to Brody [5]. Since then, many models have been proposed, including the commonly used model by Wood [6], which differ both in the mathematical form of the function and the number of model parameters that need to be estimated. More recently, random regression approaches have been proposed, which allow a model to be fitted to the yield data in which the data itself specifies the shape of the fitted model, rather than fitting to a pre-specified mathematical function [7,8]. However, random regression models require a large number of lactations with many yields recorded per lactation, which may be a problem in voluntary recording systems. Furthermore, extraction of summary information from fitted random regression models is difficult.

As well as being able to predict lactation yield at any stage of lactation, lactation curve models can be used to summarise key features of the lactation curve, such as time of peak yield and the peak flow, and characteristics related to extended lactation. One of the difficulties with some of these lactation length measures is the variety of definitions provided in the literature (e.g. Sölkner and Fuchs [2]), although there are some common threads. Lactation persistency is usually defined as the (rate of) decline of milk yield after the peak [6,9]. Lactation persistency has also been defined as the difference in yield between peak to a defined day later in lactation [10]. However, Grossman et al. [11] developed a new model for lactation curves which allows the definition of persistency as the number of days at peak lactation. However, the models developed by Grossman et al. [11] implied a lactation curve with a plateau, which is not the case for the typical lactation curve in sheep which has a relatively early but sharp peak. These derived measures of lactation persistency can then be used as traits in QTL mapping studies [12,13].

Compared with lactation persistency, extended lactation has received relatively less attention in the literature. Extended lactation deals with the ability to maintain lactation at productive levels beyond the usual "drying off' period or a fixed reference day, i.e. beyond 305 days in cattle, or 100 days in dairy sheep. Both parametric and semi-parametric models can be used to quantify extended lactation. They all consider the level of milk production beyond the reference day. Some studies fit models to empirical data from extended lactation records [14], whilst others use fitted models to predict the potential for extended lactation in the absence of actual data records beyond the reference date. Vargas et al. [15] analysed the accuracy of different models to predict daily milk yield in standard and extended lactations in cattle. Dematawewa et al. [16] compared empirical (including the Wood model) and mechanistic lactation models for their suitability to predict standard (305-day) and extended (999-day) lactations in cattle (with 0.1% of the cows having actual data on the extended lactation). They recommended the use of the model of Rook et al. [17] or the Wood model as being appropriate to describe extended lactations including milk, fat and protein yields [16]. Previously, we showed that the Wood model was appropriate to describe lactation curves in sheep [18], so here we extend the use of this model to predict lactation persistency and extended lactation in sheep.

To our knowledge, there are no reported heritability estimates for lactation persistency in sheep, nor for extended lactation in sheep or cattle. Reported heritability estimates for milk persistency in cattle were between 0.02 and 0.27 [19]. Heritability estimates for milk composition and somatic cell count persistency were between 0.03 and 0.2 [4,20]. Thus, under conventional mass selection, the rate of gain for persistency of milk yield and milk composition is expected to be low.

Growing intensification of sheep and goat production systems during the recent decades have led to increased research interests into the yield and composition of milk from small ruminants [21-23]. Milk yield and milk components (%) are negatively correlated in sheep and these correlations change during the lactation [24,25]. It has been shown that the quality of cheese processing also declines during the course of lactation in sheep with poor lactation persistency [26]. Changes in milk composition of dairy sheep throughout the lactation typically result in variable cheese yields and quality [27]. Selecting animals with improved lactation persistency should therefore also aim at more persistent milk composition and yield during the course of the lactation. This would benefit the use of milk for processing into cheese and other products.

Until now, only a few linkage studies have been published for traits related to persistency of milk yield and milk composition in cattle but none have been reported in sheep. QTL were reported for milk, fat, protein or energy persistency [10,12,28] and shape of the lactation curve [29] in cattle. Comparison of results across studies is problematic as there is no universally accepted method to calculate persistency, and results are therefore method-dependent. The objective of this study was to report on an appropriate model to predict lactation persistency and extended lactation for milk yield and milk composition traits in sheep, and to provide additional information on QTL for these traits. The application of QTL information for these traits is of potential interest in marker-assisted selection since these traits have sex-limited expression, and are potentially difficult to measure, thus limiting the rate of progress under conventional mass selection.


Resource population and phenotypes

A resource population from crosses between the improved dairy type of Awassi (A) and the apparel wool Merino (M) sheep was established as described in detail by Raadsma et al. [30]. Genotypic information from the backcross progeny of the first sire was used for linkage mapping.

Lactation and milk composition records were collected for 172 backcross ewes between 1999 and 2007 [18]. Milk yield (MY) was recorded every second day, and protein, fat, lactose percentages and somatic cell count (SCC) were analysed from milk sampled once weekly. Somatic cell score (SCS) was calculated as the logarithm to base 10 of the raw SCC. The yields of protein (PY), fat (FY), lactose (LY) and somatic cell score (SCY) were calculated by multiplication of the contents (expressed as a percentage) with MY. Useful yield (UY) was calculated as FY + 1.85 × PY [31]. Lactation curves for MY, PY, FY, LY, UY and SCY were estimated by fitting Wood models [32] to the yield data. Due to its simplicity and reliable fit to the data, the Wood model was chosen as the preferred method to summarise lactation curve characteristics. The basic form of the Wood model is defined as follows:

where W(t) is the expected yield at time t, and k = ln(a), b and c are parameters controlling the shape of the curve. Specifically, a is related to the total area under of the curve, b is related to the sharpness of the early rise, and c describes the decline rate in milk production. Multiple lactation curves were fitted simultaneously using the nlme function in R, treating each ewe by lactation as a random effect. Further description of the Wood model parameters and procedure for model fitting is in Raadsma et al. [18].

From this model, the cumulative yield up to day T (e.g. 100 days) (CumY(T)) can be calculated numerically as , where γ(·) is the lower incomplete gamma function, . Furthermore, the day of the maximum yield (tmax) and the maximum yield (maxY) were calculated as tmax = b/c and maxY = a(b/c)ae-b respectively.

Milk persistency was defined as the expected milk yield on day T, relative to that on the peak day, namely PersY(T) = W(T)/maxY. A more persistent lactation will have a flatter curve, with the persistency proportion approaching one. In the same way, derived persistency variables were also calculated for PY, FY, LY, UY and SCY; furthermore, lactation persistency is defined as a general term for milk as well as persistency of milk composition.

A previous study showed that 100 day lactation performance, measured as cumulative milk yield until day 100, is highly correlated with cumulative milk yield to day 300, whereas lactation performance taken as either cumulative milk yield to day 50 or 80 was unreliable in predicting cumulative lactation performance [18]. Furthermore, the number of observations in this data set was still very high (n = 565) at day 100 of lactation, which reinforced the choice of evaluating persistency at T = 100 days [see Additional file 1].

Additional file 1. Number of observations at different stages of the lactation. Number of observations during the lactation: the number of observations dropped from n = 565 at day 100 to n = 334 at day 150 and n = 172 at day 200 of lactation.

Format: PDF Size: 21KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Extended lactation deals with the ability of the ewe to sustain milk production beyond a certain time. In the absence of a standard lactation length, as in cattle and as outlined in paper II of this series [18], we have adopted day 100 as the standard reference for lactation length in sheep (in cattle this would nominally be day 305), even though this is shorter than most of the average ovine lactation lengths in the literature [33,34]. The definition of extended lactation adopted here is the ratio of expected production from day 100 to 300, relative to the cumulative production (of MY, PY, FY, LY, UY and SCY) up to day 100, i.e. [CumY(300)-CumY(100)]/CumY(100). The greater this ratio is, the more extended the lactation is. Note that the cut off at day 300 is arbitrarily made in order not to have an "infinite" lactation length.

QTL mapping procedure

A genome-scan using 189 polymorphic microsatellite markers covering all 26 autosomes was conducted on 172 backcross ewes. The procedures of DNA extraction, genotyping, allele calling, and the linkage map have been described previously [30].

Based on a Type I error of 0.05, the design had a predicted power of 0.72 to detect QTL with 0.4 SD effect [35]. To achieve normality, all traits were log-transformed prior to using two methods for QTL analyses. Solutions were first obtained using a maximum likelihood procedure, named QTL-MLE in R [30,36], which is suitable for the analysis of a backcross population. For QTL-MLE, a LOD of 1.75-2.0 was deemed suggestive, LOD 2.0-3.0 significant, and LOD > 3.0 highly significant. The 1-LOD drop-off method was used to calculate the respective confidence intervals for this method. QTL effects were derived from the model and normalised against a model with no QTL to express the size of QTL effects in units of phenotypic standard deviations.

The second method for QTL detection was based on the regression analysis for a half-sib design using QTL Express [37]. The half-sib model was used since no maternal genotypes were available. For this method, QTL with chromosome-wide significance thresholds (P < 0.05) were described as suggestive, chromosome-wide levels P < 0.01 as significant and experiment-wide levels (P < 0.05) as highly significant QTL. A two-QTL model was also fitted to all data using the same program [37] and this was conducted over all chromosomes. 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. The QTL heritability was calculated as the proportion of the phenotypic variance accounted for by the QTL, and this is calculated from the residual mean squares from the regression analysis of fitting a QTL or no QTL as [1-(residual mean square of full model/residual mean square of reduced model)].


To facilitate a comparative genome analysis between sheep and cattle, individual QTL locations and bovine QTL were extracted from the literature, and were loaded into the ovine genome database, which can be browsed at webcite. While a meta-assembly of all known QTL was performed as outlined in previous papers [18,38], no meta-score was calculated for lactation persistency due to the inconsistency of trait definitions of persistency and the inconsistency of QTL positions. However, as publications on persistency become available in the future, we plan to update the ovine genome database to facilitate a full meta-assembly.


Analysis of lactation persistency

The Wood model parameters were standardised against the predicted lactation curves for the following fixed effect classifications: five years of age, singleton birth type, third parity, twice milking/day and fourth season for milk yield; five years of age, second parity, and twice milking/day for protein yield; three years of age, second parity, twice milking/day and singleton birth type for fat yield; three years of age, twice milking/day, singleton birth type and second parity for useful yield; and seven years of age for lactose and somatic cell score yield. The factors used for standardisation differed across the lactation traits, as not every factor had a significant effect on all traits. Lactation persistency was highest for useful yield (0.65), while daily milk yield persistency had the lowest average value (0.23) (Table 1). Similarly, extended lactation had the highest average for useful yield, followed by fat and protein, and was lowest for milk yield (Table 1). Most of the lactation persistency characteristics for different milk traits showed low to moderate phenotypic correlations, with the exception of persistency of protein yield, which was highly correlated to persistency of fat and somatic cells (Table 2). The highest correlations between lactation persistency and extended lactation were observed for milk, lactose and somatic cell yields (0.79-0.89, P ≤ 0.01) (Table 2).

Table 1. Descriptive statistics of the persistency and extended lactation traits used in this study

Table 2. Phenotypic correlations among lactation persistency and extended lactation traits

Putative QTL for lactation persistency

Four significant (LOD > 2.0) and one highly significant (LOD > 3.0) QTL were detected on chromosomes 3, 10 and 11 using single QTL analysis across both methods. Six suggestive (1.7 ≤ LOD < 2.0), four significant (2.0 ≤ LOD < 3.0) and one highly significant (LOD ≥ 3.0) QTL were identified using QTL-MLE (Table 3, Figure 1, 2). Using QTL Express, eight suggestive (chromosome-wide P < 0.05), six significant (chromosome-wide P < 0.01) and two highly significant (experiment-wide P < 0.05) QTL were identified [see Additional file 2]. On OAR11, QTL for lactation persistency, extended lactation of milk, and extended lactation of protein were detected using both linkage methods. QTL for lactation persistency and extended lactation of milk and lactose were also identified on chromosome 17 and 21, respectively, but these were not verified by using QTL-MLE. We hypothesised that persistency of different lactation traits and extended lactation are likely to be controlled by different genes, with the possible exception of a common gene located on OAR11, which controls extended lactation for milk and protein yield, as well as milk persistency (Table 3).

Table 3. Results of the QTL analysis using QTL-MLE

thumbnailFigure 1. QTL map of the entire genome for extended lactation of milk and milk composition.

thumbnailFigure 2. QTL map of the entire genome for milk and milk composition persistency.

Additional file 2. Results of the QTL analysis using QTL Express. Information on the average QTL position (Peak) and the confidence interval (CI) in cM; F-value and respective significance threshold [*chromosome-wide significance level P < 0.05, **chromosome-wide significance level P < 0.01, ***experiment-wide significance level P < 0.05, ****experiment-wide significance level P < 0.01]; standardised QTL effect (SD) with standard error (SE); and phenotypic variance explained by the QTL are presented.

Format: PDF Size: 8KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Most of the QTL identified here showed positive (not necessarily favourable) effects for the Awassi allele, compared with the Merino allele. The largest effect was identified for extended lactation of lactose (0.69 SD) on OAR21. Persistency and extended lactation of milk showed a negative effect of the Awassi allele (-0.64 SD) on OAR11 (Table 3). The average QTL heritability, expressed as the proportion of the phenotypic variance accounted for by the QTL, and standardised QTL effects were slightly higher for extended lactation (5.48% and 0.58 SD) than for lactation persistency (5.11% and 0.56 SD) (Table 3) and [see Additional file 2]. The QTL heritability reached its highest value for fat traits and its lowest for useful yield (the QTL heritability is shown in Additional file 2). The standardised QTL effects were similar for all traits and reached values between 0.48 and 0.69 phenotypic SD. (Table 3). The greatest standardized QTL effects were found for somatic cells, and the smallest for lactose [see Additional file 2].

Using the two-QTL model for all traits across all chromosomes, evidence for two pairs of QTL was found on OAR10 (Table 4). For extended lactation of useful yield, neither of the two QTL were detected using single-QTL methods, while for extended lactation of somatic cells, one of the QTL positions was identified in both single-QTL approaches. Both pairs of QTL were in repulsion phase of opposite (unequal) effect and were located in separate flanking marker intervals, giving support for these being separate QTL but this needs to be validated in a separate study.

Table 4. Results of the two-QTL analysis using QTL Express


Two genome-wide and two partial bovine QTL studies [10,12,28,29] were summarised in the meta-assembly; results from two genome-wide association studies were not loaded into the ovine genome database [39,40]. Until now, no QTL study for lactation persistency in sheep has been reported. QTL for milk and milk composition persistency in cattle have been reported on chromosomes 1, 2, 6, 9, 14, 15, 17, 18, 21, 22, 25 and × using a granddaughter design of German Holstein dairy cattle [28]. QTL were also identified for shape of the lactation curve for milk, protein, fat and SCS in U.S. Holstein cattle, confirming the QTL on chromosomes 6, 14, 21 and 22 identified in German Holstein Friesians, and finding additional QTL regions on chromosomes 3 and 7 [29], while Weller et al. [10] could not verify a QTL for milk persistency on BTA7. Other QTL have been reported on BTA11 and 17 [12]. All QTL, except those on BTA15, 22 and 25, were confirmed by association studies, and additional associations have been identified on all other bovine chromosomes [39-41].

Four bovine QTL for lactation persistency traits have been reported on chromosomes that are orthologous to OAR3, 9, 17 and 24, where we identified QTL for these traits in sheep. The QTL on OAR9 and OAR17 aligned to QTL for the same milk persistence characteristics reported on the orthologous bovine chromosomes (BTA14 and BTA17, respectively) [12,29]. The QTL for milk persistency on BTA14 had its peak within the region of the DGAT1 gene. The QTL identified for fat yield extended lactation on OAR9 (BM4513-RJH1) in our study aligned to the bovine QTL for fat persistency at 139 cM (BM4513-BL1036).


This paper is a continuation of our previous study of lactation characteristics in a sheep population [18]. It extends the application of the Wood model using a nonlinear mixed model to describe lactation curves by deriving persistency and extended lactation measures for milk, protein, fat, lactose, useful yield and somatic cell yield. It identifies novel QTL for both persistency and extended lactation traits and compares linkage and genome-wide association studies for these traits across sheep and cattle.

Lactation persistency and extended lactation

Measures of persistency and extended lactation can be used to define the shape of lactation curves, the former focusing on shape within the conventional lactation length (say 100 days in sheep, 305 days in cattle) and the latter on sustained production after this period. Both these measures are useful indicators for practical management applications. Many other measures of lactation curve shapes have been described in the literature, some based on parameters in complex lactation curve models [11], others related to the time to reach a certain percentage of total milk production [2] or using smoothing cubic splines [39]. It has been suggested that persistency measures that are uncorrelated with total milk yield, for example using the shape of lactation curve independent of level of production, allows simultaneous selection for total lactation yield and persistency [4,42].

To our knowledge, there are no comparable studies of lactation persistency of different milk components in sheep, but in a study in cattle, the highest persistency was found for protein yield, followed by fat, and the lowest for milk yield [4]. We could confirm the lower values for persistency of milk yield in sheep, but the persistency of fat yield was higher than the persistency of protein yield. Considering that the major use of sheep milk is for producing cheese, protein and fat yields, and their combination in the form of useful yield, are important traits, and their persistency over the length of the lactation, enables a more sustained outflow of milk for cheese production [43].

A study in cattle [44] reported negative phenotypic correlations of total merit traits (namely lifetime net merit, cheese merit and fluid merit) with persistency of SCS. In the same study, it was also found that phenotypic correlations of milk, fat, and protein with SCS were negative, and that persistency of SCS had a close to zero correlation with SCS. They also found that persistency of SCS had low phenotypic correlations with milk, fat, and protein yield across breeds, but these findings were not observed for most of the traits in our study. It is difficult to resolve if this represents a difference between cattle and sheep, or if it is due to differences in the definition of lactation persistency [44].

While there is a reasonable amount of literature describing lactation persistency, particularly in cattle, little has been published in relation to extended lactation. The general assumption is that persistency and extended lactation represent similar traits. However, some studies have specifically aimed at describing extended lactation, e.g. Grossman and Koops [45], who investigated appropriate modelling techniques, and Haile-Mariam and Goddard [14], who reported on genetic parameters associated with extended lactation in cattle. While conventional lactation curve modelling studies may be extrapolated to predict extended lactation characteristics, empirical data on actual yields is needed to allow model fits to be assessed. The current study is based on observed milk yield beyond the standard lactation length in sheep (day 100), as are the two cattle studies reported above (beyond day 305).

Putative QTL

As shown in the first three papers of this series [17,29,37], the results from the maximum likelihood approaches to QTL detection (QTL-MLE) were in good agreement with those of the least-squares methodology of QTL Express [37]. However, it should be pointed out that QTL-MLE is based on fitting a finite mixture model, unlike regression-based methods, which mirrors the underlying Mendelian segregation process of the putative QTL [46].

To our knowledge, this is the first study describing a whole genome-scan for detecting QTL for lactation persistency and extended lactation of milk and milk composition in sheep. The results suggest the existence of important genomic regions for these traits but overall, the different traits are influenced by many different QTL (Tables 3, additional file 2). In a previous meta-assembly of published QTL reports, QTL for milk yield were identified on ovine chromosomes 1 to 3, 6, 9, 14, 16, 20, 22, and 24 [18]. However none of the regions identified for lactation persistency and extended lactation of milk on OAR11, 12, or 17 corresponded with these findings and only the QTL for milk persistency on OAR-17 was located in the comparative region to BTA17, where a QTL for milk persistency was identified in dairy cattle [12].

Among the two QTL identified for protein persistency and extended lactation on OAR11 and 24, only the QTL on OAR24 aligned to a QTL for milk persistency on the comparative bovine chromosomes 25 [28]. The QTL for lactation persistency and extended lactation of fat on OAR3 aligned with the QTL regions for average fat yield and energy persistency on BTA3 [28]. The QTL on OAR9 aligned with QTL for fat persistency on BTA14 at 139 cM. Rodriguez-Zas et al. [28] also reported on a QTL for the shape of the bovine lactation curve of fat at the proximal end of BTA14, which overlaps with the position of the DGAT1 gene [28,29]. The previous meta-assembly [18] also showed strong evidence of QTL for fat yield on both ovine chromosomes 3 and 9. Given the lack of QTL studies in both sheep and cattle for lactation persistency and extended lactation, the power of comparative mapping cannot be evaluated.

Among the three QTL regions for persistency and extended lactation of somatic cells identified on OAR2, 4 and 10, only the QTL on chromosome 2 mapped to a previously reported QTL in an ovine meta-assembly [18]. The QTL on OAR2 and 4 were supported by an association study for milk persistency [39]. Our QTL findings for useful yield and lactose yield are novel since there is no other published study.

From this study there does not appear to be a single major gene that regulates lactation persistency and extended lactation. Despite the reduced power to detect QTL for lactation traits since they are expressed in females only, the power of this experiment was sufficient to detect QTL with a minimum effect size of 0.4 SD. Such QTL may have utility in marker-assisted selection in the Awassi-Merino population used here. The Awassi breed has been used extensively to improve milk production through upgrading or crossbreeding programs and therefore QTL detected in the Awassi breed may also support marker-assisted introgression programs. Before such programs can be adopted, the QTL need validation in terms of effect sizes and location. The QTL findings also provide a foundation for targeted fine-mapping using SNP before positional candidate genes can be investigated. In particular, the QTL on OAR11 with effects on milk persistency and on extended lactation for milk and protein yield warrants further investigation.


This paper presents measures of lactation persistency and extended lactation in dairy ewes using the Wood lactation curve model. For the first time, QTL were identified for lactation persistency and extended lactation for milk components in a sheep resource population. It is shown that the QTL for lactation persistency and extended lactation for milk components (protein, fat, lactose, useful yield, SCS) differ from QTL for the corresponding yield traits. Furthermore, for many traits, the QTL affecting persistency are different from those affecting extended lactation. While a number of QTL of moderate effect size have been mapped, no individual QTL of major effect has been identified for persistency or extended lactation. Some QTL affecting similar traits in cattle and sheep were found to map to the same homologous regions. However, our ability to undertake this species comparison is limited by the relatively large size of support intervals in the sheep QTL study and by different trait definitions across the studies.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

EJ ran the final data and QTL analyses, and took the lead role in the preparation of the manuscript. PCT developed the statistical methodology for the lactation curve analysis and QTL methodology, implemented the QTL-MLE program, and contributed to the manuscript's preparation and the overall design. EJSH and DMcG carried out the analysis of the lactation curves, developed the models and helped drafting the manuscript. MKL ran the early stage QTL analyses, was responsible for the data assembly, and the phenotypic analysis. HWR had principal lead in the overall design, undertook the project management, and was involved in analysing the data and writing the manuscript. All authors read and approved the final manuscript.


The authors are particularly grateful to Dr Matthew Hobbs for his work on the QTL assembly and uploading the QTL onto the ovine genome browser. A special thanks also 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. The research was approved by the University of Sydney Animal Ethics Committee. The authors are grateful to the former CRC for Innovative Dairy Products, Australia for partial funding of this project. The resource flock and part of the sheep dairy was established with contributions from Mr Tom Grant and Mr Phillip Grant from Awassi Australia and a grant from the Australian Research Council.


  1. Jamrozik J, Fatehi J, Schaeffer LR: Application of robust procedures for estimation of breeding values in multiple-trait random regression test-day model.

    J Anim Breed Genet 2007, 124:3-11. PubMed Abstract | Publisher Full Text OpenURL

  2. Soelkner J, Fuchs W: A comparison of different measures of persistency with speacial respect to variation of test-day yields.

    Livest Prod Sci 1987, 16:305-319. Publisher Full Text OpenURL

  3. Appuhamy J, Cassell BG, Dechow CD, Cole JB: Phenotypic relationships of common health disorders in dairy cows to lactation persistency estimated from daily milk weights.

    J Dairy Sci 2007, 90:4424-4434. PubMed Abstract | Publisher Full Text OpenURL

  4. Cole JB, VanRaden PM: Genetic evaluation and best prediction of lactation persistency.

    J Dairy Sci 2006, 89:2722-2728. PubMed Abstract | Publisher Full Text OpenURL

  5. Brody S, Ragsdale AC, Turner CW: The rate of decline of milk secretion with the advance of the period of lactation.

    J Gen Physiol 1923, 5:441-444. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Wood PDP: Algebraic model of the lactation curve in cattle.

    Nature 1967, 216:164-165. Publisher Full Text OpenURL

  7. Kirkpatrick M, Heckman N: A quantitative genetic model for growth, shape, reaction norms, and other infinite-dimensional characters.

    J Math Biol 1989, 27:429-450. PubMed Abstract | Publisher Full Text OpenURL

  8. Jamrozik J, Schaeffer LR, Dekkers JCM: Genetic evaluation of dairy cattle using test day yields and random regression model.

    J Dairy Sci 1997, 80:1217-1226. PubMed Abstract | Publisher Full Text OpenURL

  9. Pollott G, Gootwine E: A genetic analysis of complete lactation milk production in improved Awassi sheep.

    Livest Prod Sci 2001, 71:37-47. Publisher Full Text OpenURL

  10. Weller JI, Ezra E, Leitner G: Genetic analysis of persistency in the Israeli Holstein population by the multitrait animal model.

    J Dairy Sci 2006, 89:2738-2746. PubMed Abstract | Publisher Full Text OpenURL

  11. Grossman M, Hartz SM, Koops WJ: Persistency of lactation yield: A novel approach.

    J Dairy Sci 1999, 82:2192-2197. PubMed Abstract | Publisher Full Text OpenURL

  12. Boichard D, Grohs C, Bourgeois F, Cerqueira F, Faugeras R, Neau A, Rupp R, Amigues Y, Boscher MY, Leveziel H: Detection of genes influencing economic traits in three French dairy cattle breeds.

    Genet Sel Evol 2003, 35:77-101. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  13. Weller JI, Golik M, Reikhav S, Domochovsky R, Seroussi E, Ron M: Detection and analysis of quantitative trait loci affecting production and secondary traits on chromosome 7 in Israeli Holsteins.

    J Dairy Sci 2008, 91:802-813. PubMed Abstract | Publisher Full Text OpenURL

  14. Haile-Mariam M, Goddard ME: Genetic and phenotypic parameters of lactations longer than 305 days (extended lactations).

    Animal 2008, 2:325-335. OpenURL

  15. Vargas B, Koops WJ, Herrero M, van Arendonk JAM: Modeling extended lactations of dairy cows.

    J Dairy Sci 2000, 83:1371-1380. PubMed Abstract | Publisher Full Text OpenURL

  16. Dematawewa CMB, Pearson RE, VanRaden PM: Modeling extended lactations of Holsteins.

    J Dairy Sci 2007, 90:3924-3936. PubMed Abstract | Publisher Full Text OpenURL

  17. Rook AJ, France J, Dhanoa MS: On the mathematical description of lactation curves.

    J Agric Sci 1993, 121:97-102. Publisher Full Text OpenURL

  18. 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 OpenURL

  19. Hammami H, Rekik B, Soyeurt H, Bastin C, Bay E, Stoll J, Gengler N: Accessing genotype by environment interaction using within-and across-country test-day random regression sire models.

    J Anim Breed Genet 2009, 126:366-377. PubMed Abstract | Publisher Full Text OpenURL

  20. Strabel T, Jamrozik J: Genetic analysis of milk production traits of polish black and white cattle using large-scale random regression test-day models.

    J Dairy Sci 2006, 89:3152-3163. PubMed Abstract | Publisher Full Text OpenURL

  21. Gonzalo C, Carriedo JA, Baro JA, Primitivo FS: Factors influencing variation of test day milk yield, somatic cell count, fat, and protein in dairy sheep.

    J Dairy Sci 1994, 77:1537-1542. PubMed Abstract | Publisher Full Text OpenURL

  22. Othmane MH, Carriedo JA, de la Fuente Crespo LF, San Primitivo F: An individual laboratory cheese-making method for selection in dairy ewes.

    Small Rumin Res 2002, 45:67-73. Publisher Full Text OpenURL

  23. Haenlein GFW, Park YW, Raynal-Ljutovac K, Pirisi A: Special issue: goat and sheep milk.

    Small Rumin Res 2007, 68:232. OpenURL

  24. Barillet F: Genetic improvement for dairy production in sheep and goats.

    Small Rumin Res 2007, 70:60-75. Publisher Full Text OpenURL

  25. Fuertes JA, Gonzalo C, Carriedo JA, San Primitivo F: Parameters of test day milk yield and milk components for dairy ewes.

    J Dairy Sci 1998, 81:1300-1307. PubMed Abstract | Publisher Full Text OpenURL

  26. Sorensen A, Muir DD, Knight CH: Extended lactation in dairy cows: effects of miking frequency, calving season and nutrition on lactation persistency and milk quality.

    J Dairy Res 2008, 75:90-97. PubMed Abstract | Publisher Full Text OpenURL

  27. Wendorff B: Milk composition and cheese yield.

    Proceedings of the 7th Great Lakes Dairy Sheep Symposium: 1-3 November 2001; Eau Claire: 2002104-117. OpenURL

  28. Harder B, Bennewitz J, Reinsch N, Thaller G, Thomsen H, Kuhn C, Schwerin M, Erhardt G, Forster M, Reinhardt F, Kalm E: Mapping of quantitative trait loci for lactation persistency traits in German Holstein dairy cattle.

    J Anim Breed Genet 2006, 123:89-96. PubMed Abstract | Publisher Full Text OpenURL

  29. Rodriguez-Zas SL, Southey BR, Heyen DW, Lewin HA: Detection of quantitative trait loci influencing dairy traits using a model for longitudinal data.

    J Dairy Sci 2002, 85:2681-2691. PubMed Abstract | Publisher Full Text OpenURL

  30. Raadsma HW, Thomson PC, Zenger KR, Lam MK, Jonas E, Cavanagh C, 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 OpenURL

  31. Barillet F: Genetics of milk production. In The Genetics of Sheep. Edited by Piper L, Ruvinsky A. Department of Animal Science, University of New England, Armidale, Australia: CABI; 1997:539-564. OpenURL

  32. Wood PDP: Factors Affecting Persistency of Lactation in Cattle.

    Nature 1968, 218:894-895. Publisher Full Text OpenURL

  33. Berger Y, Mikolayunas C, Thomas D: Dairy Sheep Fact Sheet. [;] webcite

    Dairy Business Innovation Center, and University of Wisconsin-Extension; 2009.

  34. The Awassi sheep with special reference to the improved dairy type-Milk and milking [] webcite

  35. Van Der Werf J[] webcite

  36. Thomson PC, Brown SC, Raadsma HW: QTL-MLE: A maximum likelihood QTL mapping program for flexible modelling using the R computing environment.

    Proceedings of the Association for the Advancement of Animal Breeding and Genetics: 28 September-1 October 2007: 2007; Armidale288-290. OpenURL

  37. 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:339-340. PubMed Abstract | Publisher Full Text OpenURL

  38. Cavanagh CR, Jonas E, Hobbs M, Thomson PC, Tammen I, Raadsma HW: 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.

    Genet Sel Evol 2010, 42:36. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  39. Verbyla KL, Verbyla AP: Estimated breeding values and association mapping for persistency and total milk yield using natural cubic smoothing splines.

    Genet Sel Evol 2009, 41:48. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  40. Pryce JE, Haile-Mariam M, Verbyla K, Bowman PJ, Goddard ME, Hayes BJ: Genetic markers for lactation persistency in primiparous Australian dairy cows.

    J Dairy Sci 2010, 93:2202-2214. PubMed Abstract | Publisher Full Text OpenURL

  41. Pant SD, Schenkel FS, Leyva-Baca I, Sharma BS, Karrow NA: Identification of single nucleotide polymorphisms in bovine CARD15 and their associations with health and production traits in Canadian Holsteins.

    BMC Genomics 2007, 8:1471-2164. OpenURL

  42. Appuhamy JADRN, Cassell BG, Cole JB: Phenotypic and genetic relationships of common health disorders with milk and fat yield persistencies from producer-recorded health data and test-day yields.

    J Dairy Sci 2009, 92:1785-1795. PubMed Abstract | Publisher Full Text OpenURL

  43. Barillet F, Marie C, Jacquin M, Lagriffoul G, Astruc JM: The French Lacaune dairy sheep breed: use in France and abroad in the last 40 years.

    Livest Prod Sci 2001, 71:17-29. Publisher Full Text OpenURL

  44. Cole JB, Null DJ: Genetic evaluation of lactation persistency for five breeds of dairy cattle.

    J Dairy Sci 2009, 92:2248-2258. PubMed Abstract | Publisher Full Text OpenURL

  45. Grossman M, Koops WJ: Modeling extended lactation curves of dairy cattle: A biological basis for the multiphasic approach.

    J Dairy Sci 2003, 86:988-998. PubMed Abstract | Publisher Full Text OpenURL

  46. Jansen RC: A general mixture model for mapping quantitative trait loci by using molecular markers.

    Theor Appl Genet 1992, 85:252-260. OpenURL