Genetic-economic comparison of New Zealand Holstein and Jersey-New Zealand Holstein cross under a pastoral scheme in southern Chile
Evaluación genético-económica de Holstein Neozelandés y cruzamiento Jersey-Holstein Neozelandés en un sistema pastoril de la X Región de Chile
a Texas A&M University, College Station, TX, USA.
ABSTRACT. Chilean bovine dairy production needs to increase milk solid content. This study evaluates and compares genetic and economic progress of New Zealand Holstein (HoNZ) and New Zealand Jersey-HoNz crossbreed (Fx) using selection indices. The traits analysed were: milk (L) with 3% fat and 3% protein (LB), fat (kg) over 3% in LB (PF), protein (kg) over 3% in LB (PP) and body weight (BW). A multivariate mixed model was used for genetic parameters estimation using AIREMLF90 software. Economic values were estimated using costs associated with dry matter (DM) intake required to supply energy requirements of each trait and income as defined by PROLESUR payment to producer scheme for LB, PF and PP. Genetic progress of PF, PP and BW were similar for the two studied genotypes, but LB progress was higher in HoNz. Economic progress was similar between genotypes, but the offspring of Fx achieves a greater economic benefit per animal and per surface of land, due to higher production of PF and PP. Genetic improvement of Fx using a selection index, allows to achieve higher average of solids production and higher economic benefit.
Keywords: genetic improvement, milk, solids.
RESUMEN. La producción de leche bovina requiere incrementar el contenido de sólidos. Este estudio evaluó y comparó el progreso genético y económico del Holstein Neozelandés (HoNz) y del cruzamiento Jersey Neozelandés-HoNz (Fx) utilizando índices de selección. Las características analizadas fueron: leche (L) con 3% de grasa y 3% de proteína (LB), grasa (kg) sobre 3% en LB (PF), proteína (kg) sobre 3% en LB (PP) y peso vivo (BW). Se utilizó un modelo mixto multivariado para realizar estimaciones de parámetros genéticos usando AIREMLF90. Se estimaron valores económicos usando los costos asociados al consumo de materia seca necesario para suplir los requerimientos de energía de cada característica y los ingresos de la pauta de pago de PROLESUR para LB, PF y PP. El progreso genético de PF, PP y BW fue similar entre los biotipos, pero el de LB fue mayor en HoNz. El progreso económico fue similar entre los biotipos, pero la descendencia de Fx permite conseguir un mayor retorno económico por animal y por superficie gracias a su mayor producción de PF y PP. El mejoramiento genético de Fx utilizando índices de selección permite conseguir mayores promedios de producción de sólidos y un mayor beneficio económico para el plantel.
Palabras clave: mejoramiento genético, leche, sólidos.
Seventy percent of industrialised milk in Chile is produced in the southern part of the country, and the most important production areas are Araucanía (IX), Los Lagos (X) and Los Ríos (XIV) (Lerdón et al 2010) where milk production is mainly based on grazing pastures where the goal is to obtain high production per unit of land (González-Verdugo et al 2004).
The main factors influencing economic return are: 1) liters (L) of milk containing 3% weight/volume (w/v) of milk fat plus 3% w/v of protein, which is known as standardized liter of milk (LB); 2) kilograms (kg) of fat above 3% w/v in each L of LB (PF) and, 3) kg of protein above 3% w/v in each L of LB (PP). The payment scheme of the industry for raw milk is influenced by these three factors (PROLESUR)1. On the other hand, body weight (BW) also represents an important factor to be accounted for in this type of pastoral milk production, because as BW per cow decreases, more cows can be carried per unit of land. Lighter cows have lower nutritional requirements compared with heavier cows (Hazard 2004).
As part of its "Strategies for Competitive Development of the Chilean Dairy Sector 2010-2020" document, the Milk Consortium S. A.2 outlined a production blueprint to be followed by the dairy sector. In relation to the improvements of economically important traits, it was pointed out that there is the need to move toward dairy breeds whose milk has higher milk solids (fat and protein) concentrations. The goal, as indicated by the Consortium, is to a reach national milk solids concentration of 7.6% in 2020 compared with the current 7.1% in 2010.
Improvement of the economically important traits can be reached by two ways; firstly, through farm management practices, such as improvement in animal nutrition and health. Change in management practices leads to immediate improvements however, are temporal and may sometimes lead to increased production costs. The second way to increase milk solids concentrations is through genetic improvement, which in the medium and long term, leads to cumulative and permanent changes (Conington et al 2001).
Selection indexes maximize economic profit by selecting breeding stock not only based on genetic merit; but also include the candidate's overall economic merit (Hazel 1943). A selection index represents an estimation of the selection objective or aggregate genotype, which is the true genetic-economic value of an animal obtained as the sum of the economic contribution of each trait included in the index. Therefore, knowledge on the economic value of each trait is needed to develop an estimation of the index value, with the economic value being defined as the additional net return obtained by increasing the output by one unit of the trait keeping the other traits unchanged (Cameron 1997).
In Chile, there are studies concerning the selection indexes for productive domestic animals, among them, dual purpose and meat sheep (Lembeye et al 2014), Corriedale sheep (Barría et al 1998) and Hereford cattle (Amthauer 1995); however, in spite of the advantages of genetic improvement using selection indexes and the economic importance of pastoral dairy production in southern Chile, there is no Chilean information regarding the use of selection indexes in pastoral dairy cattle.
Among dairy breeds and crossbreds suitable for pastoral systems, the New Zealand Holstein (HoNz) and the New Zealand Jersey-HoNz (Fx) cross are found in Chile; these genotypes are under the dairy genetic improvement program in New Zealand, which is focused on increasing milk solids content (López-Villalobos et al 2000, Livestock Improvement 2004, López-Villalobos and Garrick 2005, Sneddon 2011).
Due to some similarities between production objectives in dairy pastoral systems of New Zealand and Chile, HoNz genotype is one of the breeds that forms the dairy herd of southern Chile; however in the scientific literature reviewed here, there are no Chilean studies aimed at determining which breed or crossbred provides the greatest economic progress based on increasing milk solids output per unit of land.
The development of selection indexes, specifically for HoNz and Fx which considers milk, milk solids and BW, would provide information about which is the most suitable genotype for economic returns under southern Chilean pastoral conditions.
The objective of this study was to determine, under the current Chilean milk production economic scenario, which genotype would yield greater genetic and economic progress by jointly selecting for lower BW and higher LB, PF and PP production, using selection indexes.
MATERIAL AND METHODS
Between 2001 and 2011 data were gathered at Oromo Research Station, a farm owned by the Faculty of Agricultural Sciences of the Universidad de Chile and located in Purranque, Lake District, Chile (41°08 S, 73°09 W).
The records from 317 and 128 HoNz and Fx cows, respectively, were considered for the study, HoNz cows had a genetic makeup of at least 75% HoNz; Fx genotype were crossbred cows obtained using New Zealand Jersey semen on HoNz cows, Jersey genes fluctuated between 25 to 75%. The final data set contained 823 and 440 lactations for HoNz and Fx, respectively.
Besides genotype, the variables used in this study were cow identification, sire and dam identification, month and year of calving parturition number (NOP), days in milk, BW and lactation yield (kg) of milk, fat and protein.
The pasture at Oromo Research station has a clearly marked productive seasonality; the growth reaches its maximum in the first and second week of October, peaking up to 70 kg dry matter (DM)/day. In Summer, both productivity and nutritional quality decrease. Minimal growth occurs in winter time. The average botanical composition of the pasture is: Lolium perenne (70%); Dactilis glomerata (15%), Anthoxanthum odoratum, Bromus sp., Trifolium repens y Achillea mellefolium (5%) and other species (10%). Grazing is rotational, keeping a stocking rate of 2.45 cows/ha/year. During the non-lactating period (May and June) the cows are fed silage ad libitum. Two to three weeks before calving the cows are moved to calving paddocks where, along with silage, they are supplemented with concentrate without exceeding to 2 kg/cow/day. After calving, cows are fed with silage and the amount of concentrate is doubled to 4 kg/day. The amounts of silage offered begins to decrease as spring pasture growth replaces it; concentrate supplementation also decreases until pasture completely replaces silage and concentrate in the third week of September. Consequently the calving date of the cow determines the amount of concentrate given at the beginning of lactation. Mineral supplementation is provided around the year according to physiological condition and production level. To ensure that mineral supplementation is consumed by the cows it is mixed with the silage and, when silage is not offered, it is offered at the milking parlour. The goal of reproductive management is to concentrate calving toward the end of winter to match the largest nutritional requirement of the herd with the highest pasture growth.
Lactations began between July and November in each year. Parturition number ranged between 1 and 11 however, due to both genotypes having fewer records beyond parturition 6; records from NOP 7 to 11 were grouped in a single category. A single BW record per lactation was recorded on each cow a day after calving. Lactations had a minimum of 200 days in milk and those exceeding 305 days in milk were shortened to the nearest 305 days record.
Traits with economic relevance were: 1) L of standardised milk containing 3% w/v of fat and 3% w/v of protein (LB), this was calculated as follows: LB = LA/DC where, LA is total milk yield per lactation (kg) and DC is the milk density coefficient (1.03; FAO, 2012); 2) kg of fat above 3% w/v in each L of LB (PF), this was defined as: PF = FAT - (0.03 * LB), where FAT is total kg of fat yield per lactation; and 3) kg of protein above 3% w/v, this was defined as: PP = PROT - (0.03 * LB), where PROT is total kg of protein yield per lactation.
Least squares means (LSM) of BW were estimated to assess economic values of the selection index. Furthermore, LSM for BW, LB, PF and PP were used to assess population changes after one generation of selection using the indexes derived for both genotypes. This was done to estimate the expected cow performance, which a dairy producer would achieve, in the medium term, using this methodology, under environmental and management conditions similar to those of Oromo research farm.
Additive genetic values were estimated using a multivariate linear mixed model, in order to be used in the selection indexes. The BW model considered the fixed effect of genotype, month and year of calving, NOP, and the random effects of animal ~N (0,Aσa2) and permanent environment. For LB, PF, and PP the model included the same aforementioned variables plus fixed effect of days in milk as a covariate. The model also included relationship information among animals, and was solved using BLUP methodology by the AIREMLF90 software (Misztal et al 2002).
Aggregate genotype or total genetic merit function for each animal (Ĥ) was estimated as follows: Ĥ = v1 Ȃ1 + v2 Ȃ2 + v3 Ȃ3 + v4 Ȃ4, where vi is the economic value of the ith trait in the index, i = LB, PF, PP and BW. Ȃ1 is the estimated breeding value of a given animal for the ith trait.
Economic value estimation. Economic value was defined as the difference between income and cost associated to increase a single unit of each trait included in the selection index.
Additional income from the production of a L of LB, a kg of PF and a kg of PP were obtained from the current raw milk payment scheme of the PROLESUR company. Therefore, income was 119, 1,200 and 5,800 Chilean pesos for each L of LB, kg of PF and kg of PP, respectively. Additional income for a kg of BW was set to zero due to the fact that incremental changes in BW do not imply direct economic benefits.
To standardise costs, consumption of DM necessary to supply the metabolizable energy (ME) requirement associated with an incremental change in a unit of each trait, was used (Daza 1997, Chamberlain and Wilkinson 2002).
As diet composition varies within a given lactation, an average lactation feed composition was estimated based on the percentage of DM consumption of each diet ingredient. The ingredients were pasture, grass silage and concentrate. Total consumption per lactation was the sum of DM consumption from the three ingredients aforementioned.
Lactation pasture consumption. This was calculated using the following equation: CP = CPT - CC - CE - (TS * CC), where CP is lactation pasture consumption (kg of DM), CPT is the lactation potential consumption (kg of DM), CC is lactation concentrate consumption (kg of DM), CE is lactation silage consumption (kg of DM), TS is the pasture by concentrate substitution rate which, for this study, was defined as 0.4. The chosen TS value is within the range from 0.26 to 0.55 described for Holstein cows under low and high pasture allowance, respectively (Bargo et al 2002). To estimate CPT it was considered that DM consumption differences existed within lactation, therefore, a weighted average consumption was calculated according to stage of lactation. Consumption was calculated as a percentage of BW based on the assumption that within lactation, there are DM consumption differences as follows (Hazard 2004), a) DM consumption in the initial third of lactation = 3.5% of BW, b) DM consumption in the middle part of lactation = 3.0% of BW, c) DM consumption in the last third of lactation = 2.5% of BW.
Average body weight in each third of lactation. To account for the fact that BW varies within lactation, a simulation was used to model BW changes in each lactation third, the assumptions used were (Hazard 2004), a) from the first to second month of lactation BW decreases by 1 kg per day, b) from the third to fifth month of lactation there is no change in BW, c) from the sixth to the tenth month of lactation BW gradually increases until the beginning of dry period, reaching the same BW as it was previous to calving. LSM for BW in both genotypes were used as starting values for the simulation; then average BW of each lactation third, for both genotypes, were calculated to estimate potential consumption.
Lactation concentrate consumption. Concentrate supplementation scheme at Oromo Research Station, by default, consists of offering 4 kg of concentrate per day, from July 15 to September 20. Therefore, a consumption weighted average was calculated per lactation using records of cows that calved within the above mentioned period. Dry matter consumption was estimated from concentrate consumption, assuming concentrate had a DM of 88% (Anrique et al 2014).
Lactation silage consumption. Silage supplementation at Oromo Research Station is routinely practiced with cows calving between July 15 and September 7. Furthermore, silage supplementation changes overtime; each cow is offered 8, 5 and 3 kg in July, August and September, respectively. From this, average kg of silage consumption, per cow, was calculated according to calving date and, DM consumption was estimated assuming that silage has 42% DM (Anrique et al 2014).
Feed metabolizable energy. Estimation of metabolizable energy (ME) content per kg of feed DM was calculated based on the energetic information and consumption percentage of each feed ingredient, the following equation was used: EMA = (EMC * %C) + (EME * %E) + (EMP * %P), where EMA is mega calories (Mcal) of ME per kg of DM of the lactation average diet, EMC is Mcal of ME per kg of DM from concentrate feed, EME is Mcal of ME per kg of DM from silage, EMP is Mcal of ME per kg of DM from pasture, %C, %E and %P are consumption percentages of concentrate, grass silage and pasture, respectively.
Concentrate and grass silage energy. Information regarding energetic content per kg of DM in concentrate and silage was obtained from dairy cow supplement nutrient tables; which were specifically designed for grazing dairy cows (Parga and Lanuza 2006). Energy content used was 2.8 Mcal of ME per kg of DM for concentrate (oat grain) and 2.5 Mcal of ME per kg of grass silage.
Pasture energy. ME per kg of DM in pasture varies over time, therefore an ME weighted average was used which was based on nutritional information provided by the Institute for Agricultural Research (INIA) for permanent pastures, such as those found at Oromo Research Station (Parga and Lanuza 2006). According to INIA, winter pasture contains 2.7 Mcal of ME per kg of DM, spring pasture contains 2.9 Mcal of ME per kg of DM, summer pasture contains 2.4 Mcal of ME per kg of DM and autumn pasture contains 2.7 Mcal of ME per kg of DM.
Ingredient consumption percentage of each genotype. DM consumption of each genotype was calculated as follows:
where Cons.It is DM consumption of the ith. feed ingredient (i = pasture, silage, concentrate) of each genotype, and Cons.T is total DM consumption of each genotype.
Additional DM consumption and requirement, per output unit, of each trait according to genotype. Estimation of DM consumption, for FX and HoNz, was calculated based on the relationship between ME requirements, associated with an additional unit of LB, PF, PP and BW, and ME content of DM in the average feed diet.
ME requirements to produce and additional unit of LB, PF and PP. Requirements, in Mcal of ME per L of LB were calculated using the following formula:
where d is milk density (1.03) and kl is ME utilization efficiency expressed as net energy for lactation (EN1). In this study the kl value used was 0.64 (González-Verdugo et al 2004).
For one kg of PF the formula was:
The additional requirement to produce an additional kg of PP was:
ME requirement to maintain an additional kg of BW
DM Consumption to produce an additional unit of LB, PF and PP.
DM necessary to produce an additional L of LB was:
for an additional kg of PF was
and for each additional kg of PP was:
DM Consumption to maintain an additional kg of BW was:
Additional cost of a kg of DM in each diet ingredient
To estimate additional DM intake costs for each genotype, the DM requirement associated with an incremental change in the output unit of each trait was multiplied by the cost of a kg of DM of each feed ingredient.
Cost of a kg of DM of the diet. The cost of a kg of DM of the diet ($DMkgFeed) was calculated as follows:
where: $DMkgConc is the price of a kg of DM from concentrate.
$DMkgSil is the price of a kg of DM from grass silage.
where CM are the costs associated with pasture maintenance plus 6% unforeseen costs (Teuber 2009), PRP is the kg of DM/ha grown each year, which, was assumed to be 10,000 kg DM/ha per year (Teuber 2009). %Ofer is pasture percentage available for consumption leaving 25% for pasture regrow, and %Cons is the offered percentage of the pasture which is effectively consumed by the cow. According to historical production records of Oromo Research Farm %Cons was assumed to be 75%.
Costs associated to the additional output of 1L of LB ($LB), 1 kg of PF ($F1kg), 1 kg of PP ($P1kg) and maintenance of one additional kg of BW ($BW1kg).
The cost for LB was: $ LB = DMKg / LBL * $ DMKg Feed .
The cost of PF was: $ F1kg = DMKg / F1kg * $ DMKg Feed .
The cost of PP was: $ P1kg = DMKg / P1kg * $ DMKg Feed .
The costofBWwas: $ BW1kg = DMKg / 1KgBW * $ DMKg Feed .
Genetic-economic progress selecting using a selection index. A t test, modified by Welch (Welch 1938), was used to compare selection index average values between both genotypes.
Least square means and their standard errors for both genotypes are presented in table 1. Milk solid yield for the Fx was greater than the HoNz however, milk yield was greater for the HoNz.
1/_ milk containing 3% weight/volume (w/v) of milk fat plus 3% w/v of protein; 2/_ kg of fat
DRY MATTER CONSUMPTION
After simulation of BW across lactation for both genotypes, potential DM consumption was estimated. Consequently Fx had a potential consumption of 3,819 kg of DM per lactation (12.5 kg DM daily average). Potential consumption for HoNz was estimated in 4,150 kg of DM per lactation (13.6 kg DM daily average).
Feed ingredients ME was estimated at 2.80, 2.50 and 2.67 Mcal of ME per kg of DM for concentrate, silage and pasture, respectively; using this, consumption percentage and energy intake was estimated as shown in table 2.
Table 2. Metabolizable energy (ME) intake per kg of average feed dry mater (DM)
1/_ Jersey-New Zealand Holstein cross; 2/_ New Zealand Holstein
REQUIREMENT AND CONSUMPTION OF METABOLIZABLE ENERGY PER UNIT OF EACH TRAIT
The estimated metabolisable energy required to increase one unit of LB, PF, PP and BW was 1.04, 14.45, 8.03 and 0.02 Mcal, respectively. The result is consistent with the expectations, because the high energy concentrated in one kg of PF and PP demands a higher ME than the ME needed to produce one extra liter of LB, and to satisfy the maintenance requirements of an extra kg of BW. As a consequence, additional DM consumption associated to a unit of LB, PF, PP and BW was 0.39, 5.41, 3.01 and 0.01 kg, respectively.
According to current Chilean feedstock market prices the estimated cost of a kg of DM from each feed ingredient was 227, 48 and 43 Chilean pesos ($) for concentrate, silage and pasture, respectively.
Using the previous information, the resulting economic values needed to build the selection index, are shown in table 3 for both genotypes.
1/_ milk containing 3% weight/volume (w/v) of milk fat plus 3% w/v of protein; 2/_ kg of fat
With the economic values presented in table 3 the aggregate genotype estimations for Fx and HoNz were:
where Ȃij is the estimated additive genetic value for the ith trait of the jth animal. Aggregate genotypes were estimated for all animals within breed and averages of Fx and HoNz were compared. Average aggregate genotype for Fx was 4,047.53 Chilean pesos which is significantly greater than the HoNz, whose aggregate genotype was -2,248.81.
To estimate the economic values for the traits included in the selection indexes, it was necessary to develop a system to standardise the cost associated with each trait. Feeding is assumed to be the most important production cost (Hazard 2004, Parga and Lanuza 2006), therefore, additional DM consumption costs to produce a unit of a given trait were used. Also, changes in BW during lactation for both genotypes were simulated by using BW LSM estimates, allowing the inclusion of the variables affecting the amount and type of consumed feed to obtain results that are a close approximation to the complexity of the system (Hazard 2004).
The results show that LSM of LB, PF, PP and BW were significantly different between genotypes (table 1). For productive traits (LB, PF and PP) HoNz showed higher average LB, exceeding the Fx average by approximately 119.7 L LB. However, despite the lower Fx average for LB, the FX had a significantly greater average for PF and PP, these being traits of economic importance for dairy producers according to the current payment scheme of the raw milk buying company (PROLESUR). These findings are similar to those indicated by other authors when comparing the same breeds (González-Verdugo et al 2004, Prendiville et al 2009).
Fx had a lower average BW than HoNz (table 1), this is in agreement with the findings of González-Verdugo et al (2004) and Prendiville et al (2009). This lower BW means that maintenance ME requirement and, in turn, DM consumption of Fx is lower, therefore the number of cows that can be kept per unit area of pasture is higher compared with HoNz. These results indicate that Fx would be better adapted to pastoral productive system, maximizing economic return per unit of grazing land. Since dairy production in southern Chile is mainly pastoral-based, economic returns of the dairy enterprise could be optimised by using Fx or similar genotype.
As part of the selection index design process, a difference in pasture consumption between genotypes was found. This can be explained due to differences in potential consumption of the two genotypes and due to the supplementation strategy of the herd, where all cows were offered, on average, the same amount of concentrate and grass silage per lactation. However, the lower potential DM consumption of Fx indicates that the feed fraction corresponding to pasture be lower than that of the HoNz. In spite of these differences, when matching the energy content of the feed ingredients to consumption percentage, it was calculated that consumed ME per kg of feed is identical in both genotypes (table 2).
Metabolizabe energy required by cows for an incremental change in a unit of each trait was as expected, that is, productive traits had a greater energetic demand than an additional kg of BW; this can be explained as a consequence of the metabolic process associated and energy concentration per unit of LB, PF and PP (Nelson and Cox 2005). As ME requirements were greater for incremental changes in productive traits, DM consumption was also increased, mainly for changes in PF.
After the cost of an average kg of DM of feed intake was estimated, an estimation of the feeding costs for both genotypes was calculated. The Fx had a greater feeding cost per kg of DM in comparison with HoNz, due to the inclusion of a lower percentage of pasture in the average kg of DM of feed consumed during lactation. In spite of the fact that both genotypes were fed using the same concentrate and silage supplementation protocol, the lower percentage of pasture in the Fx's average lactation diet, caused the relative percentage of the more expensive ingredients (silage and concentrate) to be higher in the diet. As a consequence, the cost of the estimated average kg of feed eaten by Fx across the lactation was greater than HoNz.
By knowing feed cost and additional DM consumption it was possible to estimate and standardize additional cost that the herd would have due to increasing a unit of each trait. With this result, and the additional income associated to increasing a unit of each trait, it was possible to estimate economic values of LB, PF, PP and BW for both genotypes (table 3). The economic values of LB, PF, and PP for the FX were lower than those for HoNz, as a result of higher feeding costs involved to produce a unit of these traits. For BW both genotypes had the same economic value, and the selection index weighted negatively this trait when the cow estimated breeding value was above zero.
Using the aforementioned results plus economic values, aggregated genotypes for HoNz and Fx were estimated and compared; Fx had a significantly greater genetic-economic superiority of $6,296 compared with HoNz.
In these findings one should consider that part of Fx better performance must be due to heterosis effects (López-Villalobos et al 2000); therefore, within the selection strategy of this genotype a rotational crossing scheme, aimed to keep Fx heterozygosity, must be considered.
Finally, of the two genotypes used in this study, Fx was the best genetic-economic alternative for this farm. A significant portion of this better economic result can be explained due to the fact that Fx had a greater milk solids yield than HoNz; and this is a direct consequence of Jersey breed genes involved in the cross. These findings indicate that Chilean dairy farmers would, in term of profitability, benefit by preferring Fx cross over HoNz breed. If dairy farmers want to increase milk solid output and hence profitability, the recommendation would be the use a rotational crossbreeding scheme. Construction of a selection index in a breeding program, which uses local economic parameters, is crucial in a breeding program to move the herd toward a production profit goal.
1 Prolesur. 2014. Resumen de Pauta de Pago de leche para la compra que PROLESUR realiza a los actuales productores de leche de la X Región de Los Lagos, Continental. Vigencia: Desde el 1 de septiembre de 2014. www.prolesur.cl/component/docman/doc_download/249-pauta-de-precios-leche-cruda-x-continental-sept-14; accessed July 28, 2013.
2 Consorcio Lechero, 2010. Estrategia de desarrollo competitivo del sector lácteo chileno, 2010-2020. www.consorciolechero.cl/chile/docs/Estrategia-Desarrollo-Sectorial-2010-2020.pdf; accessed April 10, 2013.
Amthauer JL. 1995. Índices de selección para reproductores Hereford bajo las condiciones de la zona sur de Chile. Biblioteca Central INIA, Santiago, Chile.
Anrique R, X Molina, M Alfaro, R Saldaña. 2014. Composición de alimentos para el ganado bovino. 4a ed. Consorcio Lechero, Universidad Austral de Chile e Instituto de Investigaciones Agropecuarias Remehue, Valdivia, Chile.
Bargo F, LD Muller, JE Delahoy, TW Cassidy. 2002. Milk response to concentrate supplementation of high producing dairy cows grazing at two pasture allowances. J Dairy Sci 85, 1777-1792.
Barría N, H Montaldo, A Jara, V Martínez. 1998. Selection Index construction for the Corriedale breed in Magallanes. Proceedings 6th World Congress on Genetics Applied to Livestock Production 24, 39-42.
Cameron ND. 1997. Selection index methodology. In: Cameron ND (ed). Selection indices and prediction of genetic merit in animal breeding. CAB International, Oxford, UK, Pp 58-67.
Chamberlain AT, JM Wilkinson. 2002. Alimentación de la vaca lechera. Acribia, Zaragoza, España.
Conington J, SC Bishop, B Groudy, A Waterhouse, G Simm. 2001. Multi-trait selection indexes for sustainable UK hill sheep production. JAni Sci 73, 413-423.
Daza A. 1997. Reproducción y Sistemas de Explotación del ganado ovino. Mundi Prensa, Madrid, España.
González-Verdugo H, JC Magofke, C Mella. 2004. Productivity, intake and biological efficiency in New Zealand Friesian and F1 (Jersey-New Zealand Friesian) cows calved during late winter in the Xth Region, Chile. Arch Med Vet 37, 37-47.
Hazard S. 2004. Alimentación de vacas lecheras. Boletín INIA N° 112, Pp 52-60.
Hazel LN. 1943. The genetic basis for constructing selection indexes. Genetics 28, 476-490.
Lembeye F, G Castellaro, JC Magofke, H Uribe. 2014. Comparación de distintos índices y criterios de selección en ovinos Merino Precoz manejados en condiciones extensivas de la zona central de Chile. Arch Med Vet 46, 389-397.
Lerdón J, C Muñoz, V Moreira. 2010. Productive and economic characterization of medium and large dairy farmers of southern Chile. IDESIA 28, 41-49.
Livestock Improvement. 2004. Dairy Statistics 2003-2004. Livestock Improvement Corp. Ltd., Hamilton, New Zealand.
López-Villalobos N, D Garrick. 2005. The breeding program for the genetic improvement of New Zealand dairy cattle: an example of a systematic approach to design breeding program. Simposio sobre recursos genéticos y bienestar animal aplicados a la producción animal, Sociedad Chilena de producción animal. SOCHIPA A.G. 11, 41-72.
López-Villalobos N, D Garrick, CW Holmes, HT Blair, RJ Spelman. 2000. Profitabilities of some mating system for dairy herds in New Zealand. J Dairy Sci 83, 144-153.
Misztal I, S Tsuruta, T Strabel, B Auvray, T Druet, D Lee. 2002. BLUPF90 and related program (BGF90). Proceedings 7th World Congress on Genetics Applied to Livestock Production, Montpellier, France, CD-Rom Communication 28, Pp 1-7.
Nelson DL, MM Cox. 2005. Lehninger principios de bioquímica. 4th ed. Omega, Barcelona, España, Pp 343-368.
Parga J, F Lanuza. 2006. Suplementación de vacas lecheras a pastoreo. En: Navarro H, Siebald E, Celis S (eds). Manual de producción de leche para pequeños y medianos productores. Boletín INIA N°148, Pp 57-66.
Prendiville R, KM Pierce, F Buckley. 2009. An evaluation of production efficiencies among lactating Holstein-Friesian, Jersey, and Jersey x Holstein-Friesian cows at pasture. J Dairy Sci 92, 6176-6185.
Sneddon NW. 2011. Review on differences for production, reproduction and environment impact traits between Holstein-Friesian, Jersey and Holstein-Friesian x Jersey crossbred cows. Postgraduate Thesis, Massey University, Palmerston North, New Zealand.
Teuber N. 2009. Praderas permanentes en las zonas lecheras de Chile. Curvas de crecimiento, distribución y producción. Imprenta América, Osorno, Chile.
Welch BL. 1938. The significance of the difference between two means when the population variances are unequal. Biometrika 29, 350-362.
© 2017 •
Facultad de Ciencias Veterinarias, Universidad Austral de Chile
Teléfono/Fax: 56 63 221459 • Casilla 567 • Campus Isla Teja S/N • Valdivia • Chile