Exploring definitions of daily enteric methane emission phenotypes for genetic evaluations using a population of indoor-fed multi-breed growing cattle with feed intake data.

Genetic selection has been identified as a promising approach for reducing enteric methane (CH4) emissions; a prerequisite for genetic evaluations; however, these are estimates of the necessary genetic parameters based on a population representative of where the genetic evaluations will be used. The objective of this study was, therefore, to derive genetic parameters for a series of definitions of CH4, carbon dioxide (CO2), and dry matter intake (DMI) as well as genetic correlations between CH4, CO2, and DMI in a bid to address the paucity of studies involving methane emissions measured in beef cattle using GreenFeed systems. Lastly, estimated breeding values (EBV) were generated for nine alternative definitions of CH4 using the derived genetic parameters; the EBV were validated against both phenotypic performance (adjusted for non-genetic effects) and the Legarra and Reverter method comparing EBV generated for a subset of the dataset compared to EBV generated from the entire dataset. Individual animal CH4 and CO2 records were available from a population of 1,508 multi-breed growing beef cattle using 10 GreenFeed Emission Monitoring systems. Nine trait definitions for CH4 and CO2 were derived: individual spot measures, the average of all spot measures within a 3-h, 6-h, 12-h, 1-d, 5-d, 10-d, and 15-d period and the average of all spot measures across the full test period (20 to 114 d on test). Heritability estimates from 1,155 animals, for CH4, increased as the length of the averaging period increased and ranged from 0.09 ±â€…0.03 for the individual spot measures trait to 0.43 ±â€…0.11 for the full test average trait; a similar trend existed for CO2 with the estimated heritability ranging from 0.17 ±â€…0.04 to 0.50 ±â€…0.11. Enteric CH4 was moderately to strongly genetically correlated with DMI with a genetic correlation of 0.72 ±â€…0.02 between the spot measures of CH4 and a 1-d average DMI. Correlations, adjusted for heritability, between the adjusted phenotype and (parental average) EBV ranged from 0.56 to 1.14 across CH4 definitions and the slope between the adjusted phenotype and EBV ranged from 0.92 to 1.16 (expectation = 1). Validation results from the Legarra and Reverter regression method revealed a level bias of between -0.81 and -0.45, a dispersion bias of between 0.93 and 1.17, and ratio accuracy (ratio of the partial evaluation accuracies on whole evaluation accuracies) from 0.28 to 0.38. While EBV validation results yielded no consensus, CH4 is a moderately heritable trait, and selection for reduced CH4 is achievable.

Livestock production is a significant contributor to greenhouse gas emissions. Animal breeding programs have been proposed as a sustainable mitigation strategy to reduce enteric methane emissions in livestock production. Before creating a genetic evaluation for enteric methane production, it is important to estimate how much inter-animal genetic variability contributes to the observed differences in enteric methane production. The purpose of this study was to explore multiple enteric methane phenotypes and estimate how much phenotypic variation was due to genetic differences among 1,508 growing cattle of multiple breeds and crosses; also of interest was the extent of similarity in the genetic control of enteric methane, carbon dioxide, and feed intake (i.e., the genetic correlation) and to determine if selection of animals on the estimated genetic merit for methane emissions of their parents would manifest itself in differences in actual methane produced by those animals. Between 9% and 43% of the inter-animal differences in daily enteric methane production were due to differences in the genetic composition of those animals; the genetic control influencing methane production was similar to that of feed intake (i.e., a strong genetic correlation between methane emissions and feed intake of up to 0.72).

Adenylate kinase 9 is essential for sperm function and male fertility in mammals.

Despite passing routine laboratory tests for semen quality, bulls used in artificial insemination exhibit significant variation in fertility. Routine analysis of fertility data identified a dairy bull with extreme subfertility (10% pregnancy rate). To characterize the subfertility phenotype, a range of in vitro, in vivo, and molecular assays were carried out. Sperm from the subfertile bull exhibited reduced motility and severely reduced caffeine-induced hyperactivation compared to controls. Ability to penetrate the zona pellucida, cleavage rate, cleavage kinetics, and blastocyst yield after IVF or AI were significantly lower than in control bulls. Whole-genome sequencing from semen and RNA sequencing of testis tissue revealed a critical mutation in adenylate kinase 9 (AK9) that impaired splicing, leading to a premature termination codon and a severely truncated protein. Mice deficient in AK9 were generated to further investigate the function of the gene; knockout males were phenotypically indistinguishable from their wild-type littermates but produced immotile sperm that were incapable of normal fertilization. These sperm exhibited numerous abnormalities, including a low ATP concentration and reduced motility. RNA-seq analysis of their testis revealed differential gene expression of components of the axoneme and sperm flagellum as well as steroid metabolic processes. Sperm ultrastructural analysis showed a high percentage of sperm with abnormal flagella. Combined bovine and murine data indicate the essential metabolic role of AK9 in sperm motility and/or hyperactivation, which in turn affects sperm binding and penetration of the zona pellucida. Thus, AK9 has been found to be directly implicated in impaired male fertility in mammals.

Accelerated matrix-vector multiplications for matrices involving genotype covariates with applications in genomic prediction.

In the last decade, a number of methods have been suggested to deal with large amounts of genetic data in genomic predictions. Yet, steadily growing population sizes and the suboptimal use of computational resources are pushing the practical application of these approaches to their limits. As an extension to the C/CUDA library miraculix, we have developed tailored solutions for the computation of genotype matrix multiplications which is a critical bottleneck in the empirical evaluation of many statistical models. We demonstrate the benefits of our solutions at the example of single-step models which make repeated use of this kind of multiplication. Targeting modern Nvidia® GPUs as well as a broad range of CPU architectures, our implementation significantly reduces the time required for the estimation of breeding values in large population sizes. miraculix is released under the Apache 2.0 license and is freely available at https://github.com/alexfreudenberg/miraculix.

Efficient large-scale single-step evaluations and indirect genomic prediction of genotyped selection candidates.

BACKGROUND: Single-step genomic best linear unbiased prediction (ssGBLUP) models allow the combination of genomic, pedigree, and phenotypic data into a single model, which is computationally challenging for large genotyped populations. In practice, genotypes of animals without their own phenotype and progeny, so-called genotyped selection candidates, can become available after genomic breeding values have been estimated by ssGBLUP. In some breeding programmes, genomic estimated breeding values (GEBV) for these animals should be known shortly after obtaining genotype information but recomputing GEBV using the full ssGBLUP takes too much time. In this study, first we compare two equivalent formulations of ssGBLUP models, i.e. one that is based on the Woodbury matrix identity applied to the inverse of the genomic relationship matrix, and one that is based on marker equations. Second, we present computationally-fast approaches to indirectly compute GEBV for genotyped selection candidates, without the need to do the full ssGBLUP evaluation. RESULTS: The indirect approaches use information from the latest ssGBLUP evaluation and rely on the decomposition of GEBV into its components. The two equivalent ssGBLUP models and indirect approaches were tested on a six-trait calving difficulty model using Irish dairy and beef cattle data that include 2.6 million genotyped animals of which about 500,000 were considered as genotyped selection candidates. When using the same computational approaches, the solving phase of the two equivalent ssGBLUP models showed similar requirements for memory and time per iteration. The computational differences between them were due to the preprocessing phase of the genomic information. Regarding the indirect approaches, compared to GEBV obtained from single-step evaluations including all genotypes, indirect GEBV had correlations higher than 0.99 for all traits while showing little dispersion and level bias. CONCLUSIONS: In conclusion, ssGBLUP predictions for the genotyped selection candidates were accurately approximated using the presented indirect approaches, which are more memory efficient and computationally fast, compared to solving a full ssGBLUP evaluation. Thus, indirect approaches can be used even on a weekly basis to estimate GEBV for newly genotyped animals, while the full single-step evaluation is done only a few times within a year.

A breed-of-origin of alleles model that includes crossbred data improves predictive ability for crossbred animals in a multi-breed population.

BACKGROUND: Recently, crossbred animals have begun to be used as parents in the next generations of dairy and beef cattle systems, which has increased the interest in predicting the genetic merit of those animals. The primary objective of this study was to investigate three available methods for genomic prediction of crossbred animals. In the first two methods, SNP effects from within-breed evaluations are used by weighting them by the average breed proportions across the genome (BPM method) or by their breed-of-origin (BOM method). The third method differs from the BOM in that it estimates breed-specific SNP effects using purebred and crossbred data, considering the breed-of-origin of alleles (BOA method). For within-breed evaluations, and thus for BPM and BOM, 5948 Charolais, 6771 Limousin and 7552 Others (a combined population of other breeds) were used to estimate SNP effects separately within each breed. For the BOA, the purebreds' data were enhanced with data from ~ 4K, ~ 8K or ~ 18K crossbred animals. For each animal, its predictor of genetic merit (PGM) was estimated by considering the breed-specific SNP effects. Predictive ability and absence of bias were estimated for crossbreds and the Limousin and Charolais animals. Predictive ability was measured as the correlation between PGM and the adjusted phenotype, while the regression of the adjusted phenotype on PGM was estimated as a measure of bias. RESULTS: With BPM and BOM, the predictive abilities for crossbreds were 0.468 and 0.472, respectively, and with the BOA method, they ranged from 0.490 to 0.510. The performance of the BOA method improved as the number of crossbred animals in the reference increased and with the use of the correlated approach, in which the correlation of SNP effects across the genome of the different breeds was considered. The slopes of regression for PGM on adjusted phenotypes for crossbreds showed overdispersion of the genetic merits for all methods but this bias tended to be reduced by the use of the BOA method and by increasing the number of crossbred animals. CONCLUSIONS: For the estimation of the genetic merit of crossbred animals, the results from this study suggest that the BOA method that accommodates crossbred data can yield more accurate predictions than the methods that use SNP effects from separate within-breed evaluations.

Phenotypic relationship and repeatability of methane emissions and performance traits in beef cattle using a GreenFeed system.

Rumen methanogenesis results in the loss of 6% to 10% of gross energy intake in cattle and globally is the single most significant source of anthropogenic methane (CH4) emissions. The purpose of this study was to analyze greenhouse gas traits recorded in a commercial feedlot unit to gain an understanding into the relationships between greenhouse gas traits and production traits. Methane and carbon dioxide (CO2) data recorded via multiple GreenFeed Emission Monitoring (GEM), systems as well as feed intake, live weight, ultrasound scanning data, and slaughter data were available on 1,099 animals destined for beef production, of which 648 were steers, 361 were heifers, and 90 were bulls. Phenotypic relationships between GEM emission measurements with feed intake, weight traits, muscle ultrasound data, and carcass traits were estimated. Utilization of GEM systems, daily patterns of methane output, and repeatability of GEM system measurements across averaging periods were also assessed. Methane concentrations varied with visit number, duration, and time of day of visit to the GEM system. Mean CH4 and CO2 varied between sex, with mean CH4 of 256.1 g/day ± 64.23 for steers, 234.7 g/day ± 59.46 for heifers, and 156.9 g/day ± 55.98 for young bulls. A 10-d average period of GEM system measurements were required for steers and heifers to achieve a minimum repeatability of 0.60; however, higher levels of repeatability were observed in animals that attended the GEM system more frequently. In contrast, CO2 emissions reached repeatability estimates >0.6 for steers and heifers in all averaging periods greater than 2-d, suggesting that cattle have a moderately consistent CO2 emission pattern across time periods. Animals with heavier bodyweights were observed to have higher levels of CH4 (correlation = 0.30) and CO2 production (correlation = 0.61), and when assessing direct methane, higher levels of dry matter intake were associated with higher methane output (correlation = 0.31). Results suggest that reducing CH4 can have a negative impact on growth and body composition of cattle. Methane ratio traits, such as methane yield and intensity were also evaluated, and while easy to understand and compare across populations, ratio traits are undesirable in animal breeding, due to the unpredictable level of response. Methane adjusted for dry matter intake and liveweight (Residual CH4) should be considered as an alternative emission trait when selecting for reduced emissions within breeding goals.

Methane production from cattle digestion results in the loss of 6% to 10% of gross energy intake in cattle and globally is the single most significant source of anthropogenic methane (CH4) emissions. The purpose of this study was to analyze greenhouse gas traits recorded in a commercial feedlot unit to gain an understanding into the relationships between greenhouse gas traits and production traits of economic importance. Methane and carbon dioxide emissions recorded using Greenfeed systems were available on a total of 1,099 animals. In addition, performance indicators such as feed intake, live weight, ultrasound scanning data, and slaughter data were also available on all animals. Phenotypic repeatability of CH4 ranged from 0.13 to 0.74, with a CH4 repeatability of >0.6 achieved by both heifers and steers in 10-d measuring period. Due to the high repeatability of CH4 measures, an accurate portrayal of CH4 production can be observed from a 10-d measuring period when measures are averaged. Methane emission data were positively correlated with traits of economic importance. Phenotypically, animals with heavier body weights and greater feed intake had higher emissions.

Detection of Genomic Imprinting for Carcass Traits in Cattle Using Imputed High-Density Genotype Data.

Genomic imprinting is an epigenetic phenomenon defined as the silencing of an allele, at least partially, at a given locus based on the sex of the transmitting parent. The objective of the present study was to detect the presence of SNP-phenotype imprinting associations for carcass weight (CW), carcass conformation (CC) and carcass fat (CF) in cattle. The data used comprised carcass data, along with imputed, high-density genotype data on 618,837 single nucleotide polymorphisms (SNPs) from 23,687 cattle; all animal genotypes were phased with respect to parent of origin. Based on the phased genotypes and a series of single-locus linear models, 24, 339, and 316 SNPs demonstrated imprinting associations with CW, CC, and CF, respectively. Regardless of the trait in question, no known imprinted gene was located within 0.5 Mb of the SNPs demonstrating imprinting associations in the present study. Since all imprinting associations detected herein were at novel loci, further investigation of these regions may be warranted. Nonetheless, knowledge of these associations might be useful for improving the accuracy of genomic evaluations for these traits, as well as mate allocations systems to exploit the effects of genomic imprinting.

The Association Between Genomic Heterozygosity and Carcass Merit in Cattle.

The objective of the present study was to quantify the association between both pedigree and genome-based measures of global heterozygosity and carcass traits, and to identify single nucleotide polymorphisms (SNPs) exhibiting non-additive associations with these traits. The carcass traits of interest were carcass weight (CW), carcass conformation (CC) and carcass fat (CF). To define the genome-based measures of heterozygosity, and to quantify the non-additive associations between SNPs and the carcass traits, imputed, high-density genotype data, comprising of 619,158 SNPs, from 27,213 cattle were used. The correlations between the pedigree-based heterosis coefficient and the three defined genomic measures of heterozygosity ranged from 0.18 to 0.76. The associations between the different measures of heterozygosity and the carcass traits were biologically small, with positive associations for CW and CC, and negative associations for CF. Furthermore, even after accounting for the pedigree-based heterosis coefficient of an animal, part of the remaining variability in some of the carcass traits could be captured by a genomic heterozygosity measure. This signifies that the inclusion of both a heterosis coefficient based on pedigree information and a genome-based measure of heterozygosity could be beneficial to limiting bias in predicting additive genetic merit. Finally, one SNP located on Bos taurus (BTA) chromosome number 5 demonstrated a non-additive association with CW. Furthermore, 182 SNPs (180 SNPs on BTA 2 and two SNPs on BTA 21) demonstrated a non-additive association with CC, while 231 SNPs located on BTA 2, 5, 11, 13, 14, 18, 19 and 21 demonstrated a non-additive association with CF. Results demonstrate that heterozygosity both at a global level and at the level of individual loci contribute little to the variability in carcass merit.

Predicting carcass cut yields in cattle from digital images using artificial intelligence.

Deep Learning (DL) has proven to be a successful tool for many image classification problems but has yet to be applied to carcass images. The aim of this study was to train DL models to predict carcass cut yields and compare predictions to more standard machine learning (ML) methods. Three approaches were undertaken to predict the grouped carcass cut yields of Grilling cuts and Roasting cuts from a large dataset of 54,598 and 69,246 animals respectively. The approaches taken were (1) animal phenotypic data used as features for a range of ML algorithms, (2) carcass images used to train Convolutional Neural Networks, and (3) carcass dimensions measured directly from the carcass images, combined with the associated phenotypic data and used as feature data for ML algorithms. Results showed that DL models can be trained to predict carcass cuts yields but an approach that uses carcass dimensions in ML algorithms performs slightly better in absolute terms.

Contribution of herd characteristics to best linear unbiased estimates of slaughter traits in beef cattle.

Genetic evaluations separate phenotypes into their contributing additive genetic effects and non-(additive) genetic effects, with the former termed best linear unbiased predictions, and the latter termed best linear unbiased estimates (BLUEs). For the purpose of the present study, genetic evaluations, along with phenotypic data from 4 137 376 animals, were used to generate herd, year of slaughter and sex contemporary group BLUEs for various slaughter-related traits. These slaughter traits included carcass weight (CW), carcass conformation (CC) and carcass fat (CF). For the 4 665 herds that were consistently slaughtering ≥10 animals/year between the years 2014 and 2018, inclusive, all relevant contemporary group BLUEs were collapsed into a single herd-year value; results herein relate to these herds. The within-year herd-year BLUE correlations between CW and CC, between CW and CF, and between CC and CF were 0.51, 0.10 and -0.04, respectively. The repeatability across years of the herd-year BLUEs for CW, CC and CF was 0.66, 0.59 and 0.50, respectively. Furthermore, when the herds were stratified, within year, on the percentile rank of their herd-year BLUEs, herds had the greatest probability of remaining in the same BLUE stratum from one year to the next. In addition, results from the present study determined that various herd characteristics are associated with differences in the herd BLUEs. Results from the present study could be used to advise beef producers on the most promising strategy to improve the carcass merit of their animals.

Herd solutions from genetic evaluations can be used as a tool to rescale the expected expression of genetic potential in cattle.

The objective of this study was to determine whether response to selection for carcass weight (CW), conformation (CC) and fat (CF), and the association between heterosis and carcass performance varied by herd environment in cattle. Following edits, carcass information was available for 4,616,761 cattle, of which the majority were some crossbred combination of the following breeds: Angus, Aubrac, Belgian Blue, Blonde d'Aquitaine, Charolais, Hereford, Holstein-Friesian, Jersey, Limousin, Saler, Shorthorn and Simmental. Herd environment was defined separately for each carcass trait using herd solutions outputted from carcass trait genetic evaluations. A total of 3,859 herds were stratified, for each trait, into one of five strata based on their corresponding percentile herd solution rank, with the response to selection and the effect of heterosis then estimated within each stratum. The response in CW and CC from selection on the respective estimated breeding values (EBV) increased between the lowest (0.71 kg and 0.89 CC score increase per unit increase in the respective EBV) and highest (0.99 kg and 1.25 CC score increase per unit increase in the respective EBV) corresponding herd stratum. The response in CF from selection on CF EBV, however, reduced between the lowest and highest CF herd stratum (respective increases of 0.93 and 0.83 CF scores per unit increase in CF EBV). In addition, the effect of a unit increase in heterosis coefficient on CW, CC and CF also varied by herd stratum. Furthermore, results (i.e. the area under relative operating characteristic curves) from the present study demonstrated that the response to selection and heterosis effects estimated for the different herd stratum can be used, along with EBVs and the herd solutions themselves, to improve the accuracy of phenotypic predictions. Results from the present study could help producers to make more informed breeding decisions that are bespoke to their herd.

Identification of genomic regions that exhibit sexual dimorphism for size and muscularity in cattle.

Sexual dimorphism, the phenomenon whereby males and females of the same species are distinctive in some aspect of appearance or size, has previously been documented in cattle for traits such as growth rate and carcass merit using a quantitative genetics approach. No previous study in cattle has attempted to document sexual dimorphism at a genome level; therefore, the objective of the present study was to determine whether genomic regions associated with size and muscularity in cattle exhibited signs of sexual dimorphism. Analyses were undertaken on 10 linear-type traits that describe the muscular and skeletal characteristics of both males and females of five beef cattle breeds: 1,444 Angus (AA), 6,433 Charolais (CH), 1,129 Hereford, 8,745 Limousin (LM), and 1,698 Simmental. Genome-wide association analyses were undertaken using imputed whole-genome sequence data for each sex separately by breed. For each single-nucleotide polymorphism (SNP) that was segregating in both sexes, the difference between the allele substitution effect sizes for each sex, in each breed separately, was calculated. Suggestively (P ≤ 1 × 10-5) sexually dimorphic SNPs that were segregating in both males and females were detected for all traits in all breeds, although the location of these SNPs differed by both trait and breed. Significantly (P ≤ 1 × 10-8) dimorphic SNPs were detected in just three traits in the AA, seven traits in the CH, and three traits in the LM. The vast majority of all segregating autosomal SNPs (86% in AA to 94% in LM) had the same minor allele in both males and females. Differences (P ≤ 0.05) in allele frequencies between the sexes were observed for between 36% (LM) and 66% (AA) of the total autosomal SNPs that were segregating in both sexes. Dimorphic SNPs were located within a number of genes related to muscularity and/or size including the NAB1, COL5A2, and IWS1 genes on BTA2 that are located close to, and thought to be co-inherited with, the MSTN gene. Overall, sexual dimorphism exists in cattle at the genome level, but it is not consistent by either trait or breed.

An index framework founded on the future profit potential of female beef cattle to aid the identification of candidates for culling.

Meticulous culling decisions, coupled with careful breeding decisions, are fundamental to shifting a population distribution in the favorable direction and improving profit per cow. Nevertheless, there is a paucity of easy-to-use dynamic tools to aid in culling decisions in beef cattle. The motivation for the present study was to develop a monetary-based culling tool, here referred to as the Beef Female's Profit Potential (BFPP), to identify females for culling. The BFPP reflects the expected lifetime profitability of an individual female in a herd for the expected remainder of her lifetime; this profit included that of the beef female herself as well as her progeny. The BFPP index framework was composed of 4 subindexes reflecting the value of an animal: (1) as a nulliparae (this was voided if the cow had already calved), (2) for the remainder of her current parity, (3) summed across each of her expected remaining parities, and (4) when she is retained within the herd and not voluntarily culled. Each subindex was comprised of different components reflecting both genetic and non-genetic effects associated with each female. Transition matrices predicting the expected longevity of each female and their expected month of calving were also utilized in calculating the expected remaining lifetime profitability of each female. The BFPP index was validated on 21,102 beef cows as well as their harvested progeny from 875 herds by stratifying the cows, within herd, into 4 strata based on their BFPP. The mean of the within-herd correlation between the BFPP and the Irish national replacement (i.e., breeding) index was, on average, 0.45 indicating the shortcomings of the breeding index as a culling tool. Cows within the top BFPP stratum had a genetic expectation of accruing almost an additional €36 profit per calving, relative to cows within the worst stratum; when validated on the cow's own calving interval and survival performance as well as their progeny's carcass performance, the actual phenotypic value was estimated to be an additional €32 profit per calving. A proportion of this additional profit was due to the harvested progeny of the high BFPP cows having, on average, heavier, more conformed carcasses with less fat cover relative to their poor BFPP contemporaries. This BFPP framework is a useful and easy-to-use tool to aid in producer decision making on the choice of females to voluntarily cull but also on which replacement heifers to graduate into the mature herd.

Animal-level factors associated with the achievement of desirable specifications in Irish beef carcasses graded using the EUROP classification system.

Beef carcasses in Europe are classified on measures of carcass weight, conformation, and fat cover. These measurements provide the basis for payment to producers, with financial penalties for carcasses that do not conform to desirable characteristics. The objective of the present study was to identify animal-level factors associated with the achievement of a desirable carcass weight, conformation score, fat score, and age at harvest, as stipulated by Irish beef processors in accordance with the EUROP carcass classification system. The stipulated specifications were a EUROP conformation score ≥O=, a carcass weight between 270 and 380 kg, a EUROP fat score between 2+ and 4=, and an age at harvest ≤ 30 mo. In the present study, 59% of cattle failed to achieve at least one of these desired specifications. The logit of the probability of achieving the desired specifications was estimated using multivariable logistic regression and carcass data from 4,717,989 cattle finished and harvested in Ireland between the years 2003 and 2017. In comparison to beef-origin carcasses and after accounting for breed differences, the likelihood of dairy-origin carcasses achieving the desired age, conformation, fat, and weight specifications was 0.97, 0.88, 1.14, and 1.05, respectively. In comparison to heifer carcasses, the odds ratio (OR) of bull and steer carcasses simultaneously achieving all of the desired specifications (i.e. the overall specification) was 0.35 and 0.95, respectively. Additionally, after accounting for breed differences, heifers from the dairy herd were half as likely as heifers from the beef herd to achieve the overall specification, whereas the odds of dairy-origin bulls (OR = 3.46) and steers (OR = 2.41) achieving the overall specification was greater than that of their respective beef-origin counterparts. Finally, cattle with a greater breed proportion of Angus were most likely to achieve the overall specification. Results from the present study could provide a deeper understanding as to why animals fail to achieve desirable carcass specifications and could be implemented into the management decisions made on farm to ensure that the supply of beef carcasses that achieve the desired metrics is maximized.

Breed- and trait-specific associations define the genetic architecture of calving performance traits in cattle.

Reducing the incidence of both the degree of assistance required at calving, as well as the extent of perinatal mortality (PM) has both economic and societal benefits. The existence of heritable genetic variability in both traits signifies the presence of underlying genomic variability. The objective of the present study was to locate regions of the genome, and by extension putative genes and mutations, that are likely to be underpinning the genetic variability in direct calving difficulty (DCD), maternal calving difficulty (MCD), and PM. Imputed whole-genome single-nucleotide polymorphism (SNP) data on up to 8,304 Angus (AA), 17,175 Charolais (CH), 16,794 Limousin (LM), and 18,474 Holstein-Friesian (HF) sires representing 5,866,712 calving events from descendants were used. Several putative quantitative trait loci (QTL) regions associated with calving performance both within and across dairy and beef breeds were identified, although the majority were both breed- and trait-specific. QTL surrounding and encompassing the myostatin (MSTN) gene were associated (P < 5 × 10-8) with DCD and PM in both the CH and LM populations. The well-known Q204X mutation was the fifth strongest association with DCD in the CH population and accounted for 5.09% of the genetic variance in DCD. In contrast, none of the 259 segregating variants in MSTN were associated (P > × 10-6) with DCD in the LM population but a genomic region 617 kb downstream of MSTN was associated (P < 5 × 10-8). The genetic architecture for DCD differed in the HF population relative to the CH and LM, where two QTL encompassing ZNF613 on Bos taurus autosome (BTA)18 and PLAG1 on BTA14 were identified in the former. Pleiotropic SNP associated with all three calving performance traits were also identified in the three beef breeds; 5 SNP were pleiotropic in AA, 116 in LM, and 882 in CH but no SNP was associated with more than one trait within the HF population. The majority of these pleiotropic SNP were on BTA2 surrounding MSTN and were associated with both DCD and PM. Multiple previously reported, but also novel QTL, associated with calving performance were detected in this large study. These also included QTL regions harboring SNP with the same direction of allele substitution effect for both DCD and MCD thus contributing to a more effective simultaneous selection for both traits.

The achievement of a given carcass specification is under moderate genetic control in cattle.

The objective of the present study was to estimate the genetic parameters associated with the achievement of desirable weight, conformation, and fat specifications, represented by a series of binary traits. The desired specifications were those stipulated by Irish beef processors, in accordance with the EUROP carcass grading system, and were represented by a carcass weight between 270 and 380 kg, a fat score between 2+ and 4= (between 6 and 11 on a 15-point scale), and a conformation score of O= or better (≥5 on a 15-point scale). Using data from 58,868 beef carcasses, variance components were estimated using linear mixed models for these binary traits, as well as their underlying continuous measures. Heritability estimates for the continuous traits ranged from 0.63 to 0.73; heritability estimates for the binary traits ranged from 0.05 to 0.19. An additional trait was defined to reflect if all desired carcass specifications were met. All genetic correlations between this trait and the individual contributing binary traits were positive (0.38 to 0.87), while all genetic correlations between this trait and the continuous carcass measures were negative (-0.87 to -0.07). The genetic parameters estimated in the present study signify that potential exists to breed cattle that more consistently achieve desirable carcass metrics at harvest.

Genomic Regions Associated With Skeletal Type Traits in Beef and Dairy Cattle Are Common to Regions Associated With Carcass Traits, Feed Intake and Calving Difficulty.

Linear type traits describing the skeletal characteristics of an animal are moderately to strongly genetically correlated with a range of other performance traits in cattle including feed intake, reproduction traits and carcass merit; thus, type traits could also provide useful insights into the morphological differences among animals underpinning phenotypic differences in these complex traits. The objective of the present study was to identify genomic regions associated with five subjectively scored skeletal linear traits, to determine if these associated regions are common in multiple beef and dairy breeds, and also to determine if these regions overlap with those proposed elsewhere to be associated with correlated performance traits. Analyses were carried out using linear mixed models on imputed whole genome sequence data separately in 1,444 Angus, 1,129 Hereford, 6,433 Charolais, 8,745 Limousin, 1,698 Simmental, and 4,494 Holstein-Friesian cattle, all scored for the linear type traits. There was, on average, 18 months difference in age at assessment of the beef versus the dairy animals. While the majority of the identified quantitative trait loci (QTL), and thus genes, were both trait-specific and breed-specific, a large-effect pleiotropic QTL on BTA6 containing the NCAPG and LCORL genes was associated with all skeletal traits in the Limousin population and with wither height in the Angus. Other than that, little overlap existed in detected QTLs for the skeletal type traits in the other breeds. Only two QTLs overlapped the beef and dairy breeds; both QTLs were located on BTA5 and were associated with height in both the Angus and the Holstein-Friesian, despite the difference in age at assessment. Several detected QTLs in the present study overlapped with QTLs documented elsewhere that are associated with carcass traits, feed intake, and calving difficulty. While most breeding programs select for the macro-traits like carcass weight, carcass conformation, and feed intake, the higher degree of granularity with selection on the individual linear type traits in a multi-trait index underpinning the macro-level goal traits, presents an opportunity to help resolve genetic antagonisms among morphological traits in the pursuit of the animal with optimum performance metrics.

Genomic regions associated with muscularity in beef cattle differ in five contrasting cattle breeds.

BACKGROUND: Linear type traits, which reflect the muscular characteristics of an animal, could provide insight into how, in some cases, morphologically very different animals can yield the same carcass weight. Such variability may contribute to differences in the overall value of the carcass since primal cuts vary greatly in price; such variability may also hinder successful genome-based association studies. Therefore, the objective of our study was to identify genomic regions that are associated with five muscularity linear type traits and to determine if these significant regions are common across five different breeds. Analyses were carried out using linear mixed models on imputed whole-genome sequence data in each of the five breeds, separately. Then, the results of the within-breed analyses were used to conduct an across-breed meta-analysis per trait. RESULTS: We identified many quantitative trait loci (QTL) that are located across the whole genome and associated with each trait in each breed. The only commonality among the breeds and traits was a large-effect pleiotropic QTL on BTA2 that contained the MSTN gene, which was associated with all traits in the Charolais and Limousin breeds. Other plausible candidate genes were identified for muscularity traits including PDE1A, PPP1R1C and multiple collagen and HOXD genes. In addition, associated (gene ontology) GO terms and KEGG pathways tended to differ between breeds and between traits especially in the numerically smaller populations of Angus, Hereford, and Simmental breeds. Most of the SNPs that were associated with any of the traits were intergenic or intronic SNPs located within regulatory regions of the genome. CONCLUSIONS: The commonality between the Charolais and Limousin breeds indicates that the genetic architecture of the muscularity traits may be similar in these breeds due to their similar origins. Conversely, there were vast differences in the QTL associated with muscularity in Angus, Hereford, and Simmental. Knowledge of these differences in genetic architecture between breeds is useful to develop accurate genomic prediction equations that can operate effectively across breeds. Overall, the associated QTL differed according to trait, which suggests that breeding for a morphologically different (e.g. longer and wider versus shorter and smaller) more efficient animal may become possible in the future.

Genomic Regions Associated With Gestation Length Detected Using Whole-Genome Sequence Data Differ Between Dairy and Beef Cattle.

While many association studies exist that have attempted to relate genomic markers to phenotypic performance in cattle, very few have considered gestation length as a phenotype, and of those that did, none used whole genome sequence data from multiple breeds. The objective of the present study was therefore to relate imputed whole genome sequence data to estimated breeding values for gestation length using 22,566 sires (representing 2,262,706 progeny) of multiple breeds [Angus (AA), Charolais (CH), Holstein-Friesian (HF), and Limousin (LM)]. The associations were undertaken within breed using linear mixed models that accounted for genomic relatedness among sires; a separate association analysis was undertaken with all breeds analysed together but with breed included as a fixed effect in the model. Furthermore, the genome was divided into 500 kb segments and whether or not segments harboured a single nucleotide polymorphism (SNP) with a P ≤ 1 × 10-4 common to different combinations of breeds was determined. Putative quantitative trait loci (QTL) regions associated with gestation length were detected in all breeds; significant associations with gestation length were only detected in the HF population and in the across-breed analysis of all 22,566 sires. Twenty-five SNPs were significantly associated (P ≤ 5 × 10-8) with gestation length in the HF population. Of the 25 significant SNPs, 18 were located within three QTLs on Bos taurus autosome number (BTA) 18, six were in two QTL on BTA19, and one was located within a QTL on BTA7. The strongest association was rs381577268, a downstream variant of ZNF613 located within a QTL spanning from 58.06 to 58.19 Mb on BTA18; it accounted for 1.37% of the genetic variance in gestation length. Overall there were 11 HF animals within the edited dataset that were homozygous for the T allele at rs381577268 and these had a 3.3 day longer (P < 0.0001) estimated breeding value (EBV) for gestation length than the heterozygous animals and a 4.7 day longer (P < 0.0001) EBV for gestation length than the homozygous CC animals. The majority of the 500 kb windows harboring a SNP with a P ≤ 1 × 10-4 were unique to a single breed and no window was shared among all four breeds for gestation length, suggesting any QTLs identified are breed-specific associations.

Using the difference in actual and expected calf liveweight relative to its dam liveweight as a statistic for interherd and intraherd benchmarking and genetic evaluations1.

The importance of improving the efficiency of beef production systems using both genetic and management strategies has long been discussed. Despite the contribution of the mature beef herd to the overall cost of production in the sector as a whole, most strategies for improving (feed) efficiency have focused on the growing animal. The objective of the present study was to quantify the phenotypic and genetic variability in several novel measures that relate the weight of a calf to that of its dam and vice versa. Two novel residual traits, representing the deviation in calf weight relative to its expectation from the population based on its dam's weight (DIFFcalf) or the deviation in the weight of the dam relative to its expectation from the population based on its calf's weight (DIFFdam), were calculated while simultaneously accounting for some nuisance factors in a multiple regression model. Four supplementary traits were also calculated, namely, 1) the deviation in calf weight from its expectation expressed relative to the weight of the dam (DIFFcalf_ratio), 2) the deviation in dam weight from its expectation relative to the weight of the dam (DIFFdam_ratio), 3) DIFFcalf-DIFFdam, and 4) the simple ratio of calf weight to its dam's weight (RATIOcalfdam). Genetic and residual variance components for each of the 6 traits were estimated using animal-dam linear mixed models. The phenotypic SD for DIFFcalf was 42 kg and, when expressed relative to the weight of the dam (i.e., DIFFcalf_ratio), was 0.07. The genetic SD for DIFFcalf and DIFFcalf_ratio was 16.66 kg and 0.02, respectively. The direct and maternal heritability estimated for DIFFcalf was 0.28 (SE = 0.04) and 0.11 (SE = 0.02), respectively, and for DIFFcalf_ratio was 0.24 (SE = 0.04) and 0.17 (SE = 0.03), respectively. The genetic SD for DIFFdam was 47.09 kg; the direct heritability was 0.50 (SE = 0.03), and the dam repeatability was 0.75 (SE = 0.01). The genetic SD for RATIOcalfdam was 0.03; the direct and maternal heritability was 0.24 (SE = 0.04) and 0.24 (SE = 0.03), respectively. The suggested traits outlined in the present study provide useful metrics for benchmarking dam-calf efficiency; in addition, the genetic variability detected in these traits suggest genetic progress for more efficient dam-calf pairs is indeed possible.