Bayesian single-step genomic evaluations combining local and foreign information in Walloon Holsteins.

Most dairy cattle populations found in different countries around the world are small to medium sized and use many artificial insemination bulls imported from different foreign countries. The Walloon population in the southern part of Belgium is a good example for such a small-scale population. Wallonia has also a very active community of Holstein breeders requesting high level genetic evaluation services. Single-step Genomic BLUP (ssGBLUP) methods allow the simultaneous use of genomic, pedigree and phenotypic information and could reduce potential biases in the estimation of genomically enhanced breeding values (GEBV). Therefore, in the context of implementing a Walloon genomic evaluation system for Holsteins, it was considered as the best option. However, in contrast to multi-step genomic predictions, natively ssGBLUP will only use local phenotypic information and is unable to use directly important other sources of information coming from abroad, for example Multiple Across Country Evaluation (MACE) results as provided by the Interbull Center (Uppsala, Sweden). Therefore, we developed and implemented single-step Genomic Bayesian Prediction (ssGBayes), as an alternative method for the Walloon genomic evaluations. The ssGBayes method approximated the correct system of equations directly using estimated breeding values (EBV) and associated reliabilities (REL) without any explicit deregression step. In the Walloon genomic evaluation, local information refers to Walloon EBV and REL and foreign information refers to MACE EBV and associated REL. Combining simultaneously all available genotypes, pedigree, local and foreign information in an evaluation can be achieved but adding contributions to left-hand and right-hand sides subtracting double-counted contributions. Correct propagation of external information avoiding double counting of contributions due to relationships and due to records can be achieved. This ssGBayes method computed more accurate predictions for all types of animals. For example, for genotyped animals with low Walloon REL (<0.25) without MACE results but sired by genotyped bulls with MACE results, the average increase of REL for the studied traits was 0.38 points of which 0.08 points could be traced to the inclusion of MACE information. For other categories of genotyped animals, the contribution by MACE information was also high. The Walloon genomic evaluation system passed for the first time the Interbull GEBV tests for several traits in July 2013. Recent experiences reported here refer to its use in April 2016 for the routine genomic evaluations of milk production, udder health and type traits. Results showed that the proposed methodology should also be of interest for other, similar, populations.

Assessing the effect of pregnancy stage on milk composition of dairy cows using mid-infrared spectra.

Changes in milk production traits (i.e., milk yield, fat, and protein contents) with the pregnancy stage are well documented. To our knowledge, the effect of pregnancy on the detailed milk composition has not been studied so far. The mid-infrared (MIR) spectrum reflects the detailed composition of a milk sample and is obtained by a nonexhaustive and widely used method for milk analysis. Therefore, this study aimed to investigate the effect of pregnancy on milk MIR spectrum in addition to milk production traits (milk yield, fat, and protein contents). A model including regression on the number of days pregnant was applied on milk production traits (milk yield, fat, and protein contents) and on 212 spectral points from the MIR spectra of 9,757 primiparous Holstein cows from Walloon herds. Effects of pregnancy stage were expressed on a relative scale (effect divided by the squared root of the phenotypic variance); this allowed comparisons between effects on milk traits and on 212 spectral points. Effect of pregnancy stage on production traits were in line with previous studies indicating that the model accounted well for the pregnancy effect. Trends of the relative effect of the pregnancy stage on the 212 spectral points were consistent with known and observed effect on milk traits. The highest effect of the pregnancy was observed in the MIR spectral region from 968 to 1,577 cm-1. For some specific wavenumbers, the effect was higher than for fat and protein contents in the beginning of the pregnancy (from 30 to 90 or 120 d pregnant). In conclusion, the effect of early pregnancy can be observed in the detailed milk composition through the analysis of the MIR spectrum of bovine milk. Further analyses are warranted to explore deeply the use of MIR spectra of bovine milk for breeding and management of dairy cow pregnancy.

Genetic parameters for direct and maternal calving ease in Walloon dairy cattle based on linear and threshold models.

Calving ease scores from Holstein dairy cattle in the Walloon Region of Belgium were analysed using univariate linear and threshold animal models. Variance components and derived genetic parameters were estimated from a data set including 33,155 calving records. Included in the models were season, herd and sex of calf × age of dam classes × group of calvings interaction as fixed effects, herd × year of calving, maternal permanent environment and animal direct and maternal additive genetic as random effects. Models were fitted with the genetic correlation between direct and maternal additive genetic effects either estimated or constrained to zero. Direct heritability for calving ease was approximately 8% with linear models and approximately 12% with threshold models. Maternal heritabilities were approximately 2 and 4%, respectively. Genetic correlation between direct and maternal additive effects was found to be not significantly different from zero. Models were compared in terms of goodness of fit and predictive ability. Criteria of comparison such as mean squared error, correlation between observed and predicted calving ease scores as well as between estimated breeding values were estimated from 85,118 calving records. The results provided few differences between linear and threshold models even though correlations between estimated breeding values from subsets of data for sires with progeny from linear model were 17 and 23% greater for direct and maternal genetic effects, respectively, than from threshold model. For the purpose of genetic evaluation for calving ease in Walloon Holstein dairy cattle, the linear animal model without covariance between direct and maternal additive effects was found to be the best choice.

Progressive resistance training did not improve walking but can improve muscle performance, quality of life and fatigue in adults with multiple sclerosis: a randomized controlled trial.

BACKGROUND: Few high-quality trials have examined the effects of progressive resistance training (PRT) on people with multiple sclerosis (MS). OBJECTIVE: To determine the effectiveness of PRT for people with MS, focusing on improving the gait deficits common in this population. METHODS: Using a single blind randomized controlled trial, people with relapsing-remitting MS were randomly allocated to either a PRT program targeting the lower limb muscles twice a week for 10 weeks (n = 36), or usual care plus an attention and social program conducted once a week for 10 weeks (n = 35). Outcomes were recorded at baseline, week 10 and week 22. RESULTS: Participants attended 92% of training sessions, with no serious adverse events. At 10 weeks, no differences were detected in walking performance. However, compared with the comparison group PRT demonstrated increased leg press strength (16.8%, SD 4.5), increased reverse leg press strength (29.8%, SD 12.7), and increased muscle endurance of the reverse leg press (38.7%, SD 32.8). Improvements in favor of PRT were also found for physical fatigue (Mean difference -3.9 units, 95%CI -6.6 to -1.3), and the physical health domain of quality of life (Mean difference 1.5 units, 95%CI 0.1 to 2.9). At week 22 almost no between-group differences remained. CONCLUSION: PRT is a relatively safe intervention that can have short-term effects on reducing physical fatigue, increasing muscle endurance and can lead to small improvements in muscle strength and quality of life in people with relapsing-remitting MS. However, no improvements in walking performance were observed and benefits do not appear to persist if training is completely stopped.

Modeling milk urea of Walloon dairy cows in management perspectives.

The aim of this study was to develop an adapted random regression test-day model for milk urea (MU) and to study the possibility of using predictions and solutions given by the model for management purposes. Data included 607,416 MU test-day records of first-lactation cows from 632 dairy herds in the Walloon Region of Belgium. Several advanced features were used. First, to detect the herd influence, the classical herd x test-day effect was split into 3 new effects: a fixed herd x year effect, a fixed herd x month-period effect, and a random herd test-day effect. A fixed time period regression was added in the model to take into account the yearly oscillations of MU on a population scale. Moreover, first autoregressive processes were introduced and allowed us to consider the link between successive test-day records. The variance component estimation indicated that large variance was associated with the random herd x test-day effect (48% of the total variance), suggesting the strong influence of herd management on the MU level. The heritability estimate was 0.13. By comparing observed and predicted MU levels at both the individual and herd levels, target ranges for MU concentrations were defined to take into account features of each cow and each herd. At the cow level, an MU record was considered as deviant if it was <200 or >400 mg/L (target range used in the field) and if the prediction error was >50 mg/L (indicating a significant deviation from the expected level). Approximately 7.5% of the MU records collected between June 2007 and May 2008 were beyond these thresholds. This combination allowed for the detection of potentially suspicious cows. At the herd level, the expected MU level was considered as the sum of the solutions for specific herd effects. A herd was considered as deviant from its target range when the prediction error was greater than the standard deviation of MU averaged by herd test day. Results showed that 6.7% of the herd test-day MU levels between June 2007 and May 2008 were considered deviant. These deviations seemed to occur more often during the grazing period. Although theoretical considerations developed in this study should be validated in the field, this research showed the potential use of a test-day model for analyzing functional traits to advise dairy farmers.

Estimation of heritability and genetic correlations for the major fatty acids in bovine milk.

The current cattle selection program for dairy cattle in the Walloon region of Belgium does not consider the relative content of the different fatty acids (FA) in milk. However, interest by the local dairy industry in differentiated milk products is increasing. Therefore, farmers may be interested in selecting their animals based on the fat composition. The aim of this study was to evaluate the feasibility of genetic selection to improve the nutritional quality of bovine milk fat. The heritabilities and correlations among milk yield, fat, protein, and major FA contents in milk were estimated. Heritabilities for FA in milk and fat ranged from 5 to 38%. The genetic correlations estimated among FA reflected the common origin of several groups of FA. Given these results, an index including FA contents with the similar metabolic process of production in the mammary gland could be used, for example, to increase the monounsaturated and conjugated fatty acids in milk. Moreover, the genetic correlations between the percentage of fat and the content of C14:0, C12:0, C16:0, and C18:0 in fat were -0.06, 0.55, 0.60, and 0.84, respectively. This result demonstrates that an increase in fat content is not directly correlated with undesirable changes in FA profile in milk for human health. Based on the obtained genetic parameters, a future selection program to improve the FA composition of milk fat could be initiated.

Linear and curvilinear effects of inbreeding on production traits for walloon Holstein cows.

The nonlinear effects of inbreeding were studied by comparing linear and curvilinear regression models of phenotypic performances on inbreeding coefficients for production traits (milk, fat, and protein yields) of Holstein cows in their first lactation. Three different regression models (linear, quadratic, and cubic) were introduced separately into a single-trait, single-lactation, random regression test-day model. The significance of the different regression coefficients was studied based on a t-test after estimation of error variances and covariances associated with the different regression coefficients. All of the tested regression coefficients were significantly different from 0. The traditional regression coefficients of milk, fat, and protein yields on inbreeding were, respectively, -22.10, -1.10, and -0.72 kg for Holstein cows in their first lactation. However, the estimates of 305-d production losses for various classes of animals based on inbreeding coefficients showed that the effect of inbreeding was not a linear function of the percentage of inbreeding. The 305-d milk yield loss profiles attributable to inbreeding, obtained by the various regression models, were different. However, for inbreeding coefficients between 0 and 10%, these differences were small.

Variation in fatty acid contents of milk and milk fat within and across breeds.

The aim of this research was to study the potential for selection of cows with a higher nutritional quality of milk fat by studying the differences in fatty acid profiles within and across the following breeds: Dual Purpose Belgian Blue, Holstein-Friesian, Jersey, Montbeliarde, and non-Holstein Meuse-Rhine-Yssel type Red and White. Six hundred milk samples from 275 animals were taken from 7 herds. Several types of fatty acids in milk and milk fat were quantified using mid-infrared spectrometry and previously obtained calibration equations. Statistical analyses were made using a mixed linear model with a random animal effect. The variance components were estimated by using REML. Results showed breed differences for the fatty acid profile. The repeatability estimate obtained in the present study may suggest the existence of moderate additive genetic variance for the fatty acid profile within each breed. Results also indicated variation for each analyzed milk component in the whole cow population studied. Genetic improvement of the nutritional quality of milk fat based on fatty acid profiles might be possible, and further research and development are warranted.

Inbreeding depression for global and partial economic indexes, production, type, and functional traits.

The objective of this research was to examine the effects of inbreeding in the population of Holstein cattle in the Walloon region of Belgium. The effects of inbreeding on the global economic index and its components were studied by using data from the genetic evaluations of February 2004 for production, somatic cell score (SCS), computed from somatic cell counts and type. Inbreeding coefficients for 956,516 animals were computed using a method that allows assigning an inbreeding coefficient to individuals without known parents. These coefficients were equal to the mean inbreeding coefficient of contemporary individuals with known parents. The significance of inbreeding effects on the different evaluated traits and on the different indexes were tested using a t-test comparing estimated standard errors and effects. The inbreeding effect was significantly different from zero for the vast majority of evaluated traits and for all of the indexes. Inbreeding had the greatest deleterious effects on production traits. Inbreeding decreased yield of milk, fat, and protein during a lactation by 19.68, 0.96, and 0.69 kg, respectively, per each 1% increase in inbreeding. The regression coefficient of SCS per 1% increase in inbreeding was +0.005 SCS units. The inbreeding depression was thus relatively low for SCS, but inbred animals had higher SCS than non-inbred animals, indicating that inbred animals would be slightly more sensitive to mastitis than non-inbred animals. Estimates of inbreeding effects on evaluated type traits per 1% increase were small. The most strongly affected type traits were chest width, rear leg, and overall development on a standardized scale. For several type traits, particularly traits linked to the udder, the estimates suggested a favorable effect of inbreeding. The global economic index was depressed by around 6.13 euro of lifetime profit per 1% increase in inbreeding for the Holstein animals in the Walloon region of Belgium.

Adjustment for heterogeneous covariance due to herd milk yield by transformation of test-day random regressions.

A method of accounting for differences in covariance components of test-day milk records was developed based on transformation of regressions for random effects. Preliminary analysis indicated that genetic and nongenetic covariance structures differed by herd milk yield. Differences were found for phenotypic covariances and also for genetic, permanent environmental, and herd-time covariances. Heritabilities for test-day milk yield tended to be lower at the end and especially at the start of lactation; they also were higher (maximum of approximately 25%) for high-yield herds and lower (maximum of 15%) for low-yield herds. Permanent environmental variances were on average 10% lower in high-yield herds. Relative herd-time variances were approximately 10% at start of lactation and then began to decrease regardless of herd yield; high-yield herds increased in midlactation followed by another decrease, and medium-yield herds increased at the end of lactation. Regressors for random regression effects were transformed to adjust for heterogeneity of test-day yield covariances. Some animal reranking occurred because of this transformation of genetic and permanent environmental effects. When genetic correlations between environments were allowed to differ from 1, some additional animal re-ranking occurred. Correlations of variances of genetic and permanent-environmental regression solutions within herd, test-day, and milking frequency class with class mean milk yields were reduced with adjustment for heterogeneous covariance. The method suggests a number of innovative solutions to issues related to heterogeneous covariance structures, such as adjusted estimates in multibreed evaluation.

Estimated heterogeneity of phenotypic variance of test-day yield with a structural variance model.

First-lactation test-day milk, fat, and protein yields from New York, Wisconsin, and California herds from 1990 through 2000 were adjusted additively for age and lactation stage. A random regression model with third-order Legendre polynomials for permanent environmental and genetic effects was used. The model included a random effect with the same polynomial regressions for 2 yr of calvings within herd (herd-time effect) to provide herd-specific lactation curves that can change every 2 yr. (Co)variance components were estimated using expectation-maximization REML simultaneously with phenotypic variances that were modeled using a structural variance model. Maximum heritability for test-day milk yield was estimated to be approximately 20% around 200 to 250 d in milk; heritabilities were slightly lower for test-day fat and protein yields. Herd-time effects explained 12 to 20% of phenotypic variance and had the greatest impact at start of lactation. Variances of test-day yields increased with time, subclass size, and milking frequency. Test month had limited influence on variance. Variance increased for cows in herds with low and high milk yields and for early and late lactation stages. Repeatabilities of variances observed for a given class of herd, test-day, and milking frequency were 14 to 17% across nested variance subclasses based on lactation stage.