Early detection of infectious bovine keratoconjunctivitis with artificial intelligence.

Artificial intelligence (AI) was developed to distinguish cattle by their muzzle patterns and identify early cases of disease, including infectious bovine keratoconjunctivitis (IBK). It was tested on 870 cattle in four locations, with 170 developing IBK. The AI identified 169 of the 170 cases prior to their identification by veterinarians, and another 17 cases that remained free of IBK signs (sensitivity = 99.4%, specificity = 97.6%). These results indicate the AI can detect emerging IBK cases by muzzle images very early in the disease process and be used as an intervention tool in the prevention of IBK outbreaks.

Using pooled data for genomic prediction in a bivariate framework with missing data.

Pooling samples to derive group genotypes can enable the economically efficient use of commercial animals within genetic evaluations. To test a multivariate framework for genetic evaluations using pooled data, simulation was used to mimic a beef cattle population including two moderately heritable traits with varying genetic correlations, genotypes and pedigree data. There were 15 generations (n = 32,000; random selection and mating), and the last generation was subjected to genotyping through pooling. Missing records were induced in two ways: (a) sequential culling and (b) random missing records. Gaps in genotyping were also explored whereby genotyping occurred through generation 13 or 14. Pools of 1, 20, 50 and 100 animals were constructed randomly or by minimizing phenotypic variation. The EBV was estimated using a bivariate single-step genomic best linear unbiased prediction model. Pools of 20 animals constructed by minimizing phenotypic variation generally led to accuracies that were not different than using individual progeny data. Gaps in genotyping led to significantly different EBV accuracies (p < .05) for sires and dams born in the generation nearest the pools. Pooling of any size generally led to larger accuracies than no information from generation 15 regardless of the way missing records arose, the percentage of records available or the genetic correlation. Pooling to aid in the use of commercial data in genetic evaluations can be utilized in multivariate cases with varying relationships between the traits and in the presence of systematic and randomly missing phenotypes.

A Reference Genome Assembly of Simmental Cattle, Bos taurus taurus.

Genomics research has relied principally on the establishment and curation of a reference genome for the species. However, it is increasingly recognized that a single reference genome cannot fully describe the extent of genetic variation within many widely distributed species. Pangenome representations are based on high-quality genome assemblies of multiple individuals and intended to represent the broadest possible diversity within a species. A Bovine Pangenome Consortium (BPC) has recently been established to begin assembling genomes from more than 600 recognized breeds of cattle, together with other related species to provide information on ancestral alleles and haplotypes. Previously reported de novo genome assemblies for Angus, Brahman, Hereford, and Highland breeds of cattle are part of the initial BPC effort. The present report describes a complete single haplotype assembly at chromosome-scale for a fullblood Simmental cow from an F1 bison-cattle hybrid fetus by trio binning. Simmental cattle, also known as Fleckvieh due to their red and white spots, originated in central Europe in the 1830s as a triple-purpose breed selected for draught, meat, and dairy production. There are over 50 million Simmental cattle in the world, known today for their fast growth and beef yields. This assembly (ARS_Simm1.0) is similar in length to the other bovine assemblies at 2.86 Gb, with a scaffold N50 of 102 Mb (max scaffold 156.8 Mb) and meets or exceeds the continuity of the best Bos taurus reference assemblies to date.

A Reference Genome Assembly of American Bison, Bison bison bison.

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

Lamb wool shedding is a good predictor of ewe wool shedding.

Interest in reducing labour costs due to shearing has led to development of breed types that shed their wool naturally. Selection at young ages can facilitate response. Reliability of predictions of adult from lamb wool shedding (WS) is thus key in the design of breeding programmes to increase shedding. Our objectives were to estimate heritabilities and genetic relationships between WS measured once in lambs and repeatedly in ewes and to assess the accuracy of lamb WS EBV to predict ewe WS EBV based on a multi-trait threshold or a repeatability model. Data were 4,971 lamb and 3,335 ewe WS records on a Romanov, White Dorper and Katahdin composite flock. For the multivariate model, WS heritability ranged from 0.47 ± 0.03 in lambs to 0.59 ± 0.04 at 1 year of age. For the repeatability model, WS in adult ewes was moderately heritable (0.50 ± 0.03) and repeatable (0.60 ± 0.02). Genetic correlations were 0.72 ± 0.04, 0.65 ± 0.05, 0.50 ± 0.09 and 0.51 ± 0.09 between lamb WS and 1st through 4th record, respectively. Given the moderately high heritability and high correlations between WS performance in lambs and ewes, selecting animals early in life would effectively increase WS in crossbred flocks.

Longitudinal study of humoral immunity to bovine coronavirus, virus shedding, and treatment for bovine respiratory disease in pre-weaned beef calves.

BACKGROUND: Bovine coronavirus (BCV) is associated with respiratory infections in cattle of all ages; however, a temporal study to evaluate the effect of BCV immunity on virus shedding and bovine respiratory disease (BRD) incidence in pre-weaned beef calves has not been reported. Thus, we report here a prospective study in three herds of crossbred beef calves (n = 817) with endemic BCV. Serial blood samples for measurement of serum anti-BCV antibody titers and nasal swabs for detection of BCV and other common viral and bacterial BRD pathogens were collected from all calves or subsets of calves at predetermined times from birth through weaning. The calves were monitored for BRD and those that developed signs of respiratory disease were sampled for diagnostic testing. To discover additional risk factors that could have influenced BRD development, sequence analysis of the BCV strain(s) circulating in each herd, and the prevalence of common opportunistic bacterial pathogens in the upper respiratory tract of sick and apparently healthy cattle were also evaluated. RESULTS: Two hundred forty-eight of the 817 study calves (30.4%) were treated for BRD prior to weaning; 246 of those were from a single herd involved in two outbreaks of BRD leading to mass treatment of all calves in that group. Molecular diagnostic testing found BCV and Histophilus somni in nasal swabs taken at the time of BRD treatment. Between herd analyses revealed anti-BCV serum antibody abundance did not associate with the incidence of BRD or BCV shedding, though these measurements may have been hindered by the long periods between sample collections. Analysis of the BCV spike gene hypervariable region revealed four polymorphisms in 15 isolates from the three herds, making strain variation unlikely to account for differences in treatment rates between herds. Persistent or recurrent shedding episodes of BCV occurred in some animals treated for BRD. CONCLUSION: Co-detection of BCV and H. somni at the time of the disease outbreak suggests that these pathogens contributed to disease pathogenesis. Developing appropriate control measures for respiratory BCV infections may help decrease the incidence of pre-weaning BRD. The role of antibodies in protection must still be further defined.

RNA-Seq Meta-analysis identifies genes in skeletal muscle associated with gain and intake across a multi-season study of crossbred beef steers.

BACKGROUND: Feed intake and body weight gain are economically important inputs and outputs of beef production systems. The purpose of this study was to discover differentially expressed genes that will be robust for feed intake and gain across a large segment of the cattle industry. Transcriptomic studies often suffer from issues with reproducibility and cross-validation. One way to improve reproducibility is by integrating multiple datasets via meta-analysis. RNA sequencing (RNA-Seq) was performed on longissimus dorsi muscle from 80 steers (5 cohorts, each with 16 animals) selected from the outside fringe of a bivariate gain and feed intake distribution to understand the genes and pathways involved in feed efficiency. In each cohort, 16 steers were selected from one of four gain and feed intake phenotypes (n = 4 per phenotype) in a 2 × 2 factorial arrangement with gain and feed intake as main effect variables. Each cohort was analyzed as a single experiment using a generalized linear model and results from the 5 cohort analyses were combined in a meta-analysis to identify differentially expressed genes (DEG) across the cohorts. RESULTS: A total of 51 genes were differentially expressed for the main effect of gain, 109 genes for the intake main effect, and 11 genes for the gain x intake interaction (Pcorrected < 0.05). A jackknife sensitivity analysis showed that, in general, the meta-analysis produced robust DEGs for the two main effects and their interaction. Pathways identified from over-represented genes included mitochondrial energy production and oxidative stress pathways for the main effect of gain due to DEG including GPD1, NDUFA6, UQCRQ, ACTC1, and MGST3. For intake, metabolic pathways including amino acid biosynthesis and degradation were identified, and for the interaction analysis the pathways identified included GADD45, pyridoxal 5'phosphate salvage, and caveolar mediated endocytosis signaling. CONCLUSIONS: Variation among DEG identified by cohort suggests that environment and breed may play large roles in the expression of genes associated with feed efficiency in the muscle of beef cattle. Meta-analyses of transcriptome data from groups of animals over multiple cohorts may be necessary to elucidate the genetics contributing these types of biological phenotypes.

Association of Circulating Transfer RNA fragments with antibody response to Mycoplasma bovis in beef cattle.

BACKGROUND: High throughput sequencing allows identification of small non-coding RNAs. Transfer RNA Fragments are a class of small non-coding RNAs, and have been identified as being involved in inhibition of gene expression. Given their role, it is possible they may be involved in mediating the infection-induced defense response in the host. Therefore, the objective of this study was to identify 5' transfer RNA fragments (tRF5s) associated with a serum antibody response to M. bovis in beef cattle. RESULTS: The tRF5s encoding alanine, glutamic acid, glycine, lysine, proline, selenocysteine, threonine, and valine were associated (P < 0.05) with antibody response against M. bovis. tRF5s encoding alanine, glutamine, glutamic acid, glycine, histidine, lysine, proline, selenocysteine, threonine, and valine were associated (P < 0.05) with season, which could be attributed to calf growth. There were interactions (P < 0.05) between antibody response to M. bovis and season for tRF5 encoding selenocysteine (anticodon UGA), proline (anticodon CGG), and glutamine (anticodon TTG). Selenocysteine is a rarely used amino acid that is incorporated into proteins by the opal stop codon (UGA), and its function is not well understood. CONCLUSIONS: Differential expression of tRF5s was identified between ELISA-positive and negative animals. Production of tRF5s may be associated with a host defense mechanism triggered by bacterial infection, or it may provide some advantage to a pathogen during infection of a host. Further studies are needed to establish if tRF5s could be used as a diagnostic marker of chronic exposure.

Incorporation of genetic technologies associated with applied reproductive technologies to enhance world food production.

Animal breeding and reproductive physiology have been closely related throughout the history of animal production science. Artificial insemination provides the best method of increasing the influence of sires with superior genetics to improve production traits. Multiple ovulation embryo transfer (MOET) provides some ability to increase the genetic influence of the maternal line as well. The addition of genetic technologies to this paradigm allows for improved methods of selecting sires and dams carrying the best genes for production and yield of edible products and resistance to diseases and parasites. However, decreasing the number of influential parents within a population also increases the risk of propagating a recessive gene that could negatively impact the species (Reprod Domest Anim 44:792-796, 2009; BMC Genomics 11:337, 2010). Furthermore, antagonistic genotypic relationships between production traits and fertility (Anim Prod Sci 49:399-412, 2009; Anim Genet 43:442-446, 2012) suggest that care must be taken to ensure that increasing the frequency of genes with a positive influence on production does not negatively impact the fertility of the replacement females entering the herd.

Beef cattle body temperature during climatic stress: a genome-wide association study.

Cattle are reared in diverse environments and collecting phenotypic body temperature (BT) measurements to characterize BT variation across diverse environments is difficult and expensive. To better understand the genetic basis of BT regulation, a genome-wide association study was conducted utilizing crossbred steers and heifers totaling 239 animals of unknown pedigree and breed fraction. During predicted extreme heat and cold stress events, hourly tympanic and vaginal BT devices were placed in steers and heifers, respectively. Individuals were genotyped with the BovineSNP50K_v2 assay and data analyzed using Bayesian models for area under the curve (AUC), a measure of BT over time, using hourly BT observations summed across 5-days (AUC summer 5-day (AUCS5D) and AUC winter 5-day (AUCW5D)). Posterior heritability estimates were moderate to high and were estimated to be 0.68 and 0.21 for AUCS5D and AUCW5D, respectively. Moderately positive correlations between direct genomic values for AUCS5D and AUCW5D (0.40) were found, although a small percentage of the top 5% 1-Mb windows were in common. Different sets of genes were associated with BT during winter and summer, thus simultaneous selection for animals tolerant to both heat and cold appears possible.

Genetic parameters for carcass traits of progeny of beef bulls mated to dairy cows.

Beef × dairy crossbred cattle (n = 615) were used to evaluate the effect of preharvest indicator traits and genotypes on the accuracy of estimated breeding values (EBVs) of seedstock candidates for selection. Genotypes for 100,000 single nucleotide polymorphisms were provided by the American Simmental Association of purebred and crossbred seedstock animals (n = 2,632). Five hundred and ninety-five of the 615 beef × dairy cattle had carcass camera and ultrasound data. Phenotypes were not used for any of the seedstock animals even though some may have had performance and ultrasound data. We estimated the genomic relationship matrix among 3,247 animals including both phenotyped and unphenotyped animals. We computed genetic parameters among 37 traits using 666 bivariate restricted maximum likelihood analyses. The accuracy of EBV depends on heritability. For the sake of brevity, we report accuracy for marbling as a proxy for other traits with similar heritability. We focus on accuracy for marbling because marbling is the primary determinant of carcass value. We computed EBV for all 3,247 animals for marbling based on camera data postharvest using best linear unbiased prediction. We report evidence of overlap in causative genes among postharvest carcass traits; marbling, ribeye area, yield grade, fat thickness, and hot carcass weight (HCW) based on genetic correlations. Genetic correlations range from -0.73 to 0.89. Several live animal traits (frame size, body weight and ultrasound fat thickness and ribeye area) were genetically correlated with postharvest traits; including HCW, ribeye area, yield grade, fat thickness, and marbling. Genetic correlations between pre- and postharvest traits ranged from -0.53 to 0.95. Accuracy for marbling ranged from 0.64 to 0.80 for animals with marbling recorded, and from 0.09 to 0.60 for animals without marbling recorded. The accuracy of animals without phenotypes was related to the genomic relationship between animals with phenotype and those without. Live animal traits were useful for predicting economically important carcass traits based on genetic correlations. The accuracy of EBV for seedstock animals that were not phenotyped was low, but this is consistent with theory, and accuracy is expected to increase with the addition of genotypes and carcass data from beef × dairy animals.

Low-cost genotyping platforms and sexed-semen have enabled the production of high breeding value dairy replacement heifers from a fraction of the herd representing the most elite cows. The remainder of the cow herd can be bred to beef bulls using male-sexed-semen. Camera carcass data postharvest and ultrasound carcass estimates preharvest (live animals) on beef × dairy animals combined with genotypes and ultrasound on seedstock animals may provide an efficient scheme for selecting beef bulls to mate to dairy cows in the future to maximize carcass value of the progeny. Genotypes are needed to link carcass data from previously harvested to seedstock bull selection candidates because pedigree is typically not available for beef × dairy cattle. We report that live animal ultrasound carcass estimates are predictive of postharvest economically important carcass traits. The accuracy of genetic evaluation of selection candidates without recorded carcass traits was low but is expected to increase with more genotypes and phenotypes on beef × dairy cattle. Genotypes, ultrasound estimates, and camera carcass data on thousands of beef × dairy cattle could enable increased accuracy of selection with periodic infusion of new phenotypes from future generations.

Genes involved in the cholecystokinin receptor signaling map were differentially expressed in the jejunum of steers with variation in residual feed intake.

The jejunum is a critical site for nutrient digestion and absorption, and variation in its ability to take up nutrients within the jejunum is likely to affect feed efficiency. The purpose of this study was to determine differences in gene expression in the jejunum of beef steers divergent for residual feed intake (RFI) in one cohort of steers (Year 1), and to validate those genes in animals from a second study (Year 2). Steers from Year 1 (n = 16) were selected for high and low RFI. Jejunum mucosal tissue was obtained for RNA-seq. Thirty-two genes were differentially expressed (PFDR≤0.15), and five were over-represented in pathways including inflammatory mediator, cholecystokinin receptor (CCKR) signaling, and p38 MAPK pathways. Several differentially expressed genes (ALOX12, ALPI, FABP6, FABP7, FLT1, GSTA2, MEF2B, PDK4, SPP1, and TTF2) have been previously associated with RFI in other studies. Real-time qPCR was used to validate nine differentially expressed genes in the Year 1 steers used for RNA-seq, and in the Year 2 validation cohort. Six genes were validated as differentially expressed (P < 0.1) using RT-qPCR in the Year 1 population. In the Year 2 population, five genes displayed the same direction of expression as the Year 1 population and 3 were differentially expressed (P < 0.1). The CCKR pathway is involved in digestion, appetite control, and regulation of body weight making it a compelling candidate for feed efficiency in cattle, and the validation of these genes in a second population of cattle is suggestive of a role in feed efficiency.

Influence of therapeutic ceftiofur treatments of feedlot cattle on fecal and hide prevalences of commensal Escherichia coli resistant to expanded-spectrum cephalosporins, and molecular characterization of resistant isolates.

In the United States, the blaCMY-2 gene contained within incompatibility type A/C (IncA/C) plasmids is frequently identified in extended-spectrum-cephalosporin-resistant (ESC(r)) Escherichia coli strains from both human and cattle sources. Concerns have been raised that therapeutic use of ceftiofur in cattle may increase the prevalence of ESC(r) E. coli. We report that herd ESC(r) E. coli fecal and hide prevalences throughout the residency of cattle at a feedlot, including during the period of greatest ceftiofur use at the feedlot, were either not significantly different (P ≥ 0.05) or significantly less (P < 0.05) than the respective prevalences at arrival. Longitudinal sampling of cattle treated with ceftiofur demonstrated that once the transient increase of ESC(r) E. coli shedding that follows ceftiofur injection abated, ceftiofur-injected cattle were no more likely than untreated members of the same herd to shed ESC(r) E. coli. Pulsed-field gel electrophoresis (PFGE) genotyping, antibiotic resistance phenotyping, screening for presence of the blaCMY-2 gene, and plasmid replicon typing were performed on 312 ESC(r) E. coli isolates obtained during six sampling periods spanning the 10-month residence of cattle at the feedlot. The identification of only 26 unique PFGE genotypes, 12 of which were isolated during multiple sampling periods, suggests that clonal expansion of feedlot-adapted blaCMY-2 E. coli strains contributed more to the persistence of blaCMY-2 than horizontal transfer of IncA/C plasmids between E. coli strains at this feedlot. We conclude that therapeutic use of ceftiofur at this cattle feedlot did not significantly increase the herd prevalence of ESC(r) E. coli.

Comparison of molecular breeding values based on within- and across-breed training in beef cattle.

BACKGROUND: Although the efficacy of genomic predictors based on within-breed training looks promising, it is necessary to develop and evaluate across-breed predictors for the technology to be fully applied in the beef industry. The efficacies of genomic predictors trained in one breed and utilized to predict genetic merit in differing breeds based on simulation studies have been reported, as have the efficacies of predictors trained using data from multiple breeds to predict the genetic merit of purebreds. However, comparable studies using beef cattle field data have not been reported. METHODS: Molecular breeding values for weaning and yearling weight were derived and evaluated using a database containing BovineSNP50 genotypes for 7294 animals from 13 breeds in the training set and 2277 animals from seven breeds (Angus, Red Angus, Hereford, Charolais, Gelbvieh, Limousin, and Simmental) in the evaluation set. Six single-breed and four across-breed genomic predictors were trained using pooled data from purebred animals. Molecular breeding values were evaluated using field data, including genotypes for 2227 animals and phenotypic records of animals born in 2008 or later. Accuracies of molecular breeding values were estimated based on the genetic correlation between the molecular breeding value and trait phenotype. RESULTS: With one exception, the estimated genetic correlations of within-breed molecular breeding values with trait phenotype were greater than 0.28 when evaluated in the breed used for training. Most estimated genetic correlations for the across-breed trained molecular breeding values were moderate (> 0.30). When molecular breeding values were evaluated in breeds that were not in the training set, estimated genetic correlations clustered around zero. CONCLUSIONS: Even for closely related breeds, within- or across-breed trained molecular breeding values have limited prediction accuracy for breeds that were not in the training set. For breeds in the training set, across- and within-breed trained molecular breeding values had similar accuracies. The benefit of adding data from other breeds to a within-breed training population is the ability to produce molecular breeding values that are more robust across breeds and these can be utilized until enough training data has been accumulated to allow for a within-breed training set.

Use of overlapping DNA pools to discern genetic differences despite pooling error.

Genotyping pools of commercial cattle and individual seedstock animals may reveal hidden relationships between sectors enabling use of commercial data for genetic evaluation. However, commercial data capture may be compromised by inexact pool formation. We aimed to estimate the concordance between distances or genomic covariance among pooling allele frequencies (PAFs) of DNA pools comprised of 100 animals with 0% or 50% overlap of animals in common between pools. Cattle lung samples were collected from a commercial beef processing plant on a single day. Six pools of 100 animals each were constructed so that overlap between pools was 0% or 50%. Two pools of all 200 animals were constructed to estimate PAFs for all 200 animals. Frozen lung tissue (0.01 g) from each animal was weighed into a tube containing a pool; there were two pools of 200 animals each and six pools of 100 animals each. Every contribution of an individual animal was an independent measurement to insure independence of pooling errors. Lung samples were kept on dried ice during the pooling process to keep them from thawing. The eight pools were then assayed for approximately 100,000 single nucleotide polymorphisms (SNP). PAF for each SNP and pool was based on the relative intensity of the two dyes used to detect the alleles rather than genotype calls which are not tractable from pooling data. Euclidean distances and genomic relationships among the PAFs for the eight pools were estimated and distances were tested for concordance with pool overlap using permutation-based analysis of distance. Distances among pools were concordant with the planned overlap of animals shared between pools (P = 0.0024); pool overlap accounted for 70% of the variation and pooling error accounted for 30%. Pools containing 100 animals with no overlap were the most distant from one another and pools with 50% overlap were the least distant. This work shows that we can discern differences in distance between pairs of overlapping DNA pools sharing 0% and 50% of the animals. Genomic correlations among nonoverlapping pools indicated that nonoverlapping pool pairs did not share many related animals because genomic correlations were near zero for these pairs. On the other hand, one pair of nonoverlapping pools likely contained related animals between pools because the correlation was 0.21. Pools sharing 50% overlap ranged in genomic relationship between 0.21 and 0.39 (N = 12).

Genetic evaluation of seedstock cattle could benefit from commercial data. There are hidden relationships between commercial and seedstock sectors because many commercial producers buy bulls from the seedstock sector. Relationships are hidden because pedigree is not tracked in commercial populations. Single nucleotide polymorphism genotypes could reveal these hidden relationships; however, genotyping can be cost prohibitive. Cost of commercial data capture could be decreased by pooling DNA which is a method to genotype groups of animals to use their data in genetic evaluation; however, error from inexact pool formation can complicate interpretation. Results from pools of overlapping random unrelated animals mimic the results from pools sharing relatives with the same degree of shared genomes. For example, a pool of progeny and a pool of the dams of the pooled progeny would produce the same result as two pools sharing 50% overlap of random unrelated animals. We can estimate the relatedness between unknown pools even in the presence of pooling error if an unknown pool comparison is similar to an overlapping pool comparison. Knowing the relationship between seedstock cattle and pools of commercial cattle may allow commercial data to enhance genetic evaluation of seedstock animals.

Variance component estimates for growth traits in beef cattle using selected variants from imputed low-pass sequence data.

A beef cattle population (n = 2,343) was used to assess the impact of variants identified from the imputed low-pass sequence (LPS) on the estimation of variance components and genetic parameters of birth weight (BWT) and post-weaning gain (PWG). Variants were selected based on functional impact and were partitioned into four groups (low, modifier, moderate, high) based on predicted functional impact and re-partitioned based on the consequence of mutation, such as missense and untranslated region variants, into six groups (G1-G6). Each subset was used to construct a genomic relationship matrix (GRM) for univariate animal models. Multiple analyses were conducted to compare the proportion of additive genetic variation explained by the different subsets individually and collectively, and these estimates were benchmarked against all LPS variants in a single GRM and array (e.g., GeneSeek Genomic Profiler 100K) genotypes. When all variants were included in a single GRM, heritability estimates for BWT and PWG were 0.43 ±â€…0.05 and 0.38 ±â€…0.05, respectively. Heritability estimates for BWT ranged from 0.10 to 0.42 dependent on which variant subsets were included. Similarly, estimates for PWG ranged from 0.05 to 0.38. Results showed that variants in the subsets modifier and G1 (untranslated region) yielded the highest heritability estimates and were similar to the inclusion of all variants, while estimates from GRM containing only variants in the categories High, G4 (non-coding transcript exon), and G6 (start and stop loss/gain) were the lowest. All variants combined provided similar heritability estimates to chip genotypes and provided minimal to no additional information when combined with chip data. This suggests that the chip single nucleotide polymorphisms and the variants from LPS predicted to be less consequential are in relatively high linkage disequilibrium with the underlying causal variants as a whole and sufficiently spread throughout the genome to capture larger proportions of additive genetic variation.

Animals from a crossbred beef cattle population were sequenced at low depth (i.e., 0.5×) and different subsets of selected imputed variants were investigated relative to their ability to explain variation in birth weight (BWT) and post-weaning gain (PWG). Variants were classified by both their predicted functional impact and by the consequence of the mutation and partitioned into subsets within these two criteria. When ~ 1 million variants were included in the same genomic relationship matrix, heritability estimates were similar to a 100k chip array. Heritability estimates for BWT ranged from 0.10 to 0.42 dependent on which variant subsets were included. Similarly, estimates for PWG ranged from 0.05 to 0.38. Differences in minor allele frequency were observed among subsets and these differences likely contributed to differences in heritability estimates. Results suggest that linkage disequilibrium between the variants categorized as being less consequential and underlying causal variants is high as indicated by the high percentage of variation explained.

Heritability of beef cow metabolizable energy for maintenance.

Most of the metabolizable energy that a cow uses during a production year is for maintenance; however, less is known about the heritability of maintenance compared to other traits that can be measured directly. Feed intake is a heritable trait in the mature cow and most of the feed consumed is used for maintenance. We hypothesized that maintenance energy was a heritable trait. Individual feed intake was measured for 84 or 85 d on 5 yr old pregnant cows (N = 887) from a pedigreed population of cattle that represent prominent breeds in the United States. Phenotypic mean (± SD) values were 654 ± 68 kg for cow body weight, 0.21 ± 0.24 kg/d for average daily gain, and 175 ± 17 d for midpoint fetal age. Dry matter intake averaged (± SD) 10.84 ± 1.41 kg/d. Metabolizable energy for maintenance was estimated by subtracting the metabolizable energy used for conceptus growth and tissue accretion from metabolizable energy intake. Metabolizable energy for maintenance averaged (± SD) 139 ± 18 ME kcal/d/BW kg0.75 and had a heritability of 0.31 ± 0.11. Cows have a moderate heritability for maintenance suggesting an opportunity for selection.

Feed is one of the greatest costs of beef production. Most of the feed used annually by a cow is to maintain her body. A study was conducted measuring individual feed intake of mature pregnant cows. We have determined that the amount of energy that a cow uses to maintain her body is heritable suggesting that cows can be selected for differences in the energy required to maintain their bodies.

Genetic and phenotypic associations of mitochondrial DNA copy number, SNP, and haplogroups with growth and carcass traits in beef cattle.

Mitochondrial DNA copy number (mtDNA CN) is heritable and easily obtained from low-pass sequencing (LPS). This study investigated the genetic correlation of mtDNA CN with growth and carcass traits in a multi-breed and crossbred beef cattle population. Blood, leucocyte, and semen samples were obtained from 2,371 animals and subjected to LPS that resulted in nuclear DNA (nuDNA) and mtDNA sequence reads. Mitochondrial DNA CN was estimated as the ratio of mtDNA to nuDNA coverages. Variant calling was performed from mtDNA, and 11 single nucleotide polymorphisms (SNP) were identified in the population. Samples were classified in taurine haplogroups. Haplogroup and mtDNA type were further classified based on the 11 segregating SNP. Growth and carcass traits were available for between 7,249 and 60,989 individuals. Associations of mtDNA CN, mtDNA haplogroups, mtDNA types, and mtDNA SNP with growth and carcass traits were estimated with univariate animal models, and genetic correlations were estimated with a bivariate animal model based on pedigree. Mitochondrial DNA CN tended (P-value ≤0.08) to be associated with birth weight and weaning weight. There was no association (P-value >0.10) between mtDNA SNP, haplogroups, or types with growth and carcass traits. Genetic correlation estimates of mtDNA CN were -0.30 ± 0.16 with birth weight, -0.31 ± 0.16 with weaning weight, -0.15 ± 0.14 with post-weaning gain, -0.11 ± 0.19 with average daily dry-matter intake, -0.04 ± 0.22 with average daily gain, -0.29 ± 0.13 with mature cow weight, -0.11 ± 0.13 with slaughter weight, -0.14 ± 0.13 with carcass weight, -0.07 ± 0.14 with carcass backfat, 0.14 ± 0.14 with carcass marbling, and -0.06 ± 0.14 with ribeye area. In conclusion, mtDNA CN was negatively correlated with most traits investigated, and the genetic correlation was stronger with growth traits than with carcass traits.

This study investigated mitochondrial DNA copy number (mtDNA CN) as a potential genetic indicator of growth and carcass traits in a composite beef cattle population. Mitochondrial DNA CN was previously shown to be under genetic control. The current study estimated the genetic relationship of mtDNA CN with growth and carcass traits. Blood, leucocyte, and semen samples were obtained from 2,371 animals and subjected to whole-genome sequencing at a low depth that resulted in nuclear DNA and mtDNA sequence reads. Mitochondrial DNA CN was estimated as the ratio of mtDNA to nuclear DNA coverages. Growth and carcass traits were available for between 7,249 and 60,989 individuals. Genetic parameters were estimated from an animal model based on pedigree. Genetic correlation estimates of mtDNA CN with growth and carcass traits were low to moderate and mostly negative. These indicate that selection for lower mtDNA would be associated with an increase in birth weight, weaning weight, post-weaning gain, average daily dry-matter intake, mature cow weight, slaughter weight, and carcass weight. Therefore, the by-product of whole-genome sequencing at a low depth could be used as an indicator trait for growth and carcass traits in genetic evaluations, but the genetic relationships are not likely strong enough to prioritize mtDNA ahead of routinely used indicator traits.

Heritability and variance component estimation for feed and water intake behaviors of feedlot cattle.

Feed and water intake are two important aspects of cattle production that greatly impact the profitability, efficiency, and sustainability of producers. Feed and, to a lesser degree, water intake have been studied previously; however, there is little research on their associated animal behaviors and there is a lack of standardized phenotypes for these behaviors. Feed and water intakes obtained with an Insentec system (Hokofarm Group, The Netherlands) from 830 crossbred steers were used to compute five intake behaviors for both feed and water: daily sessions (DS), intake rate (IR), session size (SS), time per session (TS), and session interval (SI). Variance components and heritabilities were estimated for each trait. Heritabilities for feed intake behaviors were 0.50 ±â€…0.12, 0.63 ±â€…0.12, 0.40 ±â€…0.13, 0.35 ±â€…0.12, and 0.60 ±â€…0.12 for DS, IR, SS, TS, and SI, respectively. Heritabilities for water intake behaviors were 0.56 ±â€…0.11, 0.88 ±â€…0.07, 0.70 ±â€…0.11, 0.54 ±â€…0.12, and 0.80 ±â€…0.10 for NS, IR, SS, TS, and SI, respectively. Daily dry matter intake (DDMI) and daily water intake (DWI) had heritabilities of 0.57 ±â€…0.11 and 0.44 ±â€…0.11. Phenotypic correlations varied between pairs of traits (-0.83 to 0.82). Genetic correlations between DDMI and feed intake behaviors were moderate to high, while genetic correlations between DWI and water intake behaviors were low to moderate. Several significant single nucleotide polymorphisms (SNP) were identified for the feed and water intake behaviors. Genes and previously reported quantitative trait loci near significant SNPs were evaluated. The results indicated that feed and water intake behaviors are influenced by genetic factors and are heritable, providing one additional route to evaluate or manipulate feed and water intake.

Feed and water intake are important aspects of cattle production to understand because they impact producer profitability and sustainability. While feed intake and, to a lesser degree, water intake have previously been studied, the associated feeding and drinking behaviors are relatively unknown and lack standardized phenotypes. Using individual animal feed and water intake records, five behaviors were evaluated for feed and water intake from crossbred feedlot steers. The behaviors evaluated were daily sessions (no./d), session size (kg), time per session (s), intake rate (g/s), and session interval (min). The impact of season (winter vs. summer) and bunk management (ad libitum vs. slick) on feeding and drinking behaviors was evaluated. Heritability and variance components were estimated for all feeding and drinking behaviors. Pairwise phenotypic correlations between behaviors were discussed. The relationship between intake and feeding or drinking behaviors was evaluated with genetic correlations. A genome-wide association study identified several significant single nucleotide polymorphisms for feeding and drinking behaviors. The results indicate that feeding and drinking behaviors are heritable and may be one additional route to evaluate feed and water intake.