Prediction of likelihood of conception in dairy cows using milk mid-infrared spectra collected before the first insemination and machine learning algorithms.

Accurate and ex-ante prediction of cows' likelihood of conception (LC) based on milk composition information could improve reproduction management on dairy farms. Milk composition is already routinely measured by mid-infrared (MIR) spectra, which are known to change with advancing stages of pregnancy. For lactating cows, MIR spectra may also be used for predicting the LC. Our objectives were to classify the LC at first insemination using milk MIR spectra data collected from calving to first insemination and to identify the spectral regions that contribute the most to the prediction of LC at first insemination. After quality control, 4,866 MIR spectra, milk production, and reproduction records from 3,451 Holstein cows were used. The classification accuracy and area under the curve (AUC) of 6 models comprising different predictors and 3 machine learning methods were estimated and compared. The results showed that partial least square discriminant analysis (PLS-DA) and random forest had higher prediction accuracies than logistic regression. The classification accuracy of good and poor LC cows and AUC in herd-by-herd validation of the best model were 76.35 ± 10.60% and 0.77 ± 0.11, respectively. All wavenumbers with values of variable importance in the projection higher than 1.00 in PLS-DA belonged to 3 spectral regions, namely from 1,003 to 1,189, 1,794 to 2,260, and 2,300 to 2,660 cm-1. In conclusion, the model can predict LC in dairy cows from a high productive TMR system before insemination with a relatively good accuracy, allowing farmers to intervene in advance or adjust the insemination schedule for cows with a poor predicted LC.

Identification and characterisation of de novo germline structural variants in two commercial pig lines using trio-based whole genome sequencing.

BACKGROUND: De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints. RESULTS: We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data. CONCLUSIONS: Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs.

Observational study on associations between resilience indicators based on daily milk yield in first lactation and lifetime profitability.

Resilience is the ability of cows to be minimally affected by disturbances, such as pathogens, heat waves, and changes in feed quality, or to quickly recover. Obvious advantages of resilience are good animal welfare and easy and pleasant management for farmers. Furthermore, economic effects are also expected, but these remain to be determined. The goal of this study was to investigate the association between resilience and lifetime gross margin, using indicators of resilience calculated from fluctuations in daily milk yield using an observational study. Resilience indicators and lifetime gross margin were calculated for 1,325 cows from 21 herds. These cows were not alive anymore and, therefore, had complete lifetime data available for many traits. The resilience indicators were the natural log-transformed variance (LnVar) and the lag-1 autocorrelation (rauto) of daily milk yield deviations from cow-specific lactation curves in parity 1. Good resilience is indicated by low LnVar (small yield response to disturbances) and low rauto (quick yield recovery to baseline). Lifetime gross margin was calculated as the sum of all revenues minus the sum of all costs throughout life. Included revenues were from milk, calf value, and slaughter of the cow. Included costs were from feed, rearing, insemination, management around calving, disease treatments, and destruction in case of death on farm. Feed intake was unknown and, therefore, lifetime feed costs had to be estimated based on milk yield records. The association of each resilience indicator with lifetime gross margin, and also with the underlying revenues and costs, was investigated using analysis of covariance (ANCOVA) models. Mean daily milk yield in first lactation, herd, and year of birth were included as covariates and factors. Natural log-transformed variance had a significantly negative association with lifetime gross margin, which means that cows with stable milk yield (low LnVar, good resilience) in parity 1 generated on average a higher lifetime gross margin than cows that had the same milk yield level but with more fluctuations. The association with lifetime gross margin could be mainly attributed to higher lifetime milk revenues for cows with low LnVar, due to a longer lifespan. Unlike LnVar, rauto was not significantly associated with lifetime gross margin or any of the underlying lifetime costs and revenues. However, it was significantly associated with yearly treatment costs, which is important for ease of management. In conclusion, the importance of resilience for total profit generated by a cow at the end of life was confirmed by the significant association of LnVar with lifetime gross margin, although effects of differences in feed efficiency between resilient and less resilient cows remain to be studied. The economic advantage can be mainly ascribed to benefits of long lifespan.

Validation of resilience indicators by estimating genetic correlations among daughter groups and with yield responses to a heat wave and disturbances at herd level.

Resilient cows are minimally affected in their functioning by disturbances, and if affected, they quickly recover. Previously, the variance and autocorrelation of daily deviations from a lactation curve were proposed as resilience indicators. These traits were heritable and genetically associated with good health and longevity. However, it was unknown if selection for these indicators would lead to desired changes in the phenotype. The first aim of this study was to investigate if forward prediction of the resilience indicators in another environment was possible. Therefore, the resilience indicator records were split into 2 subsets, each containing half of the daughters of each sire, split within sire into cows that calved in early year-seasons and cows that calved in more recent year-seasons. Genetic correlations between the subsets were then estimated for each resilience indicator. The second aim was to estimate genetic correlations between the resilience indicators and traits describing production responses to actual disturbances. The disturbances were a heat wave in July 2015 and yield disturbances at herd level. The latter were selected by decreases in mean yield of all primiparous cows in a herd, indicating that a disturbance occurred. The data set used for calculation of the resilience indicators and the traits describing yield responses contained 62,932,794 daily milk yield records on 199,104 primiparous cows. Genetic correlations (rg) between recent and earlier daughter groups were 1 for both resilience indicators, which suggests that selection will result in changes in the phenotype in the next generation. Furthermore, low variance was genetically correlated with weak response in milk yield to both the heat wave and herd disturbances (rg 0.47 to 0.97). Low autocorrelation was genetically correlated with reduced perturbation length and quick recovery after the heat wave and herd disturbances (0.28 to 0.97). These results suggest that variance and autocorrelation cover different aspects of resilience, and should be combined in a resilience index. In conclusion, genetic selection for the resilience indicators will likely result in favorable changes in the traits themselves, and in response and recovery to actual disturbances, which confirms that they are useful resilience indicators.

Between-herd variation in resilience and relations to herd performance.

Resilient cows are minimally affected in their functioning by infections and other disturbances, and recover quickly. Herd management is expected to have an effect on disturbances and the resilience of cows, and this effect was investigated in this study. Two resilience indicators were first recorded on individual cows. The effect of herd-year on these resilience indicators was then estimated and corrected for genetic and year-season effects. The 2 resilience indicators were the variance and the lag-1 autocorrelation of daily milk yield deviations from an expected lactation curve. Low variance and autocorrelation indicate that a cow does not fluctuate much around her expected milk yield and is, thus, subject to few disturbances, or little affected by disturbances (resilient). The herd-year estimates of the resilience indicators were estimated for 9,917 herd-year classes based on records of 227,655 primiparous cows from 2,644 herds. The herd-year estimates of the resilience indicators were then related to herd performance variables. Large differences in the herd-year estimates of the 2 resilience indicators (variance and autocorrelation) were observed between herd-years, indicating an effect of management on these traits. Furthermore, herd-year classes with a high variance tended to have a high proportion of cows with a rumen acidosis indication (r = 0.31), high SCS (r = 0.19), low fat content (r = -0.18), long calving interval (r = 0.14), low survival to second lactation (r = -0.13), large herd size (r = 0.12), low lactose content (r = -0.12), and high production (r = 0.10). These correlations support that herds with high variance are not resilient. The correlation between the variance and the proportion of cows with a rumen acidosis indication suggests that feed management may have an important effect on the variance. Herd-year classes with a high autocorrelation tended to have a high proportion of cows with a ketosis indication (r = 0.14) and a high production (r = 0.13), but a low somatic cell score (r = -0.17) and a low proportion of cows with a rumen acidosis indication (r = -0.12). These correlations suggest that high autocorrelation at herd level indicates either good or poor resilience, and is thus a poor resilience indicator. However, the combination of a high variance and a high autocorrelation is expected to indicate many fluctuations with slow recovery. In conclusion, herd management, in particular feed management, seems to affect herd resilience.

Genetic analysis of resilience indicators based on milk yield records in different lactations and at different lactation stages.

Resilience is the ability of cows to cope with disturbances, such as pathogens or heat waves. To breed for improved resilience, it is important to know whether resilience genetically changes throughout life. Therefore, the aim was to perform a genetic analysis on 2 resilience indicators based on data from 3 periods of the first lactation (d 11-110, 111-210, and 211-340) and the first 3 full lactations, and to estimate genetic correlations with health traits. The resilience indicators were the natural log-transformed variance (LnVar) and lag-1 autocorrelation (rauto) of daily deviations in milk yield from an expected lactation curve. Low LnVar and rauto indicate low variability in daily milk yield and quick recovery, and were expected to indicate good resilience. Data of 200,084 first, 155,784 second, and 89,990 third lactations were used. Heritabilities were similar based on different lactation periods (0.12-0.15 for LnVar, 0.05-0.06 for rauto). However, the heritabilities of the resilience indicators based on full first lactation were higher than those based on lactation periods (0.20 for LnVar, 0.08 for rauto), due to lower residual variances. Heritabilities decreased from 0.20 in full lactation 1 to 0.19 in full lactation 3 for LnVar and from 0.08 to 0.06 for rauto. For LnVar, as well as for rauto, the strongest genetic correlation between lactation periods was between period 2 and 3 (0.97 for LnVar, 0.96 for rauto) and the weakest between period 1 and 3 (0.81 for LnVar, 0.65 for rauto). Similarly, for both traits the genetic correlation between full lactations was strongest between lactations 2 and 3 (0.99 for LnVar, 0.95 for rauto) and weakest between lactations 1 and 3 (0.91 for LnVar, 0.71 for rauto). For LnVar, genetic correlations with resilience-related traits, such as udder health, ketosis, and longevity, adjusted for correlations with milk yield, were almost always favorable (-0.59 to 0.02). In most cases these genetic correlations were stronger based on full lactations than on lactation periods. Genetic correlations were similar across full lactations, but the correlation with udder health increased substantially from -0.31 in lactation 1 to -0.51 in lactation 3. For rauto, genetic correlations with resilience-related traits were always favorable in lactation period 1 and in most full lactations, but not in the other lactation periods. However, correlations were weak (-0.27 to 0.15). Therefore, as a resilience indicator for breeding, LnVar is preferred over rauto. A multitrait index based on estimated breeding values for LnVar in lactations 1, 2, and 3 is recommended to improve resilience throughout the lifetime of a cow.

Competitive gene flow does not necessarily maximize the genetic gain of genomic breeding programs in the presence of genotype-by-environment interaction.

National and international across-population selection is often recommended and fairly common in the current breeding practice of dairy cattle, with the primary aims to increase genetic gain and genetic variability. The aim of this study was to test the hypothesis that the strategy of truncation selection of sires across populations [i.e., competitive gene flow strategy (CGF)] may not necessarily maximize genetic gain in the long term in the presence of genotype-by-environment interaction (G×E). Two alternative strategies used to be compared with CGF were forced gene flow (FGF) strategies, with 10 or 50% of domestic dams forced to be mated with foreign sires (FGF10%, FGF50%). Two equal-size populations (Ndams = 1,000) that were selected for the same breeding goal trait (h2 = 0.3) under G×E correlation (rg) of either 0.9 or 0.8 were simulated to test these 3 different strategies. Each population first experienced either 5 or 20 differentiation generations (Gd), then 15 migration generations. Discrete generations were simulated for simplicity. Each population performed a within-population conventional breeding program during differentiation generations and the 3 across-population sire selection strategies based on joint genomic prediction during migration generations. The 4 Gd_rg combinations defined 4 different levels of differentiation degree between the 2 populations at the start of migration. The true rate of inbreeding over the last 10 migration generations in each scenario was constrained at 0.01 to provide a fair basis for comparison of genetic gain across scenarios. Results showed that CGF maximized the genetic gain after 15 migration generations in 5_0.9 combination only, the case of the lowest differentiation degree, with a superiority of 0.4% (0.04 genetic SD units) over the suboptimal strategy. While in 5_0.8, 20_0.9, and 20_0.8 combinations, 2 FGF strategies had a superiority in genetic gain of 2.3 to 12.5% (0.21-1.07 genetic SD units) over CGF after 15 migration generations, especially FGF50%. The superiority of FGF strategies over CGF was that they alleviated inbreeding, introduced new genetic variance in the early migration period, and improved accuracy in the entire migration period. Therefore, we concluded that CGF does not necessarily maximize the genetic gain of across-population genomic breeding programs given moderate G×E. The across-population selection strategy remains to be optimized to maximize genetic gain.

Exploration of variance, autocorrelation, and skewness of deviations from lactation curves as resilience indicators for breeding.

The ability of a cow to cope with environmental disturbances, such as pathogens or heat waves, is called resilience. To improve resilience through breeding, we need resilience indicators, which could be based on the fluctuation patterns in milk yield resulting from disturbances. The aim of this study was to explore 3 traits that describe fluctuations in milk yield as indicators for breeding resilient cows: the variance, autocorrelation, and skewness of the deviations from individual lactation curves. We used daily milk yield records of 198,754 first-parity cows, recorded by automatic milking systems. First, we estimated a lactation curve for each cow using 4 different methods: moving average, moving median, quantile regression, and Wilmink curve. We then calculated the log-transformed variance (LnVar), lag-1 autocorrelation (rauto), and skewness (Skew) of the daily deviations from these curves as resilience indicators. A genetic analysis of the resilience indicators was performed, and genetic correlations between resilience indicators and health, longevity, fertility, metabolic, and production traits were estimated. The heritabilities differed between LnVar (0.20 to 0.24), rauto (0.08 to 0.10), and Skew (0.01 to 0.02), and the genetic correlations among the indicators were weak to moderate. For rauto and Skew, genetic correlations with health, longevity, fertility, and metabolic traits were weak or the opposite of what we expected. Therefore, rauto and Skew have limited value as resilience indicators. However, lower LnVar was genetically associated with better udder health (genetic correlations from -0.22 to -0.32), better longevity (-0.28 to -0.34), less ketosis (-0.27 to -0.33), better fertility (-0.06 to -0.17), higher BCS (-0.29 to -0.40), and greater dry matter intake (-0.53 to -0.66) at the same level of milk yield. These correlations support LnVar as an indicator of resilience. Of all 4 curve-fitting methods, LnVar based on quantile regression systematically had the strongest genetic correlations with health, longevity, and fertility traits. Thus, quantile regression is considered the best curve-fitting method. In conclusion, LnVar based on deviations from a quantile regression curve is a promising resilience indicator that can be used to breed cows that are better at coping with disturbances.

Possibilities for a specific breeding program for organic dairy production.

Organic dairy production differs from conventional dairy production in many aspects. However, breeding programs for the 2 production systems are the same in most countries. Breeding goals (BG) might be different for the 2 production systems and genotype × environment interaction may exist between organic and conventional dairy production, both of which have an effect on genetic gain in different breeding strategies. Other aspects also need to be considered, such as the application of multiple ovulation and embryo transfer (MOET), which is not allowed in organic dairy production. The general aim of this research was to assess different environment-specific breeding strategies for organic dairy production. The specific aim was to study differences in BG weights and include the effect of genotype × environment interaction, MOET, and the selection of breeding bulls from the conventional environment. Different scenarios were simulated. In the current scenario, the present-day situation for dairy production in Denmark was emulated as much as possible. The BG was based on a conventional dairy production system, MOET was applied in both environments, and conventional bulls could be selected as breeding bulls in the organic environment. Four alternative scenarios were simulated, all with a specific organic BG in the organic breeding program but differences in the usage of MOET and the selection of conventional bulls as breeding bulls. Implementation of a specific BG in organic dairy production slightly increased genetic gain in the aggregate genotype compared with the breeding program that is currently implemented in organic dairy production. Not using embryo transfer or only selecting breeding bulls from the organic environment decreased genetic gain in the aggregate genotype by as much as 24%. However, the use of embryo transfer is debatable because this is not allowed according to current regulations for organic dairy production. Assessing genetic gain on trait levels showed that a significant increase for functional traits was possible compared with the current breeding program in the organic environment without a decrease in genetic gain in the aggregate genotype. This difference on trait level was even more present when selection of conventional bulls as breeding bulls in the organic environment was not possible. This finding is very relevant when breeding for the desired cow in organic dairy production.

Genetic analysis of udder conformation traits derived from automatic milking system recording in dairy cows.

Cartesian teat coordinates measured by automatic milking systems (AMS) provide new opportunities to record udder conformation traits and to study changes in udder conformation genetically and phenotypically within and between parities. The objective of this study was to estimate heritabilities and repeatabilities of AMS-based udder conformation traits within parities, to estimate genetic correlations between parities for AMS-based udder conformation traits, and to estimate genetic correlations between AMS-based udder conformation traits and classifier-based udder conformation traits, longevity, and udder health. Data from 70 herds, including 12,663 first-parity cows, 10,206 second-parity cows, and 7,627 third-parity cows, were analyzed using univariate and bivariate mixed animal models. Heritabilities of the AMS udder conformation traits were large (0.37-0.67) and genetic correlations between the AMS udder conformation traits and classifier-based traits were strong (>0.91). Repeatabilities within parities were large as well (0.89-0.97), indicating that a single record on udder conformation per lactation reflects udder conformation well. Genetic correlations of AMS udder conformation traits between parities were strong (0.88-1.00) and were stronger than the permanent environmental correlations. This shows that udder conformation changes over parities, but this change is mostly due to nongenetic factors. Based on these results, the current herd classification system, where cows are scored on udder conformation once in first parity, is sufficient. The AMS udder conformation traits as defined in this study have limited value as replacement for classifier-based udder conformation traits because they have smaller genetic correlations with functional traits than classifier-based traits. In summary, udder conformation hardly changes genetically between parities and is highly repeatable within parities. Udder conformation traits based on AMS need fine-tuning before they can replace classifier-based traits, and AMS teat coordinates probably contain additional information about udder health that is yet to be explored.

Genomic selection improves the possibility of applying multiple breeding programs in different environments.

One joint breeding program (BP) for different dairy cattle environments can be advantageous for genetic gain depending on the genetic correlation between environments (rg). The break-even correlation (rb) refers to the specific rg where genetic gain with 1 joint BP is equal to the genetic gain of 2 environment-specific BP. One joint BP has the highest genetic gain if rg is higher than rb, whereas 2 environment-specific BP have higher genetic gain if rg is lower than rb. Genetic gain in this context is evaluated from a breeding company's perspective that aims to improve genetic gain in both environments. With the implementation of genomic selection, 2 types of collaboration can be identified: exchanging breeding animals and exchanging genomic information. The aim of this study was to study genetic gain in multiple environments with different breeding strategies with genomic selection. The specific aims were (1) to find rb when applying genomic selection; (2) to assess how much genetic gain is lost when applying a suboptimal breeding strategy; (3) to study the effect of the reliability of direct genomic values, number of genotyped animals, and environments of different size on rb and genetic gain; and (4) to find rb from each environment's point of view. Three breeding strategies were simulated: 1 joint BP for both environments, 2 environment-specific BP with selection of bulls across environments, and 2 environment-specific BP with selection of bulls within environments. The rb was 0.65 and not different from rb with progeny-testing breeding programs when compared at the same selection intensity. The maximum loss in genetic gain in a suboptimal breeding strategy was 24%. A higher direct genomic value reliability and an increased number of genotyped selection candidates increased genetic gain, and the effect on rb was not large. A different size in 2 environments decreased rb by, at most, 0.10 points. From a large environment's point of view, 1 joint BP was the optimal breeding strategy in most scenarios. From a small environment's point of view, 1 joint BP was only the optimal breeding strategy at high rg. When the exchange of breeding animals between environments was restricted, genetic gain could still increase in each environment. This was due to the exchange of genomic information between environments, even when rg between environments were as low as 0.4. Thus, genomic selection improves the possibility of applying environment-specific BP.

Fluctuations in milk yield are heritable and can be used as a resilience indicator to breed healthy cows.

Automatic milking systems record an enormous amount of data on milk yield and the cow itself. These type of big data are expected to contain indicators for health and resilience of cows. In this study, the aim was to define and estimate heritabilities for traits related with fluctuations in daily milk yield and to estimate genetic correlations with existing functional traits, such as udder health, fertility, claw health, ketosis, and longevity. We used daily milk yield records from automatic milking systems of 67,025 lactations in the first parity from 498 herds in the Netherlands. We defined 3 traits related to the number of drops in milk yield using Student t-tests based on either a rolling average (drop rolling average) or a regression (drop regression) and the natural logarithm of the within-cow variance of milk yield (LnVar). Average milk yield was added to investigate the relationships between milk yield and these new traits. ASReml was used to estimate heritabilities, breeding values (EBV), and genetic correlations among these new traits and average milk yield. Approximate genetic correlations were calculated using correlations between EBV of the new traits and existing EBV for health and functional traits correcting for nonunity reliabilities using the Calo method. Partial genetic correlations controlling for persistency and average milk yield and relative contributions to reliability were calculated to investigate whether the new traits add new information to predict fertility, health, and longevity. Heritabilities were 0.08 for drop rolling average, 0.06 for drop regression, and 0.10 for LnVar. Approximate genetic correlations between the new traits and the existing health traits differed quite a bit, with the strongest correlations (-0.29 to -0.52) between LnVar and udder health, ketosis, persistency, and longevity. This study shows that fluctuations in daily milk yield are heritable and that the variance of milk production is best among the 3 fluctuations traits tested to predict udder health, ketosis, and longevity. Using the residual variance of milk production instead of the raw variance is expected to further improve the trait to breed healthy, resilient, and long-lasting dairy cows.

Estimating variance components and breeding values for number of oocytes and number of embryos in dairy cattle using a single-step genomic evaluation.

Reproductive technologies such as multiple ovulation and embryo transfer (MOET) and ovum pick-up (OPU) accelerate genetic improvement in dairy breeding schemes. To enhance the efficiency of embryo production, breeding values for traits such as number of oocytes (NoO) and number of MOET embryos (NoM) can help in selection of donors with high MOET or OPU efficiency. The aim of this study was therefore to estimate variance components and (genomic) breeding values for NoO and NoM based on Dutch Holstein data. Furthermore, a 10-fold cross-validation was carried out to assess the accuracy of pedigree and genomic breeding values for NoO and NoM. For NoO, 40,734 OPU sessions between 1993 and 2015 were analyzed. These OPU sessions originated from 2,543 donors, from which 1,144 were genotyped. For NoM, 35,695 sessions between 1994 and 2015 were analyzed. These MOET sessions originated from 13,868 donors, from which 3,716 were genotyped. Analyses were done using only pedigree information and using a single-step genomic BLUP (ssGBLUP) approach combining genomic information and pedigree information. Heritabilities were very similar based on pedigree information or based on ssGBLUP [i.e., 0.32 (standard error = 0.03) for NoO and 0.21 (standard error = 0.01) for NoM with pedigree, 0.31 (standard error = 0.03) for NoO, and 0.22 (standard error = 0.01) for NoM with ssGBLUP]. For animals without their own information as mimicked in the cross-validation, the accuracy of pedigree-based breeding values was 0.46 for NoO and NoM. The accuracies of genomic breeding values from ssGBLUP were 0.54 for NoO and 0.52 for NoM. These results show that including genomic information increases the accuracies. These moderate accuracies in combination with a large genetic variance show good opportunities for selection of potential bull dams.

Genotype by environment interaction for production, somatic cell score, workability, and conformation traits in Dutch Holstein-Friesian cows between farms with or without grazing.

In this study, genotype by environment interaction was investigated for production traits, somatic cell score (SCS), workability traits, and conformation traits for Holstein-Friesian cows producing on farms with or without grazing in the Netherlands. Additionally, heritabilities and repeatabilities were estimated in both farm systems. Data were available for 1,019 Dutch farms, and farm type was known for those farms, 142 farms without grazing and 877 farms with grazing. The data set consisted of 428,600 test-day records for production from 49,412 cows, and from this data set a subset for SCS was created, consisting of 374,734 test-day records from 45,955 cows. For workability and conformation traits, the data set consisted of 30,180 cows. Bivariate mixed models with multiple fixed effects and random sire and random permanent environment effects were applied. The majority of sires had daughters in both farm types. The heritabilities for milk yield (0.27), fat yield (0.19), and protein yield (0.20) were higher in farms with grazing than in farms without grazing with heritabilities of 0.24 for milk yield, 0.18 for fat yield, and 0.18 for protein yield. Repeatability was lower in the grazing farms for milk yield, fat yield, and protein yield, probably because of alternating quality of dry matter intake during grazing. Genetic correlations between grazing and no grazing were 0.99, 0.98, 0.97, and 1.00 for milk yield, fat yield, protein yield, and SCS, respectively. Genetic correlations for workability traits and conformation traits between grazing and no grazing varied between 0.93 and 1.00. For all traits, genetic correlations were close to unity, indicating no genotype by environment interaction between farms with or without grazing for production traits, SCS, workability traits, and conformation traits in Dutch Holstein-Friesians. Therefore, the same sires can be used for farms both with and without grazing.

Genome-wide association study reveals novel loci for litter size and its variability in a Large White pig population.

BACKGROUND: In many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method. RESULTS: In total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB. CONCLUSIONS: To the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci.

Genetic variance in micro-environmental sensitivity for milk and milk quality in Walloon Holstein cattle.

Animals that are robust to environmental changes are desirable in the current dairy industry. Genetic differences in micro-environmental sensitivity can be studied through heterogeneity of residual variance between animals. However, residual variance between animals is usually assumed to be homogeneous in traditional genetic evaluations. The aim of this study was to investigate genetic heterogeneity of residual variance by estimating variance components in residual variance for milk yield, somatic cell score, contents in milk (g/dL) of 2 groups of milk fatty acids (i.e., saturated and unsaturated fatty acids), and the content in milk of one individual fatty acid (i.e., oleic acid, C18:1 cis-9), for first-parity Holstein cows in the Walloon Region of Belgium. A total of 146,027 test-day records from 26,887 cows in 747 herds were available. All cows had at least 3 records and a known sire. These sires had at least 10 cows with records and each herd × test-day had at least 5 cows. The 5 traits were analyzed separately based on fixed lactation curve and random regression test-day models for the mean. Estimation of variance components was performed by running iteratively expectation maximization-REML algorithm by the implementation of double hierarchical generalized linear models. Based on fixed lactation curve test-day mean models, heritability for residual variances ranged between 1.01×10(-3) and 4.17×10(-3) for all traits. The genetic standard deviation in residual variance (i.e., approximately the genetic coefficient of variation of residual variance) ranged between 0.12 and 0.17. Therefore, some genetic variance in micro-environmental sensitivity existed in the Walloon Holstein dairy cattle for the 5 studied traits. The standard deviations due to herd × test-day and permanent environment in residual variance ranged between 0.36 and 0.45 for herd × test-day effect and between 0.55 and 0.97 for permanent environmental effect. Therefore, nongenetic effects also contributed substantially to micro-environmental sensitivity. Addition of random regressions to the mean model did not reduce heterogeneity in residual variance and that genetic heterogeneity of residual variance was not simply an effect of an incomplete mean model.

Variance component and breeding value estimation for genetic heterogeneity of residual variance in Swedish Holstein dairy cattle.

Trait uniformity, or micro-environmental sensitivity, may be studied through individual differences in residual variance. These differences appear to be heritable, and the need exists, therefore, to fit models to predict breeding values explaining differences in residual variance. The aim of this paper is to estimate breeding values for micro-environmental sensitivity (vEBV) in milk yield and somatic cell score, and their associated variance components, on a large dairy cattle data set having more than 1.6 million records. Estimation of variance components, ordinary breeding values, and vEBV was performed using standard variance component estimation software (ASReml), applying the methodology for double hierarchical generalized linear models. Estimation using ASReml took less than 7 d on a Linux server. The genetic standard deviations for residual variance were 0.21 and 0.22 for somatic cell score and milk yield, respectively, which indicate moderate genetic variance for residual variance and imply that a standard deviation change in vEBV for one of these traits would alter the residual variance by 20%. This study shows that estimation of variance components, estimated breeding values and vEBV, is feasible for large dairy cattle data sets using standard variance component estimation software. The possibility to select for uniformity in Holstein dairy cattle based on these estimates is discussed.

Unraveling the genetic architecture of environmental variance of somatic cell score using high-density single nucleotide polymorphism and cow data from experimental farms.

In recent years, it has been shown that not only is the phenotype under genetic control, but also the environmental variance. Very little, however, is known about the genetic architecture of environmental variance. The main objective of this study was to unravel the genetic architecture of the mean and environmental variance of somatic cell score (SCS) by identifying genome-wide associations for mean and environmental variance of SCS in dairy cows and by quantifying the accuracy of genome-wide breeding values. Somatic cell score was used because previous research has shown that the environmental variance of SCS is partly under genetic control and reduction of the variance of SCS by selection is desirable. In this study, we used 37,590 single nucleotide polymorphism (SNP) genotypes and 46,353 test-day records of 1,642 cows at experimental research farms in 4 countries in Europe. We used a genomic relationship matrix in a double hierarchical generalized linear model to estimate genome-wide breeding values and genetic parameters. The estimated mean and environmental variance per cow was used in a Bayesian multi-locus model to identify SNP associated with either the mean or the environmental variance of SCS. Based on the obtained accuracy of genome-wide breeding values, 985 and 541 independent chromosome segments affecting the mean and environmental variance of SCS, respectively, were identified. Using a genomic relationship matrix increased the accuracy of breeding values relative to using a pedigree relationship matrix. In total, 43 SNP were significantly associated with either the mean (22) or the environmental variance of SCS (21). The SNP with the highest Bayes factor was on chromosome 9 (Hapmap31053-BTA-111664) explaining approximately 3% of the genetic variance of the environmental variance of SCS. Other significant SNP explained less than 1% of the genetic variance. It can be concluded that fewer genomic regions affect the environmental variance of SCS than the mean of SCS, but genes with large effects seem to be absent for both traits.

Imputation of genotypes with low-density chips and its effect on reliability of direct genomic values in Dutch Holstein cattle.

Genomic selection using 50,000 single nucleotide polymorphism (50k SNP) chips has been implemented in many dairy cattle breeding programs. Cheap, low-density chips make genotyping of a larger number of animals cost effective. A commonly proposed strategy is to impute low-density genotypes up to 50,000 genotypes before predicting direct genomic values (DGV). The objectives of this study were to investigate the accuracy of imputation for animals genotyped with a low-density chip and to investigate the effect of imputation on reliability of DGV. Low-density chips contained 384, 3,000, or 6,000 SNP. The SNP were selected based either on the highest minor allele frequency in a bin or the middle SNP in a bin, and DAGPHASE, CHROMIBD, and multivariate BLUP were used for imputation. Genotypes of 9,378 animals were used, from which approximately 2,350 animals had deregressed proofs. Bayesian stochastic search variable selection was used for estimating SNP effects of the 50k chip. Imputation accuracies and imputation error rates were poor for low-density chips with 384 SNP. Imputation accuracies were higher with 3,000 and 6,000 SNP. Performance of DAGPHASE and CHROMIBD was very similar and much better than that of multivariate BLUP for both imputation accuracy and reliability of DGV. With 3,000 SNP and using CHROMIBD or DAGPHASE for imputation, 84 to 90% of the increase in DGV reliability using the 50k chip, compared with a pedigree index, was obtained. With multivariate BLUP, the increase in reliability was only 40%. With 384 SNP, the reliability of DGV was lower than for a pedigree index, whereas with 6,000 SNP, about 93% of the increase in reliability of DGV based on the 50k chip was obtained when using DAGPHASE for imputation. Using genotype probabilities to predict gene content increased imputation accuracy and the reliability of DGV and is therefore recommended for applications of imputation for genomic prediction. A deterministic equation was derived to predict accuracy of DGV based on imputation accuracy, which fitted closely with the observed relationship. The deterministic equation can be used to evaluate the effect of differences in imputation accuracy on accuracy and reliability of DGV.

Reliability of direct genomic values for animals with different relationships within and to the reference population.

Accuracy of genomic selection depends on the accuracy of prediction of single nucleotide polymorphism effects and the proportion of genetic variance explained by markers. Design of the reference population with respect to its family structure may influence the accuracy of genomic selection. The objective of this study was to investigate the effect of various relationship levels within the reference population and different level of relationship of evaluated animals to the reference population on the reliability of direct genomic breeding values (DGV). The DGV reliabilities, expressed as squared correlation between estimated and true breeding value, were calculated for evaluated animals at 3 heritability levels. To emulate a trait that is difficult or expensive to measure, such as methane emission, reference populations were kept small and consisted of females with own performance records. A population reflecting a dairy cattle population structure was simulated. Four chosen reference populations consisted of all females available in the first genotyped generation. They consisted of highly (HR), moderately (MR), or lowly (LR) related animals, by selecting paternal half-sib families of decreasing size, or consisted of randomly chosen animals (RND). Of those 4 reference populations, RND had the lowest average relationship. Three sets of evaluated animals were chosen from 3 consecutive generations of genotyped animals, starting from the same generation as the reference population. Reliabilities of DGV predictions were calculated deterministically using selection index theory. The randomly chosen reference population had the lowest average relationship within the reference population. Average reliabilities increased when average relationship within the reference population decreased and the highest average reliabilities were achieved for RND (e.g., from 0.53 in HR to 0.61 in RND for a heritability of 0.30). A higher relationship to the reference population resulted in higher reliability values. At the average squared relationship of evaluated animals to the reference population of 0.005, reliabilities were, on average, 0.49 (HR) and 0.63 (RND) for a heritability of 0.30; 0.20 (HR) and 0.27 (RND) for a heritability of 0.05; and 0.07 (HR) and 0.09 (RND) for a heritability of 0.01. Substantial decrease in the reliability was observed when the number of generations to the reference population increased [e.g., for heritability of 0.30, the decrease from evaluated set I (chosen from the same generation as the reference population) to II (one generation younger than the reference population) was 0.04 for HR, and 0.07 for RND]. In this study, the importance of the design of a reference population consisting of cows was shown and optimal designs of the reference population for genomic prediction were suggested.