Genome-enabled prediction of indicator traits of resistance to gastrointestinal nematodes in sheep using parametric models and artificial neural networks.

This study aimed to assess the predictive ability of parametric models and artificial neural network method for genomic prediction of the following indicator traits of resistance to gastrointestinal nematodes in Santa Inês sheep: packed cell volume (PCV), fecal egg count (FEC), and Famacha© method (FAM). After quality control, the number of genotyped animals was 551 (PCV), 548 (FEC), and 565 (FAM), and 41,676 SNP. The average prediction accuracy (ACC) calculated by Pearson correlation between observed and predicted values and mean squared errors (MSE) were obtained using genomic best unbiased linear predictor (GBLUP), BayesA, BayesB, Bayesian least absolute shrinkage and selection operator (BLASSO), and Bayesian regularized artificial neural network (three and four hidden neurons, BRANN_3 and BRANN_4, respectively) in a 5-fold cross-validation technique. The average ACC varied from moderate to high according to the trait and models, ranging between 0.418 and 0.546 (PCV), between 0.646 and 0.793 (FEC), and between 0.414 and 0.519 (FAM). Parametric models presented nearly the same ACC and MSE for the studied traits and provided better accuracies than BRANN. The GBLUP, BayesA, BayesB and BLASSO models provided better accuracies than the BRANN_3 method, increasing by around 23% for PCV, and 18.5% for FEC. In conclusion, parametric models are suitable for genome-enabled prediction of indicator traits of resistance to gastrointestinal nematodes in sheep. Due to the small differences in accuracy found between them, the use of the GBLUP model is recommended due to its lower computational costs.

Estimating genetic parameters of reproductive, carcass, and meat quality traits in Polled Nellore cattle.

Considering the economic and commercial efficiency of the beef production chain, the yield and quality of the meat produced must also be included in breeding programs. For the Nellore breed, including the polled herd, these aspects have not been much studied. The aim of this study was to estimate genetic parameters for scrotal circumference adjusted to 365 (SC365) and 450 (SC450) days of age, age at first calving (AFC), accumulated productivity (AP), stayability (STAY), longissimus muscle area (LMA), thickness of subcutaneous fat over the 12th-13th ribs (BF), thickness of subcutaneous fat over the rump (RF), and shear force measured by Warner-Bratzler shear force (WBSF) of polled Nellore cattle. Bayesian analyses were performed by adopting a linear animal model, whereas STAY analyses used the linear threshold model. Heritability estimates were 0.31 (SC365), 0.37 (SC450), 0.16 (AFC), 0.25 (AP), 0.16 (STAY), 0.30 (LMA), 0.13 (BF), 0.24 (RF), and 0.15 (WBSF), indicating moderate response to selection. Genetic and residual correlations between SC365 and SC450 were high (0.91 and 0.74, respectively), as well as the genetic correlations of AP with SC365, SC450, AFC, and STAY (0.61, 0.62, - 0.69, and 0.83, respectively). Genetic and residual correlations of WBSF with reproductive and carcass characteristics exhibited high standard deviations, however favorable. Based on the results, it is expected that in the medium term, animals with greater sexual precocity will also have greater accumulated productivity and longer permanence of females in the herd, along with superior carcass traits. However, due to the low heritabilities and small genetic associations with reproductive traits, fat thickness characteristics (BF and RF) will still require direct selection.

Combined effect of purulent vaginal discharge and anovulation on pregnancy status in a large multi-state population of Holstein cows.

The objective of this observational prospective cohort study was to evaluate the combined effect of purulent vaginal discharge (PVD) and anovulation (ANOV) on the reproductive performance of a large multi-state population of Holstein cows. Data were prospectively collected from 11,729 cows in 16 herds located in 4 regions in the United States [Northeast (4 herds), Midwest (6), Southeast (1), and Southwest (5)]. Cows were enrolled at calving and monitored weekly for disease occurrence, reproductive events, and survival. Prevalence of PVD was evaluated at 28 ± 3 d in milk and defined by the presence of mucopurulent to fetid vaginal discharge. Resumption of ovarian cyclicity was determined via transrectal ultrasonography at 40 ± 3 and 54 ± 3 d postpartum. Pregnancy diagnosis was performed by ultrasonography on d 32 ± 3 after artificial insemination (AI) and reconfirmed at d 60 ± 3 of gestation. Pregnancy loss (PL) was defined as a cow diagnosed pregnant at 32 ± 3 but nonpregnant at 60 ± 3 d after AI. The association of PVD and ANOV with pregnancy traits was analyzed using 4 PVD-cyclicity categories that considered the following combinations: NPVD-CYC = absence of PVD and cycling; PVD-CYC = presence of PVD and cycling; NPVD-ANOV = absence of PVD and anovular; and PVD-ANOV = presence of PVD and anovular. Multiple logistic regression and Cox proportional regression were used for the analysis of potential associations between PVD and cyclicity categories and pregnancy at first AI (PAI1), days from calving to pregnancy, and PL at first AI. The odds (95% confidence intervals) of pregnancy increased from cows in the PVD-ANOV category (reference category) to cows in NPVD-ANOV [2.09 (1.62-2.50)], PVD-CYC [2.52 (2.02-3.14)], and NPVD-CYC [3.46 (2.84-4.23)]. Similarly, days from calving to pregnancy were less for NPVD-CYC, followed by PVD-CYC, NPVD-ANOV, and PVD-ANOV (121.4, 137.2, 137.3, and 157.4 d, respectively). On the contrary, no clear association was identified between groups and PL. The results indicate that both PVD and ANOV had a negative impact on PAI1 and days from calving to pregnancy. The results indicated a variable magnitude in the negative impact on the reproductive traits analyzed when both conditions were combined.

Investigating the effect of temporal, geographic, and management factors on US Holstein lactation curve parameters.

We fit the Wood's lactation model to an extensive database of test-day milk production records of US Holstein cows to obtain lactation-specific parameter estimates and investigated the effects of temporal, spatial, and management factors on lactation curve parameters and 305-d milk yield. Our approach included 2 steps as follows: (1) individual animal-parity parameter estimation with nonlinear least-squares optimization of the Wood's lactation curve parameters, and (2) mixed-effects model analysis of 8,595,413 sets of parameter estimates from individual lactation curves. Further, we conducted an analysis that included all parities and a separate analysis for first lactation heifers. Results showed that parity had the most significant effect on the scale (parameter a), the rate of decay (parameter c), and the 305-d milk yield. The month of calving had the largest effect on the rate of increase (parameter b) for models fit with data from all lactations. The calving month had the most significant effect on all lactation curve parameters for first lactation models. However, age at first calving, year, and milking frequency accounted for a higher proportion of the variance than month for first lactation 305-d milk yield. All parameter estimates and 305-d milk yield increased as parity increased; parameter a and 305-d milk yield rose, and parameters b and c decreased as year and milking frequency increased. Calving month estimates parameters a, b, c, and 305-d milk yield were the lowest values for September, May, June, and July, respectively. The results also indicated the random effects of herd and cow improved model fit. Lactation curve parameter estimates from the mixed-model analysis of individual lactation curve fits describe well US Holstein lactation curves according to temporal, spatial, and management factors.

Effect of body condition change and health status during early lactation on performance and survival of Holstein cows.

Body condition score (BCS) and disease records are commonly available in dairy operations. However, the effect of BCS changes (ΔBCS) considering specific health profiles has not been investigated extensively. The objective of this study was to assess the effects of different levels of ΔBCS on fertility, milk yield, and survival of Holstein cows diagnosed with reproductive disorders (REP; dystocia, twins, retained fetal membranes, metritis, and clinical endometritis), other health disorders (OTH; subclinical ketosis, left displaced abomasum, lameness, clinical mastitis, and respiratory disease), or with no disease events (HLT) within 40 days in milk (DIM). Data included lactation information from 11,733 cows calving between November 2012 and October 2014 in 16 herds across 4 geographical regions in the United States (Northeast, Midwest, Southwest, Southeast). Cows were evaluated for BCS at 5 ± 3 DIM (BCS5) and at 40 ± 3 DIM (BCS40) and the difference between BCS40 and BCS5 was classified as excessive loss of BCS (EL; ΔBCS ≤-0.75), moderate loss (ML; ΔBCS = -0.5 to -0.25), no change (NC; ΔBCS = 0), or gain of BCS (GN; ΔBCS ≥0.25). Multivariable logistic regression was used for assessing potential associations between the outcomes of interest and ΔBCS and health. The effect of the interaction term ΔBCS by health group was not statistically significant for any of the study outcomes. The odds of resumption of ovarian cyclicity (ROC), in GN, NC, and ML cows were 1.94 (95% CI: 1.57-2.40), 1.59 (1.28-1.97), and 1.27 (1.10-1.47) times greater than the odds of ROC in EL cows, respectively. The odds of pregnancy at 150 DIM (P150) in GN cows were 1.61 (1.20-2.17) times greater than the odds of P150 in EL cows. Cows with REP or OTH disorders had smaller odds of ROC compared with HLT cows [REP: OR = 0.65 (0.56-0.76) and OTH: OR = 0.79 (0.68-0.92)]. For pregnancy outcomes, REP cows had smaller odds of pregnancy at the first artificial insemination compared with HLT cows [0.70 (0.58-0.84)]. Similarly, REP cows had smaller odds of being diagnosed pregnant by 150 and 305 DIM compared with HLT cows [P150: 0.73 (0.59-0.87), P305: 0.58 (0.49-0.69)]. Overall, average daily milk within the first 90 DIM was greater in EL (39.5 ± 1.13 kg/d) and ML (38.9 ± 1.11 kg/d) cows than in NC (37.8 ± 1.12 kg/d) and GN (36.2 ± 1.12 kg/d) cows. On the other hand, average daily milk within the first 90 DIM was lower in REP (37.0 ± 1.11 kg/d) cows compared with OTH (38.7 ± 1.12 kg/d) and HLT cows (38.6 ± 1.11 kg/d). The magnitude of ΔBCS and the health status of early lactation cows should be considered when assessing subsequent cow performance and survival.

Genome-enabled prediction of meat and carcass traits using Bayesian regression, single-step genomic best linear unbiased prediction and blending methods in Nelore cattle.

Several methods have been used for genome-enabled prediction (or genomic selection) of complex traits, for example, multiple regression models describing a target trait with a linear function of a set of genetic markers. Genomic selection studies have been focused mostly on single-trait analyses. However, most profitability traits are genetically correlated, and an increase in prediction accuracy of genomic breeding values for genetically correlated traits is expected when using multiple-trait models. Thus, this study was carried out to assess the accuracy of genomic prediction for carcass and meat quality traits in Nelore cattle, using single- and multiple-trait approaches. The study considered 15 780, 15 784, 15 742 and 526 records of rib eye area (REA, cm2), back fat thickness (BF, mm), rump fat (RF, mm) and Warner-Bratzler shear force (WBSF, kg), respectively, in Nelore cattle, from the Nelore Brazil Breeding Program. Animals were genotyped with a low-density single nucleotide polymorphism (SNP) panel and subsequently imputed to arrays with 54 and 777 k SNPs. Four Bayesian specifications of genomic regression models, namely, Bayes A, Bayes B, Bayes Cπ and Bayesian Ridge Regression; blending methods, BLUP; and single-step genomic best linear unbiased prediction (ssGBLUP) methods were compared in terms of prediction accuracy using a fivefold cross-validation. Estimates of heritability ranged from 0.20 to 0.35 and from 0.21 to 0.46 for RF and WBSF on single- and multiple-trait analyses, respectively. Prediction accuracies for REA, BF, RF and WBSF were all similar using the different specifications of regression models. In addition, this study has shown the impact of genomic information upon genetic evaluations in beef cattle using the multiple-trait model, which was also advantageous compared to the single-trait model because it accounted for the selection process using multiple traits at the same time. The advantage of multi-trait analyses is attributed to the consideration of correlations and genetic influences between the traits, in addition to the non-random association of alleles.

Genome-wide association study of milk production traits in a crossbred dairy sheep population using three statistical models.

Milk production is one of the most important characteristics of dairy sheep, and the identification of genes affecting milk production traits is critical to understanding the genetics and improve milk production in future generations. Three statistical techniques, namely GWAS, ridge-regression BLUP and BayesC π , were used to identify SNPs in significant association with three milk production traits (milk yield, fat yield and protein yield) in a crossbred dairy sheep population. The results suggested that chromosomes 1, 3, 4, 5, 7 and 11 were likely to harbor genes important to milk production because these chromosomes had the greatest top-100-SNP variance contributions on the three milk production traits. The GWAS analysis identified between 74 and 288 genome-wide significant SNP (P < 0.05) whereas the BayesCπ model revealed between six and 63 SNPs, each with >95% posterior probability of inclusion as having a non-zero association effect on at least one of the three milk production traits. Positional candidate genes for milk production in sheep were searched, based on the sheep genomic assembly OAR version 3.1, such as those which map position coincided with or was located within 0.1 Mbp of a genome-wide suggestive or significant SNP. These identified SNPs and candidate genes supported some previous findings and also added new information about genetic markers for genetic improvement of lactation in dairy sheep, but keeping in mind that the majority of these positional candidate genes are not necessarily true causative loci for these traits and future validations are thus necessary.

Estimation of genomic breed composition of individual animals in composite beef cattle.

Three statistical models (an admixture model, linear regression, and ridge-regression BLUP) and two strategies for selecting SNP panels (uniformly spaced vs. maximum Euclidean distance of SNP allele frequencies between ancestral breeds) were compared for estimating genomic-estimated breed composition (GBC) in Brangus and Santa Gertrudis cattle, respectively. Animals were genotyped with a GeneSeek Genomic Profiler bovine low-density version 4 SNP chip. The estimated GBC was consistent among the uniformly spaced SNP panels, and values were similar between the three models. However, estimated GBC varied considerably between the three methods when using fewer than 10 000 SNPs that maximized the Euclidean distance of allele frequencies between the ancestral breeds. The admixture model performed most consistently across various SNP panel sizes. For the other two models, stabilized estimates were obtained with an SNP panel size of 20 000 SNPs or more. Based on the uniformly spaced 20K SNP panel, the estimated GBC was 69.8-70.5% Angus and 29.5-30.2% Brahman for Brangus, and 63.9-65.3% Shorthorn and 34.7-36.1% Brahman in Santa Gertrudis. The estimated GBC of ancestries for Santa Gertrudis roughly agreed with the pedigree-expected values. However, the estimated GBC in Brangus showed a considerably larger Angus composition than the pedigree-expected value (62.5%). The elevated Angus composition in the Brangus could be due to the mixture of some 1/2 Ultrablack animals (Brangus × Angus). Another reason could be the consequences of selection in Brangus cattle for phenotypes where the Angus breed has advantages.

Mining data from milk infrared spectroscopy to improve feed intake predictions in lactating dairy cows.

Feed intake is one of the most important components of feed efficiency in dairy systems. However, it is a difficult trait to measure in commercial operations for individual cows. Milk spectrum from mid-infrared spectroscopy has been previously used to predict milk traits, and could be an alternative to predict dry matter intake (DMI). The objectives of this study were (1) to evaluate if milk spectra can improve DMI predictions based only on cow variables; (2) to compare artificial neural network (ANN) and partial least squares (PLS) predictions; and (3) to evaluate if wavelength (WL) selection through Bayesian network (BN) improves prediction quality. Milk samples (n = 1,279) from 308 mid-lactation dairy cows [127 ± 27 d in milk (DIM)] were collected between 2014 and 2016. For each milk spectra time point, DMI (kg/d), body weight (BW, kg), milk yield (MY, kg/d), fat (%), protein (%), lactose (%), and actual DIM were recorded. The DMI was predicted with ANN and PLS using different combinations of explanatory variables. Such combinations, called covariate sets, were as follows: set 1 (MY, BW0.75, DIM, and 361 WL); set 2 [MY, BW0.75, DIM, and 33 WL (WL selected by BN)]; set 3 (MY, BW0.75, DIM, and fat, protein, and lactose concentrations); set 4 (MY, BW0.75, DIM, 33 WL, fat, protein, and lactose); set 5 (MY, BW0.75, DIM, 33 WL, and visit duration in the feed bunk); set 6 (MY, DIM, and 33 WL); set 7 (MY, BW0.75, and DIM); set-WL (included 361 WL); and set-BN (included just 33 selected WL). All models (i.e., each combination of covariate set and fitting approach, ANN or PLS) were validated with an external data set. The use of ANN improved the performance of models 2, 5, 6, and BN. The use of BN combined with ANN yielded the highest accuracy and precision. The addition of individual WL compared with milk components (set 2 vs. set 3) did not improve prediction quality when using PLS. However, when ANN was employed, the model prediction with the inclusion of 33 WL was improved over the model containing only milk components (set 2 vs. set 3; concordance correlation coefficient = 0.80 vs. 0.72; coefficient of determination = 0.67 vs. 0.53; root mean square error of prediction 2.36 vs. 2.81 kg/d). The use of ANN and the inclusion of a behavior parameter, set 5, resulted in the best predictions compared with all other models (coefficient of determination = 0.70, concordance correlation coefficient = 0.83, root mean square error of prediction = 2.15 kg/d). The addition of milk spectra information to models containing cow variables improved the accuracy and precision of DMI predictions in lactating dairy cows when ANN was used. The use of BN to select more informative WL improved the model prediction when combined with cow variables, with further improvement when combined with ANN.

Improving accuracy of genomic prediction in Brangus cattle by adding animals with imputed low-density SNP genotypes.

Reliable genomic prediction of breeding values for quantitative traits requires the availability of sufficient number of animals with genotypes and phenotypes in the training set. As of 31 October 2016, there were 3,797 Brangus animals with genotypes and phenotypes. These Brangus animals were genotyped using different commercial SNP chips. Of them, the largest group consisted of 1,535 animals genotyped by the GGP-LDV4 SNP chip. The remaining 2,262 genotypes were imputed to the SNP content of the GGP-LDV4 chip, so that the number of animals available for training the genomic prediction models was more than doubled. The present study showed that the pooling of animals with both original or imputed 40K SNP genotypes substantially increased genomic prediction accuracies on the ten traits. By supplementing imputed genotypes, the relative gains in genomic prediction accuracies on estimated breeding values (EBV) were from 12.60% to 31.27%, and the relative gain in genomic prediction accuracies on de-regressed EBV was slightly small (i.e. 0.87%-18.75%). The present study also compared the performance of five genomic prediction models and two cross-validation methods. The five genomic models predicted EBV and de-regressed EBV of the ten traits similarly well. Of the two cross-validation methods, leave-one-out cross-validation maximized the number of animals at the stage of training for genomic prediction. Genomic prediction accuracy (GPA) on the ten quantitative traits was validated in 1,106 newly genotyped Brangus animals based on the SNP effects estimated in the previous set of 3,797 Brangus animals, and they were slightly lower than GPA in the original data. The present study was the first to leverage currently available genotype and phenotype resources in order to harness genomic prediction in Brangus beef cattle.

Bayesian analyses of genetic parameters for growth traits in Nellore cattle raised on pasture.

This study was carried out to investigate (co)variance components and genetic parameters for growth traits in beef cattle using a multi-trait model by Bayesian methods. Genetic and residual (co)variances and parameters were estimated for weights at standard ages of 120 (W120), 210 (W210), 365 (W365), and 450 days (W450), and for pre- and post-weaning daily weight gain (preWWG and postWWG) in Nellore cattle. Data were collected over 16 years (1993-2009), and all animals were raised on pasture in eight farms in the North of Brazil that participate in the National Association of Breeders and Researchers. Analyses were run by the Bayesian approach using Gibbs sampler. Additive direct heritabilities for W120, W210, W365, and W450 and for preWWG and postWWG were 0.28 ± 0.013, 0.32 ± 0.002, 0.31 ± 0.002, 0.50 ± 0.026, 0.61 ± 0.047, and 0.79 ± 0.055, respectively. The estimates of maternal heritability were 0.32 ± 0.012, 0.29 ± 0.004, 0.30 ± 0.005, 0.25 ± 0.015, 0.23 ± 0.017, and 0.22 ± 0.016, respectively, for W120, W210, W365, and W450 and for preWWG and postWWG. The estimates of genetic direct additive correlation among all traits were positive and ranged from 0.25 ± 0.03 (preWWG and postWWG) to 0.99 ± 0.00 (W210 and preWWG). The moderate to high estimates of heritability and genetic correlation for weights and daily weight gains at different ages is suggestive of genetic improvement in these traits by selection at an appropriate age. Maternal genetic effects seemed to be significant across the traits. When the focus is on direct and maternal effects, W210 seems to be a good criterium for the selection of Nellore cattle considering the importance of this breed as a major breed of beef cattle not only in Northern Brazil but all regions covered by tropical pastures. As in this study the genetic correlations among all traits were high, the selection based on weaning weight might be a good choice because at this age there are two important effects (maternal and direct genetic effects). In contrast, W120 should be preferred when the objective is improving the maternal ability of the dams. Furthermore, selection for postWWG can be used if the animals show both heavier weaning weights and high growth rate after weaning because it is possible to shorten the time between weaning and slaughter based on weaning weight, postWWG, and desired weight at the time of slaughter.

Genomewide association mapping and pathway analysis of meat tenderness in Polled Nellore cattle.

Brazil is one of the world's largest beef exporters, although the product has a low price due to quality issues. The meat exported by Brazil is considered medium and low quality by international buyers, mainly due to lack of tenderness. The predominant Zebu breeds (80% Nellore) are known for producing tougher beef than taurine breeds. Nonetheless, some studies have shown that there is substantial genetic variability for tenderness within the Nellore breed, although it is a difficult trait to improve by conventional selection methods. Therefore, the aim of this study was to perform a genomewide association study (GWAS) and a gene set enrichment analysis to identify genomic regions and biologically relevant pathways associated with meat tenderness in Polled Nellore cattle. Data consisted of Warner-Bratzler shear force values of LM from 427 Polled Nellore animals divided into 3 experimental slaughters (years 2005, 2008, and 2010). The animals were genotyped with either the Illumina BovineHD BeadChip (777k, on 61 samples) or the GGP Indicus HD chip (77k, on 366 samples). Single nucleotide polymorphisms were excluded when the call rate was <90%, the Hardy-Weinberg proportions -value was <1% (Fisher exact test, Bonferroni adjusted), and the minor allele frequency was <1%. Imputation from the GGP Indicus HD chip to the Illumina BovineHD BeadChip was performed using the FImput program. Genomewide association analysis was performed using the Efficient Mixed Model Association eXpedited (EMMAx) and the population parameters previously determined (P3D) methods. The GWAS was complemented with a gene set enrichment analysis performed using the FatiGO procedure. Significant markers ( < 0.0001) explaining a larger proportion of variation than other significant SNPs were located on chromosomes 3, 13, 17, 20, 21, and 25, indicating QTL associated with meat tenderness throughout the genome. Additionally, gene set analysis identified 22 Gene Ontology functional terms and 2 InterPro entries that showed significant enrichment of genes associated with tenderness. The functional categories included protein tyrosine and serine/threonine kinase activity, calcium ion binding, lipid metabolic process, and growth factors, among others. These results help to elucidate the genetic architecture and metabolic pathways underlying this trait, which is of extreme economic and social importance to Brazil, because Nellore is the dominant beef cattle breed in the country.

Evaluation of longevity modeling censored records in Nellore.

The aim of the present study was to evaluate the prediction ability of models that cope with longevity phenotypic expression as uncensored and censored in Nellore cattle. Longevity was defined as the difference between the dates of last weaned calf and cow birth. There were information of 77 353 females, being 61 097 cows with uncensored phenotypic information and 16 256 cows with censored records. These data were analyzed considering three different models: (1) Gaussian linear model (LM), in which only uncensored records were considered; and two models that consider both uncensored and censored records: (2) Censored Gaussian linear model (CLM); and (3) Weibull frailty hazard model (WM). For the model prediction ability comparisons, the data set was randomly divided into training and validation sets, containing 80% and 20% of the records, respectively. There were considered 10 repetitions applying the following restrictions: (a) at least three animals per contemporary group in the training set; and (b) sires with more than 10 progenies with uncensored records (352 sires) should have daughters in the training and validation sets. The variance components estimated using the whole data set in each model were used as true values in the prediction of breeding values of the animals in the training set. The WM model showed the best prediction ability, providing the lowest χ 2 average and the highest number of sets in which a model had the smallest value of χ 2 statistics. The CLM and LM models showed prediction abilities 2.6% and 3.7% less efficient than WM, respectively. In addition, the accuracies of sire breeding values for LM and CLM were lower than those obtained for WM. The percentages of bulls in common, considering only 10% of sires with the highest breeding values, were around 75% and 54%, respectively, between LM-CLM and LM-WM models, considering all sires, and 75% between LM-CLM and LM-WM, when only sires with more than 10 progenies with uncensored records were taken into account. These results are indicative of reranking of animals in terms of genetic merit between LM, CLM and WM. The model in which censored records of longevity were excluded from the analysis showed the lowest prediction ability. The WM provides the best predictive performance, therefore this model would be recommended to perform genetic evaluation of longevity in this population.

Pathway-based genome-wide association analysis of milk coagulation properties, curd firmness, cheese yield, and curd nutrient recovery in dairy cattle.

It is becoming common to complement genome-wide association studies (GWAS) with gene-set enrichment analysis to deepen the understanding of the biological pathways affecting quantitative traits. Our objective was to conduct a gene ontology and pathway-based analysis to identify possible biological mechanisms involved in the regulation of bovine milk technological traits: coagulation properties, curd firmness modeling, individual cheese yield (CY), and milk nutrient recovery into the curd (REC) or whey loss traits. Results from 2 previous GWAS studies using 1,011 cows genotyped for 50k single nucleotide polymorphisms were used. Overall, the phenotypes analyzed consisted of 3 traditional milk coagulation property measures [RCT: rennet coagulation time defined as the time (min) from addition of enzyme to the beginning of coagulation; k20: the interval (min) from RCT to the time at which a curd firmness of 20 mm is attained; a30: a measure of the extent of curd firmness (mm) 30 min after coagulant addition], 6 curd firmness modeling traits [RCTeq: RCT estimated through the CF equation (min); CFP: potential asymptotic curd firmness (mm); kCF: curd-firming rate constant (% × min-1); kSR: syneresis rate constant (% × min-1); CFmax: maximum curd firmness (mm); and tmax: time to CFmax (min)], 3 individual CY-related traits expressing the weight of fresh curd (%CYCURD), curd solids (%CYSOLIDS), and curd moisture (%CYWATER) as a percentage of weight of milk processed and 4 milk nutrient and energy recoveries in the curd (RECFAT, RECPROTEIN, RECSOLIDS, and RECENERGY calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk), milk pH, and protein percentage. Each trait was analyzed separately. In total, 13,269 annotated genes were used in the analysis. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway databases were queried for enrichment analyses. Overall, 21 Gene Ontology and 17 Kyoto Encyclopedia of Genes and Genomes categories were significantly associated (false discovery rate at 5%) with 7 traits (RCT, RCTeq, kCF, %CYSOLIDS, RECFAT, RECSOLIDS, and RECENERGY), with some being in common between traits. The significantly enriched categories included calcium signaling pathway, salivary secretion, metabolic pathways, carbohydrate digestion and absorption, the tight junction and the phosphatidylinositol pathways, as well as pathways related to the bovine mammary gland health status, and contained a total of 150 genes spanning all chromosomes but 9, 20, and 27. This study provided new insights into the regulation of bovine milk coagulation and cheese ability that were not captured by the GWAS.

Genome-wide association study for cheese yield and curd nutrient recovery in dairy cows.

Cheese production and consumption are increasing in many countries worldwide. As a result, interest has increased in strategies for genetic selection of individuals for technological traits of milk related to cheese yield (CY) in dairy cattle breeding. However, little is known about the genetic background of a cow's ability to produce cheese. Recently, a relatively large panel (1,264 cows) of different measures of individual cow CY and milk nutrient and energy recoveries in the cheese (REC) became available. Genetic analyses showed considerable variation for CY and for aptitude to retain high proportions of fat, protein, and water in the coagulum. For the dairy industry, these characteristics are of major economic importance. Nevertheless, use of this knowledge in dairy breeding is hampered by high costs, intense labor requirement, and lack of appropriate technology. However, in the era of genomics, new possibilities are available for animal breeding and genetic improvement. For example, identification of genomic regions involved in cow CY might provide potential for marker-assisted selection. The objective of this study was to perform genome-wide association studies on different CY and REC measures. Milk and DNA samples from 1,152 Italian Brown Swiss cows were used. Three CY traits expressing the weight (wt) of fresh curd (%CYCURD), curd solids (%CYSOLIDS), and curd moisture (%CYWATER) as a percentage of weight of milk processed, and 4 REC (RECFAT, RECPROTEIN, RECSOLIDS, and RECENERGY, calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk) were analyzed. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2. Single marker regressions were fitted using the GenABEL R package (genome-wide association using mixed model and regression-genomic control). In total, 103 significant associations (88 single nucleotide polymorphisms) were identified in 10 chromosomes (2, 6, 9, 11, 12, 14, 18, 19, 27, 28). For RECFAT and RECPROTEIN, high significance peaks were identified in Bos taurus autosome (BTA) 6 and BTA11, respectively. Marker ARS-BFGL-NGS-104610 (∼104.3 Mbp) was highly associated with RECPROTEIN and Hapmap52348-rs29024684 (∼87.4 Mbp), closely located to the casein genes on BTA6, with RECFAT. Genomic regions identified may enhance marker-assisted selection in bovine cheese breeding beyond the use of protein (casein) and fat contents, whereas new knowledge will help to unravel the genomic background of a cow's ability for cheese production.

Inferring phenotypic causal structures among meat quality traits and the application of a structural equation model in Japanese Black cattle.

Meat quality is one of the most important traits determining carcass price in the Japanese beef market. Optimized breeding goals and management practices for the improvement of meat quality traits requires knowledge regarding any potential functional relationships between them. In this context, the objective of this research was to infer phenotypic causal networks involving beef marbling score (BMS), beef color score (BCL), firmness of beef (FIR), texture of beef (TEX), beef fat color score (BFS), and the ratio of MUFA to SFA (MUS) from 11,855 Japanese Black cattle. The inductive causation (IC) algorithm was implemented to search for causal links among these traits and was conditionally applied to their joint distribution on genetic effects. This information was obtained from the posterior distribution of the residual (co)variance matrix of a standard Bayesian multiple trait model (MTM). Apart from BFS, the IC algorithm implemented with 95% highest posterior density (HPD) intervals detected only undirected links among the traits. However, as a result of the application of 80% HPD intervals, more links were recovered and the undirected links were changed into directed ones, except between FIR and TEX. Therefore, 2 competing causal networks resulting from the IC algorithm, with either the arrow FIR → TEX or the arrow FIR ← TEX, were fitted using a structural equation model () to infer causal structure coefficients between the selected traits. Results indicated similar genetic and residual variances as well as genetic correlation estimates from both structural equation models. The genetic variances in BMS, FIR, and TEX from the structural equation models were smaller than those obtained from the MTM. In contrast, the variances in BCL, BFS, and MUS, which were not conditioned on any of the other traits in the causal structures, had no significant differences between the structural equation model and MTM. The structural coefficient for the path from MUS (BCL) to BMS showed that a 1-unit improvement in MUS (BCL) resulted in an increase of 0.85 or 1.45 (an decrease of 0.52 or 0.54) in BMS in the causal structures. The analysis revealed some interesting functional relationships, direct genetic effects, and the magnitude of the causal effects between these traits, for example, indicating that BMS would be affected by interventions on MUS and BCL. In addition, if interventions existed in this scenario, a breeding strategy based only on the MTM would lead to a mistaken selection for BMS.

Genome scan for postmortem carcass traits in Nellore cattle.

Carcass traits measured after slaughter are economically relevant traits in beef cattle. In general, the slaughter house payment system is based on HCW. Ribeye area (REA) is associated with the amount of the meat in the carcass, and a minimum of backfat thickness (BFT) is necessary to protect the carcass during cooling. The aim of this study was to identify potential genomic regions harboring candidate genes affecting those traits in Nellore cattle. The data set used in the present study consisted of 1,756 Nellore males with phenotype records. A subset of 1,604 animals had both genotypic and phenotypic information. Genotypes were generated based on a panel with 777,962 SNPs from the Illumina Bovine HD chip. The SNP effects were calculated based on the genomic breeding values obtained by using the single-step GBLUP approach and a genomic matrix re-weighting procedure. The proportion of the variance explained by moving windows of 100 consecutive SNPs was used to assess potential genomic regions harboring genes with major effects on each trait. The top 10 non-overlapping SNP-windows explained 8.72%, 11.38%, and 9.31% of the genetic variance for REA, BFT, and HCW, respectively. These windows are located on chromosomes 5, 7, 8, 10, 12, 20, and 29 for REA; chromosomes 6, 8, 10, 13, 16, 17, 18, and 24 for BFT; and chromosomes 4, 6, 7, 8, 14, 16, 17, and 21 for HCW. For REA, there were identified genes ( and ) involved in the cell cycle biological process which affects many aspects of animal growth and development. The and genes, both from AA transporter family, was also associated with REA. The AA transporters are essential for cell growth and proliferation, acting as carriers of tissue nutrient supplies. Various genes identified for BFT (, , , , , and ) have been associated with lipid metabolism in different mammal species. One of the most promising genes identified for HCW was the . There is evidence, in the literature, that this gene is located in putative QTL affecting carcass weight in beef cattle. Our results showed several genomic regions containing plausible candidate genes that may be associated with carcass traits in Nellore cattle. Besides contributing to a better understanding of the genetic control of carcass traits, the identified genes can also be helpful for further functional genomic studies.

Accuracy of genomic breeding values for meat tenderness in Polled Nellore cattle.

Zebu () cattle, mostly of the Nellore breed, comprise more than 80% of the beef cattle in Brazil, given their tolerance of the tropical climate and high resistance to ectoparasites. Despite their advantages for production in tropical environments, zebu cattle tend to produce tougher meat than Bos taurus breeds. Traditional genetic selection to improve meat tenderness is constrained by the difficulty and cost of phenotypic evaluation for meat quality. Therefore, genomic selection may be the best strategy to improve meat quality traits. This study was performed to compare the accuracies of different Bayesian regression models in predicting molecular breeding values for meat tenderness in Polled Nellore cattle. The data set was composed of Warner-Bratzler shear force (WBSF) of longissimus muscle from 205, 141, and 81 animals slaughtered in 2005, 2010, and 2012, respectively, which were selected and mated so as to create extreme segregation for WBSF. The animals were genotyped with either the Illumina BovineHD (HD; 777,000 from 90 samples) chip or the GeneSeek Genomic Profiler (GGP Indicus HD; 77,000 from 337 samples). The quality controls of SNP were Hard-Weinberg Proportion -value ≥ 0.1%, minor allele frequency > 1%, and call rate > 90%. The FImpute program was used for imputation from the GGP Indicus HD chip to the HD chip. The effect of each SNP was estimated using ridge regression, least absolute shrinkage and selection operator (LASSO), Bayes A, Bayes B, and Bayes Cπ methods. Different numbers of SNP were used, with 1, 2, 3, 4, 5, 7, 10, 20, 40, 60, 80, or 100% of the markers preselected based on their significance test (-value from genomewide association studies [GWAS]) or randomly sampled. The prediction accuracy was assessed by the correlation between genomic breeding value and the observed WBSF phenotype, using a leave-one-out cross-validation methodology. The prediction accuracies using all markers were all very similar for all models, ranging from 0.22 (Bayes Cπ) to 0.25 (Bayes B). When preselecting SNP based on GWAS results, the highest correlation (0.27) between WBSF and the genomic breeding value was achieved using the Bayesian LASSO model with 15,030 (3%) markers. Although this study used relatively few animals, the design of the segregating population ensured wide genetic variability for meat tenderness, which was important to achieve acceptable accuracy of genomic prediction. Although all models showed similar levels of prediction accuracy, some small advantages were observed with the Bayes B approach when higher numbers of markers were preselected based on their -values resulting from a GWAS analysis.

Comparison of models for the genetic evaluation of reproductive traits with censored data in Nellore cattle.

In typical genetic evaluation, often some females have missing records due to reproductive failure and due to voluntary and involuntary culling before the breeding season. These partially or unobserved phenotypes are known as censored records and their inclusion into genetic evaluations might lead to better inferences and breeding value predictions. Then, the objective was to compare prediction ability of models in which the phenotypic expression of age at the first calving (AFC) and days to calving (DC) were considered to be censored and uncensored in a Nellore cattle population. Age at first calving and days to calving were analyzed as following: uncensored animals (LM); penalization of 21 d (PLM); censored records simulated from truncated normal distributions (CLM); threshold-linear model in which censored records were handled as missing (TLM) or coded as the upper AFC/DC value within contemporary group (PTLM); and Weibull frailty hazard model (WM). Pearson correlations (PC), the percentage of the 10% best bulls in common (pTOP10%), accuracy of estimated breeding values (), and a cross-validation scheme were performed. Heritability estimates for AFC were 0.18, 0.12, 0.12, 0.17, 0.14, and 0.07 for LM, PLM, CLM, TLM, PTLM, and WM, respectively. PC and pTOP10% were higher among linear models and smaller between these models and WM. The models provided similar r of sire breeding values. Heritability estimates for DC were 0.03, 0.08, 0.06, 0.02, 0.07, and 0.10 for LM, PLM, CLM, TLM, PTLM, and WM, respectively. Strongly associated predictions were observed in CLM, PLM, PTLM, and WM. The highest coincidence levels of sires in the TOP10% were between CLM, PLM, and PTLM. The r of sire breeding values obtained applying CLM, PLM, PTLM, and WM were similar and higher than those obtained with LM and TLM. In terms of prediction ability, WM, PLM, TLM, and PTLM showed similar prediction performance for AFC. On the other hand, CLM, PLM, PTLM, and WM showed the similar prediction ability for DC Therefore, these models would be recommended to perform genetic evaluation of age at first calving and days to calving in this Nellore population.