Search

The objective of this study was to investigate genetic variation and genotype by environment (G × E) interactions for fertility (including age at first calving and calving interval), somatic cell score (SCS), and milk production traits for Iranian Holsteins. Different environments were defined based on the climatic zones (cold, semi-cold, and moderate) and considering the herd location. Data were collected between 2003 and 2018 by the National Animal Breeding Center of Iran (Karaj). Variance and covariance components and genetic correlations were estimated using 2 different models, which were analyzed using Bayesian methods. For both models, performance of traits in each climate were considered as different traits. Fertility traits were analyzed using a trivariate model. Furthermore, SCS and production traits were analyzed using trivariate random regression models (records in different climate zones considered as different traits). For the fertility traits, the largest estimates of heritability were observed in cold climate. Fertility performance was always better in cold environment. Genetic correlations between climatic zones ranged from 0.85 to 0.94. For daily measurements of SCS and production traits, heritability ranged from 0.031 to 0.037 and 0.069 to 0.209, respectively. Genetic variances were the highest in the semi-cold and moderate climates for the SCS and production traits, respectively. Furthermore, across the studied climates, 305-d genetic correlation ranged from 0.756 to 0.884 for SCS and from 0.925 to 0.957 for the production traits. The structure of genetic correlation within each climate indicated a negative correlation between early and late lactation for SCS, especially in the cold climate and for milk production in the moderate climate. For fat percentage, in all climatic zones, the lowest genetic correlations were observed between early and mid-lactation. In addition, for protein production in the cold climate, a negative correlation was observed between early and late lactation. Results indicated heterogeneous variance components for all the studied traits across various climatic zones. Estimated genetic correlations for SCS revealed that the genetic expression of animals may vary by climatic zone. Results indicated the existence of G × E interaction due to the climatic condition, only for SCS. Therefore, in Iranian Holsteins, the effect of G × E interactions should not be neglected, especially for SCS, as different sires might be optimal for use in different climatic zones.

In this study, computationally efficient methods to approximate the reliabilities of genomic estimated breeding values (GEBV) in a single-step genomic prediction model including a residual polygenic (RPG) effect are described. In order to calculate the reliabilities of the genotyped animals, a single nucleotide polymorphism best linear unbiased prediction (SNPBLUP) or a genomic BLUP (GBLUP), was used, where two alternatives to account for the RPG effect were tested. In the direct approach, the genomic model included the RPG effect, while in the blended method, it did not but an index was used to weight the genomic and pedigree-based BLUP (PBLUP) reliabilities. In order to calculate the single-step GBLUP reliabilities for the breeding values for the non-genotyped animals, a simplified weighted-PBLUP model that included a general mean and additive genetic effects with weights accounting for the non-genomic and genomic information was used. We compared five schemes for the weights. Two datasets, i.e., a small (Data 1) one and a large (Data 2) one were used.For the genotyped animals in Data 1, correlations between approximate reliabilities using the blended method and exact reliabilities ranged from 0.993 to 0.996 across three lactations. The slopes observed by regressing the reliabilities of GEBV from the exact method on those from the blended method were 1.0 for all three lactations. For Data 2, the correlations and slopes ranged, respectively, from 0.980 to 0.986 and from 0.91 to 0.96, and for the non-genotyped animals in Data 1, they ranged, respectively, from 0.987 to 0.994 and from 0.987 to 1, which indicate that the approximations were in line with the exact results. The best approach achieved correlations of 0.992 to 0.994 across lactations.Our results demonstrate that the approximated reliabilities calculated using our proposed approach are in good agreement with the exact reliabilities. The blended method for the genotyped animals is computationally more feasible than the direct method when RPG effects are included, particularly for large-scale datasets. The approach can serve as an effective strategy to estimate the reliabilities of GEBV in large-scale single-step genomic predictions.

Fat-tailed sheep have a unique characteristic of depositing fat in their tails. In the present study, we conducted genome-wide association studies (GWAS) on traits related to tail fat deposition and body size in the Hulun Buir sheep. A total number of 300 individuals belonging to two fat-tailed lines of the Hulun Buir sheep breed genotyped with the Ovine Infinium HD SNP BeadChip were included in the current study. Two mixed models, one for continuous and one for binary phenotypic traits, were employed to analyse ten traits, that is, body length (BL), body height (BH), chest girth (CG), tail length (TL), tail width (TW), tail circumference (TC), carcass weight (CW), tail fat weight (TF), ratio of CW to TF (RCT) and tail type (TT). We identified 7, 6, 7, 2, 10 and 1 SNPs significantly associated with traits TF, CW, RCT, TW, TT and CG, respectively. Their associated genomic regions harboured 42 positional candidate genes. Out of them, 13 candidate genes including SMURF2, FBF1, DTNBP1, SETD7 and RBM11 have been associated with fat metabolism in sheep. The RBM11 gene has already been identified in a previous study on signatures of selection in this specific sheep population. Two more genes, that is, SMARCA5 and GAB1 were associated with body size in sheep. The present study has identified candidate genes that might be implicated in tail fat deposition and body size in sheep.

This work focuses on the effects of variable amount of genomic information in the Bayesian estimation of unknown variance components associated with single-step genomic prediction. We propose a quantitative criterion for the amount of genomic information included in the model and use it to study the relative effect of genomic data on efficiency of sampling from the posterior distribution of parameters of the single-step model when conducting a Bayesian analysis with estimating unknown variances. The rate of change of estimated variances was dependent on the amount of genomic information involved in the analysis, but did not depend on the Gibbs updating schemes applied for sampling realizations of the posterior distribution. Simulation revealed a gradual deterioration of convergence rates for the locations parameters when new genomic data were gradually added into the analysis. In contrast, the convergence of variance components showed continuous improvement under the same conditions. The sampling efficiency increased proportionally to the amount of genomic information. In addition, an optimal amount of genomic information in variance-covariance matrix that guaranty the most (computationally) efficient analysis was found to correspond a proportion of animals genotyped ***0.8. The proposed criterion yield a characterization of expected performance of the Gibbs sampler if the analysis is subject to adjustment of the amount of genomic data and can be used to guide researchers on how large a proportion of animals should be genotyped in order to attain an efficient analysis.

Records on groups of individuals could be valuable for predicting breeding values when a trait is difficult or costly to measure on single individuals, such as feed intake and egg production. Adding genomic information has shown improvement in the accuracy of genetic evaluation of quantitative traits with individual records. Here, we investigated the value of genomic information for traits with group records. Besides, we investigated the improvement in accuracy of genetic evaluation for group-recorded traits when including information on a correlated trait with individual records. The study was based on a simulated pig population, including three scenarios of group structure and size. The results showed that both the genomic information and a correlated trait increased the accuracy of estimated breeding values (EBVs) for traits with group records. The accuracies of EBV obtained from group records with a size 24 were much lower than those with a size 12. Random assignment of animals to pens led to lower accuracy due to the weaker relationship between individuals within each group. It suggests that group records are valuable for genetic evaluation of a trait that is difficult to record on individuals, and the accuracy of genetic evaluation can be considerably increased using genomic information. Moreover, the genetic evaluation for a trait with group records can be greatly improved using a bivariate model, including correlated traits that are recorded individually. For efficient use of group records in genetic evaluation, relatively small group size and close relationships between individuals within one group are recommended.

Random regression test-day model has become the most commonly adopted model for routine genetic evaluations for different dairy populations, which allows accurately accounting for genetic and environmental effects at different periods during lactation. The objective of this study was to explore appropriate random regression test-day model for genetic evaluation of milk yield in Chinese Holstein population. Data included 419,567 test-day records from 54,417 cows in the first lactation. Variance components and breeding values were estimated using random regression test-day model with different order (first order to fifth order) of Legendre polynomials, and accounted for homogeneous or heterogeneous residual variance across the lactation. The goodness of fit of the models was evaluated by total residual variance (TRV) and − 2logL. Further, the predictive ability of the models was assessed by Spearman’s rank correlation between estimated breeding values for 305d milk yield (EBV305) from the full data set and reduced data set in which the records from the last calving year were masked. The results showed that random regression models using third order Legendre polynomials (LP3) with heterogeneous residual variance achieved the lower TRV and − 2logL value and the highest correlation for EBV305 between full data and reduced data. Heritability estimated by this model was 0.250 for 305d milk yield and ranged from 0.163 to 0.304 for test-day milk yield. We suggest random regression model with Legendre polynomial of order 3 and accounting for heterogeneous residual variances could be an appropriate model to be used for genetic evaluation of milk yield for Chinese Holstein population.

It is of practical importance to ensure the data quality from a milk-recording system before use for genetic evaluation. A procedure was developed for detection of multivariate outliers based on an approximation for Mahalanobis distance and was implemented in the Nordic Holstein and Red population. The general target of this procedure is based on the Nordic Cattle Genetic Evaluation yield model, which is a 9-trait model for milk, protein, and fat in the first 3 lactations. The procedure is based on the phenotypic correlation structure as a function of days in milk (DIM) and on computation of trait means and standard deviations within a production year, lactation, and DIM. For each record in the data, a Mahalanobis distance value was computed based on the trait mean and the covariance matrix for the actual production year, lactation, and DIM. A set of cutoff values, ranging from 10 to 100 with steps of 10, for discarding multivariate outliers was investigated. Prediction accuracy was calculated as the Pearson correlations between estimated breeding values predicted by full data set and estimated breeding values predicted by reduced data set for cows without records in the reduced data set and with 1 or more records deleted due to the editing rules on Mahalanobis distance. The results showed that, averaged over all scenarios, gains of 0.005 to 0.048 on prediction accuracy have been obtained by deleting the multivariate outliers. The improvements were more profound for progeny of young bulls compared with progeny of proven bulls. It is easy to implement this multivariate outlier-detection procedure in the routine genetic evaluation for different dairy cattle breeds; however, an optimal cutoff value for Mahalanobis distance needs to be defined to achieve an acceptable compromise between genetic evaluation accuracy and data deletion.

Inbreeding, which has several causes including genetic drift, population bottlenecks, mating of close relatives and selection, can leave tracts of runs of homozygosity (ROH) along genomes. Recently, decreasing performance on reproductive traits, which might have resulted from inbreeding, has been observed in Chinese pigs. In this study, 830 individuals from Western and Chinese pig breeds were genotyped using the reduced-representation sequencing method. After imputation and quality control, 60 850 high-confidence SNPs were retained for ROH detection. A simulation was performed to explore the reliability of ROH detection with imputed data. Different ROH-related variables were compared between imputed and non-missing genotypes used in ROH detection. Furthermore, ROH islands were evaluated and annotated to find genes influenced by inbreeding in these pigs. The simulation results showed that imputed data with 0.7 as the average missing genotype rate and three heterozygotes allowed in a sliding window have comparable ROH detected compared with data with no missing genotypes. Compared with Western pig breeds, Chinese pigs had more autosomes covered by ROH longer than 5 Mb, indicating higher inbreeding in Chinese pigs in recent times. Genes related to reproduction, immunity, meat quality and adaptability in Chinese pigs and several genes related to growth speed and immunity in Western pigs were observed in short ROH islands. The reproduction-related gene PRM1 was found to be located in the most frequent long ROH island in Chinese pigs, which might explain the decreasing fertility in Chinese pig breeds.

Single-step genomic prediction models utilizing both genotyped and nongenotyped animals are likely to become the prevailing tool in genetic evaluations of livestock. Various single-step prediction models have been proposed, based either on estimation of individual marker effects or on direct prediction via a genomic relationship matrix. In this study, a classical pedigree-based animal model, a regular single-step genomic BLUP (ssGBLUP) model, algorithm for proven and young (APY) with 2 strategies for choosing core animals, and a single-step Bayesian regression (ssBR) model were compared for 305-d production traits (milk, fat, protein) in the Finnish red dairy cattle population. A residual polygenic effect with 10% of total genetic variance was included in the single-step models to reduce inflation of genomic predictions. Validation reliability was calculated as the squared Pearson correlation coefficient between genomically enhanced breeding value (GEBV) and yield deviation for masked records for 2,056 validation cows from the last year in the data set investigated. The results showed that gains of 0.02 to 0.04 on validation reliability were achieved by using single-step methods compared with the classical animal model. The regular ssGBLUP model and ssBR model with an extra polygenic effect yielded the same results. The APY methods yielded similar reliabilities as the regular ssGBLUP and ssBR. Exact prediction error variance of GEBV could be obtained by ssBR to avoid any approximation methods used for ssGBLUP when inversion left-hand side of mixed model equations is computationally infeasible for large data sets.

Including genotyped females in a reference population (RP) is an obvious way to increase the RP in genomic selection, especially for dairy breeds of limited population size. However, the incorporation of these females must be conducted cautiously because of the potential preferential treatment of the genotyped cows and lower reliabilities of phenotypes compared with the proven pseudo-phenotypes of bulls. Breeding organizations in Denmark, Finland, and Sweden have implemented a female-genotyping project with the possibility of genotyping entire herds using the low-density (LD) chip. In the present study, 5 scenarios for building an RP were investigated in the Nordic Jersey population: (1) bulls only, (2) bulls with females from the LD project, (3) bulls with females from the LD project plus non-LD project females genotyped before their first calving, (4) bulls with females from the LD project plus non-LD project females genotyped after their first calving, and (5) bulls with all genotyped females. The genomically enhanced breeding value (GEBV) was predicted for 8 traits in the Nordic total merit index through a genomic BLUP model using deregressed proof (DRP) as the response variable in all scenarios. In addition, (daughter) yield deviation and raw phenotypic data were studied as response variables for comparison with the DRP, using stature as a model trait. The validation population was formed using a cut-off birth year of 2005 based on the genotyped Nordic Jersey bulls with DRP. The average increment in reliability of the GEBV across the 8 traits investigated was 1.9 to 4.5 percentage points compared with using only bulls in the RP (scenario 1). The addition of all the genotyped females to the RP resulted in the highest gain in reliability (scenario 5), followed by scenario 3, scenario 2, and scenario 4. All scenarios led to inflated GEBV because the regression coefficients are less than 1. However, scenario 2 and scenario 3 led to less bias of genomic predictions than scenario 5, with regression coefficients showing less deviation from scenario 1. For the study on stature, the daughter yield deviation/daughter yield deviation performed slightly better than the DRP as the response variable in the genomic BLUP (GBLUP) model. Therefore, adding unselected females in the RP could significantly improve the reliabilities and tended to reduce the prediction bias compared with adding selectively genotyped females. Although the DRP has performed robustly so far, the use of raw data is recommended with a single-step model as an optimal solution for future genomic evaluations.

Background To compare the efficacy of in situ replacement with rifampicin-soaked silver-coated polyester (RSSCP) to that of expanded polytetrafluoroethylene (ePTFE) graft replacements in a porcine model for early aortic prosthetic vascular graft infection (PVGI). Material and methods Sixty pigs received polyester or silver-coated grafts with an 8 mm diameter implanted end-to-end in the infrarenal aorta, and the grafts were inoculated with approximately 10 6 Staphylococcus aureus . All developed S. aureus PVGI. Two weeks later, the 52 surviving pigs were randomised to undergo in situ graft replacement with ePTFE or RSSCP grafts followed by oral administration of 300 mg rifampicin and 750 mg ciprofloxacin twice a day, postoperatively. After three weeks, all pigs were sacrificed. In situ perigraft swabs and graft material were analysed for S. aureus quantitatively. Results Only one out of 25 RSSCP grafts were infected with S. aureus , whereas 15 of 27 ePTFE grafts were infected after 3 weeks (OR = 0.022, 95% CI: 0.002, 0.219, P = 0.001). Conclusions In situ replacement with RSSCP grafts and oral rifampicin plus ciprofloxacin is more efficiency in eradicating S. aureus PVGI than ePTFE grafts treated with same oral antibiotics in a porcine aortic PVGI model.

Objectives To evaluate the efficacy of a silver-coated vascular polyester graft in the prevention of graft infection after inoculation with Staphylococcus aureus in a porcine model. Material and methods Eighty-four pigs were randomly selected 1:1 to receive a silver-coated or non-silver-coated 8-mm-wide polyester graft implanted end-to-end in the infrarenal aorta. At the end of implantation, 106 colony forming units (CFUs) S. aureus in 0.3 ml suspension were inoculated directly on the graft surface. Blood samples assayed for white blood corpuscles (WBCs) and C-reactive protein (CRP) were taken before implantation and on the postoperative days 2, 5, 7, 11 and 14. Two weeks after implantation, the perigraft swabs were analysed for S. aureus or contaminants. CFUs of S. aureus were quantified and logarithmised. Student's t-tests, repeated measurement analysis of variance (ANOVA) and chi-square test were employed to compare the two grafts. Results All pigs developed graft infection. There were no statistically significant differences between the silver-coated and non-silver-coated grafts in the quantity of S. aureus, macroscopic signs of infection and postoperative changes in the temperature, WBC and CRP. Conclusions Silver–coated polyester grafts failed to prevent graft infections after inoculation with 106S. aureus in 0.3 ml suspension in a porcine model.

A procedure was developed for detection of multivariate outliers based on an approximation for Mahanalobis Distance (MD) and was implemented in the Nordic Jersey population. Evaluations are carried out by Nordic Cattle Genetic Evaluation (NAV), who uses a 9 trait model for milk, protein and fat in the first 3 lactations. It is based on the phenotypic correlation structure as a function of days in milk (DIM) and on computation of trait means and standard deviations (SDs) within classes of production year (PY), lactation and days in milk (DIM). For each record in the data, MD is computed based on trait means and co-variance matrix for the actual PY, lactation and DIM. Accuracy of EBV’s is improved for animals having extreme outlier record(s) deleted compared to EBV’s based on data not filtered for MD.

Background A single-step blending approach allows genomic prediction using information of genotyped and non-genotyped animals simultaneously. However, the combined relationship matrix in a single-step method may need to be adjusted because marker-based and pedigree-based relationship matrices may not be on the same scale. The same may apply when a GBLUP model includes both genomic breeding values and residual polygenic effects. The objective of this study was to compare single-step blending methods and GBLUP methods with and without adjustment of the genomic relationship matrix for genomic prediction of 16 traits in the Nordic Holstein population. Methods The data consisted of de-regressed proofs (DRP) for 5 214 genotyped and 9 374 non-genotyped bulls. The bulls were divided into a training and a validation population by birth date, October 1, 2001. Five approaches for genomic prediction were used: 1) a simple GBLUP method, 2) a GBLUP method with a polygenic effect, 3) an adjusted GBLUP method with a polygenic effect, 4) a single-step blending method, and 5) an adjusted single-step blending method. In the adjusted GBLUP and single-step methods, the genomic relationship matrix was adjusted for the difference of scale between the genomic and the pedigree relationship matrices. A set of weights on the pedigree relationship matrix (ranging from 0.05 to 0.40) was used to build the combined relationship matrix in the single-step blending method and the GBLUP method with a polygenetic effect. Results Averaged over the 16 traits, reliabilities of genomic breeding values predicted using the GBLUP method with a polygenic effect (relative weight of 0.20) were 0.3% higher than reliabilities from the simple GBLUP method (without a polygenic effect). The adjusted single-step blending and original single-step blending methods (relative weight of 0.20) had average reliabilities that were 2.1% and 1.8% higher than the simple GBLUP method, respectively. In addition, the GBLUP method with a polygenic effect led to less bias of genomic predictions than the simple GBLUP method, and both single-step blending methods yielded less bias of predictions than all GBLUP methods. Conclusions The single-step blending method is an appealing approach for practical genomic prediction in dairy cattle. Genomic prediction from the single-step blending method can be improved by adjusting the scale of the genomic relationship matrix.

Objectives To study the feasibility and efficacy of experimental laparoscopy in the diagnosis of aortic graft infection in pigs. Material and methods Eight pigs had an aortic tube graft implanted and inoculated with either 5 × 10 4 or 10 6 CFU of Staphylococcus aureus ATCC 29213. Laparoscopy was performed after a median of 20 days with debridement and sampling for bacterial culture. Thereafter, the grafts were locally soaked in rifampicin and postoperatively, the pigs received rifampicin and ciprofloxacin orally for two weeks and were then sacrificed. Results All pigs developed graft infection. One pig died from severe clostridial septicaemia before laparoscopy could be performed. The remaining pigs had all samples for bacterial culture taken by laparoscopy from the inflamed tissue. The temperature dropped significantly after laparoscopy, and no macroscopic signs of infection presented at autopsy. However, only culture from one pig was without S. aureus at autopsy. Conclusions Laparoscopy is a potential diagnostic tool for aortic graft infection and also affords the opportunity to carry out bacteriological sampling and local antibiotic treatment. The efficacy of laparoscopic treatment needs further evaluation.

This study compared genomic predictions based on imputed high-density markers (~777,000) in the Nordic Holstein population using a genomic BLUP (GBLUP) model, 4 Bayesian exponential power models with different shape parameters (0.3, 0.5, 0.8, and 1.0) for the exponential power distribution, and a Bayesian mixture model (a mixture of 4 normal distributions). Direct genomic values (DGV) were estimated for milk yield, fat yield, protein yield, fertility, and mastitis, using deregressed proofs (DRP) as response variable. The validation animals were split into 4 groups according to their genetic relationship with the training population. Groupsmgs had both the sire and the maternal grandsire (MGS), Groupsire only had the sire, Groupmgs only had the MGS, and Groupnon had neither the sire nor the MGS in the training population. Reliability of DGV was measured as the squared correlation between DGV and DRP divided by the reliability of DRP for the bulls in validation data set. Unbiasedness of DGV was measured as the regression of DRP on DGV. The results indicated that DGV were more accurate and less biased for animals that were more related to the training population. In general, the Bayesian mixture model and the exponential power model with shape parameter of 0.30 led to higher reliability of DGV than did the other models. The differences between reliabilities of DGV from the Bayesian models and the GBLUP model were statistically significant for some traits. We observed a tendency that the superiority of the Bayesian models over the GBLUP model was more profound for the groups having weaker relationships with training population. Averaged over the 5 traits, the Bayesian mixture model improved the reliability of DGV by 2.0 percentage points for Groupsmgs, 2.7 percentage points for Groupsire, 3.3 percentage points for Groupmgs, and 4.3 percentage points for Groupnon compared with GBLUP. The results showed that a Bayesian model with intense shrinkage of the explanatory variable, such as the Bayesian mixture model and the Bayesian exponential power model with shape parameter of 0.30, can improve genomic predictions using high-density markers.