Genomic imprinting variances of beef carcass traits and physiochemical characteristics in Japanese Black cattle

The objective of this study was to estimate variance components related to imprinting for carcass traits and physiochemical characteristics in Japanese Black cattle. The carcass records obtained from 4,220 Japanese Black feedlot cattle included carcass weight (CW), rib eye area (REA), rib thickness, subcutaneous fat thickness, and beef marbling score (BMS), and the physiochemical characteristics were fat, moisture, glycogen per proportion of moisture content, oleic acid, and monounsaturated fatty acids (MUFA). To detect gametic effects, an imprinting model was fitted. High additive heritabilities were estimated for all traits (from 0.516 for glycogen to 0.853 for fat) and were reduced in Mendelian heritability. The range of the differences was from 0.002 (CW) to 0.331 (fat and moisture), and the reductions were due to their imprinting variances. The ratio of the imprinting variance to the total additive genetic variance for REA (0.374), BMS (0.291), fat (0.387), moisture (0.388), and MUFA (0.337) were large (p < 0.05). These imprinting variances were due to the maternal contribution and suggested the existence of maternally expressed genomic imprinting effects on the traits in Japanese Black cattle. Therefore, maternal gametic effects should be considered in breeding programs for Japanese Black cattle.     

Genetic analysis of calf market weight and carcass traits in Japanese Black cattle

Heritabilities of and genetic correlations between additive direct and maternal genetic effects for calf market weight, and additive direct genetic effects for carcass traits, were estimated for Japanese Black cattle by REML procedures under 2-trait animal models. Data were collected from calf and carcass markets in Hyogo and Tottori prefectures and analyzed separately by prefecture. Calf market weight was measured on 42,745 and 23,566 calves in Hyogo and Tottori, respectively. Only the fattening animals with calf market weight were extracted from the carcass database and used for estimation. The carcass traits analyzed were carcass weight, ribeye area, rib thickness, subcutaneous fat thickness, yield estimate, beef marbling score, and 4 meat characters (color, brightness, firmness, and texture). Direct and maternal heritabilities for calf market weight were estimated to be 0.22 and 0.07 in Hyogo, and 0.37 and 0.15 in Tottori, respectively. The estimates of heritabilities for carcass traits were moderate to high in both prefectures. The estimates of direct-maternal genetic correlations for calf market weight were positive (0.17) in Hyogo and negative (-0.63) in Tottori. The direct effect for calf market weight was positively correlated with the direct effect for carcass weight (0.87 and 0.56 in Hyogo and Tottori, respectively) but negatively correlated with the direct effect for beef marbling score (-0.10 in both prefectures). The estimates of genetic correlations between the maternal effect for calf market weight and the direct effects for carcass traits varied from -0.13 to 0.34 in Hyogo and from -0.14 to 0.15 in Tottori. Because direct and maternal genetic effects for early growth traits can be evaluated from calf market weight data in the production system of Japanese Black cattle, this information should be incorporated into selection and mating schemes of the breed.     

Comparison of genetic gains per year for carcass traits among breeding programs in the Japanese Brown and the Japanese Black cattle

The breeding program for beef cattle in Japan has changed dramatically over 4 decades. Visual judging was done initially, but progeny testing in test stations began in 1968. In the 1980s, the genetic evaluation program using field records, so-called on-farm progeny testing, was first adopted in Oita, Hyogo, and Kumamoto prefectures. In this study, genetic trends for carcass traits in these 3 Wagyu populations were estimated, and genetic gains per year were compared among the 3 different beef cattle breeding programs. The field carcass records used were collected between 1988 and 2003. The traits analyzed were carcass weight, LM area, rib thickness, s.c. fat thickness, and beef marbling standard number. The average breeding values of reproducing dams born the same year were used to estimate the genetic trends for the carcass traits. For comparison of the 3 breeding programs, birth years of the dams were divided into 3 periods reflecting each program. Positive genetic trends for beef marbling standard number were clearly shown in all populations. The genetic gains per year for all carcass traits were significantly enhanced by adopting the on-farm progeny testing program. These results indicate that the on-farm progeny testing program with BLUP is a very powerful approach for genetic improvement of carcass traits in Japanese Wagyu beef cattle.     

Genetic relationships between temperament of calves at auction and carcass traits in Japanese Black cattle

Correlations of calves’ temperament with carcass traits were estimated to clarify the genetic relationships between them in Japanese Black cattle. The temperament records for 3128 calves during auction at a calf market were scored on a scale of 1 (calm) to 5 (nervous) as temperament score (TS ), and the TS were divided into two groups (TSG ): TS 1 and 2 comprised TSG 1, and 3 to 5 constituted TSG 2. Carcass data were obtained from 33 552 fattened cattle. A threshold animal model was used for analyzing the underlying liability for TSG , whereas a linear one was used for TS and carcass traits. The heritability estimates for TS and TSG were 0.12 and 0.11, respectively. On the other hand, moderate to high heritability estimates were obtained for carcass traits (0.40 to 0.68). The temperament scores were negatively correlated with carcass weight, rib thickness and subcutaneous fat thickness (?0.13 to ?0.59). In contrast, weak to moderate positive correlations were found between the temperament scores and rib eye area or yield estimate (0.16 to 0.45). The temperament scores and beef marbling score had no correlation. These results showed that it is possible to improve temperament and carcass traits simultaneously.     

Estimation of direct and maternal genetic parameters for growth and carcass traits in a herd of Japanese Black cattle in Miyagi prefecture, using a multitrait animal model

Data from 1170 records of fattening calves were collected on growth and carcass traits from a Japanese Black cattle herd located in Miyagi prefecture, Japan. The objective was to determine direct and maternal heritabilities, direct and maternal genetic correlations and phenotypic correlations between bodyweight at the beginning of the fattening period (BWS), bodyweight at the end of the fattening period (BWF), carcass weight (CW), average daily gain during the fattening period (ADG), rib eye area (REA), rib thickness, subcutaneous backfat thickness (SFT), yield estimate (YE) and beef marbling score (BMS). Direct heritability estimates of 0.16 (SFT) and 0.07 (BMS) were low, whereas estimates of the other traits were medium to high and ranged between 0.44 (REA) and 0.78 (CW). Direct genetic correlations were all positive, except those that were between BWS and SFT, and between BWS and YE (?0.49 and ?0.14, respectively). The lowest positive genetic correlation was between BWS and BMS (0.04) and the highest was between BWF and CW (0.99). The phenotypic correlation coefficients ranged between ?0.41 (between SFT and YE) and 0.96 (between BWF and CW). Maternal heritability estimates were generally low and ranged between 0.00 for BMS and 0.08 for BWS, CW and ADG. Selection programs comprising information on growth and carcass traits of calves and maternal traits of dams were suggested.     

Parent‐of‐origin effects on carcass traits in Japanese Black cattle

Variances caused by the differential expression of paternally and maternally imprinted genes controlling carcass traits in Japanese Black cattle were estimated in this study. Data on marbling score (BMS), carcass weight, rib thickness, rib‐eye area (REA) and subcutaneous fat thickness (SFT) were collected from a total of 13,115 feedlot steers and heifers in a commercial population. A sire–maternal grandsire model was used to analyse the data, and then, imprinting parameters were derived by replacing the genetic effect of the dam with the effect of the maternal grandsire in the imprinting model to calculate the genetic parameter estimates. The proportions of the total genetic variance attributable to imprinted genes ranged from 8.7% (SFT) to 35.2% (BMS). The remarkably large imprinting variance of BMS was mainly contributed by maternally expressed inheritance because the maternal contribution of the trait was much larger than that of the paternal trait. The parent‐of‐origin effect originating from maternal gene expression was also observed for REA. The results suggested the existence of genomic imprinting effects on the traits of the Japanese Black cattle. Hence, the parent‐of‐origin effect should be considered for the genetic evaluation of Japanese Black cattle in breeding programmes.     

Genetic parameters of serum vitamin A and total cholesterol concentrations and the genetic relationships with carcass traits in an F1 cross between Japanese Black sires and Holstein dams

Blood sample and carcass trait data were collected from 841 and 3,219 fattening animals of an F(1) cross between Japanese Black sires and Holstein dams, respectively. Data on serum vitamin A and total cholesterol concentrations containing 582 to 739 records at 4 stages of beef-fattening production were used to estimate variance components and heritabilities for the 2 traits, based on fattening periods: less than 13 mo of age (stage 1), 14 to 18 mo of age (stage 2), 19 to 21 mo of age (stage 3), and greater than 22 mo of age (stage 4). Furthermore, the genetic correlations of serum vitamin A and total cholesterol concentrations with beef marbling standards, carcass weight, ribeye area, rib thickness, and subcutaneous fat thickness were estimated. The heritability estimates of serum vitamin A concentration were consistently and drastically decreased from 0.37 ± 0.15 to 0.07 ± 0.07 from stages 1 to 4 because of considerable decreases in sire variances, whereas the residual variances remained large and stable throughout all 4 stages. Serum total cholesterol concentration was moderately heritable (approximately 0.35 to 0.64) throughout all 4 stages. The genetic correlations among serum vitamin A concentrations and beef marbling standards were high and negative (-0.94) in stage 4. The genetic correlations between serum vitamin A concentration and carcass weight in stages 3 and 4 were moderate and positive (approximately 0.26 to 0.36). Moderate to high positive genetic correlations between serum vitamin A concentration and subcutaneous fat thickness were obtained throughout the stages (approximately 0.40 to 0.75). Genetic correlations of serum total cholesterol concentration with carcass weight and rib thickness were moderate and positive (approximately 0.29 to 0.46) in stages 2, 3, and 4. These results indicate serum vitamin A and total cholesterol concentrations could be effective physiological indicators for improving carcass traits.     

Effects of parent‐of‐origin models with different pedigree information on beef carcass traits and fatty acid composition in Japanese Black cattle

Genomic imprinting should be considered in animal breeding systems to avoid lead in bias in genetic parameter estimation. The objective of this study was to clarify the effects of pedigree information on imprinting variances for carcass traits and fatty acid composition in Japanese Black cattle. Carcass records [carcass weight, rib eye area, rib thickness (RT), subcutaneous fat thickness and beef marbling score (BMS)] and fatty acid composition were obtained for 11,855 Japanese Black feedlot cattle. To estimate and compare the imprinting variances for the traits, two imprinting models with different pedigree information [the sire–dam gametic relationship matrix (Model 1) and the sire–maternal grandsire (MGS) numerator relationship matrix (Model 2)] were fitted. The ratio of the imprinting variance to the total additive genetic variance for RT (6.33%) and BMS (19.00%) was significant in Model 1, but only that for BMS (21.09%) was significant in Model 2. This study revealed that fitting the sire–MGS model could be useful in estimating imprinting variance under certain conditions, such as when restricted pedigree information is available. Furthermore, the present result suggested that the maternal gametic effects on BMS should be included in breeding programmes for Japanese Black cattle to avoid selection bias caused by imprinting effects.     

Genetic parameters of postnatal mortality and birth weight in Japanese Black calves

Postnatal mortality (PM) of calves is an economically important trait, because a high mortality rate of calves decreases farm income and significantly increases production costs. In spite of the economic importance of PM, genetic studies on PM have been scarce for Japanese Black cattle. The objectives of this study were to estimate the direct and maternal heritability for PM using a Bayesian threshold model with Gibbs sampling in Japanese Black cattle. The original data consisted of 43 723 records of animals born from January 2005 to June 2006. Four binary traits of mortality were considered: D1‐14, D15‐60, D61‐180 and D1‐180 with numbers indicating the period of risk in days after birth. Direct heritabilities for PM ranged from 0.12 to 0.31 and maternal heritabilities for PM were lower than the direct heritabilities and ranged from 0.02 to 0.16. Direct genetic correlations between D1‐14 and D15‐60, between D15‐60 and D61‐180, and between D1‐14 and D61‐180 were 0.501, 0.689 and 0.062, respectively. Direct genetic correlations of birth weight with D1‐14, D15‐60, D61‐180 and D1‐180 were 0.089, 0.482, 0.434 and 0.394, respectively. These results indicated that different genes were responsible for early and late PM and early PM (D1‐14) can be genetically reduced without influencing birth weight.     

Estimation of variance components for carcass traits in Japanese Black cattle using 50K SNP genotype data

Genomic selection using high‐density single nucleotide polymorphism (SNP) genotype data may accelerate genetic improvements in livestock animals. In this study, we attempted to estimate the variance components of six carcass traits in fattened Japanese Black steers using SNP genotype data. Six hundred and seventy‐three steers were genotyped using an Illumina Bovine SNP50 BeadChip and phenotyped for cold carcass weight, ribeye area, rib thickness, subcutaneous fat thickness, estimated yield percent and marbling score. Additive polygenic variance and the variance attributable to a set of SNPs that had statistically significant effects on the trait were estimated via Gibbs sampling with two models: (i) a model with the chosen SNPs and the additive polygenic effects; and (ii) a model with the polygenic effects alone. The proportion of the estimated variance attributable to the SNPs became higher as the number of SNP effects that fit increased. High correlations between breeding values estimated with the model containing the polygenic effect alone and those estimated by chosen SNPs were obtained. No fraction of the total genetic variance was explained by SNPs associated with the trait at P ≥ 0.1. Our results suggest that for the carcass traits of Japanese Black cattle, a maximum of half of the total additive genetic variance may be explained by SNPs between 100 several tens to several 100s.     

A genetic analysis of meat compositions in Japanese Black cattle: Genetic parameters and sex influence

Meat composition in beef is related to eating quality and food functionality. Genetic parameters for several meat compositions including free amino acid, peptide and sugar, however, remain poorly described. In this study, we estimated genetic parameters for 51 meat components, including free amino acids, peptides, sugars and fatty acid compositions, and two carcase traits in 1,354 heifers and 1,797 steers of Japanese Black cattle. Heritability estimates were generally equivalent to or moderately greater than those in previous studies of this breed. Genetic correlations between free amino acids, peptides and sugars and carcase traits were often negative, suggesting a trade‐off between traits. Using two‐trait animal models that treat records from the two sexes as different traits, we estimated sex‐specific heritabilities and cross‐sex genetic correlations which indicate the sex differences in genetic architecture. In these analyses, 12 traits showed significant heritability differences between sexes and cross‐sex genetic correlations occasionally deviated from unity. These results could be used to inform future breeding schemes and investigations of the genetic architecture of meat compositions in beef.     

Genetic relationship of body measurement traits at early age with carcass traits in Japanese black cattle

Estimates of genetic parameters were obtained for body measurement traits of 648 animals at 4 months of age, of 545 at 8 months and carcass traits of 14 972 animals with the use of an animal model by the restricted maximum likelihood procedure. The estimated heritabilities for carcass traits were high (0.41 to 0.54). At 4 months the estimated direct heritabilities for body measurement traits were moderate to high (0.28 to 0.64), except for chest width (0.19); at 8 months they were also moderate to high (0.23 to 0.49), except for chest depth and chest width (0.18 and 0.06, respectively). Maternal heritabilities for all body measurement traits were low at both ages. The results indicate that because of their moderate direct genetic correlations with body measurement traits, carcass weight, rib thickness and subcutaneous fat thickness can be improved; however, rib eye area and beef marbling standard show little such possibility considering their correlation with body measurement traits.     

Estimation of genetic parameters of defective appearances in Japanese Black heifer calves in Kagoshima

Defective appearances, including white spotting (WS), tongue defect (TD), and nipple defect (ND), in Japanese Black cattle potentially lead to economic losses to farmers in Japan. We estimated genetic parameters of defective appearances using 553,433 records of Japanese Black heifer calves housed in the Kagoshima Prefecture. Variance and covariance were estimated using the Gibbs sampling algorithm. The estimated heritability ranged from 0.29 for TD to 0.76 for WS. Percent breeding value (%BV) estimates indicated high variation in WS and ND among sires, reflecting higher heritability. Furthermore, there was a positive linear relationship between the %BV estimate of a sire and the mean incidence rate of each defect in his female offsprings. TD was positively associated with other defects. Therefore, genetic factors strongly affect the incidence of defective appearances in Japanese Black cattle.     

Genetic analysis of residual feed intakes and other performance test traits of Japanese Black cattle from revised protocol

The performance test protocol for Japanese Black cattle was revised in April 2002. This resulted in restriction of access to concentrate (based on body weight) and modification of the concentrate's ingredients. Genetic parameters of growth and feed utilization traits of the performance test were estimated using 1304 records using the revised protocol. Residual feed intakes (RFIs) as alternative indicators for feed utilization efficiency were included. (Co)variance components were estimated by EM-REML. Heritabilities for growth traits were between 0.26 and 0.47. Heritabilities for feed intakes and conversions ranged from 0.25 to 0.37 and from 0.03 to 0.29, respectively. Genetic variances and heritabilities were lower for the revised protocol. Highly positive genetic correlations of daily gain (DG) with feed intakes indicated selection on DG is expected to increase feed intake. Selection on feed conversion may lead to higher DG. The heritability estimates for RFIs ranged from 0.10 to 0.33 and were generally higher than corresponding estimates for feed conversion ratios. RFI of TDN showed positive genetic correlations with all feed intakes. The reduction of feed intakes could be expected through selection on the RFI without changing body size. RFIs were considered to be alternative indicators to improve feed utilization efficiency under the new performance test.     

Genotype x environment interaction effects on carcass traits in Japanese Black cattle

The importance of genotype x environment (region or management system) interactions for carcass traits in Japanese Black cattle was investigated using both univariate and multivariate animal models. The univariate approach was used mainly to test the significance of interaction effects. The multivariate approach was used to estimate genetic correlations, which indicated the magnitude of genotype x environment (GE) interactions. The more a genetic correlation deviates from 1, the larger the interaction. From the univariate approach, the addition of genotype x environment (region or management system) interaction (co)variance components resulted in an improved fit of the model for all traits in both cases (P < 0.001). However, estimates of genetic correlation between regions obtained from the multivariate approach for hot carcass weight, LM area, rib thickness, s.c. fat thickness, and marbling score were 0.97, 0.95, 0.93, 0.97, and 0.93, respectively. The corresponding estimates between management systems were 0.84, 0.92, 0.84, 0.90, and 0.97, respectively. These results indicate that GE interaction effects on carcass traits of Japanese Black cattle may be biologically unimportant. Therefore, breeding values obtained using the multivariate method would rank sires similarly in all environments. Consequently, carcass traits measured in these two different regions or management systems can be treated as the same traits.     

Estimated genetic parameters for carcass traits of Brahman cattle

Heritabilities and genetic and phenotypic correlations were estimated from feedlot and carcass data collected from Brahman calves (n = 504) in central Florida from 1996 to 2000. Data were analyzed using animal models in MTDFREML. Models included contemporary group (n = 44; groups of calves of the same sex, fed in the same pen, slaughtered on the same day) as a fixed effect and calf age in days at slaughter as a continuous variable. Estimated feedlot trait heritabilities were 0.64, 0.67, 0.47, and 0.26 for ADG, hip height at slaughter, slaughter weight, and shrink. The USDA yield grade estimated heritability was 0.71; heritabilities for component traits of yield grade, including hot carcass weight, adjusted 12th rib backfat thickness, loin muscle area, and percentage kidney, pelvic, and heart fat were 0.55, 0.63, 0.44, and 0.46, respectively. Heritability estimates for dressing percentage, marbling score, USDA quality grade, cutability, retail yield, and carcass hump height were 0.77, 0.44, 0.47, 0.71, 0.5, and 0.54, respectively. Estimated genetic correlations of adjusted 12th rib backfat thickness with ADG, slaughter weight, marbling score, percentage kidney, pelvic, and heart fat, and yield grade (0.49, 0.46, 0.56, 0.63, and 0.93, respectively) were generally larger than most literature estimates. Estimated genetic correlations of marbling score with ADG, percentage shrink, loin muscle area, percentage kidney, pelvic, and heart fat, USDA yield grade, cutability, retail yield, and carcass hump height were 0.28, 0.49, 0.44, 0.27, 0.45, -0.43, 0.27, and 0.43, respectively. Results indicate that sufficient genetic variation exists within the Brahman breed for design and implementation of effective selection programs for important carcass quality and yield traits.     

Genetic parameters for gestation length and the relationship with birth weight and carcass traits in Japanese Black cattle

The objectives of this study were to estimate genetic parameters for gestation length (GL), including estimation of maternal effects, and to investigate the genetic relationships of GL with birth weight and carcass traits in a Japanese Black cattle population. The original data comprised 34 775 records of animals born from October 1999 to August 2003. Two different models were used to analyze the data for GL. The first model (M1) included direct genetic effect of the calf and maternal genetic effect as random effects. The second model (M2) treated GL as a trait of the dam and included direct genetic effects only. M1 was used in bi-variate analysis. The direct and maternal heritabilities for GL estimated from M1 were 0.53 and 0.14, respectively. This result shows that GL is moderately inherited and can be controlled genetically. The direct × maternal genetic correlation for GL was −0.73. Direct genetic correlations of GL with carcass traits were close to zero. However, genetic correlation of maternal GL with carcass weight was moderate (0.25).     

荷斯坦牛皮肤褶皱厚度及体况评分性状遗传分析

本研究通过对荷斯坦牛颈侧部、肋部和后乳房基部皮肤褶皱厚度（分别简称为颈部、肋部和乳房皮褶厚）及体况评分（body condition score,BCS）进行大群测定,旨在探究荷斯坦牛不同部位皮肤褶皱厚度和BCS的群体特征,并建模估计荷斯坦牛不同部位皮肤褶皱厚度及BCS的遗传参数。本研究于2015—2020年夏季测定了北京地区12个规模化牧场10 915头泌乳荷斯坦牛的皮肤褶皱厚度,同时对测定牛群进行了体况评分,使用单性状和四性状动物模型分别对颈部皮褶厚、肋部皮褶厚、乳房皮褶厚和BCS进行遗传分析并估计遗传参数,并计算了上述4个性状与奶牛寿命和繁殖等重要功能性状之间的近似遗传相关。结果显示,荷斯坦牛乳房皮褶厚、颈部皮褶厚、肋部皮褶厚和体况评分分别为（7.49±1.45） mm、（7.27±1.34） mm、（11.74±1.90） mm和（2.94±0.79）;乳房皮褶厚和颈部皮褶厚为中等遗传力性状,遗传力分别为0.11和0.15;肋部皮褶厚和体况评分为中高遗传力性状,遗传力分别为0.28和0.22;不同部位的皮褶厚之间存在中高遗传相关,颈部和肋部皮褶厚之间的遗传相关最高（0.69）,颈部和乳房皮褶厚之间的遗传相关最低（0.21）;BCS与不同部位皮褶厚之间遗传相关的差异较大,BCS与颈部、肋部皮褶厚之间为正遗传相关（0.28、0.13）,BCS与乳房皮褶厚之间为负遗传相关（-0.25）;体况评分和皮褶厚性状与部分重要功能性状之间存在中等遗传相关。本研究对奶牛3个部位的皮肤褶皱厚度和BCS性状进行了遗传分析,研究获得的遗传参数有助于通过皮肤褶皱厚度和BCS理解奶牛体脂沉积的遗传基础,助力我国奶牛的平衡育种。     

基于系谱和基因组信息估计荷斯坦青年母牛体重性状遗传参数

旨在设计利用不同信息来源的模型估计荷斯坦后备牛不同月龄体重性状的遗传参数。本研究于2014—2020年测定并收集了7 122头荷斯坦牛32 338条0～12月龄体重数据,分别利用系谱信息(linear mixed model with pedigree relationship matrix, LM＿A)和系谱-基因组信息构建亲缘关系矩阵(linear mixed model with genotype-pedigree joint relationship matrix, LM＿H),基于母体效应动物模型估计初生重,基于是否考虑初生重作为协变量的单性状动物模型估计2～12月龄各月龄体重遗传力,并利用双性状动物模型估计初生重与其它月龄体重的遗传相关。结果显示,对于初生重,根据赤池信息量准则(Akaike information criterion, AIC),LM＿H方法的拟合程度显著优于LM＿A方法,但两种方法估计的遗传参数相差不大：直接遗传力分别为0.30和0.32,母体遗传力分别为0.08和0.09,个体直接遗传效应和母体遗传效应遗传相关系数分别为-0.65和-0.64;对于2～...     

Threshold-linear versus linear-linear analysis of birth weight and calving ease using an animal model: I. Variance component estimation.

Birth weight and calving difficulty were analyzed with Bayesian methodology using univariate linear models, a bivariate linear model, a threshold model for calving difficulty, and a joint threshold-linear model using a probit approach. Field data included 26,006 records of Gelbvieh cattle. Simulated populations were generated using parameters estimated from the field data. The Gibbs sampler was used to obtain estimates of the marginal posterior mean and standard deviation of the (co)variance components, heritabilities, and correlations. In the univariate analyses, the posterior mean of direct heritability for calving difficulty was .23 with the threshold model and .18 with the linear model. Maternal heritabilities were .10 and .08, respectively. In the bivariate analysis, posterior means of direct heritability for calving difficulty were .21 and .18 for the bivariate linear-threshold and linear-linear model, respectively. Maternal heritabilities were .09 and .06, respectively. Direct heritability for birth weight was .25 for the univariate model and .26 for bivariate models. Maternal heritability was .05 for the linear-threshold model and the univariate model and .06 for the bivariate linear model. Genetic correlation between direct genetic effects in both traits was .81 for the linear-threshold model and .79 for the bivariate linear. Residual correlation was .35 for the bivariate linear model and .50 for the bivariate linear-threshold. A simulation study confirmed that the posterior mean of the marginal distribution was suitable as a point estimate for univariate threshold and bivariate linear-threshold models.