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

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

Genome-wide association study for claw disorders and trimming status in dairy cattle.

Performing a genome-wide association study (GWAS) might add to a better understanding of the development of claw disorders and the need for trimming. Therefore, the aim of the current study was to perform a GWAS on claw disorders and trimming status and to validate the results for claw disorders based on an independent data set. Data consisted of 20,474 cows with phenotypes for claw disorders and 50,238 cows with phenotypes for trimming status. Recorded claw disorders used in the current study were double sole (DS), interdigital hyperplasia (IH), sole hemorrhage (SH), sole ulcer (SU), white line separation (WLS), a combination of infectious claw disorders consisting of (inter-)digital dermatitis and heel erosion, and a combination of laminitis-related claw disorders (DS, SH, SU, and WLS). Of the cows with phenotypes for claw disorders, 1,771 cows were genotyped and these cow data were used for the GWAS on claw disorders. A SNP was considered significant when the false discovery rate≤0.05 and suggestive when the false discovery rate≤0.20. An independent data set of 185 genotyped bulls having at least 5 daughters with phenotypes (6,824 daughters in total) for claw disorders was used to validate significant and suggestive SNP detected based on the cow data. To analyze the trait "trimming status" (i.e., the need for claw trimming), a data set with 327 genotyped bulls having at least 5 daughters with phenotypes (18,525 daughters in total) was used. Based on the cow data, in total 10 significant and 45 suggestive SNP were detected for claw disorders. The 10 significant SNP were associated with SU, and mainly located on BTA8. The suggestive SNP were associated with DS, IH, SU, and laminitis-related claw disorders. Three of the suggestive SNP were validated in the data set of 185 bulls, and were located on BTA13, BTA14, and BTA17. For infectious claw disorders, SH, and WLS, no significant or suggestive SNP associations were detected. For trimming status, 1 significant and 1 suggestive SNP were detected, both located close to each other on BTA15. Some significant and suggestive SNP were located close to SNP detected in studies on feet and leg conformation traits. Genes with major effects could not be detected and SNP associations were spread across the genome, indicating that many SNP, each explaining a small proportion of the genetic variance, influence claw disorders. Therefore, to reduce the incidence of claw disorders by breeding, genomic selection is a promising approach.

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.

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

Phenotypic and genetic relationships of bovine natural antibodies binding keyhole limpet hemocyanin in plasma and milk.

To improve the health status (resilience) of dairy cows, levels of natural antibodies (NAb) might be useful. The objective of the present study was to compare levels and to estimate genetic parameters for NAb measured in milk and plasma samples. Titers of NAb IgM and IgG isotype-binding keyhole limpet hemocyanin of 2,919 cows, in both plasma and milk, were measured using ELISA. Analysis revealed that NAb levels in milk significantly increased with parity, whereas they remained constant in plasma. Moderate positive phenotypic correlations were found between NAb levels in milk and in plasma: 0.18 for IgG and 0.40 for IgM. This indicates that NAb from milk and plasma might reflect different aspects of dairy cow health status. However, high genetic correlations were found for NAb in milk and plasma: 0.81 for IgG and 0.79 for IgM. Heritabilities (SE in parentheses) for NAb measured in plasma [0.15 (0.05) for IgG and 0.25 (0.06) for IgM] were higher than heritabilities of NAb measured in milk [0.08 (0.03) for IgG and 0.23 (0.05) for IgM]. Our results indicate that NAb measured in milk and plasma are heritable and likely have a common genetic background, suggesting that NAb levels measured in milk might be useful for genetic improvement of disease resistance.

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.

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.

A quantitative trait locus on Bos taurus autosome 17 explains a large proportion of the genetic variation in de novo synthesized milk fatty acids.

A genomic region associated with milk fatty acid (FA) composition has been detected on Bos taurus autosome (BTA)17 based on 50,000 (50K) single nucleotide polymorphism (SNP) genotypes. The aim of our study was to fine-map BTA17 with imputed 777,000 (777 K) SNP genotypes to identify candidate genes associated with milk FA composition. Phenotypes consisted of gas chromatography measurements of 14 FA based on winter and summer milk samples. Phenotypes and genotypes were available on 1,640 animals in winter milk, and on 1,581 animals in summer milk samples. Single-SNP analyses showed that several SNP in a region located between 29.0 and 34.0 Mbp were in strong association with C6:0, C8:0, and C10:0. This region was further characterized based on haplotypes. In summer milk samples, for example, these haplotypes explained almost 10% of the genetic variance in C6:0, 9% in C8:0, 3.5% in C10:0, 1.8% in C12:0, and 0.9% in C14:0. Two groups of haplotypes with distinct predicted effects could be defined, suggesting the presence of one causal variant. Predicted haplotype effects tended to increase from C6:0 to C14:0; however, the proportion of genetic variance explained by the haplotypes tended to decrease from C6:0 to C14:0. This is an indication that the quantitative trait locus (QTL) region is involved either in the elongation process or in early termination of de novo synthesized FA. Although many genes are present in this QTL region, most of these genes on BTA17 have not been characterized yet. The strongest association was found close to the progesterone receptor membrane component 2 (PGRMC2) gene, which has not yet been associated with milk FA composition. Therefore, no clear candidate gene associated with milk FA composition could be identified for this QTL.

Optimizing the design of small-sized nucleus breeding programs for dairy cattle with minimal performance recording.

Dairy cattle breeding programs in developing countries are constrained by minimal and erratic pedigree and performance recording on cows on commercial farms. Small-sized nucleus breeding programs offer a viable alternative. Deterministic simulations using selection index theory were performed to determine the optimum design for small-sized nucleus schemes for dairy cattle. The nucleus was made up of 197 bulls and 243 cows distributed in 8 non-overlapping age classes. Each year 10 sires and 100 dams were selected to produce the next generation of male and female selection candidates. Conception rates and sex ratio were fixed at 0.90 and 0.50, respectively, translating to 45 male and 45 female candidates joining the nucleus per year. Commercial recorded dams provided information for genetic evaluation of selection candidates (bulls) in the nucleus. Five strategies were defined: nucleus records only [within-nucleus dam performance (DP)], progeny records in addition to nucleus records [progeny testing (PT)], genomic information only [genomic selection (GS)], dam performance records in addition to genomic information (GS+DP), and progeny records in addition to genomic information (GS+PT). Alternative PT, GS, GS+DP, and GS+PT schemes differed in the number of progeny per sire and size of reference population. The maximum number of progeny records per sire was 30, and the maximum size of the reference population was 5,000. Results show that GS schemes had higher responses and lower accuracies compared with other strategies, with the higher response being due to shorter generation intervals. Compared with similar sized progeny-testing schemes, genomic-selection schemes would have lower accuracies but these are offset by higher responses per year, which might provide additional incentive for farmers to participate in recording.