Genetic origin, admixture and population history of aurochs (Bos primigenius) and primitive European cattle.

This corrects the article DOI: 10.1038/hdy.2016.79.

Fine-mapping of BTA17 using imputed sequences for associations with de novo synthesized fatty acids in bovine milk.

The aim of this study was to fine-map a genomic region associated with milk fatty acids (FA) on Bos taurus autosome (BTA) 17. This genomic region has been discovered with 50,000 (50k) single nucleotide polymorphisms (SNP) imputed to 777,000 (777k) SNP. In this study, high-density genotypes were imputed to whole-genome sequences level to identify candidate gene(s) associated with milk FA composition on BTA17. Phenotypes and genotypes were available for 1,640 cows sampled in winter, and for 1,581 cows sampled in summer. Phenotypes consisted of gas chromatography measurements in winter and in summer milk samples of 6 individual FA and the indicator of de novo synthesis, C6:0-C14:0. Genotypes consisted of imputed 777k SNP, and 89 sequenced ancestors of the population of genotyped cows. In addition, 450 whole-genome sequences from the 1,000 Bull Genome Consortium were available. Using 495 Holstein-Friesian sequences as a reference population, the 777k SNP genotypes of the cows were imputed to sequence level. We then applied single-variant analyses with an animal model, and identified thousands of significant associations with C6:0, C8:0, C10:0, C12:0, C14:0, and C6:0-C14:0. For C8:0 in summer milk samples, the genomic region located between 29 and 34 Mbp on BTA17 revealed a total of 646 significant associations. The most significant associations [-log10(P-value) = 7.82] were 8 SNP in perfect linkage disequilibrium. After fitting one of these 8 SNP as a fixed effect in the model, and re-running the single-variant analyses, no further significant associations were found for any of the 6 FA or C6:0-C14:0. These findings suggest that one polymorphism underlying this QTL on BTA17 influences multiple de novo synthesized milk FA. Thirteen genes in the QTL region were identified and analyzed carefully. Six out of the 8 SNP that showed the strongest associations were located in the La ribonucleoprotein domain family, member 1B (LARP1B) gene, and we suggest LARP1B as a primary candidate gene. Another gene of interest for this QTL region might be PKL4. None of these suggested candidate genes have previously been associated with milk fat synthesis or milk FA composition.

Multiple criteria decision-making process to derive consensus desired genetic gains for a dairy cattle breeding objective for diverse production systems.

Dairy cattle industries contribute to food and nutrition security and are a source of income for numerous households in many developing countries. Selective breeding can enhance efficiency in these industries. Developing dairy industries are characterized by diverse production and marketing systems. In this paper, we use weighted goal aggregating procedure to derive consensus trait preferences for different producer categories and processors. We based the study on the dairy industry in Kenya. The analytic hierarchy process was used to derive individual preferences for milk yield (MY), calving interval (CIN), production lifetime (PLT), mature body weight (MBW), and fat yield (FY). Results show that classical classification of production systems into large-scale and smallholder systems does not capture all differences in trait preferences. These differences became apparent when classification was based on productivity at the individual animal level, with high and low intensity producers and processors as the most important groups. High intensity producers had highest preferences for PLT and MY, whereas low intensity producers had highest preference for CIN and PLT; processors preferred MY and FY the most. The highest disagreements between the groups were observed for FY, PLT, and MY. Individual and group preferences were aggregated into consensus preferences using weighted goal programming. Desired gains were obtained as a product of consensus preferences and percentage genetic gains (G%). These were 2.42, 0.22, 2.51, 0.15, and 0.87 for MY, CIN, PLT, MBW, and FY, respectively. Consensus preferences can be used to derive a single compromise breeding objective for situations where the same genetic resources are used in diverse production and marketing circumstances.

Effect of production quotas on economic and environmental values of growth rate and feed efficiency in sea cage fish farming.

In sea cage fish farming, production quotas aim to constrain the impact of fish farming on the surrounding ecosystem. It is unknown how these quotas affect economic profitability and environmental impact of genetic improvement. We combined bioeconomic modelling with life cycle assessment (LCA) to calculate the economic (EV) and environmental (ENV) values of thermal growth coefficient (TGC) and feed conversion ratio (FCR) of sea bass reared in sea cages, given four types of quota commonly used in Europe: annual production (Qprod), annual feed distributed (Qannual_feed), standing stock (Qstock), and daily feed distributed (Qdaily_feed). ENV were calculated for LCA impact categories climate change, eutrophication and acidification. ENV were expressed per ton of fish produced per year (ENV(fish)) and per farm per year (ENV(farm)). Results show that irrespective of quota used, EV of FCR as well as ENV(fish) and ENV(farm) were always positive, meaning that improving FCR increased profit and decreased environmental impacts. However, the EV and the ENV(fish) of TGC were positive only when quota was Qstock or Qdaily_feed. Moreover, the ENV(farm) of TGC was negative in Qstock and Qdaily_feed quotas, meaning that improving TGC increased the environmental impact of the farm. We conclude that Qstock quota and Qdaily_feed quota are economically favorable to a genetic improvement of TGC, a major trait for farmers. However, improving TGC increases the environmental impact of the farm. Improving FCR represents a good opportunity to balance out this increase but more information on its genetic background is needed to develop breeding programs improving FCR.

Economic evaluation of progeny-testing and genomic selection schemes for small-sized nucleus dairy cattle breeding programs in developing countries.

In developing countries minimal and erratic performance and pedigree recording impede implementation of large-sized breeding programs. Small-sized nucleus programs offer an alternative but rely on their economic performance for their viability. We investigated the economic performance of 2 alternative small-sized dairy nucleus programs [i.e., progeny testing (PT) and genomic selection (GS)] over a 20-yr investment period. The nucleus was made up of 453 male and 360 female animals distributed in 8 non-overlapping age classes. Each year 10 active sires and 100 elite dams were selected. Populations of commercial recorded cows (CRC) of sizes 12,592 and 25,184 were used to produce test daughters in PT or to create a reference population in GS, respectively. Economic performance was defined as gross margins, calculated as discounted revenues minus discounted costs following a single generation of selection. Revenues were calculated as cumulative discounted expressions (CDE, kg) × 0.32 (€/kg of milk) × 100,000 (size commercial population). Genetic superiorities, deterministically simulated using pseudo-BLUP index and CDE, were determined using gene flow. Costs were for one generation of selection. Results show that GS schemes had higher cumulated genetic gain in the commercial cow population and higher gross margins compared with PT schemes. Gross margins were between 3.2- and 5.2-fold higher for GS, depending on size of the CRC population. The increase in gross margin was mostly due to a decreased generation interval and lower running costs in GS schemes. In PT schemes many bulls are culled before selection. We therefore also compared 2 schemes in which semen was stored instead of keeping live bulls. As expected, semen storage resulted in an increase in gross margins in PT schemes, but gross margins remained lower than those of GS schemes. We conclude that implementation of small-sized GS breeding schemes can be economically viable for developing countries.

Genetic correlations between dressage, show jumping and studbook-entry inspection traits in a process of specialization in Dutch Warmblood horses.

Sport performance in dressage and show jumping are two important traits in the breeding goals of many studbooks. To determine the optimum selection scheme for jumping and dressage, knowledge is needed on the genetic correlation between both disciplines and between traits measured early in life and performance in competition in each discipline. This study aimed to estimate genetic parameters to support decision-making on specialization of breeding horses for dressage and show jumping in Dutch warmblood horses. Genetic correlations between performance of horses in dressage and show jumping were estimated as well as the genetic correlation between traits recorded during studbook-entry inspections and performance in dressage and show jumping competitions. The information on competition comprised the performance of 82 694 horses in dressage and 62 072 horses in show jumping, recorded in the period 1993-2012. For 26 056 horses, information was available for both disciplines. The information on traits recorded at studbook-entry inspections comprised 62 628 horses, recorded in the period 1992-2013. Genetic parameters were estimated from the whole dataset and from a subset without horses recorded in both disciplines. Additionally, the genetic parameters were estimated in three different time periods defined by horses' birth year. The genetic correlation between dressage and show jumping in the whole dataset was -0.23, and it was -0.03 when it was estimated from horses recorded in only one discipline. The genetic correlation between dressage and show jumping was more negative in the most recent time period in all the cases. The more negative correlation between disciplines in more recent time periods was not reflected in changes in the correlations between competitions traits and the traits recorded in the studbook-first inspection. These results suggest that a breeding programme under specialization might be most effective defining two separate aggregate breeding goals for each of the disciplines.

Genetic origin, admixture and population history of aurochs (Bos primigenius) and primitive European cattle.

The domestication of taurine cattle initiated ~10 000 years ago in the Near East from a wild aurochs (Bos primigenius) population followed by their dispersal through migration of agriculturalists to Europe. Although gene flow from wild aurochs still present at the time of this early dispersion is still debated, some of the extant primitive cattle populations are believed to possess the aurochs-like primitive features. In this study, we use genome-wide single nucleotide polymorphisms to assess relationship, admixture patterns and demographic history of an ancient aurochs sample and European cattle populations, several of which have primitive features and are suitable for extensive management. The principal component analysis, the model-based clustering and a distance-based network analysis support previous works suggesting different histories for north-western and southern European cattle. Population admixture analysis indicates a zebu gene flow in the Balkan and Italian Podolic cattle populations. Our analysis supports the previous report of gene flow between British and Irish primitive cattle populations and local aurochs. In addition, we show evidence of aurochs gene flow in the Iberian cattle populations indicating wide geographical distribution of the aurochs. Runs of homozygosity (ROH) reveal that demographic processes like genetic isolation and breed formation have contributed to genomic variations of European cattle populations. The ROH also indicate recent inbreeding in southern European cattle populations. We conclude that in addition to factors such as ancient human migrations, isolation by distance and cross-breeding, gene flow between domestic and wild-cattle populations also has shaped genomic composition of European cattle populations.

Genome-wide association study for behavior, type traits, and muscular development in Charolais beef cattle.

Behavior, type traits, and muscular development are of interest for beef cattle breeding. Genome-wide association studies (GWAS) enable the identification of candidate genes, which enables gene-based selection and provides insight in the genetic architecture of these traits. The objective of the current study was to perform a GWAS for 3 behavior traits, 12 type traits, and muscular development in Charolais cattle. Behavior traits, including aggressiveness at parturition, aggressiveness during gestation period, and maternal care, were scored by farmers. Type traits, including udder conformation, teat, feet and legs, and locomotion, were scored by trained classifiers. Data used in the GWAS consisted of 3,274 cows with phenotypic records and genotyping information for 44,930 SNP. When SNP had a false discovery rate (FDR) smaller than 0.05, they were referred to as significant. When SNP had a FDR between 0.05 and 0.20, they were referred to as suggestive. Four significant and 12 suggestive regions were detected for aggressiveness during gestation, maternal care, udder balance, teat thinness, teat length, foot angle, foot depth, and locomotion. These 4 significant and 12 suggestive regions were not supported by other significant SNP in close proximity. No SNP with major effects were detected for behavior and type traits, and SNP associations for these traits were spread across the genome, suggesting that behavior and type traits were influenced by many genes, each explaining a small part of genetic variance. The GWAS identified 1 region on chromosome 2 significantly associated with muscular development, which included the myostatin gene (), which is known to affect muscularity. No other regions associated with muscular development were found. Results showed that the myostatin region associated with muscular development had pleiotropic effects on udder volume, teat thinness, rear leg, and leg angle.

Accounting for genetic architecture in single- and multipopulation genomic prediction using weights from genomewide association studies in pigs.

We studied the effect of including GWAS results on the accuracy of single- and multipopulation genomic predictions. Phenotypes (backfat thickness) and genotypes of animals from two sire lines (SL1, n = 1146 and SL3, n = 1264) were used in the analyses. First, GWAS were conducted for each line and for a combined data set (both lines together) to estimate the genetic variance explained by each SNP. These estimates were used to build matrices of weights (D), which was incorporated into a GBLUP method. Single population evaluated with traditional GBLUP had accuracies of 0.30 for SL1 and 0.31 for SL3. When weights were employed in GBLUP, the accuracies for both lines increased (0.32 for SL1 and 0.34 for SL3). When a multipopulation reference set was used in GBLUP, the accuracies were higher (0.36 for SL1 and 0.32 for SL3) than in single-population prediction. In addition, putting together the multipopulation reference set and the weights from the combined GWAS provided even higher accuracies (0.37 for SL1, and 0.34 for SL3). The use of multipopulation predictions and weights estimated from a combined GWAS increased the accuracy of genomic predictions.

Analysis of competition performance in dressage and show jumping of Dutch Warmblood horses.

Most Warmblood horse studbooks aim to improve the performance in dressage and show jumping. The Dutch Royal Warmblood Studbook (KWPN) includes the highest score achieved in competition by a horse to evaluate its genetic ability of performance. However, the records collected during competition are associated with some aspects that might affect the quality of the genetic evaluation based on these records. These aspects include the influence of rider, censoring and preselection of the data. The aim of this study was to quantify the impact of rider effect, censoring and preselection on the genetic analysis of competition data of dressage and show jumping of KWPN. Different models including rider effect were evaluated. To assess the impact of censoring, genetic parameters were estimated in data sets that differed in the degree of censoring. The effect of preselection on variance components was analysed by defining a binary trait (sport-status) depending on whether the horse has a competition record or not. This trait was included in a bivariate model with the competition trait and used all horses registered by KWPN since 1984. Results showed that performance in competition for dressage and show jumping is a heritable trait (h2 ~ 0.11-0.13) and that it is important to account for the effect of rider in the genetic analysis. Censoring had a small effect on the genetic parameter for highest performance achieved by the horse. A moderate heritability obtained for sport-status indicates that preselection has a genetic basis, but the effect on genetic parameters was relatively small.

Defining a breeding objective for Nile tilapia that takes into account the diversity of smallholder production systems.

In general, livestock and fish farming systems in developing countries tend to be highly diverse in terms of agro-ecological conditions and market orientation. There are no studies that have investigated if and how this diversity translates to varying preferences for breeding objective traits. This is particularly important for breeding programmes that are organized on a national level (e.g. government-supported nucleus breeding programmes). The aim of this study was to investigate whether Nile tilapia farmers with diverse production systems and economic constraints have different preferences for breeding objective traits. The second objective was to derive a consensus breeding goal, using weighted goal programming that could be used for a national breeding programme for Nile tilapia. A survey was conducted among 100 smallholder Nile tilapia farmers in Kenya to obtain preference values for traits of economic importance, by using multiple pairwise comparisons. Individual and group preference values were estimated using analytical hierarchy process. Low-income farmers preferred harvest weight, while medium- and high-income farmers preferred growth rate and survival. Grouping farmers according to market objective (fingerling production or fattening) showed that fingerling producers preferred growth rate and survival, while fattening farmers preferred harvest weight, height and thickness. Weighted goal programming was used to obtain consensus preference values, and these were used to derive desired gains for a breeding goal of a national breeding programme that takes into account the diversity of smallholder production systems.

Effects of the diacylglycerol o-acyltransferase 1 (DGAT1) K232A polymorphism on fatty acid, protein, and mineral composition of dairy cattle milk.

Several studies have described associations between the diacylglycerol o-acyltransferase 1 (DGAT1) K232A polymorphism and routinely collected milk production traits but not much is known about effects of the DGAT1 polymorphism on detailed milk composition. The aim of this study was to estimate effects of the DGAT1 polymorphism on milk fatty acid, protein, and mineral composition. We looked for effects that were significant and consistent in Danish Holstein Friesian (HF), Danish Jersey, and Dutch HF as these are likely to be true effects of the DGAT1 K232A polymorphism rather than being effects of linked loci. For fatty acid composition, significant and consistent effects of the DGAT1 polymorphism were detected on C14:0, C16:0, C15:0, C16:1, C18:1 cis-9, conjugated linoleic acid (CLA) cis-9,trans-11, C18:2 cis-9,cis-12, and C18:3 cis-9,cis-12,cis-15 content (percent by weight, wt/wt %). For C16:0, C16:1, and C18:1 cis-9, the DGAT1 polymorphism explained more than 10% of the phenotypic variation. Significant effects on milk protein composition in Dutch HF could not be confirmed in Danish Jersey or Danish HF. For mineral content, significant and consistent effects of the DGAT1 polymorphism on calcium, phosphorus, and zinc were detected. In the Dutch HF population, the contribution of the DGAT1 K232A polymorphism to phenotypic variance was 12.0% for calcium, 8.3% for phosphorus, and 6.1% for zinc. Different from effects on fatty acid composition, effects of the DGAT1 polymorphism on yields of long-chain fatty acids C18:1 cis-9, CLA cis-9,trans-11, C18:2 cis-9,cis-12, and C18:3 cis-9,cis-12,cis-15 were not significant. This indicates that effects of DGAT1 on these fatty acids are indirect, not direct, effects: DGAT1 affects de novo synthesis of fatty acids and, consequently, the contribution of the long-chain fatty acids to total fat is decreased. In addition, effects of the DGAT1 polymorphism on yields of Ca, P, and Zn were not significant, which indicates that effects on these minerals are the result of indirect rather than direct effects of DGAT1: effects on calcium, phosphorus, and zinc content can be explained by effects of DGAT1 on milk volume. The reported effects of the DGAT1 polymorphism on fatty acid and mineral composition of milk are substantial and therefore relevant for milk quality.

Accuracy of genome-enabled prediction exploring purebred and crossbred pig populations.

Pig breeding companies keep relatively small populations of pure sire and dam lines that are selected to improve the performance of crossbred animals. This design of the pig breeding industry presents challenges to the implementation of genomic selection, which requires large data sets to obtain highly accurate genomic breeding values. The objective of this study was to evaluate the impact of different reference sets (across population and multipopulation) on the accuracy of genomic breeding values in 3 purebred pig populations and to assess the potential of using crossbreed performance in genomic prediction. Data consisted of phenotypes and genotypes on animals from 3 purebred populations (sire line [SL] 1, = 1,146; SL2, = 682; and SL3, = 1,264) and 3 crossbred pig populations (Terminal cross [TER] 1, = 183; TER2, = 106; and TER3, = 177). Animals were genotyped using the Illumina Porcine SNP60 Beadchip. For each purebred population, within-, across-, and multipopulation predictions were considered. In addition, data from the paternal purebred populations were used as a reference set to predict the performance of crossbred animals. Backfat thickness phenotypes were precorrected for fixed effects and subsequently included in the genomic BLUP model. A genomic relationship matrix that accounted for the differences in allele frequencies between lines was implemented. Accuracies of genomic EBV obtained within the 3 different sire lines varied considerably. For within-population prediction, SL1 showed higher values (0.80) than SL2 (0.61) and SL3 (0.67). Multipopulation predictions had accuracies similar to within-population accuracies for the validation in SL1. For SL2 and SL3, the accuracies of multipopulation prediction were similar to the within-population prediction when the reference set was composed by 900 animals (600 of the target line plus 300 of another line). For across-population predictions, the accuracy was mostly close to zero. The accuracies of predicting crossbreed performance were similar for the 3 different crossbred populations (ranging from 0.25 to 0.29). In summary, the differences in accuracy of the within-population scenarios may be due to line divergences in heritability and genetic architecture of the trait. Within- and multipopulation predictions yield similar accuracies. Across-population prediction accuracy was negligible. The moderate accuracy of prediction of crossbreed performance appears to be a result of the relationship between the crossbreed and its parental lines.

Genetic parameters for large-scale behavior traits and type traits in Charolais beef cows.

In the last decades, beef cattle breeding mainly focused on improving production and reproduction traits. Nowadays, there is a growing interest to include behavior and type traits in the breeding goal. There is an interest in behavior as it is associated with human safety and workability and in type traits as they might be associated with longevity of cows. The objective of the current study was to estimate the heritability for behavior and type traits in Charolais and to estimate the genetic correlations among these traits. Behavior traits, including aggressiveness at parturition, aggressiveness during gestation period, and maternal care, were scored by farmers using an on-farm recording system to enable large-scale collection of phenotypes. Type traits, including udder traits ( = 3), teat traits (3), feet and leg traits (5), and locomotion (1), were scored by 10 trained classifiers. Data was available on 6,649 cows in parity 1 to 12 and located in 380 herds. Results showed that differences between herds explained up to 23% of the total phenotypic variance in behavior traits. This might be due to differences in management or to consistent differences in scoring between farmers. Aggressiveness at parturition had higher heritability (0.19) and higher genetic coefficient of variation (CV = 11%) than aggressiveness during gestation ( = 0.06 and CV = 4%) and maternal care ( = 0.02 and CV = 2%). Heritabilities for udder traits (0.14 to 0.20) and teat traits (0.17 to 0.35) were higher than for feet and leg traits (0.02 to 0.19). Genetic coefficients of variation for udder and teat traits were also higher (up to 21%) than for feet and leg traits (up to 11%). Strong genetic correlations were found between behavior traits (with absolute values from 0.71 to 0.98). The genetic correlations indicate that it is difficulty to simultaneously improve maternal care and reduce aggressiveness. We concluded that there are good opportunities to implement selection for improved udder and teat traits and against aggressiveness at parturition using a simple on-farm recording system of behavior.

Short communication: Genetic study of methane production predicted from milk fat composition in dairy cows.

Dairy cows produce enteric methane, a greenhouse gas with 25 times the global warming potential of CO2. Breeding could make a permanent, cumulative, and long-term contribution to methane reduction. Due to a lack of accurate, repeatable, individual methane measurements needed for breeding, indicators of methane production based on milk fatty acids have been proposed. The aim of the present study was to quantify the genetic variation for predicted methane yields. The milk fat composition of 1,905 first-lactation Dutch Holstein-Friesian cows was used to investigate 3 different predicted methane yields (g/kg of DMI): Methane1, Methane2, and Methane3. Methane1 was based on the milk fat proportions of C17:0anteiso, C18:1 rans-10+11, C18:1 cis-11, and C18:1 cis-13 (R(2)=0.73). Methane2 was based on C4:0, C18:0, C18:1 trans-10+11, and C18:1 cis-11 (R(2)=0.70). Methane3 was based on C4:0, C6:0, and C18:1 trans-10+11 (R(2)=0.63). Predicted methane yields were demonstrated to be heritable traits, with heritabilities between 0.12 and 0.44. Breeding can, thus, be used to decrease methane production predicted based on milk fatty acids.

Genetic relationships between claw health traits of dairy cows in different parities, lactation stages, and herds with different claw disorder frequencies.

Claw disorders affect cow welfare and profitability of farms. Not only claw disorders but also the need for trimming, known as trimming status, has been shown to be heritable. Limited knowledge is available on whether claw health traits (claw disorders and trimming status) are genetically the same trait in different parities, lactation stages, or in herds with low or high frequency of claw disorders. The aim of the current study was to estimate frequencies, heritabilities, and genetic correlations of claw health traits measured in different parities (first vs. later parities), in different lactation stages (early vs. late lactation), and in herds with different frequencies of claw disorders (low vs. high frequency). Analyses revealed that heritabilities measured in different parities, lactation stages, or herds with different trait frequencies are similar for most claw health traits. Also, genetic correlations (rg) for most claw health traits were not found to be different from unity for traits in different parities, lactation stages, or herds with different trait frequencies. Sole hemorrhage and infectious lesions were genetically different traits in first or later parities (rg=0.29±0.31 and 0.66±0.15, respectively). White line separation and infectious lesions were genetically different in early and late lactation (rg=0.53±0.20 and 0.69±0.13, respectively), and sole ulcer was genetically different in herds with low or high frequency of sole ulcer (rg=0.75±0.14). In our analysis, we did not find convincing evidence to supports the use of multiple trait models for the analysis of claw health traits treating them as different traits in different parities, lactation stages, and herds with different claw disorder frequency.

Effects of the DGAT1 polymorphism on test-day milk production traits throughout lactation.

Several studies have shown that the diacylglycerol O-acyltransferase 1 (DGAT1) K232A polymorphism has a major effect on milk production traits. It is less clear how effects of DGAT1 on milk production traits change throughout lactation, if dominance effects of DGAT1 are relevant, and whether DGAT1 also affects lactose content, lactose yield, and total energy output in milk. Results from this study, using test-day records of 3 subsequent parities of around 1,800 cows, confirm previously reported effects of the DGAT1 polymorphism on milk, fat, and protein yield, as well as fat and protein content. In addition, we found significant effects of the DGAT1 polymorphism on lactose content and lactose yield. No significant effects on somatic cell score were detected. The effect of DGAT1 on total energy excreted in milk was only significant in parity 1 and is mainly due to a higher energy output in milk of heterozygous cows. Significant but relatively small dominance effects of DGAT1 on fat content and yield were detected, which are of little practical relevance. Significant DGAT1 by lactation stage interaction was detected for milk yield, lactose yield, fat content, and protein content, indicating that the effect of the DGAT1 polymorphism changes during lactation. In general, the DGAT1 effect shows a large increase during early lactation (from the start of lactation to d 50 to 150) and tends to decrease later in lactation. No DGAT1 by lactation stage interaction for fat yield was observed. Similar to DGAT1, effects of other genes also might vary throughout lactation and, therefore, using longitudinal models is recommended.

Effect of specialization on genetic parameters of studbook-entry inspection in Dutch Warmblood horses.

Recent studies on data from the Dutch Warmblood Studbook (KWPN) have shown that the ongoing specialization of horses for either dressage (DH) or show jumping (JH) has led to a decreasing genetic relationship between the two subpopulations. The aim of this study was to analyse the effect of the specialization process on the genetic parameters of traits measured in the studbook-entry inspection of KWPN during the last fifteen years. Data from 18,125 DH and 23,800 JH recorded from 1998 until 2013 were used to analyse 13 traits scored in both DH and JH. Analyses were performed in a Bayesian framework. Firstly, variance components were estimated based on the whole data set. Secondly, genetic correlations between traits measured in DH or JH were estimated using bivariate analyses. Thirdly, three time periods were defined and genetic correlations between subpopulations were estimated within each period. Heritability was moderate (0.17-0.39) for both DH and JH. Genetic correlations between traits measured in DH or JH were not different from one considering the posterior standard deviation of the estimation; however, in most of the traits, a clear trend in reduction of the genetic correlation for traits expressed in DH and JH and an increase in their posterior standard deviation for recent years was observed. These results suggest that specialization could lead to differences in traits measured in DH and JH in the recent years.

Effect of feed-related farm characteristics on relative values of genetic traits in dairy cows to reduce greenhouse gas emissions along the chain.

Breeding has the potential to reduce greenhouse gas (GHG) emissions from dairy farming. Evaluating the effect of a 1-unit change (i.e., 1 genetic standard deviation improvement) in genetic traits on GHG emissions along the chain provides insight into the relative importance of genetic traits to reduce GHG emissions. Relative GHG values of genetic traits, however, might depend on feed-related farm characteristics. The objective of this study was to evaluate the effect of feed-related farm characteristics on GHG values by comparing the values of milk yield and longevity for an efficient farm and a less efficient farm. The less efficient farm did not apply precision feeding and had lower feed production per hectare than the efficient farm. Greenhouse gas values of milk yield and longevity were calculated by using a whole-farm model and 2 different optimization methods. Method 1 optimized farm management before and after a change in genetic trait by maximizing labor income; the effect on GHG emissions (i.e., from production of farm inputs up to the farm gate) was considered a side effect. Method 2 optimized farm management after a change in genetic trait by minimizing GHG emissions per kilogram of milk while maintaining labor income and milk production at least at the level before the change in trait; the effect on labor income was considered a side effect. Based on maximizing labor income (method 1), GHG values of milk yield and longevity were, respectively, 279 and 143kg of CO2 equivalents (CO2e)/unit change per cow per year on the less efficient farm, and 247 and 210kg of CO2e/unit change per cow per year on the efficient farm. Based on minimizing GHG emissions (method 2), GHG values of milk yield and longevity were, respectively, 538 and 563kg of CO2e/unit change per cow per year on the less efficient farm, and 453 and 441kg of CO2e/unit change per cow per year on the efficient farm. Sensitivity analysis showed that, for both methods, the absolute effect of a change in genetic trait depends on model inputs, including prices and emission factors. Substantial changes in relative importance between traits due to a change in model inputs occurred only in case of maximizing labor income. We concluded that assumptions regarding feed-related farm characteristics affect the absolute level of GHG values, as well as the relative importance of traits to reduce emissions when using a method based on maximizing labor income. This is because optimizing farm management based on maximizing labor income does not give any incentive for lowering GHG emissions. When using a method based on minimizing GHG emissions, feed-related farm characteristics affected the absolute level of the GHG values, but the relative importance of the traits scarcely changed: at each level of efficiency, milk yield and longevity were equally important.

High natural antibody titers of indigenous chickens are related with increased hazard in confinement.

Natural antibody (NAb) levels and survival rates were evaluated in 4 breeds of laying hens in Ethiopia: indigenous, improved indigenous, exotic layer, and crossbred. Titers of NAb isotypes IgG and IgM binding keyhole limpet hemocyanin (KLH) in serum were measured at 20, 26, 35, and 45 wk age. Repeated-measure ANOVA showed that IgG and IgM levels vary with time within each breed (P < 0.05). Indigenous chickens had significantly (P < 0.05) higher NAb levels at all ages. The Cox proportional hazard analysis showed increased hazard with increased levels of NAbs in the exotic layers (P < 0.05). However, the reduced hazards with increased levels of NAbs were not significant in the improved indigenous and crossbred chickens. Indigenous chickens showed increased hazard with increasing levels of NAb (P > 0.05). We concluded that not only the NAb levels but also the effect of Nabs on survival vary between indigenous and improved breeds. The results indicate that NAb levels are associated with survival in elite (improved) breeds, but are associated with increased hazard in indigenous chickens.