Modeling relationships between calving traits: a comparison between standard and recursive mixed models

Background

The use of structural equation models for the analysis of recursive and simultaneous relationships between phenotypes has become more popular recently. The aim of this paper is to illustrate how these models can be applied in animal breeding to achieve parameterizations of different levels of complexity and, more specifically, to model phenotypic recursion between three calving traits: gestation length (GL), calving difficulty (CD) and stillbirth (SB). All recursive models considered here postulate heterogeneous recursive relationships between GL and liabilities to CD and SB, and between liability to CD and liability to SB, depending on categories of GL phenotype.

Methods

Four models were compared in terms of goodness of fit and predictive ability: 1) standard mixed model (SMM), a model with unstructured (co)variance matrices; 2) recursive mixed model 1 (RMM1), assuming that residual correlations are due to the recursive relationships between phenotypes; 3) RMM2, assuming that correlations between residuals and contemporary groups are due to recursive relationships between phenotypes; and 4) RMM3, postulating that the correlations between genetic effects, contemporary groups and residuals are due to recursive relationships between phenotypes.

Results

For all the RMM considered, the estimates of the structural coefficients were similar. Results revealed a nonlinear relationship between GL and the liabilities both to CD and to SB, and a linear relationship between the liabilities to CD and SB.Differences in terms of goodness of fit and predictive ability of the models considered were negligible, suggesting that RMM3 is plausible.

Conclusions

The applications examined in this study suggest the plausibility of a nonlinear recursive effect from GL onto CD and SB. Also, the fact that the most restrictive model RMM3, which assumes that the only cause of correlation is phenotypic recursion, performs as well as the others indicates that the phenotypic recursion may be an important cause of the observed patterns of genetic and environmental correlations.     

Inferring causal structures and comparing the causal effects among calving difficulty,gestation length and calf size in Japanese Black cattle

The objectives of this study were to infer phenotypic causal networks involving gestation length (GL) and calving difficulty (CD) for the primiparity of 1850 Japanese Black heifers, and the birth weight (BWT), withers height (WH) and chest girth (CHG) of their full blood calves, and to compare the causal effects among them. The inductive causation (IC) algorithm was employed to search for causal links among these traits; it was applied to the posterior distribution of the residual (co)variance matrix of a multiple-trait sire-maternal grand sire (MGS) model. The IC algorithm implemented with 95% and 90% highest posterior density intervals detected only one structure with links between GL and BWT (WH or CHG) and between BWT (WH or CHG) and CD, although their directions were not resolved. Therefore, a possible causal structure based on the networks obtained from the IC algorithm [GL→BWT (WH or CHG)→CD] was fitted using a structural equation model to infer causal structure coefficients between the traits. The structural coefficients of GL on BWT and of BWT on GL on the observable scale showed that an extra day of GL led to a 270-g gain in BWT, and a 1-kg increase in BWT increased the risk for dystocia by 1.1%, in the causal structure. Similarly, an increase in GL by 1 day resulted in a 2.1 (2.0)-mm growth in WH (CHG), and a 1-cm increase in WH (CHG) increased the risk of dystocia by 1.2% (0.9%). The structural equation model was also fitted to alternative causal structures, which involved the addition of a directed link from GL to CD, or GL→CD to the structures described above. The inferred structural coefficients with the alternative structures were almost the same as the corresponding ones that had GL→BWT (WH or CHG)→CD. However, the direct causal effect of the extra link from GL on CD was similar to the indirect causal effect of GL through the mediating effect of BWT (WH or CHG) on CD and significant (P<0.05). This suggest that maternal genetic effects might not be removed completely from the residual variance components in the sire-MGS model, and the application of the IC algorithm to the variances from the model could detect an incorrect structure. Nonetheless, fitting the structural equation model to the causal structure provided useful information such as the magnitude of the causal effects between the traits.     

Bayesian estimates of genetic parameters for pre-conception traits, gestation length and calving interval in beef cattle

A total of 5253 records obtained from 2081 Rubia Gallega beef cows managed using artificial insemination as the only reproduction system were analysed to estimate genetic parameters for days to first insemination (DFI), days from first insemination to conception (FIC), number of inseminations per conception (IN), days open (DO), gestation length (GL) and calving interval (CI) via multitrait Bayesian procedures. Estimates of the mean of posterior distribution of the heritability of DFI, FIC, IN, DO, GL and CI were, respectively, 0.050, 0.078, 0.071, 0.053, 0.037 and 0.085 and the corresponding estimates for repeatability of these traits were 0.116, 0.129, 0.147, 0.138, 0.082 and 0.132, respectively. No significant genetic correlations associated to DFI or GL were found. However, genetic correlations between the other four analysed traits were high and significant. Genetic correlations between FIC and IN, DO and CI were similar and higher than 0.85. Genetic correlations of IN-DO and IN-CI were over 0.65. The highest genetic correlation was estimated for the pair DO-CI (0.992) that can be considered the same trait in genetic terms. Results indicated that DFI can be highly affected by non-genetic factors thus limiting its usefulness to be used as an earlier indicator of reproductive performance in beef cattle. Moreover, GL could not be associated to the reproductive performance of the cow before conception. The other four analysed traits, FIC, IN, DO and CI, have close genetic relationships. The inclusion of IN as an earlier indicator of fertility in beef cattle improvement programs using artificial insemination as the main reproductive system can be advisable due to the low additional recording effort needed.     

A whole-genome association study for pig reproductive traits

A whole-genome association study was performed for reproductive traits in commercial sows using the PorcineSNP60 BeadChip and Bayesian statistical methods. The traits included total number born (TNB), number born alive (NBA), number of stillborn (SB), number of mummified foetuses at birth (MUM) and gestation length (GL) in each of the first three parities. We report the associations of informative QTL and the genes within the QTL for each reproductive trait in different parities. These results provide evidence of gene effects having temporal impacts on reproductive traits in different parities. Many QTL identified in this study are new for pig reproductive traits. Around 48% of total genes located in the identified QTL regions were predicted to be involved in placental functions. The genomic regions containing genes important for foetal developmental (e.g. MEF2C) and uterine functions (e.g. PLSCR4) were associated with TNB and NBA in the first two parities. Similarly, QTL in other foetal developmental (e.g. HNRNPD and AHR) and placental (e.g. RELL1 and CD96) genes were associated with SB and MUM in different parities. The QTL with genes related to utero-placental blood flow (e.g. VEGFA) and hematopoiesis (e.g. MAFB) were associated with GL differences among sows in this population. Pathway analyses using genes within QTL identified some modest underlying biological pathways, which are interesting candidates (e.g. the nucleotide metabolism pathway for SB) for pig reproductive traits in different parities. Further validation studies on large populations are warranted to improve our understanding of the complex genetic architecture for pig reproductive traits.     

Genetic parameters for growth, reproductive and maternal traits in a multibreed meat sheep population

The genetic parameters for growth, reproductive and maternal traits in a multibreed meat sheep population were estimated by applying the Average Information Restricted Maximum Likelihood method to an animal model. Data from a flock supported by the Programa de Melhoramento Genético de Caprinos e Ovinos de Corte (GENECOC) were used. The traits studied included birth weight (BW), weaning weight (WW), slaughter weight (SW), yearling weight (YW), weight gain from birth to weaning (GBW), weight gain from weaning to slaughter (GWS), weight gain from weaning to yearling (GWY), age at first lambing (AFL), lambing interval (LI), gestation length (GL), lambing date (LD - number of days between the start of breeding season and lambing), litter weight at birth (LWB) and litter weight at weaning (LWW). The direct heritabilities were 0.35, 0.81, 0.65, 0.49, 0.20, 0.15 and 0.39 for BW, WW, SW, YW, GBW, GWS and GWY, respectively, and 0.04, 0.06, 0.10, 0.05, 0.15 and 0.11 for AFL, LI, GL, LD, LWB and LWW, respectively. Positive genetic correlations were observed among body weights. In contrast, there was a negative genetic correlation between GBW and GWS (-0.49) and GBW and GWY (-0.56). Positive genetic correlations were observed between AFL and LI, LI and GL, and LWB and LWW. These results indicate a strong maternal influence in this herd and the presence of sufficient genetic variation to allow mass selection for growth traits. Additive effects were of little importance for reproductive traits, and other strategies are necessary to improve the performance of these animals.     

粳稻穗角与谷粒性状的相关性及谷粒性状遗传分析

粳稻穗角小的品种产量高而品质差。谷粒性状与产量和品质都有关联。为培育产量与品质相协调的品种提供相关遗传信息, 测定了粳稻直立穗品种丙8979和弯曲穗品种C堡及其重组自交系349个株系的穗角和7个谷粒性状, 分析了穗角与谷粒性状之间的相关性, 并运用主基因+多基因混合遗传模型, 对7个谷粒性状进行了遗传分析。结果表明, 穗角与粒厚、长厚比和宽厚比均无显著相关, 而与千粒重、粒长、粒宽和长宽比均呈极显著正相关。7个谷粒性状均受2对主基因+多基因控制, 2对主基因的作用方式因性状而异。千粒重、粒长、长厚比和宽厚比4个性状以主基因遗传为主; 粒宽、粒厚和长宽比3个性状以多基因遗传为主。     

A bivariate quantitative genetic model for a threshold trait and a survival trait

Many of the functional traits considered in animal breeding can be analyzed as threshold traits or survival traits with examples including disease traits, conformation scores, calving difficulty and longevity. In this paper we derive and implement a bivariate quantitative genetic model for a threshold character and a survival trait that are genetically and environmentally correlated. For the survival trait, we considered the Weibull log-normal animal frailty model. A Bayesian approach using Gibbs sampling was adopted in which model parameters were augmented with unobserved liabilities associated with the threshold trait. The fully conditional posterior distributions associated with parameters of the threshold trait reduced to well known distributions. For the survival trait the two baseline Weibull parameters were updated jointly by a Metropolis-Hastings step. The remaining model parameters with non-normalized fully conditional distributions were updated univariately using adaptive rejection sampling. The Gibbs sampler was tested in a simulation study and illustrated in a joint analysis of calving difficulty and longevity of dairy cattle. The simulation study showed that the estimated marginal posterior distributions covered well and placed high density to the true values used in the simulation of data. The data analysis of calving difficulty and longevity showed that genetic variation exists for both traits. The additive genetic correlation was moderately favorable with marginal posterior mean equal to 0.37 and 95% central posterior credibility interval ranging between 0.11 and 0.61. Therefore, this study suggests that selection for improving one of the two traits will be beneficial for the other trait as well.     

Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein-Friesian dairy cattle

Background

The focus in dairy cattle breeding is gradually shifting from production to functional traits and genetic parameters of calving traits are estimated more frequently. However, across countries, various statistical models are used to estimate these parameters. This study evaluates different models for calving ease and stillbirth in United Kingdom Holstein-Friesian cattle.

Methods

Data from first and later parity records were used. Genetic parameters for calving ease, stillbirth and gestation length were estimated using the restricted maximum likelihood method, considering different models i.e. sire (−maternal grandsire), animal, univariate and bivariate models. Gestation length was fitted as a correlated indicator trait and, for all three traits, genetic correlations between first and later parities were estimated. Potential bias in estimates was avoided by acknowledging a possible environmental direct-maternal covariance. The total heritable variance was estimated for each trait to discuss its theoretical importance and practical value. Prediction error variances and accuracies were calculated to compare the models.

Results and discussion

On average, direct and maternal heritabilities for calving traits were low, except for direct gestation length. Calving ease in first parity had a significant and negative direct-maternal genetic correlation. Gestation length was maternally correlated to stillbirth in first parity and directly correlated to calving ease in later parities. Multi-trait models had a slightly greater predictive ability than univariate models, especially for the lowly heritable traits. The computation time needed for sire (−maternal grandsire) models was much smaller than for animal models with only small differences in accuracy. The sire (−maternal grandsire) model was robust when additional genetic components were estimated, while the equivalent animal model had difficulties reaching convergence.

Conclusions

For the evaluation of calving traits, multi-trait models show a slight advantage over univariate models. Extended sire models (−maternal grandsire) are more practical and robust than animal models. Estimated genetic parameters for calving traits of UK Holstein cattle are consistent with literature. Calculating an aggregate estimated breeding value including direct and maternal values should encourage breeders to consider both direct and maternal effects in selection decisions.     

Genetic evaluation of ovulatory disorders in Austrian Fleckvieh cows: a comparison between linear models and survival analysis

The objective of this study was to compare linear models and survival analysis for genetic evaluation of ovulatory disorders, which included veterinary treatments of silent heat/anestrus and cystic ovaries. Data of 23 450 daughters of 274 Austrian Fleckvieh sires were analyzed. For linear model analyses, ovulatory disorders were defined as a binary response (presence or absence) in the time periods from calving to 150 days after calving and from calving to 300 days after calving. For survival analysis, ovulatory disorders were defined either as the number of days from calving to the day of the first treatment for an ovulatory disorder (uncensored record) or from calving to the day of culling, or the last day of the period under investigation (until 150 or 300 days after calving; censored record). Estimates of heritability were very similar (0.016 to 0.020) across methods and periods. Correlations between sire estimated breeding value from linear model and survival analysis were 0.98, whereas correlations between different time periods were somewhat lower (0.95 and 0.96). The results showed that the length of time period had a larger effect on genetic evaluation than methodology.     

National evaluation for calving ease,gestation length and birth weight by linear and threshold model methodologies

Data included 393,097 calving ease, 129,520 gestation length, and 412,484 birth weight records on 412,484 Gelbvieh cattle. Additionally, pedigrees were available on 72,123 animals. Included in the models were effects of sex and age of dam, treated as fixed, as well as direct, maternal genetic and permanent environmental effects and effects of contemporary group (herd-year-season), treated as random. In all analyses, birth weight and gestation length were treated as continuous traits. Calving ease (CE) was treated either as a continuous trait in a mixed linear model (LM), or as a categorical trait in linear-threshold models (LTM). Solutions in TM obtained by empirical Bayes (TMEB) and Monte Carlo (TMMC) methodologies were compared with those by LM. Due to the computational cost, only 10,000 samples were obtained for TMMC. For calving ease, correlations between LM and TMEB were 0.86 and 0.78 for direct and maternal genetic effects, respectively. The same correlations but between TMEB and TMMC were 1.00 and 0.98, respectively. The correlations between LM and TMMC were 0.85 and 0.75, respectively. The correlations for the linear traits were above.97 between LM and TMEB but as low as 0.91 between LM and TMMC, suggesting insufficient convergence of TMMC. Computing time required was about 2 hrs, 5 hrs, and 6 days for LM, TMEB and TMMC, respectively, and memory requirements were 169, 171, and 445 megabytes, respectively. Bayesian implementation of threshold model is simple, can be extended to multiple categorical traits, and allows easy calculation of accuracies; however, computing time is prohibitively long for large models.     

A genome‐wide association study for gestation length in swine

Selection for increased litter size in swine has potentially resulted in a correlated increase in preweaning mortality. Additional selection criteria should be considered when selecting for increased litter size to account for associated decreases in piglet quality, specifically piglet survival, initial weight and growth. Traits such as gestation length (GL), which have been associated with piglet performance, could be utilized to improve piglet development and survivability. The objective of this study was to conduct a genome‐wide association study to identify genomic regions associated with GL in differing parities in swine (n = 831) from the University of Nebraska–Lincoln reproductive longevity project. Gestation length was calculated as the number of days between last insemination administered and farrowing. Sows were genotyped with the Illumina SNP60 BeadArray, and the data were analyzed using Bayesian mixture models for GL at parity 1, 2, 3 and 4 (GL1, GL2, GL3 and GL4 respectively). Means (SD) for GL1–GL4 were 113 (1.4), 114 (1.2), 114 (1.3) and 115 (1.2) respectively. Posterior mean heritability estimates (PSD) for GL1, GL2, GL3 and GL4 were 0.33 (0.06), 0.34 (0.07), 0.32 (0.08) and 0.20 (0.08) respectively. Rank correlations between genomic estimated breeding values between GL1 and GL2, GL3 and GL4 respectively were moderate: 0.67, 0.65 and 0.60. The top SNP (ASGA0017859, SSC4, 7.8 Mb), located in the top common genomic region associated with GL1, GL2 and GL3, was associated with a difference of 1.1 days in GL1 between homozygote genotypes (P < 0.0001). The results of this study suggest that GL is a largely polygenic trait with relatively minor contributions from multiple genomic regions.     

Genetic and non-genetic effects on productive and reproductive traits of cows in dual-purpose herds in southeastern Mexico

Data on productive and reproductive performance of cows in dual-purpose herds were analyzed to determine the effect of some environmental and genetic factors on saleable milk yield (SMY), lactation length (LL), daily saleable milk per lactation (DMYL), calving interval (CI), and daily saleable milk per calving interval (MYCI) in dual-purpose herds in Yucatan, Mexico. Repeatabilities (re) for these traits were also estimated. Data were obtained from monthly visits to 162 herds from 1996 to 2000. The fixed factors studied were: region (central, eastern and southern), parity number (1 to 6), calving year (1996 to 2000) and calving season (dry, rainy and windy and rainy), genetic group (low- (< 0.50%), medium- (50%) and high- (> 50%) grade cows with European genes). About 2700 to 7700 cows were evaluated for each trait. All factors had significant effects (P < 0.05) on the traits except for region on CI and calving season on DMYL. The overall means for SMY, LL, DMYL, CI, and MYCI were: 1322.3 +/- 80.5 kg, 224.8 +/- 1.3 days, 5.8 +/- 0.1 kg, 555.1 +/- 16.5 days, and 3.0 +/- 0.1 kg, respectively. The re values for SMY, LL, DMYL, CI, and MYCI traits were: 0.19 +/- 0.03, 0.08 +/- 0.04, 0.16 +/- 0.04, 0.00 +/- 0.08, and 0.08 +/- 0.07, respectively. First parity cows had lower SMY, shorter LL, longer CI, and lower MYCI means than cows with more than one parity. Medium grade cows produced more SMY, DMYL and MYCI and had shorter CI than low- and high-grade cows. Therefore, under Yucatan conditions medium-grade cows should be exploited, and more attention should be given to first parity cows in order to improve the productivity in the herd. The relatively high re estimates for SMY and DMYL can be used to calculate most probable producing abilities, in order to identify which cows should be culled.     

The effects of crossbreeding with Norwegian Red dairy cattle on common postpartum diseases,fertility and body condition score

Norwegian Red bulls, selected in Norway, have been used for crossbreeding with Israeli Holstein on commercial farms. The aim of this project was to investigate Norwegian Red×Israeli Holstein (NRX) performance to see how the daughters perform in a different environment than the one their sires were selected in. This was done by comparing health and fertility of NRX with their Israeli Holstein (HO) counterparts. The data consisted of 71 911 HO records and 10 595 NRX records from 33 855 cows in 23 Israeli dairy herds. Calving events took place between 2006 and 2016. Five postpartum disorders (mean frequency in HO v. NRX, %) recorded by veterinarians were analyzed: anestrus (37.4 v. 41.2), metritis (40.1 v. 28.6), ketosis (11.9 v. 7.1), lameness (7.1 v. 3.1) and retained placenta (6.2 v. 4.0). The incidence of abortions was also analyzed; HO had a mean frequency of 9.9% and NRX 8.2%. These traits were defined as binary traits, with ‘1’ indicating that the disorder was present and a treatment took place at least once, or ‘0’ if the cow did not show signs of that disorder. Days open (i.e. the number of days from calving to conception), body condition score (BCS) recorded on a 1 to 5 scale and changes in BCS from calving to peak lactation were also analyzed. A logistic model was used for the health traits, while days open and BCS were analyzed with linear models. The model included breed group, herd-year of calving, birth year and parity as fixed effects. There was a significantly higher risk (odds ratio for HO v. NRX in parentheses) of ketosis (1.46), metritis (1.78), lameness (2.07), retained placenta (1.41) and abortion (1.13) in HO compared with NRX. Israeli Holstein heifers and cows in parity 3 to 6 had fewer cases of anestrus than NRX but no differences were found between the groups in parities 1 and 2. Body condition score was higher for NRX than HO and there was less change in BCS from calving to peak lactation in NRX compared with HO. Likewise, NRX had fewer days open than HO. Results indicate that crossbreeding can produce cows with better fertility that are less susceptible to postpartum disorders.     

Multi-Environmental Genetic Analysis of Grain Size Traits Based on Chromosome Segment Substitution Line in Rice (<i>Oryza sativa</i> L.)

Grain size traits are critical agronomic traits which directly determine grain yield, but the genetic bases of these traits are still not well understood. In this study, a total of 154 chromosome segment substitution lines (CSSLs) population derived from a cross between a japonica variety Koshihikari and an indica variety Nona Bokra was used to investigate grain length (GL), grain width (GW), length-width ratio (LWR), grain perimeter (GP), grain area (GA), and thousand grain weight (TGW) under four environments. QTL mapping analysis of six grain size traits was performed by QTL IciMapping 4.2 with an inclusive composite interval mapping (ICIM) model. A total of 64 QTLs were identified for these traits, which mapped to chromosomes 1, 2, 3, 4, 6, 7, 8, 10, 11, and 12 and accounted for 1.6%–27.1% of the total phenotypic variations. Among these QTLs, thirty-six loci were novel and seven QTLs were identified under four environments. One locus containing the known grain size gene, qGL3/GL3.1/OsPPKL1, also have been found. Moreover, five pairs of digenic epistatic interactions were identified except for GL and GP. These findings will facilitate fine mapping of the candidate gene and QTL pyramiding to genetically improve grain yield in rice.     

Genetic parameters for endocrine and traditional fertility traits,hyperketonemia and milk yield in dairy cattle

High-yielding cows may suffer from negative energy balance during early lactation, which can lead to ketosis and delayed ability of returning to cyclicity after calving. Fast recovery after calving is essential when breeding for improved fertility. Traditionally used fertility traits, such as the interval from calving to first insemination (CFI), have low heritabilities and are highly influenced by management decisions. Herd Navigator™ management program samples and analyses milk progesterone and β-hydroxybutyrate (BHB) automatically during milking. In this study, the genetic parameters of endocrine fertility traits (measured from milk progesterone) and hyperketonemia (measured from milk BHB) in early lactation were evaluated and compared with traditional fertility traits (CFI, interval from calving to the last insemination and interval from first to last insemination) and the milk yield in red dairy cattle herds in Finland. Data included observations from 14 farms from 2014 to 2017. Data were analyzed with linear animal models using DMU software and analyses were done for first parity cows. Heritability estimates for traditional fertility traits were low and varied between 0.03 and 0.07. Estimated heritabilities for endocrine fertility traits (interval from calving to the first heat (CFH) and commencement of luteal activity (C-LA)) were higher than for traditional fertility traits (0.19 to 0.33). Five slightly different hyperketonemia traits divided into two or three classes were studied. Linear model heritability estimates for hyperketonemia traits were low, however, when the threshold model was used for binary traits the estimates became slightly higher (0.07 to 0.15). Genetic correlation between CFH and C-LA for first parity cows was high (0.97) as expected since traits are quite similar. Moderate genetic correlations (0.47 to 0.52) were found between the endocrine fertility traits and early lactation milk yield. Results suggest that the data on endocrine fertility traits measured by automatic systems is a promising tool for improving fertility, specifically when more data is available. For hyperketonemia traits, dividing values into three classes instead of two seemed to work better. Based on the current study and previous studies, where higher heritabilities have been found for milk BHB traits than for clinical ketosis, milk BHB traits are a promising indicator trait for resistance to ketosis and should be studied more. It is important that this kind of data from automatic devices is made available to recording and breeding organizations in the future.     

Random regression analysis for body weights and main morphological traits in genetically improved farmed tilapia (<Emphasis Type="Italic">Oreochromis niloticus</Emphasis>)

To genetically analyse growth traits in genetically improved farmed tilapia (GIFT), the body weight (BWE) and main morphological traits, including body length (BL), body depth (BD), body width (BWI), head length (HL) and length of the caudal peduncle (CPL), were measured six times in growth duration on 1451 fish from 45 mixed families of full and half sibs. A random regression model (RRM) was used to model genetic changes of the growth traits with days of age and estimate the heritability for any growth point and genetic correlations between pairwise growth points. Using the covariance function based on optimal RRMs, the heritabilities were estimated to be from 0.102 to 0.662 for BWE, 0.157 to 0.591 for BL, 0.047 to 0.621 for BD, 0.018 to 0.577 for BWI, 0.075 to 0.597 for HL and 0.032 to 0.610 for CPL between 60 and 140 days of age. All genetic correlations exceeded 0.5 between pairwise growth points. Moreover, the traits at initial days of age showed less correlation with those at later days of age. With phenotypes observed repeatedly, the model choice showed that the optimal RRMs could more precisely predict breeding values at a specific growth time than repeatability models or multiple trait animal models, which enhanced the efficiency of selection for the BWE and main morphological traits.     

Bos indicus or Bos taurus mitochondrial DNA - comparison of productive and reproductive breeding values in a Guzerat dairy herd

The observation of bovine mitochondrial DNA (mtDNA) polymorphisms allows the separation of American zebu cattle, according to its maternal lineage ancestry, into two groups: one with Bos indicus mtDNA and other with Bos taurus mtDNA. The aim of the present study was to determine the productive and reproductive differences between these two groups, in a Guzerat dairy herd. The genotyping of a sample of 56 animals allowed the categorization of most of the 3835 animals in the pedigree file. The production file included 3528 calving and 3198 lactation records from 729 cows, born during the years 1947 to 2007. The traits considered were: lactation milk yield (LMY); days in milk (DIM); age at first calving (AFC), and calving interval (CI). Heritabilities and breeding values were estimated using an animal model. The regression of the average breeding values per year of birth indicated the genetic trends of the herd. The heritability coefficients estimated for LMY, DIM, AFC, and CI were 0.42, 0.43, 0.20, and 0.10, respectively. The genetic trends were similar for both groups, pointing to an improvement in the productive and a worsening in the reproductive traits. The two groups differed significantly regarding the average estimated breeding values for LMY, DIM and AFC, in the starting period, until 1970, but no differences were observed in the more recent years, after 1970. The segregation between the groups existed in the starting period, probably because the Bos taurus contributions to the herd had occurred more recently at that moment. The conclusion is that mtDNA has no significant effect on these traits.     

Agro-Morphological Characterization and Genetic Dissection of Linseed (<i>Linum usitatissimum</i> L.) Genotypes

Linseed is a multipurpose crop and the crop needs further improvement to increase production and yield due to its high value and demand. This study aimed to assess the extent and pattern of genetic variability of forty linseed genotypes based on diverse agro–morphological and yield attributes. The field experiment was conducted following a Randomized Complete Block Design with three replications. Linseed germplasm showed a wide range of phenotypic expression, genetic variability and heritability for 30 studied traits. A low to high phenotypic coeffi- cient of variation (PCV) and genotypic coefficient of variation (GCV) were observed. The lowest genotypic (σ² g) and phenotypic variances (σ² p) were found in capsule diameter (CD), length of calyx (LC), capsule length (CL), seed length (SL), and seed breadth (SB). High broad-sense heritability (h²_b) with high genetic advance as a percentage of mean (GAM) were observed in days to germination started (DGS), days to 80% emergence (DE), plant height at 28 and 40 DAS, number of flowers (NFPP), filled capsules (NFCPP) and yield per plant (YPP) indicating additive gene action exists for these characters. Hierarchical cluster analysis separated 40 genotypes into five clusters, where Clusters I to V assembled with 13, 4, 4, 5 and 14 genotypes, respectively. Considering yield and yield attributes, Cluster-IV (G3, G4, G6, G10 and G31) genotypes showed promising while, Cluster-II (G2, G16, G35, G36) and Cluster-III (G1, G33, G39 and G40) genotypes were dominant on plant morphological traits. Based on principal component analysis (PCA), few characters such as YPP, NFPP, NFCPP, days to first flowering and capsule formation, early emergence, days to branch initiation and plant heights at different growth stages revealed important and effective traits for consideration in the selection of linseed breeding programs.     

Two-trait random regression model to estimate the genetic association of scrotal circumference with female reproductive performance in Nelore cattle

In an attempt to determine when scrotal circumference (SC) could be a reasonable indicator of female reproductive performance, a series of two-trait random regression model (regression for SC on age at measurement) using Gibbs sampling was applied to field data of Nelore cattle raised in a tropical environment. The female traits evaluated were age at first calving (AFC), first calving interval (FCI), heifer pregnancy (HP), and stayability (STAY). The posterior means of heritability of female traits ranged from 0.15 for AFC to 0.46 for HP and were about 0.50 for SC. The posterior means of genetic correlations between SC and AFC, FCI, HP, STAY were up to −0.70, −0.25, 0.48, and 0.29, respectively. Genetically, SC could be a reasonable indicator of female puberty (e.g., HP) as long as it is measured at a young age (400–440 days). However, for female reproductive traits measured at an older age (e.g., STAY), SC is merely a modest or even poor indicator. The use of sire expected progeny differences for female reproductive traits will be more effective than the use of expected progeny differences for SC to improve the reproductive performance of female cattle.     

Genetic relationship between milk urea nitrogen and reproductive performance in Holstein dairy cows

The objective of this study was to describe the genetic and phenotypic relationship between milk urea nitrogen (MUN) and reproductive traits in Iranian Holstein dairy cows. Test-day MUN data obtained from 57 301 dairy cows on 20 large dairy herds in Iran between January 2005 and June 2009. Genetic parameters for MUN and reproductive traits were estimated with a five-trait model using ASREML program. Random regression test-day models were used to estimate heritabilities separately for MUN from first, second and third lactations. Regression curves were modeled using Legendre polynomials of order 3. Herd-year-season along with age at calving was included as fixed effects in all models for reproductive traits. Heritabilities for MUN and reproductive traits were estimated separately for first lactation, second lactation and third lactation. The estimated heritabilities for MUN varied from 0.18 to 0.22. The heritability estimate was low for reproductive traits, which ranged from 0.02 to 0.06 for different traits and across parities. Except for days open, phenotypic and genetic correlations of MUN with reproductive performance traits were close to zero. Genetic correlations between MUN and days open were 0.23, 0.35 and 0.45 in first, second and third lactation, respectively. However, the phenotypic correlation between MUN at different parities was moderate (0.28 to 0.35), but the genetic correlation between MUN at different parities was high and ranged from 0.84 to 0.97. This study shows a limited application of MUN for use in selection programs to improve reproductive performance.