Genomic heritability of bovine growth using a mixed model

Genomic heritability and correlation between carcass traits in Japanese Black cattle evaluated under different ceilings of relatedness among individuals

The investigation of carcass traits to produce meat with high efficiency has been in focus on Japanese Black cattle since 1972. To implement a successful breeding program in carcass production, a comprehensive understanding of genetic characteristics and relationships between the traits is of paramount importance. In this study, genomic heritability and genomic correlation between carcass traits, including carcass weight (CW), rib eye area (REA), rib thickness (RT), subcutaneous fat thickness (SFT), yield rate (YI), and beef marbling score (BMS) were estimated using the genomic data of 9,850 Japanese Black cattle (4,142 heifers and 5,708 steers). In addition, we investigated the effect of genetic relatedness degree on the estimation of genetic parameters of carcass traits in sub-populations created based on different GRM-cutoff values. Genome-based restricted maximum likelihood (GREML) analysis was applied to estimate genetic parameters. Using all animal data, the heritability values for carcass traits were estimated as moderate to relatively high magnitude, ranging from 0.338 to 0.509 with standard errors, ranging from 0.014 to 0.015. The genetic correlations were obtained low and negative between SFT and REA [-0.198 (0.034)] and between SFT and BMS [-0.096 (0.033)] traits, and high and negative between SFT and YI [-0.634 (0.022)]. REA trait was genetically highly correlated with YI and BMS [0.811 (0.012) and 0.625 (0.022), respectively]. In sub-populations created based on the genetic-relatedness ceiling, the heritability estimates ranged from 0.212 (0.131) to 0.647 (0.066). At the genetic-relatedness ceiling of 0.15, the correlation values between most traits with low genomic correlation were overestimated while the correlations between the traits with relatively moderate to high correlations, ranging from 0.380 to 0.811, were underestimated. The values were steady at the ceilings of 0.30-0.95 (sample size of 5,443-9,850) for most of the highly correlated traits. The results demonstrated that there is considerable genetic variation and also favorable genomic correlations between carcass traits. Therefore, the genetic improvement for the traits can be simultaneously attained through genomic selection. In addition, we observed that depending on the degree of relationship between individuals and sample size, the genomic heritability and correlation estimates for carcass traits may be different.

Genome-wide association and genotype by environment interactions for growth traits in U.S. Red Angus cattle

Background

Genotypic information produced from single nucleotide polymorphism (SNP) arrays has routinely been used to identify genomic regions associated with complex traits in beef and dairy cattle. Herein, we assembled a dataset consisting of 15,815 Red Angus beef cattle distributed across the continental U.S. and a union set of 836,118 imputed SNPs to conduct genome-wide association analyses (GWAA) for growth traits using univariate linear mixed models (LMM); including birth weight, weaning weight, and yearling weight. Genomic relationship matrix heritability estimates were produced for all growth traits, and genotype-by-environment (GxE) interactions were investigated.

Results

Moderate to high heritabilities with small standard errors were estimated for birth weight (0.51 ± 0.01), weaning weight (0.25 ± 0.01), and yearling weight (0.42 ± 0.01). GWAA revealed 12 pleiotropic QTL (BTA6, BTA14, BTA20) influencing Red Angus birth weight, weaning weight, and yearling weight which met a nominal significance threshold (P ≤ 1e-05) for polygenic traits using 836K imputed SNPs. Moreover, positional candidate genes associated with Red Angus growth traits in this study (i.e., LCORL, LOC782905, NCAPG, HERC6, FAM184B, SLIT2, MMRN1, KCNIP4, CCSER1, GRID2, ARRDC3, PLAG1, IMPAD1, NSMAF, PENK, LOC112449660, MOS, SH3PXD2B, STC2, CPEB4) were also previously associated with feed efficiency, growth, and carcass traits in beef cattle. Collectively, 14 significant GxE interactions were also detected, but were less consistent among the investigated traits at a nominal significance threshold (P ≤ 1e-05); with one pleiotropic GxE interaction detected on BTA28 (24 Mb) for Red Angus weaning weight and yearling weight.

Conclusions

Sixteen well-supported QTL regions detected from the GWAA and GxE GWAA for growth traits (birth weight, weaning weight, yearling weight) in U.S. Red Angus cattle were found to be pleiotropic. Twelve of these pleiotropic QTL were also identified in previous studies focusing on feed efficiency and growth traits in multiple beef breeds and/or their composites. In agreement with other beef cattle GxE studies our results implicate the role of vasodilation, metabolism, and the nervous system in the genetic sensitivity to environmental stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08667-6.

Comparison of Pure and Crossbred Japanese Black Steers in Growth Performance and Metabolic Features from Birth to Slaughter at a Spanish Fattening Farm

Simple Summary

Cattle growth performance is a determinant of beef production. Nowadays, customers demand specialized, high-quality beef products produced according to stringent health and welfare standards. Intramuscular fat, or marbling, improves beef quality, and the Japanese Black (Wagyu) is the breed with the highest rates of marbling. Wagyu steers are reared under specific conditions in Japan, which may differ from the conditions in other countries, and these differences may affect animal well-being and, therefore, growth rates and beef quality. The current study shows that purebred Wagyu and crossbred Wagyu-by-Angus steers that were raised at a cow–calf operation and fattening system in Spain with no exercise restriction, high welfare, and a local diet high in olein content showed appropriate growth and fattening rates, health status, and metabolic development. Wagyu crossbred steers did not show substantially faster growth than purebred Wagyu animals, so they may not be as profitable as purebred Japanese Black in this type of production system.

Abstract

Japanese Black (Wagyu) cattle produce high-quality beef. However, whether Wagyu steers can be profitably raised under conditions different than the traditional Japanese ones remains unclear. From 2018 to 2020, we raised 262 Wagyu purebred steers, 103 Wagyu-by-Angus (Wangus) crossbred steers, and 43 Angus-by-European (ACL) crossbred steers on a Spanish farm with high welfare standards and a locally sourced, high-olein diet. Factors and factors’ interactions impacting steer growth were analyzed using generalized linear models. ACL steers grew faster than the other two groups, with Wangus showing intermediate fattening and muscle development. Average daily weight gains (kg/day) were 0.916 for Wagyu, 1.046 for Wangus, and 1.293 for ACL during the weaning to growing period, and 0.628 for Wagyu, 0.64 for Wangus, and 0.802 for ACL during the growing to fattening phase. ACL showed the lowest marbling rates. Wagyu and Wangus usually showed higher cholesterol, triglycerides, and high-density lipoprotein than ACL. ACL calves may experience greater stress at weaning, as suggested by higher glucose, lactate, and β-hydroxybutyrate than the other groups. The results suggest that Wagyu and Wangus steers showed adequate growth, health, and metabolic development in this type of production system, with Wagyu purebreds probably being more profitable than Wangus crossbreeds.

Evaluation of Genetic Parameters of Growth of Pelibuey and Blackbelly Sheep through Pedigree in Mexico

Simple Summary

The importance of sheep production in Mexico has been increasing in recent years. Animal performance can be improved through continuous selection and cumulative genetic management. Studies on non-genetic factors that influence the growth characteristics of sheep, due to their positive genetic correlation with other live weights, can present data of great interest to small producers with limited resources. In the present study, genealogical and functional information from the historical archive of the Mexican Association of Sheep Breeders (AMCO) was used. The objective was to estimate the heritability and genetic correlations of the growth of ewes born and weaned at different times of the same year from different herds with pedigree registration. The expected results could be incorporated into genetic selection programs with a favorable impact on the local economy of small sheep producers in Mexico.

Abstract

Birth weight (BW) and weaning weight (WW) data from Pelibuey and Blackbelly lambs belonging to the Asociación Mexicana de Criadores de Ovinos (AMCO) were used with the objective of estimating genetic parameters (heritability and genetic correlations) and analyzing the growth characteristics of ewes born and weaned at different times of the same year from different herds with pedigree registration. In the case of Pelibuey lambs, the animal model included the weaning weight at 75 days of age, considering the direct additive genetic effect, maternal additive genetic effect, covariance between direct and maternal effects, as well as the permanent environmental effect of the mother. The direct estimators of heritability for Pelibuey were BW = 0.01 ± 0.021 and WW = 0.31 ± 0.074 and for Blackbelly they were BW = 0.05 ± 0.042 and WW = 0.41 ± 0.146. In the case of the maternal heritability estimators in Pelibuey they were BW = 0.02 ± 0.040 and WW = 0.21 ± 0.121 and for Blackbelly they were BW = 0.12 ± 0.054 and WW = 0.28 ± 0.121. The magnitude of the estimates of genetic correlations between direct and maternal effects for adjusted weaning weight at 75 days of age indicate that genetic progress may be slow in a breeding program. However, these selection results could be included in the short term in the breeding programs for the Pelibuey and Blackbelly breeds in Mexico, for livestock development in low-income rural areas.

Performance of pedigree and various forms of marker-derived relationship coefficients in genomic prediction and their correlations

In recent years, with development and validation of different genotyping panels, several methods have been proposed to build efficient similarity matrices among individuals to be used for genomic selection. Consequently, the estimated genetic parameters from such information may deviate from their counterpart using traditional family information. In this study, we used a pedigree-based numerator relationship matrix (A) and three types of marker-based relationship matrices ( G ) including two identical by descent, that is G K and G M and one identical by state, G V as well as four Gaussian kernel ( GK ) similarity kernels with different smoothing parameters to predict yet to be observed phenotypes. Also, we used different kinship matrices that are a linear combination of marker-derived IBD or IBS matrices with A, constructed as K = λ G + 1 - λ A , where the weight ( λ ) assigned to each source of information varied over a grid of values. A Bayesian multiple-trait Gaussian model was fitted to estimate the genetic parameters and compare the prediction accuracy in terms of predictive correlation, mean square error and unbiasedness. Results show that the estimated genetic parameters (heritability and correlations) are affected by the source of the information used to create kinship or the weight placed on the sources of genomic and pedigree information. The superiority of GK-based model depends on the smoothing parameters (θ) so that with an optimum θ value, the GK-based model statistically yielded better performance (higher predictive correlation, lowest MSE and unbiased estimates) and more stable correlations and heritability than the model with IBD, IBS or A kinship matrices or any of the linear combinations.

Genome-wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome-wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed- and species-specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.

Genome-wide associations and detection of potential candidate genes for direct genetic and maternal genetic effects influencing dairy cattle body weight at different ages

Background

Body weight (BW) at different ages are of increasing importance in dairy cattle breeding schemes, because of their strong correlation with energy efficiency traits, and their impact on cow health, longevity and farm economy. In total, 15,921 dairy cattle from 56 large-scale test-herds with BW records were genotyped for 45,613 single nucleotide polymorphisms (SNPs). This dataset was used for genome-wide association studies (GWAS), in order to localize potential candidate genes for direct and maternal genetic effects on BW recorded at birth (BW0), at 2 to 3 months of age (BW23), and at 13 to 14 months of age (BW1314).

Results

The first 20 principal components (PC) of the genomic relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G) grouped the genotyped cattle into three clusters. In the statistical models used for GWAS, correction for population structure was done by including polygenic effects with various genetic similarity matrices, such as the pedigree-based relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A), the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G-matrix, the reduced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G-matrix LOCO (i.e. exclusion of the chromosome on which the candidate SNP is located), and LOCO plus chromosome-wide PC. Inflation factors for direct genetic effects using \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A and LOCO were larger than 1.17. For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G and LOCO plus chromosome-wide PC, inflation factors were very close to 1.0. According to Bonferroni correction, ten, two and seven significant SNPs were detected for the direct genetic effect on BW0, BW23, and BW1314, respectively. Seventy-six candidate genes contributed to direct genetic effects on BW with four involved in growth and developmental processes: FGF6, FGF23, TNNT3, and OMD. For maternal genetic effects on BW0, only three significant SNPs (according to Bonferroni correction), and four potential candidate genes, were identified. The most significant SNP on chromosome 19 explained only 0.14% of the maternal de-regressed proof variance for BW0.

Conclusions

For correction of population structure in GWAS, we suggest a statistical model that considers LOCO plus chromosome-wide PC. Regarding direct genetic effects, several SNPs had a significant effect on BW at different ages, and only two SNPs on chromosome 5 had a significant effect on all three BW traits. Thus, different potential candidate genes regulate BW at different ages. Maternal genetic effects followed an infinitesimal model.

Electronic supplementary material

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