Linking variation in the casein fraction and salt composition to casein micelle size in milk of Dutch dairy goats.

The casein (CN) composition, salt composition and micelle size varies largely between milk samples of individual animals. In goats, the link between those casein characteristics are unknown and could provide useful insights into goat casein micelle structure. In this study, the casein- and salt composition of 42 individual Dutch goats from 17 farms was studied and linked to casein micelle size. Micelle size, proportions of individual caseins, and protein content were associated with each other. Milk with smaller casein micelles was higher in protein content, salt content, and proportion of αs1-CN, but lower in αs2-CN and ß-CN. The higher salt content in milk with small casein micelles was mainly attributed to a higher protein content, but changes in casein composition might additionally contribute to differences in mineralization. The non-sedimentable casein content in goat milk correlated with non-sedimentable fractions of ß-CN and κ-CN and was independent of micelle size. Between large and small casein micelles, goat casein micelles showed more differences in casein and salt composition than bovine micelles, indicating differences in internal structure. Nevertheless, the casein mineralization in goat milk was similar to casein mineralization in bovine milk, indicating that mineralization of casein micelles follows a general principle. These results can help to better understand how composition and micelle structure in goat milk are related to each other, which may be useful to improve processing and product properties of goat milk in the future.

Genetic parameters, reciprocal cross differences, and age-related heterosis of egg-laying performance in chickens.

BACKGROUND: Egg-laying performance is economically important in poultry breeding programs. Crossbreeding between indigenous and elite commercial lines to exploit heterosis has been an upward trend in traditional layer breeding for niche markets. The objective of this study was to analyse the genetic background and to estimate the heterosis of longitudinal egg-laying traits in reciprocal crosses between an indigenous Beijing-You and an elite commercial White Leghorn layer line. Egg weights were measured for the first three eggs, monthly from 28 to 76 weeks of age, and at 86 and 100 weeks of age. Egg quality traits were measured at 32, 54, 72, 86, and 100 weeks of age. Egg production traits were measured from the start of lay until 43, 72, and 100 weeks of age. Heritabilities and phenotypic and genetic correlations were estimated. Heterosis was estimated as the percentage difference of performance of a crossbred from that of the parental average. Reciprocal cross differences were estimated as the difference between the reciprocal crossbreds as a percentage of the parental average. RESULTS: Estimates of heritability of egg weights ranged from 0.29 to 0.75. Estimates of genetic correlations between egg weights at different ages ranged from 0.72 to 1.00. Estimates of heritability for cumulative egg numbers until 43, 72, and 100 weeks of age were around 0.15. Estimates of heterosis for egg weight and cumulative egg number increased with age, ranging from 1.0 to 9.0% and from 1.4 to 11.6%, respectively. From 72 to 100 weeks of age, crossbreds produced more eggs per week than the superior parent White Leghorn (3.5 eggs for White Leghorn, 3.8 and 3.9 eggs for crossbreds). Heterosis for eggshell thickness ranged from 2.7 to 6.6% when using Beijing-You as the sire breed. No significant difference between reciprocal crosses was observed for the investigated traits, except for eggshell strength at 54 weeks of age. CONCLUSIONS: The heterosis was substantial for egg weight and cumulative egg number, and increased with age, suggesting that non-additive genetic effects are important in crossbreds between the indigenous and elite breeds. Generally, the crossbreds performed similar to or even outperformed the commercial White Leghorns for egg production persistency.

Inbreeding depression and its effect on sperm quality traits in Pietrain pigs.

In most cases, inbreeding is expected to have unfavourable effects on traits in livestock. The consequences of inbreeding depression could be substantial, primarily in reproductive and sperm quality traits, and thus lead to decreased fertility. Therefore, the objectives of this study were (i) to compute inbreeding coefficients using pedigree (FPED ) and genomic data based on runs of homozygosity (ROH) in the genome (FROH ) of Austrian Pietrain pigs, and (ii) to assess inbreeding depression on four sperm quality traits. In total, 74,734 ejaculate records from 1034 Pietrain boars were used for inbreeding depression analyses. Traits were regressed on inbreeding coefficients using repeatability animal models. Pedigree-based inbreeding coefficients were lower than ROH-based inbreeding values. The correlations between pedigree and ROH-based inbreeding coefficients ranged from 0.186 to 0.357. Pedigree-based inbreeding affected only sperm motility while ROH-based inbreeding affected semen volume, number of spermatozoa, and motility. For example, a 1% increase in pedigree inbreeding considering 10 ancestor generations (FPED10 ) was significantly (p < 0.05) associated with a 0.231% decrease in sperm motility. Almost all estimated effects of inbreeding on the traits studied were unfavourable. It is advisable to properly manage the level of inbreeding to avoid high inbreeding depression in the future. Further, analysis of effects of inbreeding depression for other traits, including growth and litter size for the Austrian Pietrain population is strongly advised.

Divergent selection for natural antibodies in poultry in the presence of a major gene.

BACKGROUND: Natural antibodies (NAb) are antibodies that are present in a healthy individual without requiring previous exposure to an exogenous antigen. Selection for high NAb levels might contribute to improved general disease resistance. Our aim was to analyse the genetic background of NAb based on a divergent selection experiment in poultry, and in particular the effect of a polymorphism in the TLR1A gene. METHODS: The study population consisted of a base population from a commercial pure-bred elite white leghorn layer line and seven generations of birds from a High and Low selection line. Birds were selected for total KLH-binding NAb titer (IgTotal). An enzyme-linked immunosorbent assay was performed to determine NAb titers in blood plasma for IgTotal and the antibody isotypes IgM and IgG. NAb titers were available for 10,878 birds. Genotypes for a polymorphism in TLR1A were determined for chickens in generations 5, 6 and 7. The data were analysed using mixed linear animal models. RESULTS: The heritability estimate for IgM was 0.30 and higher than that for IgG and IgTotal (0.12). Maternal environmental effects explained 2 to 3% of the phenotypic variation in NAb. Selection for IgTotal resulted in a genetic difference between the High and Low line of 2.4 titer points (5.1 genetic standard deviation) in generation 7. For IgM, the selection response was asymmetrical and higher in the Low than the High line. The frequency of the TLR1A C allele was 0.45 in the base population and 0.66 and 0.04 in generation 7 of the High and Low line, respectively. The TLR1A polymorphism had large and significant effects on IgTotal and IgM. Estimated genotypic effects suggest full dominance of the TLR1A C allele. Significant TLR1A by generation interactions were detected for IgM and IgTotal. CONCLUSIONS: The effect of a polymorphism in the TLR1A gene on IgTotal and IgM NAb was confirmed. Furthermore, we provide experimental verification of changes in allele frequencies at a major gene with dominant gene action on a quantitative trait that is subjected to mass selection. TLR1A by generation interactions indicate sensitivity to environmental factors.

Genetic parameters of colostrum and calf serum antibodies in Swedish dairy cattle.

BACKGROUND: A sufficient IgG content in the colostrum is essential for the newborn calf, as it provides passive immunity which substantially affects the probability of survival during rearing. Failure of passive transfer (FPT) occurs when a calf does not absorb enough antibodies from the colostrum and is defined by an IgG concentration in calf serum lower than 10 g/L. Apart from delayed access to colostrum, FPT can be due to a low production of IgG in the mother or poor IgG absorption by the calf. The aim of this study was to estimate the genetic background of antibody levels and indicator traits for antibodies in the colostrum and calf serum, and their correlation with milk production. RESULTS: Colostrum data were available for 1340 dairy cows with at least one calving and calf serum data were available for 886 calves from these cows. Indicator traits for antibody concentrations were estimated using refractometry (a digital Brix refractometer for colostrum and an optical refractometer for serum), and enzyme-linked immunosorbent assays (ELISA) were used to determine the levels of total IgG and natural antibodies (NAb) of various antibody isotypes in the colostrum and calf serum. Colostrum traits had heritabilities ranging from 0.16 to 0.31 with repeatabilities ranging from 0.21 to 0.55. Brix percentages had positive genetic correlations with all colostrum antibody traits including total IgG (0.68). Calf serum antibody concentrations had heritabilities ranging from 0.25 to 0.59, with a significant maternal effect accounting for 17 to 27% of the variance. When later in life calves produced their first lactation, the lactation average somatic cell score was found to be negatively correlated with NAb levels in calf serum. CONCLUSIONS: Our results suggest that antibody levels in the colostrum and calf serum can be increased by means of selection.

Phenotypic and genetic effects of season on milk production traits in dairy cattle in the Netherlands.

Milk production systems in several countries show considerable differences between seasons. For example, in the Netherlands, cows are kept inside and fed silage in winter, whereas they are on pasture in summer. The differences between seasons affect milk yield and composition and might influence the genetic background of milk production traits. The objective of this study was to estimate phenotypic and genetic effects of season on milk production traits. For this purpose, 19,286 test-day milk production records of 1,800 first-parity Dutch Holstein-Frisian cows were available, and these cows were genotyped using a 50K SNP panel. Phenotypic effects of season were significant for all milk production traits. Effects of season were large for milk fat yield, fat content, and protein content. Genetic correlations between milk production traits in different seasons showed that genotype by season interaction effects were relatively small for most milk production traits. The genetic background of protein content and lactose content seems to be sensitive to seasonal effects. Furthermore, the genetic correlations between spring and autumn differed significantly from unity for almost all milk production traits. A genome-wide association study for genotype by season interaction identified chromosomal regions on BTA3, BTA14, BTA20, and BTA25 that showed genotype by season interaction effects, including a region containing DGAT1, which showed interaction effects for fat content and protein content.

Phenotypic and genetic effects of pregnancy on milk production traits in Holstein-Friesian cattle.

Pregnancy is a prerequisite for the initiation of lactation and for maintaining the milk production cycle. Pregnancy affects milk production and therefore should be accounted for in the genetic evaluation. Furthermore, there might be genetic differences in pregnancy effects on milk composition. The objective of this study was to estimate phenotypic and genetic effects of pregnancy on milk production traits. For this purpose, test-day records and conception dates of 1,359 first-parity Holstein-Friesian cows were analyzed. Significant effects of pregnancy on all milk production traits were detected except somatic cell score (e.g., the cumulative effects of pregnancy on milk yield were -247 kg). The pregnancy effects on milk yield, lactose yield, protein yield, fat yield, and fat content were small during early gestation (<150 d) and substantially increased in late gestation. The effects of pregnancy on milk protein yield were relatively stronger than those on fat yield. The effects of pregnancy on milk production traits differed for DGAT1 genotypes. Milk yield, lactose yield, protein yield, and fat yield of DGAT1 AA cows were more affected by pregnancy than that of DGAT1 KK cows (e.g., the cumulative effects of pregnancy on milk yield were negligible for DGAT1 KK cows and were -443 kg for DGAT1 AA cows). These results suggest that DGAT1 KK cows may be more suitable for shortening or omitting the dry period than DGAT1 AA cows.

Combined use of milk infrared spectra and genotypes can improve prediction of milk fat composition.

It has been shown that milk infrared (IR) spectroscopy can be used to predict detailed milk fat composition. In addition, polymorphisms with substantial effects on milk fat composition have been identified. In this study, we investigated the combined use of milk IR spectroscopy and genotypes of dairy cows on the accuracy of predicting milk fat composition. Milk fat composition data based on gas chromatography and milk IR spectra were available for 1,456 Dutch Holstein Friesian cows. In addition, genotypes for the diacylglycerol acyltransferase 1 (DGAT1) K232A and stearoyl-CoA desaturase 1 (SCD1) A293V polymorphisms and a SNP located in an intron of the fatty acid synthase (FASN) gene were available. Adding SCD1 genotypes to the milk IR spectra resulted in a considerable improvement of the prediction accuracy for the unsaturated fatty acids C10:1, C12:1, C14:1 cis-9, and C16:1 cis-9 and their corresponding unsaturation indices. Adding DGAT1 genotypes to the milk IR spectra resulted in an improvement of the prediction accuracy for C16:1 cis-9 and C16 index. Adding genotypes of the FASN SNP to the IR spectra did not improve prediction of milk fat composition. This study demonstrated the potential of combining milk IR spectra with genotypic information from 3 polymorphisms to predict milk fat composition. We hypothesize that prediction accuracy of milk fat composition can be further improved by combining milk IR spectra with genomic breeding values.

Genome-wide association study for genotype by lactation stage interaction of milk production traits in dairy cattle.

Substantial evidence demonstrates that the genetic background of milk production traits changes during lactation. However, most GWAS for milk production traits assume that genetic effects are constant during lactation and therefore might miss those quantitative trait loci (QTL) whose effects change during lactation. The GWAS for genotype by lactation stage interaction are aimed at explicitly detecting the QTL whose effects change during lactation. The purpose of this study was to perform GWAS for genotype by lactation stage interaction for milk yield, lactose yield, lactose content, fat yield, fat content, protein yield, and somatic cell score to detect QTL with changing effects during lactation. For this study, 19,286 test-day records of 1,800 first-parity Dutch Holstein cows were available and cows were genotyped using a 50K SNP panel. A total of 7 genomic regions with effects that change during lactation were detected in the GWAS for genotype by lactation stage interaction. Two regions on Bos taurus autosome (BTA)14 and BTA19 were also significant based on a GWAS that assumed constant genetic effects during lactation. Five regions on BTA4, BTA10, BTA11, BTA16, and BTA23 were only significant in the GWAS for genotype by lactation stage interaction. The biological mechanisms that cause these changes in genetic effects are still unknown, but negative energy balance and effects of pregnancy may play a role. These findings increase our understanding of the genetic background of lactation and may contribute to the development of better management indicators based on milk composition.

Reliability of genomic prediction for milk fatty acid composition by using a multi-population reference and incorporating GWAS results.

BACKGROUND: Large-scale phenotyping for detailed milk fatty acid (FA) composition is difficult due to expensive and time-consuming analytical techniques. Reliability of genomic prediction is often low for traits that are expensive/difficult to measure and for breeds with a small reference population size. An effective method to increase reference population size could be to combine datasets from different populations. Prediction models might also benefit from incorporation of information on the biological underpinnings of quantitative traits. Genome-wide association studies (GWAS) show that genomic regions on Bos taurus chromosomes (BTA) 14, 19 and 26 underlie substantial proportions of the genetic variation in milk FA traits. Genomic prediction models that incorporate such results could enable improved prediction accuracy in spite of limited reference population sizes. In this study, we combine gas chromatography quantified FA samples from the Chinese, Danish and Dutch Holstein populations and implement a genomic feature best linear unbiased prediction (GFBLUP) model that incorporates variants on BTA14, 19 and 26 as genomic features for which random genetic effects are estimated separately. Prediction reliabilities were compared to those estimated with traditional GBLUP models. RESULTS: Predictions using a multi-population reference and a traditional GBLUP model resulted in average gains in prediction reliability of 10% points in the Dutch, 8% points in the Danish and 1% point in the Chinese populations compared to predictions based on population-specific references. Compared to the traditional GBLUP, implementation of the GFBLUP model with a multi-population reference led to further increases in prediction reliability of up to 38% points in the Dutch, 23% points in the Danish and 13% points in the Chinese populations. Prediction reliabilities from the GFBLUP model were moderate to high across the FA traits analyzed. CONCLUSIONS: Our study shows that it is possible to predict genetic merits for milk FA traits with reasonable accuracy by combining related populations of a breed and using models that incorporate GWAS results. Our findings indicate that international collaborations that facilitate access to multi-population datasets could be highly beneficial to the implementation of genomic selection for detailed milk composition traits.

Genetic parameters for body weight and different definitions of residual feed intake in broiler chickens.

BACKGROUND: The objectives of this study were to (1) simultaneously estimate genetic parameters for BW, feed intake (FI), and body weight gain (Gain) during a FI test in broiler chickens using multi-trait Bayesian analysis; (2) derive phenotypic and genetic residual feed intake (RFI) and estimate genetic parameters of the resulting traits; and (3) compute a Bayesian measure of direct and correlated superiority of a group selected on phenotypic or genetic residual feed intake. A total of 56,649 male and female broiler chickens were measured at one of two ages ([Formula: see text] or [Formula: see text] days). BW, FI, and Gain of males and females at the two ages were considered as separate traits, resulting in a 12-trait model. Phenotypic RFI ([Formula: see text]) and genetic RFI ([Formula: see text]) were estimated from a conditional distribution of FI given BW and Gain using partial phenotypic and partial genetic regression coefficients, respectively. RESULTS: Posterior means of heritability for BW, FI and Gain were moderately high and estimates were significantly different between males and females at the same age for all traits. In addition, the genetic correlations between male and female traits at the same age were significantly different from 1, which suggests a sex-by-genotype interaction. Genetic correlations between [Formula: see text] and [Formula: see text] were significantly different from 1 at an older age but not at a younger age. CONCLUSIONS: The results of the multivariate Bayesian analyses in this study showed that genetic evaluation for production and feed efficiency traits should take sex and age differences into account to increase accuracy of selection and genetic gain. Moreover, for communicating with stakeholders, it is easier to explain results from selection on [Formula: see text] than selection on [Formula: see text], since [Formula: see text] is genetically independent of production traits and it explains the efficiency of birds in nutrient utilization independently of energy requirements for production and maintenance.

Validation strategy can result in an overoptimistic view of the ability of milk infrared spectra to predict methane emission of dairy cattle.

Because of the environmental impact of methane (CH4), it is of great interest to reduce CH4 emission of dairy cattle and selective breeding might contribute to this. However, this approach requires a rapid and inexpensive measurement technique that can be used to quantify CH4 emission for a large number of individual dairy cows. Milk infrared (IR) spectroscopy has been proposed as a predictor for CH4 emission. In this study, we investigated the feasibility of milk IR spectra to predict breath sensor-measured CH4 of 801 dairy cows on 10 commercial farms. To evaluate the prediction equation, we used random and block cross validation. Using random cross validation, we found a validation coefficient of determination (R2val) of 0.49, which suggests that milk IR spectra are informative in predicting CH4 emission. However, based on block cross validation, with farms as blocks, a negligible R2val of 0.01 was obtained, indicating that milk IR spectra cannot be used to predict CH4 emission. Random cross validation thus results in an overoptimistic view of the ability of milk IR spectra to predict CH4 emission of dairy cows. The difference between the validation strategies could be due to the confounding of farm and date of milk IR analysis, which introduces a correlation between batch effects on the IR analyses and farm-average CH4. Breath sensor-measured CH4 is strongly influenced by farm-specific conditions, which magnifies the problem. Milk IR wavenumbers from water absorption regions, which are generally considered uninformative, showed moderate accuracy (R2val = 0.25) when based on random cross validation, but not when based on block cross validation (R2val = 0.03). These results indicate, therefore, that in the current study, random cross validation results in an overoptimistic view on the ability of milk IR spectra to predict CH4 emission. We suggest prediction based on wavenumbers from water absorption regions as a negative control to identify potential dependence structures in the data.

Genome-wide association studies for genetic effects that change during lactation in dairy cattle.

Genetic effects on milk production traits in dairy cattle might change during lactation. However, most genome-wide association studies (GWAS) for milk production traits assume that genetic effects are constant during lactation. This assumption might lead to missing these quantitative trait loci (QTL) whose effects change during lactation. This study aimed to screen the whole genome specifically for QTL whose effects change during lactation. For this purpose, 4 different GWAS approaches were performed using test-day milk protein content records: (1) separate GWAS for specific lactation stages, (2) GWAS for estimated Wilmink lactation curve parameters, (3) a GWAS using a repeatability model where SNP effects are assumed constant during lactation, and (4) a GWAS for genotype by lactation stage interaction using a repeatability model and accounting for changing genetic effects during lactation. Separate GWAS for specific lactation stages suggested that the detection power greatly differs between lactation stages and that genetic effects of some QTL change during lactation. The GWAS for estimated Wilmink lactation curve parameters detected many chromosomal regions for Wilmink parameter a (protein content level), whereas 2 regions for Wilmink parameter b (decrease in protein content toward nadir) and no regions for Wilmink parameter c (increase in protein content after nadir) were detected. Twenty chromosomal regions were detected with effects on milk protein content; however, there was no evidence that their effects changed during lactation. For 5 chromosomal regions located on chromosomes 3, 9, 10, 14, and 27, significant evidence was observed for a genotype by lactation stage interaction and thus their effects on milk protein content changed during lactation. Three of these 5 regions were only identified using a GWAS for genotype by lactation stage interaction. Our study demonstrated that GWAS for genotype by lactation stage interaction offers new possibilities to identify QTL involved in milk protein content. The performed approaches can be applied to other milk production traits. Identification of QTL whose genetic effects change during lactation will help elucidate the genetic and biological background of milk production.

Genome-wide association study identifies loci influencing natural antibody titers in milk of Dutch Holstein-Friesian cattle.

Natural antibodies (NAb) are produced without any antigenic stimulation as a part of the innate immune system and provide a first line of defense against pathogens. Hence, they may be a useful trait when estimating an animal's potential immune competence and in selection for disease resistance. The aim of this study was to identify genomic regions associated with different NAb traits in milk and potentially describe candidate genes. Milk samples from 1,695 first-lactation Holstein Friesian cows with titer measurements for keyhole limpet hemocyanin, lipopolysaccharide, lipoteichoic acid, and peptidoglycan-binding total NAb and isotypes IgG1, IgM, and IgA were used. Genome-wide association study analyses were performed using imputed 777K SNP genotypes, accounting for relationships using pedigree information. Functional enrichment analysis was performed on the significantly associated genomic regions to look for candidate genes. For IgM NAb, significant associations (false discovery rate <0.05) were found on Bos taurus autosome (BTA) 17, 18, and 21 with candidate genes related to immunoglobulin structure and early B cell development. For IgG1, associations were found on BTA3, and we confirmed a quantitative trait loci on BTA21 previously reported for IgG NAb in serum. Our results provide new insights into the regulation of milk NAb that will help unravel the complex relationship between milk immunoglobulins and disease resistance in dairy cattle.

The DGAT1 K232A polymorphism and feeding modify milk fat triacylglycerol composition.

In the present study, we aimed to investigate the changes in triacylglycerol (TAG) composition as affected by alterations in the cows' diet due to seasonal variations and genetic factors. For this study, 50 milk fat samples in winter and 50 in summer were used from 25 cows with the DGAT1 KK genotype and 25 cows with the DGAT1 AA genotype. The samples were analyzed for milk fat content (%), fat composition, and TAG composition. We found that the content of TAG species CN54 was higher and that of CN34 and CN36 lower in summer than in winter. This seasonal variation in TAG profile was related to seasonal changes in the fatty acids C14:0, C16:0, C18:0, C18:1 cis-9, total unsaturated fatty acids, and total long-chain fatty acids, most likely resulting from dietary differences between seasons. Furthermore, we quantified the effect of DGAT1 K232A polymorphism on TAG profile and detected a significant effect on TAG species CN36, with higher values for the DGAT1 KK genotype. When adjusting for differences in fat content, we found no significant effects of the DGAT1 K232A polymorphism on TAG profile. We detected a significant interaction between DGAT1 K232A polymorphism and season for TAG species CN42 and CN52; in summer, the KK genotype was associated with higher levels for CN42 than the AA genotype, whereas in winter, the difference between the genotypes was small. For CN52, in summer the AA genotype was associated with higher levels than the KK genotype. In winter, the difference between the genotypes was also small. We show that, regardless of preference for DGAT1 genotype (AA or KK) and depending on the availability of FA according to season, UFA (C18:1 cis-9), short-chain FA (C6:0 and C10:0), and medium-chain FA might be esterified on the glycerol backbone of the TAG, keeping the structure characteristics of each TAG species. To our knowledge, this is the first report on the interaction effect of DGAT1 K232A polymorphism and season on the TAG composition in milk fat.

The impact of genome editing on the introduction of monogenic traits in livestock.

BACKGROUND: Genome editing technologies provide new tools for genetic improvement and have the potential to become the next game changer in animal and plant breeding. The aim of this study was to investigate how genome editing in combination with genomic selection can accelerate the introduction of a monogenic trait in a livestock population as compared to genomic selection alone. METHODS: A breeding population was simulated under genomic selection for a polygenic trait. After reaching Bulmer equilibrium, the selection objective was to increase the allele frequency of a monogenic trait, with or without genome editing, in addition to improving the polygenic trait. Scenarios were compared for time to fixation of the desired allele, selection response for the polygenic trait, and level of inbreeding. The costs, in terms of number of editing procedures, were compared to the benefits of having more animals with the desired phenotype of the monogenic trait. Effects of reduced editing efficiency were investigated. RESULTS: In a population of 20,000 selection candidates per generation, the total number of edited zygotes needed to reach fixation of the desired allele was 22,118, 7072, or 3912 with, no, moderate, or high selection emphasis on the monogenic trait, respectively. Genome editing resulted in up to four-fold faster fixation of the desired allele when efficiency was 100%, while the loss in long-term selection response for the polygenic trait was up to seven-fold less compared to genomic selection alone. With moderate selection emphasis on the monogenic trait, introduction of genome editing led to a four-fold reduction in the total number of animals showing the undesired phenotype before fixation. However, with a currently realistic editing efficiency of 4%, the number of required editing procedures increased by 72% and loss in selection response increased eight-fold compared to 100% efficiency. With low efficiency, loss in selection response was 29% more compared to genomic selection alone. CONCLUSIONS: Genome editing strongly decreased the time to fixation for a desired allele compared to genomic selection alone. Reduced editing efficiency had a major impact on the number of editing procedures and on the loss in selection response. In addition to ethical and welfare considerations of genome editing, a careful assessment of its technical costs and benefits is required.

Genome-wide association study for milk infrared wavenumbers.

Individual wavenumbers of the infrared (IR) spectra of bovine milk have been shown to be moderately to highly heritable. The objective of this study was to identify genomic regions associated with individual milk IR wavenumbers. This is expected to provide information about the genetic background of milk composition and give insight in the relation between IR wavenumbers and milk components. For this purpose, a genome-wide association study was performed for a selected set of 50 individual IR wavenumbers measured on 1,748 Dutch Holstein cows. Significant associations were detected for 28 of the 50 wavenumbers. In total, 24 genomic regions distributed over 16 bovine chromosomes were identified. Major genomic regions associated with milk IR wavenumbers were identified on chromosomes 1, 5, 6, 14, 19, and 20. Most of these regions also showed significant associations with fat, protein, or lactose percentage. However, we also identified some new regions that were not associated with any one of these routinely collected milk composition traits. On chromosome 1, we identified 2 new genomic regions and hypothesized that they are related to variation in milk phosphorus content and orotic acid, respectively. On chromosome 20, we identified a new genomic region that seems to be related to citric acid. Identification of genomic regions associated with milk phosphorus content, orotic acid, and citric acid suggest that the milk IR spectra contain direct information on these milk components. Consequently milk IR analyses probably can be used to predict these milk components, which have low concentrations in milk; this can lead to novel applications of milk IR spectroscopy for dairy cattle breeding and herd management.

Characterizing a region on BTA11 affecting ß-lactoglobulin content of milk using high-density genotyping and haplotype grouping.

BACKGROUND: Milk ß-lactoglobulin (ß-LG) content is of interest as it is associated with nutritional and manufacturing properties. It is known that milk ß-LG content is strongly affected by genetic factors. In cattle, most of the genetic differences are associated with a chromosomal region on BTA11, which contains the ß-LG gene. The aim of this study was to characterize this region using 777 k SNP data (BovineHDbeadChip) and perform a haplotype-based association study. A statistical approach was developed to build haplotypes that capture the genetic variation associated with this genomic region. RESULTS: The SNP with the most significant effect on ß-lactoglobulin content was one of the 2 causal mutations responsible for the ß-lactoglobulin protein variants A/B. Haplotypes based on 2 to 5 selected lead SNP were clustered in groups with different effects on ß-lactoglobulin content. Four different groups were identified suggesting that ß-lactoglobulin variant A and B can be further refined in A1, A2, B1 and B2. CONCLUSIONS: This study showed that ß-lactoglobulin protein variants A/B do not explain all genetic variation associated with the tail part of BTA11 but this region contains more than one mutation with an effect on ß-lactoglobulin content. These findings can be used for selection of cows with higher cheese yield, which is desirable for the dairy industry.

Modeling heterogeneous (co)variances from adjacent-SNP groups improves genomic prediction for milk protein composition traits.

BACKGROUND: Accurate genomic prediction requires a large reference population, which is problematic for traits that are expensive to measure. Traits related to milk protein composition are not routinely recorded due to costly procedures and are considered to be controlled by a few quantitative trait loci of large effect. The amount of variation explained may vary between regions leading to heterogeneous (co)variance patterns across the genome. Genomic prediction models that can efficiently take such heterogeneity of (co)variances into account can result in improved prediction reliability. In this study, we developed and implemented novel univariate and bivariate Bayesian prediction models, based on estimates of heterogeneous (co)variances for genome segments (BayesAS). Available data consisted of milk protein composition traits measured on cows and de-regressed proofs of total protein yield derived for bulls. Single-nucleotide polymorphisms (SNPs), from 50K SNP arrays, were grouped into non-overlapping genome segments. A segment was defined as one SNP, or a group of 50, 100, or 200 adjacent SNPs, or one chromosome, or the whole genome. Traditional univariate and bivariate genomic best linear unbiased prediction (GBLUP) models were also run for comparison. Reliabilities were calculated through a resampling strategy and using deterministic formula. RESULTS: BayesAS models improved prediction reliability for most of the traits compared to GBLUP models and this gain depended on segment size and genetic architecture of the traits. The gain in prediction reliability was especially marked for the protein composition traits ß-CN, κ-CN and ß-LG, for which prediction reliabilities were improved by 49 percentage points on average using the MT-BayesAS model with a 100-SNP segment size compared to the bivariate GBLUP. Prediction reliabilities were highest with the BayesAS model that uses a 100-SNP segment size. The bivariate versions of our BayesAS models resulted in extra gains of up to 6% in prediction reliability compared to the univariate versions. CONCLUSIONS: Substantial improvement in prediction reliability was possible for most of the traits related to milk protein composition using our novel BayesAS models. Grouping adjacent SNPs into segments provided enhanced information to estimate parameters and allowing the segments to have different (co)variances helped disentangle heterogeneous (co)variances across the genome.

Genomic selection for crossbred performance accounting for breed-specific effects.

BACKGROUND: Breed-specific effects are observed when the same allele of a given genetic marker has a different effect depending on its breed origin, which results in different allele substitution effects across breeds. In such a case, single-breed breeding values may not be the most accurate predictors of crossbred performance. Our aim was to estimate the contribution of alleles from each parental breed to the genetic variance of traits that are measured in crossbred offspring, and to compare the prediction accuracies of estimated direct genomic values (DGV) from a traditional genomic selection model (GS) that are trained on purebred or crossbred data, with accuracies of DGV from a model that accounts for breed-specific effects (BS), trained on purebred or crossbred data. The final dataset was composed of 924 Large White, 924 Landrace and 924 two-way cross (F1) genotyped and phenotyped animals. The traits evaluated were litter size (LS) and gestation length (GL) in pigs. RESULTS: The genetic correlation between purebred and crossbred performance was higher than 0.88 for both LS and GL. For both traits, the additive genetic variance was larger for alleles inherited from the Large White breed compared to alleles inherited from the Landrace breed (0.74 and 0.56 for LS, and 0.42 and 0.40 for GL, respectively). The highest prediction accuracies of crossbred performance were obtained when training was done on crossbred data. For LS, prediction accuracies were the same for GS and BS DGV (0.23), while for GL, prediction accuracy for BS DGV was similar to the accuracy of GS DGV (0.53 and 0.52, respectively). CONCLUSIONS: In this study, training on crossbred data resulted in higher prediction accuracy than training on purebred data and evidence of breed-specific effects for LS and GL was demonstrated. However, when training was done on crossbred data, both GS and BS models resulted in similar prediction accuracies. In future studies, traits with a lower genetic correlation between purebred and crossbred performance should be included to further assess the value of the BS model in genomic predictions.