Genetic evaluation for stillbirth and pre-weaning mortality in Australian dairy cattle

The welfare of calves is important to both farmers and consumers. Practices that increase the proportion of calves born alive and enable them to thrive through to weaning contribute to improved sustainability. Stillbirths (SB) are calvings where the calf dies at birth or within 24 h after birth. Pre-weaning mortality (PWM) refers to calves that die after the first day of life but before weaning based on termination data. Both SB and PWM are binary traits characterized by low heritability. Data collection for these traits is incomplete, compared with traits like milk yield in cows. Despite these challenges, genetic variation can be measured and used to produce breeding tools, such as estimated breeding values (EBV), to reduce calf mortality over time. The aim of this study was to compare the performance of various linear models to predict SB and PWM traits in Holstein and Jersey cattle and evaluate their applicability for industry-wide use in the Australian dairy industry. Calving records from around 2.25 million Holstein and Jersey dams were obtained from DataGene's Central Data Repository from 2000 onwards to calculate genetic parameters. About 7% of calves were recorded as stillborn in the period 2000-2021 (n = 1.48 million calvings). The prevalence of PWM was much lower than stillbirth during the same period at 2% (n = 0.89 million calves). Genetic parameters were estimated for SB direct, SB maternal and PWM using bivariate linear models with calving ease (CE) as the second trait in the model. The heritability of these calf traits was low and varied between 1 to 5% depending on the breed, trait and model. In Holstein cattle, heritabilities were 2% for PWM and SB direct and 1% for SB maternal while in Jersey cattle heritabilities were 5% for PWM, 2% for SB direct and 1% for SB maternal. The genetic trends for both SB direct and maternal in Holstein cattle indicate improvement in both traits whereas there was no apparent increase or decrease in PWM in the past 2 decades. The coefficient of genetic variation for SB direct and PWM was between 11.7 and 23.0% in Holstein and Jersey cattle demonstrating that there was considerable genetic variation in calf survival traits as a first step to using genetic selection to increase the proportion of calves born alive and calves weaned. A focus on improved calf and calving recording practices is expected to increase the reliability of genetic predictions.

Genetic parameters and evaluation of mortality and slaughter rate in Holstein and Jersey cows

The longevity of dairy cattle has economic, animal welfare, and health implications and is influenced by the frequency of mortality on the farm and sale for slaughter. In this study cows removed from the herd due to death or slaughter during the lactation were coded 1 and cows that were not terminated were coded 0. Genetic parameters for mortality rates (MR) and slaughter rates (SR) were estimated for Holstein (H) and Jersey (J) breeds by applying both linear (LM) and threshold (TM) sire models using about 1.2 million H and 286,000 J cows. Estimated breeding values (EBV) for MR and SR were predicted using animal models to assess the opportunity for selection and genetic trends. Cow termination data, recorded between 1990 and 2020 on a voluntary basis by Australian dairy farmers, were analyzed. Cow MR has increased from below 1% in the 1990s to 4.1% and 3.6% in recent years in H and J cows, respectively. Most dead cows (∼36%) left the herd before 120 d of lactation, while cows that were slaughtered left the herd toward the end of the lactation. Using the LM, heritability (h2) estimates for MR were lower (1%) than those for SR (2%-3.5%). When h2 were estimated using a TM, the estimates for both traits varied between 4% and 20%, suggesting that the difference in incidence level is one of the reasons for the difference in the h2 values between MR and SR. Early test-day milk yield (MY) and 305-d MY (305-d MY) have unfavorable genetic correlations (0.32-0.41) with MR in both breeds. The genetic correlations of calving interval with MR were stronger (0.54-0.68) than with SR (0.28-0.45) suggesting that poor fertility can serve as an early indicator of poor cow health that may lead to increased risk of death. High early test-day somatic cell count is genetically associated with increased likelihood of slaughter (0.24-0.46), but not with increased likelihood of death. In H, 305-d protein yield (PY) had the strongest genetic correlation (-0.34 to -0.40) with SR whereas in J, both 305-d PY and fat yield showed high genetic (-0.64 to -0.70) and moderate environmental (-0.35 to -0.37) correlations with SR. The genetic correlation of removal from the herd due to death and slaughter was negative (-0.3) in J and zero in H. Strong selection for improved fertility and survival and less selection emphasis for MY, has led to an improvement in the genetic trend for cow MR in H and the trend in J has stabilized. Although genetic evaluations for cow MR are feasible, the reliabilities of the EBV are low and the level of cow MR in Australia are relatively low compared with similar countries. Therefore, genetic evaluation for survival based on mortality and slaughter data could be sufficient in the current selection circumstances where breeding objectives are broadly defined. Nevertheless, all Australian farmers should be encouraged to continue recording mortality and slaughter data for monitoring of the trends and for future development of genetic evaluations.

Genome-wide assessment and mapping of inbreeding depression identifies candidate genes associated with semen traits in Holstein bulls

BACKGROUND

The reduction in phenotypic performance of a population due to mating between close relatives is called inbreeding depression. The genetic background of inbreeding depression for semen traits is poorly understood. Thus, the objectives were to estimate the effect of inbreeding and to identify genomic regions underlying inbreeding depression of semen traits including ejaculate volume (EV), sperm concentration (SC), and sperm motility (SM). The dataset comprised ~ 330 K semen records from ~ 1.5 K Holstein bulls genotyped with 50 K single nucleotide polymorphism (SNP) BeadChip. Genomic inbreeding coefficients were estimated using runs of homozygosity (i.e., F_ROH > 1 Mb) and excess of SNP homozygosity (F_SNP). The effect of inbreeding was estimated by regressing phenotypes of semen traits on inbreeding coefficients. Associated variants with inbreeding depression were also detected by regressing phenotypes on ROH state of the variants.

RESULTS

Significant inbreeding depression was observed for SC and SM (p < 0.01). A 1% increase in F_ROH reduced SM and SC by 0.28% and 0.42% of the population mean, respectively. By splitting F_ROH into different lengths, we found significant reduction in SC and SM due to longer ROH, which is indicative of more recent inbreeding. A genome-wide association study revealed two signals positioned on BTA 8 associated with inbreeding depression of SC (p < 0.00001; FDR < 0.02). Three candidate genes of GALNTL6, HMGB2, and ADAM29, located in these regions, have established and conserved connections with reproduction and/or male fertility. Moreover, six genomic regions on BTA 3, 9, 21 and 28 were associated with SM (p < 0.0001; FDR < 0.08). These genomic regions contained genes including PRMT6, SCAPER, EDC3, and LIN28B with established connections to spermatogenesis or fertility.

CONCLUSIONS

Inbreeding depression adversely affects SC and SM, with evidence that longer ROH, or more recent inbreeding, being especially detrimental. There are genomic regions associated with semen traits that seems to be especially sensitive to homozygosity, and evidence to support some from other studies. Breeding companies may wish to consider avoiding homozygosity in these regions for potential artificial insemination sires.

Genomic prediction and selection response for grain yield in safflower

In plant breeding programs, multiple traits are recorded in each trial, and the traits are often correlated. Correlated traits can be incorporated into genomic selection models, especially for traits with low heritability, to improve prediction accuracy. In this study, we investigated the genetic correlation between important agronomic traits in safflower. We observed the moderate genetic correlations between grain yield (GY) and plant height (PH, 0.272–0.531), and low correlations between grain yield and days to flowering (DF, −0.157–0.201). A 4%–20% prediction accuracy improvement for grain yield was achieved when plant height was included in both training and validation sets with multivariate models. We further explored the selection responses for grain yield by selecting the top 20% of lines based on different selection indices. Selection responses for grain yield varied across sites. Simultaneous selection for grain yield and seed oil content (OL) showed positive gains across all sites with equal weights for both grain yield and oil content. Combining g×E interaction into genomic selection (GS) led to more balanced selection responses across sites. In conclusion, genomic selection is a valuable breeding tool for breeding high grain yield, oil content, and highly adaptable safflower varieties.

Meta-analysis of genetic parameters for growth traits in meat, wool and dual-purpose sheep breeds in the world using a random-effects model

BACKGROUND

There is large variation in genetic parameters in literature for growth traits in sheep. Reliable estimation of genetic parameters is required for developing breeding programmes.

OBJECTIVES

The aim of this study was to aggregate results of different studies by meta-analysis to improve reliability of estimated parameters.

METHODS

In the current study, 221 papers that have been published between 1995 and 2021 were reviewed. Using a random-effects model in the Comprehensive Meta-Analysis software, direct and maternal heritabilities, as well as, genetic and phenotypic correlations between growth traits were estimated in meat (M), wool (W) and dual-purpose (D) sheep breeds. The growth traits in this study were birth weight, 3-month weight, 6-month weight, 9-month weight and yearling weight.

RESULTS

The combined direct heritability was the lowest for birth weight (0.190 ± 0.004, 0.198 ± 0.003 and 0.196 ± 0.004 for M, W and D breeds, respectively) and the highest for yearling weight (0.264 ± 0.010, 0.304 ± 0.005 and 0.285 ± 0.020 for M, W and D breeds, respectively). The maternal heritability was the lowest for yearling weight (0.085 ± 0.003, 0.055 ± 0.002 and 0.052 ± 0.005 for M, W and D breeds, respectively) and the highest for 6-month weight (0.240 ± 0.088, 0.164 ± 0.001 and 0.162 ± 0.006 for M, W and D breeds, respectively). The phenotypic and genetic correlations were lower between the weights measured at more distant intervals. The lowest genetic correlation was observed between birth weight and yearling weight (0.290 ± 0.051 for W breeds).

CONCLUSIONS

The small standard errors could indicate that the aggregation of results from different studies improved the reliability of estimated parameters and reduced range of 95% confidence intervals. Hence, the results could be used with greater level of confidence in sheep breeding programmes.

Sharing of either phenotypes or genetic variants can increase the accuracy of genomic prediction of feed efficiency

Background

Sharing individual phenotype and genotype data between countries is complex and fraught with potential errors, while sharing summary statistics of genome-wide association studies (GWAS) is relatively straightforward, and thus would be especially useful for traits that are expensive or difficult-to-measure, such as feed efficiency. Here we examined: (1) the sharing of individual cow data from international partners; and (2) the use of sequence variants selected from GWAS of international cow data to evaluate the accuracy of genomic estimated breeding values (GEBV) for residual feed intake (RFI) in Australian cows.

Results

GEBV for RFI were estimated using genomic best linear unbiased prediction (GBLUP) with 50k or high-density single nucleotide polymorphisms (SNPs), from a training population of 3797 individuals in univariate to trivariate analyses where the three traits were RFI phenotypes calculated using 584 Australian lactating cows (AUSc), 824 growing heifers (AUSh), and 2526 international lactating cows (OVE). Accuracies of GEBV in AUSc were evaluated by either cohort-by-birth-year or fourfold random cross-validations. GEBV of AUSc were also predicted using only the AUS training population with a weighted genomic relationship matrix constructed with SNPs from the 50k array and sequence variants selected from a meta-GWAS that included only international datasets. The genomic heritabilities estimated using the AUSc, OVE and AUSh datasets were moderate, ranging from 0.20 to 0.36. The genetic correlations (r_g) of traits between heifers and cows ranged from 0.30 to 0.95 but were associated with large standard errors. The mean accuracies of GEBV in Australian cows were up to 0.32 and almost doubled when either overseas cows, or both overseas cows and AUS heifers were included in the training population. They also increased when selected sequence variants were combined with 50k SNPs, but with a smaller relative increase.

Conclusions

The accuracy of RFI GEBV increased when international data were used or when selected sequence variants were combined with 50k SNP array data. This suggests that if direct sharing of data is not feasible, a meta-analysis of summary GWAS statistics could provide selected SNPs for custom panels to use in genomic selection programs. However, since this finding is based on a small cross-validation study, confirmation through a larger study is recommended.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00749-z.

Erratum to "Evaluation of updated Feed Saved breeding values developed in Australian Holstein dairy cattle" (JDS Commun. 3:114-119)

[This corrects the article DOI: 10.3168/jdsc.2021-0150.].

Combining NDVI and Bacterial Blight Score to Predict Grain Yield in Field Pea

Field pea is the most commonly grown temperate pulse crop, with close to 15 million tons produced globally in 2020. Varieties improved through breeding are important to ensure ongoing improvements in yield and disease resistance. Genomic selection (GS) is a modern breeding approach that could substantially improve the rate of genetic gain for grain yield, and its deployment depends on the prediction accuracy (PA) that can be achieved. In our study, four yield trials representing breeding lines' advancement stages of the breeding program (S0, S1, S2, and S3) were assessed with grain yield, aerial high-throughput phenotyping (normalized difference vegetation index, NDVI), and bacterial blight disease scores (BBSC). Low-to-moderate broad-sense heritability (0.31-0.71) and narrow-sense heritability (0.13-0.71) were observed, as the estimated additive and non-additive genetic components for the three traits varied with the different models fitted. The genetic correlations among the three traits were high, particularly in the S0-S2 stages. NDVI and BBSC were combined to investigate the PA for grain yield by univariate and multivariate GS models, and multivariate models showed higher PA than univariate models in both cross-validation and forward prediction methods. A 6-50% improvement in PA was achieved when multivariate models were deployed. The highest PA was indicated in the forward prediction scenario when the training population consisted of early generation breeding stages with the multivariate models. Both NDVI and BBSC are commonly used traits that could be measured in the early growth stage; however, our study suggested that NDVI is a more useful trait to predict grain yield with high accuracy in the field pea breeding program, especially in diseased trials, through its incorporation into multivariate models.

Evaluation of updated Feed Saved breeding values developed in Australian Holstein dairy cattle

•

The Feed Saved (FS) estimated breeding value (EBV) was updated by doubling the number of Australian and overseas cows.

•

The reliability of the residual feed intake component of FS has increased from 11% (2015 model) to 20% (current model).

•

The mean reliability of FS EBV in Holstein bulls that were born in the last 10 years has improved by 10%.

•

The genetic trend of FS EBV has been stabilizing since 2015.

Although selection for increased milk production traits has led to a genetic increase in body weight (BW), the genetic gain in milk production has exceeded the gain in BW, so gross feed efficiency has improved. Nonetheless, greater gains may be possible by directly selecting for a measure of feed efficiency. Australia first introduced Feed Saved (FS) estimated breeding value (EBV) in 2015. Feed Saved combines residual feed intake (RFI) genomic EBV and maintenance requirements calculated from mature BW EBV. The FS EBV was designed to enable the selection of cows for reduced energy requirements with similar milk production. In this study, we used a reference population of 3,711 animals in a multivariate analysis including Australian heifers (AUSh), Australian cows (AUSc), and overseas cows (OVEc) to update the Australian EBV for lifetime RFI (i.e., a breeding value that incorporated RFI in growing and lactating cows) and to recalculate the FS EBV in Australian Holstein bulls (AUSb). The estimates of genomic heritabilities using univariate (only AUSc or AUSh) to trivariate (including the OVEc) analyses were similar. Genomic heritabilities for RFI were estimated as 0.18 for AUSc, 0.27 for OVEc, and 0.36 for AUSh. The genomic correlation for RFI between AUSc and AUSh was 0.47 and that between AUSc and OVEc was 0.94, but these estimates were associated with large standard errors (range: 0.18–0.28). The reliability of lifetime RFI (a component of FS) in the trivariate analysis (i.e., including OVEc) increased from 11% to 20% compared with the 2015 model and was greater, by 12%, than in a bivariate analysis in which the reference population included only AUSc and AUSh. By applying the prediction equation of the 2020 model, the average reliability of the FS EBV in 20,816 AUSb that were born between 2010 and 2020 improved from 33% to 43%. Previous selection strategies—that is, using the predecessor of the Balanced Performance Index (Australian Profit Ranking index) that did not include FS—have resulted in an unfavorable genetic trend in FS. However, this unfavorable trend has stabilized since 2015, when FS was included in the Balanced Performance Index, and is expected to move in a favorable direction with selection on Balanced Performance Index or the Health Weighted Index. Doubling the reference population, particularly by incorporating international data for feed efficiency, has improved the reliability of the FS EBV. This could lead to increased genetic gain for feed efficiency in the Australian industry.

Machine learning approaches for the prediction of lameness in dairy cows

Lameness is one of the costliest health problems, as well as a welfare concern in dairy cows. However, it is difficult to detect cows with possible lameness, or the ones that are at risk of becoming lame e.g. in the next week or so. In this study, we investigated the ability of three machine learning algorithms, Naïve Bayes (NB), Random Forest (RF) and Multilayer Perceptron (MLP), to predict cases of lameness using milk production and conformation traits. The performance of these algorithms was compared with logistic regression (LR) as the gold standard approach for binary classification. We had a total of 2 535 lameness scores (2 248 sound and 287 unsound) and 29 predictor features from nine dairy herds in Australia to predict lameness incidence. Training was done on 80% of the data within each herd with the remainder used as validation set. Our results indicated that in terms of area under curve of receiver operating characteristics, there were negligible differences between LR (0.67) and NB (0.66) while MLP (0.62) and RF (0.61) underperformed compared to the other two methods. However, the F1-score in NB (27%) outperformed LR (1%), suggesting that NB could potentially be a more reliable method for the prediction of lameness in practice, given enough relevant data are available for proper training, which was a limitation in this study. Considering the small size of our dataset, lack of information about environmental conditions prior to the incidence of lameness, management practices, short time gap between production records and lameness scoring, and farm information, this study proved the concept of using machine learning predictive models to predict the incidence of lameness a priori to its occurrence and thus may become a valuable decision support system for better lameness management in precision dairy farming.

Identification of sequence variation in the oocyte-derived bone morphogenetic protein 15 (BMP15) gene (BMP15) associated with litter size in New Zealand sheep (Ovis aries) breeds

Advances in the study of reproductive traits indicate that functional variation in fertility genes may be useful for improving sheep fertility. The aim of this study was to search for variation in the bone morphogenetic protein 15 gene (BMP15) and ascertain any association with litter size in purebred Finnish Landrace sheep (n = 148), Finnish Landrace × Texel-cross sheep (n = 45), and composite sheep (of varying breed background; n = 58) from New Zealand (NZ). A 482 bp and 312 bp fragment of exon 1 and 2, respectively, of BMP15 were analysed using polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). The additive and dominance effect of BMP15 variation on litter size were estimated using animal and sire models. Two variants (A and B) were detected in exon 1; no sequence variation was detected in exon 2. Variant A had the nucleotide sequence CTT between positions c.31 and c.33, while variant B had a deletion (c.31_33del). The observed frequency for variant A in the Finnish Landrace sheep, Finnish Landrace × Texel-cross sheep and the composite sheep, was 0.77, 0.92, and 0.68, respectively while the frequency of variant B (c.31_33del) was 0.23, 0.08, and 0.32, respectively. An association between litter size and c.31_33del (P < 0.001) was observed in composite sheep. Analysis of more sheep will be required to confirm these results. Litter size did not differ significantly between sheep breeds regardless of the presence/absence of c.31_33del. Results suggested that c.31_33del might be a genetic marker for improving fecundity in some NZ sheep.

Including milk production, conformation, and functional traits in multivariate models for genetic evaluation of lameness

Lameness is a serious health and welfare issue that can negatively affect the economic performance of cows, especially on pasture-based dairy farms. However, most genetic predictions (GP) of lameness have low accuracy because lameness data are often incomplete as data are collected voluntarily by farmers in countries such as Australia. The objective of this study was to find routinely measured traits that are correlated with lameness and use them in multivariate evaluation models to improve the accuracy of GP for lameness. We used health events and treatments associated with lameness recorded by Australian farmers from 2002 to early 2019. The lameness incidence rates in Holstein and Jersey cows were 3.3% and 4.6%, respectively. We analyzed the records of 36 other traits (milk production, conformation, fertility, and survival traits) to estimate genetic correlations with lameness. The estimated heritability ± standard error (and repeatability ± standard error) for lameness in both Holstein and Jersey breeds were very low: 0.007 ± 0.002 (and 0.029 ± 0.002) and 0.005 ± 0.003 (and 0.027 ± 0.006), respectively, in univariate sire models. For the GP models, we tested including measurements of overall type to prediction models for Holsteins, stature and body length for Jersey, and milk yield and fertility traits for both breeds. The average accuracy of GP, calculated from prediction error variances, were 0.38 and 0.24 for Holstein and Jersey sires, respectively, when estimated using univariate sire models and both increased to 0.43 using multivariate sire models. In conclusion, we found that the accuracy of GP for lameness could be improved by including genetically correlated traits in a multivariate model. However, to further improve the accuracy of predictions of lameness, precise identification and recording incidences of hoof or leg disorder, or large-scale recording of locomotion and claw scores by trained personnel should be considered.

A conditional multi-trait sequence GWAS discovers pleiotropic candidate genes and variants for sheep wool, skin wrinkle and breech cover traits

Background

Imputation to whole-genome sequence is now possible in large sheep populations. It is therefore of interest to use this data in genome-wide association studies (GWAS) to investigate putative causal variants and genes that underpin economically important traits. Merino wool is globally sought after for luxury fabrics, but some key wool quality attributes are unfavourably correlated with the characteristic skin wrinkle of Merinos. In turn, skin wrinkle is strongly linked to susceptibility to “fly strike” (Cutaneous myiasis), which is a major welfare issue. Here, we use whole-genome sequence data in a multi-trait GWAS to identify pleiotropic putative causal variants and genes associated with changes in key wool traits and skin wrinkle.

Results

A stepwise conditional multi-trait GWAS (CM-GWAS) identified putative causal variants and related genes from 178 independent quantitative trait loci (QTL) of 16 wool and skin wrinkle traits, measured on up to 7218 Merino sheep with 31 million imputed whole-genome sequence (WGS) genotypes. Novel candidate gene findings included the MAT1A gene that encodes an enzyme involved in the sulphur metabolism pathway critical to production of wool proteins, and the ESRP1 gene. We also discovered a significant wrinkle variant upstream of the HAS2 gene, which in dogs is associated with the exaggerated skin folds in the Shar-Pei breed.

Conclusions

The wool and skin wrinkle traits studied here appear to be highly polygenic with many putative candidate variants showing considerable pleiotropy. Our CM-GWAS identified many highly plausible candidate genes for wool traits as well as breech wrinkle and breech area wool cover.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00651-0.

Improving the accuracy of predictions for cow survival by multivariate evaluation model

Context Cow survival measures the ability of cows to survive from the current to subsequent lactation. In addition to economic gain, genetic selection for survival could improve animal welfare by increasing the adaptability and resilience of the cows to both environmental and health challenges. However, survival is a complex trait because it results from a diverse range of reasons for culling of cows from the herd. Consequently, the accuracy of genetic predictions of direct survival are often low. Aims Our aim was to increase the accuracy of predictions of survival in Holstein and Jersey sires by including important predictor traits in multi-trait evaluation models. Methods Phenotypic and genetic correlations between survival trait deviations (TDs) and 35 routinely measured traits (including milk yield, fertility and type traits) were estimated using bivariate sire models. Survival TDs for 538 394 Holstein and 63 839 Jersey cows were used in our study; these cows or their close relatives also had milk, fertility and type traits records between 2002 and 2019. These genetic parameters were required to assess the potential usefulness of predictor traits for the prediction of survival. Key results Survival was genetically correlated with milk, fat and protein yields, overall type, composite mammary system and fertility TDs in both Holstein and Jersey. Further, most of the type traits related to feet and legs, and rump, were also correlated with survival TDs in Jersey. For sires, the accuracy of predictions for survival increased by 0.05 for Holsteins (from 0.54 to 0.59) and for Jerseys (from 0.48 to 0.53) through the use of multivariate models compared with univariate models. Conclusions Survival was genetically associated with traits affecting voluntary and involuntary culling and when included in multi-trait genetic evaluation models, they moderately improved the accuracy of genetic prediction of survival. Implications Predictor traits can be used to increase the accuracy of predictions of survival through the use of multi-trait models. The inclusion of breed-specific predictor traits should be considered, especially for Jerseys in genetic evaluations of survival.

Improving Genomic Prediction of Crossbred and Purebred Dairy Cattle

This study assessed the accuracy and bias of genomic prediction (GP) in purebred Holstein (H) and Jersey (J) as well as crossbred (H and J) validation cows using different reference sets and prediction strategies. The reference sets were made up of different combinations of 36,695 H and J purebreds and crossbreds. Additionally, the effect of using different sets of marker genotypes on GP was studied (conventional panel: 50k, custom panel enriched with, or close to, causal mutations: XT_50k, and conventional high-density with a limited custom set: pruned HDnGBS). We also compared the use of genomic best linear unbiased prediction (GBLUP) and Bayesian (emBayesR) models, and the traits tested were milk, fat, and protein yields. On average, by including crossbred cows in the reference population, the prediction accuracies increased by 0.01–0.08 and were less biased (regression coefficient closer to 1 by 0.02–0.16), and the benefit was greater for crossbreds compared to purebreds. The accuracy of prediction increased by 0.02 using XT_50k compared to 50k genotypes without affecting the bias. Although using pruned HDnGBS instead of 50k also increased the prediction accuracy by about 0.02, it increased the bias for purebred predictions in emBayesR models. Generally, emBayesR outperformed GBLUP for prediction accuracy when using 50k or pruned HDnGBS genotypes, but the benefits diminished with XT_50k genotypes. Crossbred predictions derived from a joint pure H and J reference were similar in accuracy to crossbred predictions derived from the two separate purebred reference sets and combined proportional to breed composition. However, the latter approach was less biased by 0.13. Most interestingly, using an equalized breed reference instead of an H-dominated reference, on average, reduced the bias of prediction by 0.16–0.19 and increased the accuracy by 0.04 for crossbred and J cows, with a little change in the H accuracy. In conclusion, we observed improved genomic predictions for both crossbreds and purebreds by equalizing breed contributions in a mixed breed reference that included crossbred cows. Furthermore, we demonstrate, that compared to the conventional 50k or high-density panels, our customized set of 50k sequence markers improved or matched the prediction accuracy and reduced bias with both GBLUP and Bayesian models.

Genome-wide association studies for methane emission and ruminal volatile fatty acids using Holstein cattle sequence data

Background

Methane emission by ruminants has contributed considerably to the global warming and understanding the genomic architecture of methane production may help livestock producers to reduce the methane emission from the livestock production system. The goal of our study was to identify genomic regions affecting the predicted methane emission (PME) from volatile fatty acids (VFAs) indicators and VFA traits using imputed whole-genome sequence data in Iranian Holstein cattle.

Results

Based on the significant-association threshold (p < 5 × 10^− 8), 33 single nucleotide polymorphisms (SNPs) were detected for PME per kg milk (n = 2), PME per kg fat (n = 14), and valeric acid (n = 17). Besides, 69 genes were identified for valeric acid (n = 18), PME per kg milk (n = 4) and PME per kg fat (n = 47) that were located within 1 Mb of significant SNPs. Based on the gene ontology (GO) term analysis, six promising candidate genes were significantly clustered in organelle organization (GO:0004984, p = 3.9 × 10^− 2) for valeric acid, and 17 candidate genes significantly clustered in olfactory receptors activity (GO:0004984, p = 4 × 10^− 10) for PME traits. Annotation results revealed 31 quantitative trait loci (QTLs) for milk yield and its components, body weight, and residual feed intake within 1 Mb of significant SNPs.

Conclusions

Our results identified 33 SNPs associated with PME and valeric acid traits, as well as 17 olfactory receptors activity genes for PME traits related to feed intake and preference. Identified SNPs were close to 31 QTLs for milk yield and its components, body weight, and residual feed intake traits. In addition, these traits had high correlations with PME trait. Overall, our findings suggest that marker-assisted and genomic selection could be used to improve the difficult and expensive-to-measure phenotypes such as PME. Moreover, prediction of methane emission by VFA indicators could be useful for increasing the size of reference population required in genome-wide association studies and genomic selection.

Mitochondrial protein gene expression and the oxidative phosphorylation pathway associated with feed efficiency and energy balance in dairy cattle

Feed efficiency and energy balance are important traits underpinning profitability and environmental sustainability in animal production. They are complex traits, and our understanding of their underlying biology is currently limited. One measure of feed efficiency is residual feed intake (RFI), which is the difference between actual and predicted intake. Variation in RFI among individuals is attributable to the metabolic efficiency of energy utilization. High RFI (H_RFI) animals require more energy per unit of weight gain or milk produced compared with low RFI (L_RFI) animals. Energy balance (EB) is a closely related trait calculated very similarly to RFI. Cellular energy metabolism in mitochondria involves mitochondrial protein (MiP) encoded by both nuclear (NuMiP) and mitochondrial (MtMiP) genomes. We hypothesized that MiP genes are differentially expressed (DE) between H_RFI and L_RFI animal groups and similarly between negative and positive EB groups. Our study aimed to characterize MiP gene expression in white blood cells of H_RFI and L_RFI cows using RNA sequencing to identify genes and biological pathways associated with feed efficiency in dairy cattle. We used the top and bottom 14 cows ranked for RFI and EB out of 109 animals as H_RFI and L_RFI, and positive and negative EB groups, respectively. The gene expression counts across all nuclear and mitochondrial genes for animals in each group were used for differential gene expression analyses, weighted gene correlation network analysis, functional enrichment, and identification of hub genes. Out of 244 DE genes between RFI groups, 38 were MiP genes. The DE genes were enriched for the oxidative phosphorylation (OXPHOS) and ribosome pathways. The DE MiP genes were underexpressed in L_RFI (and negative EB) compared with the H_RFI (and positive EB) groups, suggestive of reduced mitochondrial activity in the L_RFI group. None of the MtMiP genes were among the DE MiP genes between the groups, which suggests a non-rate limiting role of MtMiP genes in feed efficiency and warrants further investigation. The role of MiP, particularly the NuMiP and OXPHOS pathways in RFI, was also supported by our gene correlation network analysis and the hub gene identification. We validated the findings in an independent data set. Overall, our study suggested that differences in feed efficiency in dairy cows may be linked to differences in cellular energy demand. This study broadens our knowledge of the biology of feed efficiency in dairy cattle.

Identification of polymorphisms in the oocyte-derived growth differentiation growth factor 9 (GDF9) gene associated with litter size in New Zealand sheep (Ovis aries) breeds

Having the ability to control litter size is important for sheep farmers and breeders worldwide. However, making genetic gain in key livestock traits like reproductive performance needs typically a lot of time, and both the fecundity and fertility traits have a great economic importance. Attention has therefore turned to better understanding the genes that control reproductive performance. Of these genes, research has focussed on the growth differentiation growth factor 9 (GDF9) gene (GDF9). In this study, a PCR-single strand conformation polymorphism (PCR-SSCP) approach was used to investigate variation in this gene in separate groups of purebred Finnish Landrace sheep, Finnish Landrace × Texel-cross sheep and composite sheep of undefined breed background, but based on New Zealand Romney-type genetics. Three GDF9 variants (named A, B and C) were found, and upon DNA sequencing, the nucleotide substitutions c.978A>G, c.994G>A and c.1111G>A were revealed. The frequency of variant A (containing nucleotides c.978A, c.994G and c.1111G) in the Finnish Landrace, Finnish Landrace × Texel-cross and composite sheep was 0.86, 0.78 and 0.76, respectively. In these three sheep groups, the frequency of B (defined by the presence of nucleotides c.978G and c.994A) was 0.01, 0.03 and 0.23 and for C (containing c.1111A) was 0.13, 0.18 and 0.01, respectively. An animal model was used to estimate the additive effect of fertility data for Finnish Landrace × Texel-cross sheep and revealed an association between litter size and the c.1111G>A variation (p = .036), but this was not observed for the Finnish Landrace sheep (p = .27) or the composite sheep (p = .17). When all the sheep were analysed together, the presence of c.1111A was associated (p < .05) with increased litter size, when compared to ewes that had c.1111G. Litter size did not differ between sheep with and without c.994A in all three groups of sheep investigated. This study suggests that c.1111A could be a useful genetic marker for improving fecundity in New Zealand sheep breeds and that it could be introgressed into other breeds, but analysis of more sheep will be required to confirm the associations that have been observed here.

Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds

Background

Thousands of years of natural and artificial selection have resulted in indigenous cattle breeds that are well-adapted to the environmental challenges of their local habitat and thereby are considered as valuable genetic resources. Understanding the genetic background of such adaptation processes can help us design effective breeding objectives to preserve local breeds and improve commercial cattle. To identify regions under putative selection, GGP HD 150 K single nucleotide polymorphism (SNP) arrays were used to genotype 106 individuals representing five Swedish breeds i.e. native to different regions and covering areas with a subarctic cold climate in the north and mountainous west, to those with a continental climate in the more densely populated south regions.

Results

Five statistics were incorporated within a framework, known as de-correlated composite of multiple signals (DCMS) to detect signatures of selection. The obtained p-values were adjusted for multiple testing (FDR < 5%), and significant genomic regions were identified. Annotation of genes in these regions revealed various verified and novel candidate genes that are associated with a diverse range of traits, including e.g. high altitude adaptation and response to hypoxia (DCAF8, PPP1R12A, SLC16A3, UCP2, UCP3, TIGAR), cold acclimation (AQP3, AQP7, HSPB8), body size and stature (PLAG1, KCNA6, NDUFA9, AKAP3, C5H12orf4, RAD51AP1, FGF6, TIGAR, CCND2, CSMD3), resistance to disease and bacterial infection (CHI3L2, GBP6, PPFIBP1, REP15, CYP4F2, TIGD2, PYURF, SLC10A2, FCHSD2, ARHGEF17, RELT, PRDM2, KDM5B), reproduction (PPP1R12A, ZFP36L2, CSPP1), milk yield and components (NPC1L1, NUDCD3, ACSS1, FCHSD2), growth and feed efficiency (TMEM68, TGS1, LYN, XKR4, FOXA2, GBP2, GBP5, FGD6), and polled phenotype (URB1, EVA1C).

Conclusions

We identified genomic regions that may provide background knowledge to understand the mechanisms that are involved in economic traits and adaptation to cold climate in cattle. Incorporating p-values of different statistics in a single DCMS framework may help select and prioritize candidate genes for further analyses.

Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis

Background

Consecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~ 65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (F_ROH), excess of homozygosity (F_HOM), correlation between uniting gametes (F_UNI), and diagonal elements of the genomic relationship matrix (F_GRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies.

Results

In this study, 9102 ROH were identified, with an average number of 21.2 ± 13.1 and 33.2 ± 15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8 ± 120.3 Mb), and in KHZ, 5.96% (149.1 ± 107.7 Mb) of the genome was autozygous. The estimated inbreeding values based on F_HOM, F_UNI and F_GRM were higher in AZ than they were in KHZ, which was in contrast to the F_ROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ≤ 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX).

Conclusion

The calculated F_ROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that F_ROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.

Linkage disequilibrium and within-breed genetic diversity in Iranian Zandi sheep

This research aimed to measure the extent of linkage disequilibrium (LD), effective population size (Ne), and runs of homozygosity (ROHs) in one of the major Iranian sheep breeds (Zandi) using 96 samples genotyped with Illumina Ovine SNP50 BeadChip. The amount of LD (r2) for single-nucleotide polymorphism (SNP) pairs in short distances (10–20 kb) was 0.21±0.25 but rapidly decreased to 0.10±0.16 by increasing the distance between SNP pairs (40–60 kb). The Ne of Zandi sheep in past (approximately 3500 generations ago) and recent (five generations ago) populations was estimated to be 6475 and 122, respectively. The ROH-based inbreeding was 0.023. We found 558 ROH regions, of which 37 % were relatively long (>10 Mb). Compared with the rate of LD reduction in other species (e.g., cattle and pigs), in Zandi, it was reduced more rapidly by increasing the distance between SNP pairs. According to the LD pattern and high genetic diversity of Zandi sheep, we need to use an SNP panel with a higher density than Illumina Ovine SNP50 BeadChip for genomic selection and genome-wide association studies in this breed.

Accuracy of imputation to whole-genome sequence in sheep

Background

The use of whole-genome sequence (WGS) data for genomic prediction and association studies is highly desirable because the causal mutations should be present in the data. The sequencing of 935 sheep from a range of breeds provides the opportunity to impute sheep genotyped with single nucleotide polymorphism (SNP) arrays to WGS. This study evaluated the accuracy of imputation from SNP genotypes to WGS using this reference population of 935 sequenced sheep.

Results

The accuracy of imputation from the Ovine Infinium^® HD BeadChip SNP (~ 500 k) to WGS was assessed for three target breeds: Merino, Poll Dorset and F1 Border Leicester × Merino. Imputation accuracy was highest for the Poll Dorset breed, although there were more Merino individuals in the sequenced reference population than Poll Dorset individuals. In addition, empirical imputation accuracies were higher (by up to 1.7%) when using larger multi-breed reference populations compared to using a smaller single-breed reference population. The mean accuracy of imputation across target breeds using the Minimac3 or the FImpute software was 0.94. The empirical imputation accuracy varied considerably across the genome; six chromosomes carried regions of one or more Mb with a mean imputation accuracy of < 0.7. Imputation accuracy in five variant annotation classes ranged from 0.87 (missense) up to 0.94 (intronic variants), where lower accuracy corresponded to higher proportions of rare alleles. The imputation quality statistic reported from Minimac3 (R²) had a clear positive relationship with the empirical imputation accuracy. Therefore, by first discarding imputed variants with an R² below 0.4, the mean empirical accuracy across target breeds increased to 0.97. Although accuracy of genomic prediction was less affected by filtering on R² in a multi-breed population of sheep with imputed WGS, the genomic heritability clearly tended to be lower when using variants with an R² ≤ 0.4.

Conclusions

The mean imputation accuracy was high for all target breeds and was increased by combining smaller breed sets into a multi-breed reference. We found that the Minimac3 software imputation quality statistic (R²) was a useful indicator of empirical imputation accuracy, enabling removal of very poorly imputed variants before downstream analyses.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0443-5) contains supplementary material, which is available to authorized users.

Comparing allele specific expression and local expression quantitative trait loci and the influence of gene expression on complex trait variation in cattle

BACKGROUND

The mutations changing the expression level of a gene, or expression quantitative trait loci (eQTL), can be identified by testing the association between genetic variants and gene expression in multiple individuals (eQTL mapping), or by comparing the expression of the alleles in a heterozygous individual (allele specific expression or ASE analysis). The aims of the study were to find and compare ASE and local eQTL in 4 bovine RNA-sequencing (RNA-Seq) datasets, validate them in an independent ASE study and investigate if they are associated with complex trait variation.

RESULTS

We present a novel method for distinguishing between ASE driven by polymorphisms in cis and parent of origin effects. We found that single nucleotide polymorphisms (SNPs) driving ASE are also often local eQTL and therefore presumably cis eQTL. These SNPs often, but not always, affect gene expression in multiple tissues and, when they do, the allele increasing expression is usually the same. However, there were systematic differences between ASE and local eQTL and between tissues and breeds. We also found that SNPs significantly associated with gene expression (p < 0.001) were likely to influence some complex traits (p < 0.001), which means that some mutations influence variation in complex traits by changing the expression level of genes.

CONCLUSION

We conclude that ASE detects phenomenon that overlap with local eQTL, but there are also systematic differences between the SNPs discovered by the two methods. Some mutations influencing complex traits are actually eQTL and can be discovered using RNA-Seq including eQTL in the genes CAST, CAPN1, LCORL and LEPROTL1.

Genome variants associated with RNA splicing variations in bovine are extensively shared between tissues

Background

Mammalian phenotypes are shaped by numerous genome variants, many of which may regulate gene transcription or RNA splicing. To identify variants with regulatory functions in cattle, an important economic and model species, we used sequence variants to map a type of expression quantitative trait loci (expression QTLs) that are associated with variations in the RNA splicing, i.e., sQTLs. To further the understanding of regulatory variants, sQTLs were compare with other two types of expression QTLs, 1) variants associated with variations in gene expression, i.e., geQTLs and 2) variants associated with variations in exon expression, i.e., eeQTLs, in different tissues.

Results

Using whole genome and RNA sequence data from four tissues of over 200 cattle, sQTLs identified using exon inclusion ratios were verified by matching their effects on adjacent intron excision ratios. sQTLs contained the highest percentage of variants that are within the intronic region of genes and contained the lowest percentage of variants that are within intergenic regions, compared to eeQTLs and geQTLs. Many geQTLs and sQTLs are also detected as eeQTLs. Many expression QTLs, including sQTLs, were significant in all four tissues and had a similar effect in each tissue. To verify such expression QTL sharing between tissues, variants surrounding (±1 Mb) the exon or gene were used to build local genomic relationship matrices (LGRM) and estimated genetic correlations between tissues. For many exons, the splicing and expression level was determined by the same cis additive genetic variance in different tissues. Thus, an effective but simple-to-implement meta-analysis combining information from three tissues is introduced to increase power to detect and validate sQTLs. sQTLs and eeQTLs together were more enriched for variants associated with cattle complex traits, compared to geQTLs. Several putative causal mutations were identified, including an sQTL at Chr6:87392580 within the 5th exon of kappa casein (CSN3) associated with milk production traits.

Conclusions

Using novel analytical approaches, we report the first identification of numerous bovine sQTLs which are extensively shared between multiple tissue types. The significant overlaps between bovine sQTLs and complex traits QTL highlight the contribution of regulatory mutations to phenotypic variations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4902-8) contains supplementary material, which is available to authorized users.

Gene expression analysis of blood, liver, and muscle in cattle divergently selected for high and low residual feed intake

Improving feed efficiency in cattle is important because it increases profitability by reducing costs, and it also shrinks the environmental footprint of cattle production by decreasing manure and greenhouse gas emissions. Residual feed intake (RFI) is 1 measurement of feed efficiency and is the difference between actual and predicted feed intake. Residual feed intake is a complex trait with moderate heritability, but the genes and biological processes associated with its variation still need to be found. We explored the variation in expression of genes using RNA sequencing to find genes whose expression was associated with RFI and then investigated the pathways that are enriched for these genes. In this study, we used samples from growing Angus bulls (muscle and liver tissues) and lactating Holstein cows (liver tissue and white blood cells) divergently selected for low and high RFI. Within each breed-tissue combination, the correlation between the expression of genes and RFI phenotypes, as well as GEBV, was calculated to determine the genes whose expression was correlated with RFI. There were 16,039 genes expressed in more than 25% of samples in 1 or more tissues. The expression of 6,143 genes was significantly associated with RFI phenotypes, and expression of 2,343 genes was significantly associated with GEBV for RFI ( < 0.05) in at least 1 tissue. The genes whose expression was correlated with RFI phenotype (or GEBV) within each breed-tissue combination were enriched for 158 (78) biological processes (Fisher Exact Statistics for gene-enrichment analysis, EASE score < 0.1) and associated with 13 (13) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways ( < 0.05 and fold enrichment > 2). These biological processes were related to regulation of transcription, translation, energy generation, cell cycling, apoptosis, and proteolysis. However, the direction of the correlation between RFI and gene expression in some cases reversed between tissues. For instance, low levels of proteolysis in muscle were associated with high efficiency in growing bulls, but high levels of proteolysis in white blood cells were associated with efficiency of milk production in lactating cows.

Extensive variation between tissues in allele specific expression in an outbred mammal

BACKGROUND

Allele specific gene expression (ASE), with the paternal allele more expressed than the maternal allele or vice versa, appears to be a common phenomenon in humans and mice. In other species the extent of ASE is unknown, and even in humans and mice there are several outstanding questions. These include; to what extent is ASE tissue specific? how often does the direction of allele expression imbalance reverse between tissues? how often is only one of the two alleles expressed? is there a genome wide bias towards expression of the paternal or maternal allele; and finally do genes that are nearby on a chromosome share the same direction of ASE? Here we use gene expression data (RNASeq) from 18 tissues from a single cow to investigate each of these questions in turn, and then validate some of these findings in two tissues from 20 cows.

RESULTS

Between 40 and 100 million sequence reads were generated per tissue across three replicate samples for each of the eighteen tissues from the single cow (the discovery dataset). A bovine gene expression atlas was created (the first from RNASeq data), and differentially expressed genes in each tissue were identified. To analyse ASE, we had access to unambiguously phased genotypes for all heterozygous variants in the cow's whole genome sequence, where these variants were homozygous in the whole genome sequence of her sire, and as a result we were able to map reads to parental genomes, to determine SNP and genes showing ASE in each tissue. In total 25,251 heterozygous SNP within 7985 genes were tested for ASE in at least one tissue. ASE was pervasive, 89 % of genes tested had significant ASE in at least one tissue. This large proportion of genes displaying ASE was confirmed in the two tissues in a validation dataset. For individual tissues the proportion of genes showing significant ASE varied from as low as 8-16 % of those tested in thymus to as high as 71-82 % of those tested in lung. There were a number of cases where the direction of allele expression imbalance reversed between tissues. For example the gene SPTY2D1 showed almost complete paternal allele expression in kidney and thymus, and almost complete maternal allele expression in the brain caudal lobe and brain cerebellum. Mono allelic expression (MAE) was common, with 1349 of 4856 genes (28 %) tested with more than one heterozygous SNP showing MAE. Across all tissues, 54.17 % of all genes with ASE favoured the paternal allele. Genes that are closely linked on the chromosome were more likely to show higher expression of the same allele (paternal or maternal) than expected by chance. We identified several long runs of neighbouring genes that showed either paternal or maternal ASE, one example was five adjacent genes (GIMAP8, GIMAP7 copy1, GIMAP4, GIMAP7 copy 2 and GIMAP5) that showed almost exclusive paternal expression in brain caudal lobe.

CONCLUSIONS

Investigating the extent of ASE across 18 bovine tissues in one cow and two tissues in 20 cows demonstrated 1) ASE is pervasive in cattle, 2) the ASE is often MAE but ranges from MAE to slight overexpression of the major allele, 3) the ASE is most often tissue specific and that more than half the time displays divergent allele specific expression patterns across tissues, 4) across all genes there is a slight bias towards expression of the paternal allele and 5) genes expressing the same parental allele are clustered together more than expected by chance, and there are several runs of large numbers of genes expressing the same parental allele.

Estimation of genomic breeding values for residual feed intake in a multibreed cattle population

Residual feed intake (RFI) is a measure of the efficiency of animals in feed utilization. The accuracies of GEBV for RFI could be improved by increasing the size of the reference population. Combining RFI records of different breeds is a way to do that. The aims of this study were to 1) develop a method for calculating GEBV in a multibreed population and 2) improve the accuracies of GEBV by using SNP associated with RFI. An alternative method for calculating accuracies of GEBV using genomic BLUP (GBLUP) equations is also described and compared to cross-validation tests. The dataset included RFI records and 606,096 SNP genotypes for 5,614 Bos taurus animals including 842 Holstein heifers and 2,009 Australian and 2,763 Canadian beef cattle. A range of models were tested for combining genotype and phenotype information from different breeds and the best model included an overall effect of each SNP, an effect of each SNP specific to a breed, and a small residual polygenic effect defined by the pedigree. In this model, the Holsteins and some Angus cattle were combined into 1 "breed class" because they were the only cattle measured for RFI at an early age (6-9 mo of age) and were fed a similar diet. The average empirical accuracy (0.31), estimated by calculating the correlation between GEBV and actual phenotypes divided by the square root of estimated heritability in 5-fold cross-validation tests, was near to that expected using the GBLUP equations (0.34). The average empirical and expected accuracies were 0.30 and 0.31, respectively, when the GEBV were estimated for each breed separately. Therefore, the across-breed reference population increased the accuracy of GEBV slightly, although the gain was greater for breeds with smaller number of individuals in the reference population (0.08 in Murray Grey and 0.11 in Hereford for empirical accuracy). In a second approach, SNP that were significantly (P < 0.001) associated with RFI in the beef cattle genomewide association studies were used to create an auxiliary genomic relationship matrix for estimating GEBV in Holstein heifers. The empirical (and expected) accuracy of GEBV within Holsteins increased from 0.33 (0.35) to 0.39 (0.36) and improved even more to 0.43 (0.50) when using a multibreed reference population. Therefore, a multibreed reference population is a useful resource to find SNP with a greater than average association with RFI in 1 breed and use them to estimate GEBV in another breed.

The Italian National Registry of severe acquired brain injury: epidemiological, clinical and functional data of 1469 patients

BACKGROUND

The lack of knowledge about epidemiological and clinical data of patients with Acquired Brain Injury (ABI) admitted to Rehabilitation facilities in Italy led to the creation, in June 2008, of a data collection on-line registry.

AIM

To collect epidemiological and clinical data and to evaluate functional outcome of patients with severe traumatic and non-traumatic ABI admitted to Rehabilitation facilities in Italy between June 2008 and December 2011 and to compare data of patients with ABI of different aetiologies.

DESIGN

Observational retrospective study.

SETTING

The study involved 29 Italian Rehabilitation facilities.

POPULATION

The study enrolled 1469 patients with severe traumatic (TBI) and non-traumatic ABI (NTBI).

METHODS

Data collected included demographic (number of patients with TBI and NTBI, gender, age) and clinical characteristics (provenience, number of days elapsed between onset and rehabilitation admission, rehabilitation length of stay, discharge destination, death and vegetative state diagnosis, presence of percutaneous endoscopic gastrostomy, tracheostomy, pressure sores and paraosteoarthropathies). Functional outcome was evaluated using the Disability Rating Scale.

RESULTS

Of the whole population studied, 44.31% and 55.69% patients had suffered a TBI and a NTBI, respectively. In the NTBI group 40.09% had a cerebrovascular injury, 12.04% an anoxic brain damage, 3.6% had a brain injury of other causes. The mean age was 43.67 and 56.68 for subjects with TBI and NTBI, respectively. Patients with TBI showed a lower onset-admission interval (OAI), compared with NTBI group; no difference in rehabilitation length of stay (LOS) was recorded between groups. Patients with TBI presented a lower DRS score at admission and discharge and returned home more frequently than NTBI group.

CONCLUSIONS

The creation of a National registry allows the collection of data about patients with ABI in order to study the clinical course, the functional outcome and to establish a basis for comparison with other data sources. Clinical Rehabilitation Impact. Data collection could be useful in the evaluation and planning of rehabilitation pathways, and to assess the allocation of healthcare and rehabilitative resources.

Use of early indicators in rehabilitation process to predict functional outcomes in subjects with acquired brain injury

BACKGROUND

The use of Evaluation Scales in ABI is necessary for measure of outcome, but not always they are used as predictor factors for rehabilitation processes and organization.

AIM

The aim of this study was to evaluate the effectiveness and efficiency of an inpatient rehabilitation program for patients with traumatic brain injury through the use of selected indicators and to identify predictive factors for functional outcome.

DESIGN

This was a retrospective database analysis.

SETTING

Patients admitted to an Intensive Rehabilitation Unit as inpatient (Sacro Cuore-Don Calabria Hospital, Negrar-Verona).

POPULATION

The population included patients with traumatic brain injury.

METHODS

The study enrolled 175 patients admitted to an Intensive Rehabilitation Unit between 2004 and 2007. Data collected included demographic characteristics, first 24-hours worst GCS, length of acute and rehabilitative stay at admission and discharge FIM, DRS, LCF and GOS.

RESULTS

There was a statistically significant recovery over the course of admission for all assessment tools (P<0.000). When patients were subdivided on the basis of admission DRS categories a linear correlation among variables could be observed, with most disabled patients showing the longest acute and rehabilitation stays and the lowest functional gains. Within each DRS category age appeared to affect improvement (P<0.005) while final outcome was influenced not only by age but also by initial functional status (P<0.000) and time from injury to admission (P<0.004). CONCLUSIONS AND CLINICAL REHABILITATION IMPACT: Systematic data collection in intensive rehabilitation is of great importance to monitor recovery and plan appropriate programs on the basis of admission functional status.