Relationship of molecular breeding value for beef tenderness with heifer traits through weaning of their first calf

Polymorphisms in μ-calpain (CAPN1) that beneficially associate with beef tenderness are reported to antagonistically associate with calving day in beef heifers and post-partum interval to estrus in beef cows. We, therefore, hypothesized that a molecular breeding value for slice shear force, calculated based on CAPN1 and calpastatin (CAST) genotypes, would demonstrate an antagonistic relationship between genomically predicted slice shear force and ordinal calving date in replacement beef heifers. A secondary objective of this study was to evaluate the association of a polymorphism in diacylglycerol O-acyltransferase (DGAT1) with reproductive traits in beef heifers. One hundred eighty-seven MARC III heifers (¼ Angus, ¼ Hereford, ¼ Red Poll, and ¼ Pinzgauer) that had been selectively bred to increase the frequency of these polymorphisms were submitted for monthly ultrasound exams beginning at 333 d of age and continuing until the start of breeding to determine pubertal status. At the last exam before breeding, all antral follicles were counted, and the length and height of each ovary was measured to determine if genomic selection for slice shear force associated with ovarian follicle number. Calving date, calf gender, and calf birth weight were recorded at parturition. Regression analysis of the molecular breeding value for slice shear force of the heifers on ordinal calving date indicated no association between genomic prediction of tenderness and calving date (P = 0.16); however, there was a tendency for age at puberty to be delayed in heifers as genetic merit for tenderness improved (P = 0.09). The results of the present study indicate that within experimental precision, selecting for tenderness using genomic predictions had minimal or no antagonistic association with reproductive performance in heifers. Further analysis of reproductive performance as cows is needed within this population but applying these genetic markers to select for tenderness in steers does not antagonize reproductive traits influencing conception or first calf birth date and birth weight in replacement beef heifers.

Production performance of cows raised with different postweaning growth patterns

The period of heifer development is a relatively small fraction of a cow's life; however, her pattern of growth may have permanent effects on her productivity as a cow. We hypothesized that altering the growth pattern during the peri-pubertal period would increase life-time productivity across genetic types of Bos taurus cows. The objective was to determine the stayability, calf production, and weight of calf weaned across six calf crops. Heifers (n = 685) were placed on one of two developmental programs at 256 ± 1 d of age. Control heifers received a diet that provided 228 kcal ME·(body weight [BW], kg) ^-0.75 daily, and stair-step heifers were allocated 157 kcal ME·(BW, kg)^-0.75 daily for 84 or 85 d, and then the daily allocation was increased to 277 kcal ME·(BW, kg)^-0.75. Stair-step heifers (0.33 ± 0.02 kg/d) had a lower average daily gain (ADG) than control heifers (0.78 ± 0.02 kg/d; P < 0.001) during Period 1, and stair-step heifers (0.93 ± 0.03 kg/d) had a greater ADG than controls (0.70 ± 0.03 kg/d; P < 0.001) during Period 2. There were no treatment (P = 0.28) or breed type differences (P = 0.42) for the proportion of cows weaning a calf; however, the proportion of cows weaning a calf decreased with cow age (P < 0.001). Calves from stair-step dams had heavier weaning weights (193 ± 1 kg) compared to control calves (191 ± 1 kg; P = 0.007). There was not a treatment (P = 0.25) or breed type differences in cumulative BW weaned (P = 0.59). A diverse genetic population of cattle within B. taurus was tested and responses in calf production did not differ between stair-step growth pattern and a more constant nonobese growth pattern.

A large interactive visual database of copy number variants discovered in taurine cattle

BACKGROUND

Copy number variants (CNVs) contribute to genetic diversity and phenotypic variation. We aimed to discover CNVs in taurine cattle using a large collection of whole-genome sequences and to provide an interactive database of the identified CNV regions (CNVRs) that includes visualizations of sequence read alignments, CNV boundaries, and genome annotations.

RESULTS

CNVs were identified in each of 4 whole-genome sequencing datasets, which together represent >500 bulls from 17 breeds, using a popular multi-sample read-depth-based algorithm, cn.MOPS. Quality control and CNVR construction, performed dataset-wise to avoid batch effects, resulted in 26,223 CNVRs covering 107.75 unique Mb (4.05%) of the bovine genome. Hierarchical clustering of samples by CNVR genotypes indicated clear separation by breeds. An interactive HTML database was created that allows data filtering options, provides graphical and tabular data summaries including Hardy-Weinberg equilibrium tests on genotype proportions, and displays genes and quantitative trait loci at each CNVR. Notably, the database provides sequence read alignments at each CNVR genotype and the boundaries of constituent CNVs in individual samples. Besides numerous novel discoveries, we corroborated the genotypes reported for a CNVR at the KIT locus known to be associated with the piebald coat colour phenotype in Hereford and some Simmental cattle.

CONCLUSIONS

We present a large comprehensive collection of taurine cattle CNVs in a novel interactive visual database that displays CNV boundaries, read depths, and genome features for individual CNVRs, thus providing users with a powerful means to explore and scrutinize CNVRs of interest more thoroughly.

Likelihood of obtaining tender meat from confined calf

Abstract The objective was to evaluate the main factors that influence the shear strength of meat from confined steers and the probability of obtaining soft meat. For this purpose, we evaluated the literature on carcass and/or meat of beef steers in Brazil published between January 1999 and April 2019 and extracted the content from three sections for analysis: materials and methods, results, and discussion. Pearson’s correlation was used to analyze the data, and the stepwise statistic was used to determine the proportion of the synchronized effect of variables on shear force. For determining the probability of tenderness, meat with a shear force lower than 4.6 kgf/cm3 was classified as soft; meat with a higher sheer force was classified as hard. Following the classification, logistic regression analysis and odds ratio test were performed. The factors of study location, the proportion of zebu background in the genome, finishing weight, the percentage of concentrate in the diet, and finishing period and meat marbling explained 62.45% of the variability in the shear strength of beef. The following strategies were found to increase the chances of effectively obtaining soft meat from confined steers: starting the termination phase early even in animals with lower weights, prolonging the confinement time, increasing concentrate percentage in the diet, and a higher marbling degree. It is possible to estimate a large proportion of shear force variability using the production variables (ante-mortem), and the process can be adjusted accordingly to considerable increase the possibility of obtaining soft meat.

Enhanced estimates of carcass and meat quality effects for polymorphisms in myostatin and µ-calpain genes

The objective of this study was to enhance estimates of additive, dominance, and epistatic effects of marker polymorphisms on beef carcass and quality traits. Myostatin (MSTN) F94L SNP and the µ-calpain (CAPN1) 316 and 4751 SNP haplotype have previously been associated with fat and muscle traits in beef cattle. Multiyear selection in a composite population segregating these polymorphisms increased minor allele (F94L L) and chosen haplotype (CAPN1 CC and GT) frequencies to intermediate levels resulting in more precise estimates of additive and nonadditive genetic effects. During the 3 yr after selection, 176 steers were evaluated for growth, carcass, meat quality, tenderness (n = 103), and meat color traits. The statistical model included year, age of dam, age of the steer, and genotype in a random animal model. The 9 genotypes (3 CAPN1 diplotypes × 3 F94L genotypes) affected marbling score, ribeye area, adjusted fat thickness, vision yield grade (all P < 0.001), slice shear force (P = 0.03), and CIE L* reflectance (P = 0.01). Linear contrasts of the 9 genotypes estimated additive, recessive, and epistatic genetic effects. Significant additive effects of the F94L L allele decreased marbling score, adjusted fat thickness, vision yield grade, and slice shear force; and increased ribeye area and CIE L* reflectance. The homozygous F94L FF and LL genotypes differed by 1.3 to 1.9 phenotypic SD for most carcass traits and by 0.8 to 0.9 SD for slice shear force and CIE L* reflectance but carcass weight differed by only 3 kg (0.1 SD). The L allele was partially recessive to F for ribeye area (P = 0.02) and the heterozygous FL means tended to be closer to the FF genotype than the LL genotype for other carcass traits but differences from additive were not significant. The CAPN1 additive × F94L additive effect on slice shear force was the only significant epistatic estimate. The F94L L allele is prevalent in Limousin but nearly absent in other U.S. purebreds. This allele had about half of the effects on birth weight, muscle, and fat traits reported for severe MSTN mutations in Belgian Blue and Piedmontese breeds. The interaction between MSTN and CAPN1 genotypes may reflect the strong additive effects of MSTN F94L L allele on fat and muscle traits interfering with the phenotypic effect of CAPN1 genotype on meat tenderness.

Comparing allele specific expression and local expression quantitative trait loci and the influence of gene expression on complex trait variation in cattle

BACKGROUND

The mutations changing the expression level of a gene, or expression quantitative trait loci (eQTL), can be identified by testing the association between genetic variants and gene expression in multiple individuals (eQTL mapping), or by comparing the expression of the alleles in a heterozygous individual (allele specific expression or ASE analysis). The aims of the study were to find and compare ASE and local eQTL in 4 bovine RNA-sequencing (RNA-Seq) datasets, validate them in an independent ASE study and investigate if they are associated with complex trait variation.

RESULTS

We present a novel method for distinguishing between ASE driven by polymorphisms in cis and parent of origin effects. We found that single nucleotide polymorphisms (SNPs) driving ASE are also often local eQTL and therefore presumably cis eQTL. These SNPs often, but not always, affect gene expression in multiple tissues and, when they do, the allele increasing expression is usually the same. However, there were systematic differences between ASE and local eQTL and between tissues and breeds. We also found that SNPs significantly associated with gene expression (p < 0.001) were likely to influence some complex traits (p < 0.001), which means that some mutations influence variation in complex traits by changing the expression level of genes.

CONCLUSION

We conclude that ASE detects phenomenon that overlap with local eQTL, but there are also systematic differences between the SNPs discovered by the two methods. Some mutations influencing complex traits are actually eQTL and can be discovered using RNA-Seq including eQTL in the genes CAST, CAPN1, LCORL and LEPROTL1.

Effects of polymorphisms in CAPN1 and CAST genes on meat tenderness of Chinese Simmental cattle

Considerable evidence has demonstrated that the μ -calpain (CAPN1) gene and its inhibitor calpastatin (CAST) gene are major factors affecting meat quality. Marker-assisted selection (MAS) has been widely used to improve beef quality traits. Therefore, the objective of the present study was to investigate the single nucleotide polymorphisms (SNPs) of bovine CAPN1 and CAST genes using 367 animals representing the four main Chinese cattle breeds and to explore the effects of these SNPs on meat quality traits. Two SNPs within CAPN1 and one SNP in CAST were successfully identified in cattle. Genetic diversity analyses suggested that most SNPs in the four breeds exhibited a moderate genetic diversity. Moreover, associations between individual markers and meat quality traits were analyzed in Chinese Simmental cattle. The CAPN1 4558 A > G locus was found to be significantly associated with shear force value (SFV) and marbling score (BMS), and CAPN1 4684 C > T exerted a significant effect on SFV, while the CAST genotype was not significantly associated with any of the measured traits. SFV, commonly used to measure meat tenderness, represents an important quality trait as it contributes to the flavor of cooked meat. This work confirms the effect of CAPN1 on beef tenderness and lays an important foundation for future cattle breeding.

Effect of Brahman genetics on myofibrillar protein degradation, collagen crosslinking, and tenderness of the longissimus lumborum

The objective of this study was to examine the effect of percent Brahman genetics on Warner-Bratzler shear force (WBSF), desmin and troponin-T (TnT) degradation, hydroxylysyl pyridinoline (HP) crosslink content, and perimysial collagen melting temperature. Steers ( = 131) produced in 2012 and 2013 were harvested at 1.27 cm of visual s.c. back fat thickness. Steers were divided into 4 genetic categories consisting of steers that contained 6/32 or less Brahman genetics, 12/32 Brahman genetics, 14/32 to 18/32 Brahman genetics, and 23/32 to 32/32 Brahman genetics. Twenty-four hours after harvest, a 7.62-cm piece of the longissimus lumborum beginning at the 13th rib was collected and aged for 14 d. Following aging, three 2.54-cm steaks were cut for WBSF, trained sensory panel, and laboratory analyses. Laboratory analyses steaks were used to determine protein degradation, HP crosslink analysis, and perimysial collagen melting temperature. Data were analyzed using a polynomial regression for unequally spaced treatments. As the percent Brahman genetics increased, WBSF increased (linear, = 0.01). As percent Brahman genetics increased, tenderness score decreased (less tender) and connective tissue score increased (more connective tissue; linear, = 0.01). As the percentage of Brahman genetics increased, the amount of degraded desmin (38 kDa) and TnT (34 and 30 kDa) decreased (linear, < 0.03) whereas the amount of immunoreactive 36 kDa TnT increased (linear, = 0.04). Percent Brahman genetics had no effect ( = 0.14) on HP crosslink content but did tend to increase ( = 0.07) perimysial collagen melting temperature as the percent Brahman increased. The percentage of Brahman genetic influence was positively correlated to WBSF ( = 0.25), 36 kDa immunoreactive TnT ( = 0.26), and perimysial collagen melting temperature ( = 0.25, = 0.01). Sensory panel tenderness ( = -0.44), juiciness ( = -0.26), and connective tissue scores ( = -0.63); 38 kDa degraded desmin ( = -0.34), 34 ( = -0.36) and 30 kDa degraded TnT ( = -0.29); and HP collagen crosslinks ( = -0.20) were negatively correlated to percent Brahman genetic influence ( < 0.03). Increasing Brahman genetic influence in steers negatively affects tenderness, partially through a reduction in degradation of desmin and TnT. Although HP collagen crosslinks are unaffected by Brahman genetics, a tendency for increased perimysium melting temperature indicates that other collagen-stabilizing crosslinks may be affected.

A Review of Selected Genes with Known Effects on Performance and Health of Cattle

There are genetic conditions that influence production in dairy and beef cattle. The objective of this review was to describe relevant genetic conditions that have been associated with productivity and health in cattle. Genes or genomic regions that have been identified as a candidate for the condition will be included, and the genetic basis of the condition will be defined. Genes and genetic conditions included in this review are bovine leukocyte adhesion deficiency, deficiency of the uridine monophosphate synthase, bovine chronic interstitial nephritis, horn development, myostatin, complex vertebral malformation, leptin, osteopetrosis, apoptosis peptide activating factor 1, chondrodysplastic dwarfism, caseins, calpastatin, umbilical hernia, lactoglobulin, citrullinemia, cholesterol deficiency, prions, thyroglobulin, diacylglycerol acyltransferase, syndactyly, maple syrup urine disease, slick hair, Factor XI deficiency, and μ-Calpain. This review is not meant to be comprehensive, and relevant information is provided to ascertain genetic markers associated with the conditions.

Detection of genomic signatures of recent selection in commercial broiler chickens

Background

Identification of the genomic signatures of recent selection may help uncover causal polymorphisms controlling traits relevant to recent decades of selective breeding in livestock. In this study, we aimed at detecting signatures of recent selection in commercial broiler chickens using genotype information from single nucleotide polymorphisms (SNPs). A total of 565 chickens from five commercial purebred lines, including three broiler sire (male) lines and two broiler dam (female) lines, were genotyped using the 60K SNP Illumina iSelect chicken array. To detect genomic signatures of recent selection, we applied two methods based on population comparison, cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio (XP-CLR), and further analyzed the results to find genomic regions under recent selection in multiple purebred lines.

Results

A total of 321 candidate selection regions spanning approximately 1.45 % of the chicken genome in each line were detected by consensus of results of both XP-EHH and XP-CLR methods. To minimize false discovery due to genetic drift, only 42 of the candidate selection regions that were shared by 2 or more purebred lines were considered as high-confidence selection regions in the study. Of these 42 regions, 20 were 50 kb or less while 4 regions were larger than 0.5 Mb. In total, 91 genes could be found in the 42 regions, among which 19 regions contained only 1 or 2 genes, and 9 regions were located at gene deserts.

Conclusions

Our results provide a genome-wide scan of recent selection signatures in five purebred lines of commercial broiler chickens. We found several candidate genes for recent selection in multiple lines, such as SOX6 (Sex Determining Region Y-Box 6) and cTR (Thyroid hormone receptor beta). These genes may have been under recent selection due to their essential roles in growth, development and reproduction in chickens. Furthermore, our results suggest that in some candidate regions, the same or opposite alleles have been under recent selection in multiple lines. Most of the candidate genes in the selection regions are novel, and as such they should be of great interest for future research into the genetic architecture of traits relevant to modern broiler breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0430-1) contains supplementary material, which is available to authorized users.

Implementing meta-analysis from genome-wide association studies for pork quality traits

Pork quality plays an important role in the meat processing industry. Thus, different methodologies have been implemented to elucidate the genetic architecture of traits affecting meat quality. One of the most common and widely used approaches is to perform genome-wide association (GWA) studies. However, a limitation of many GWA in animal breeding is the limited power due to small sample sizes in animal populations. One alternative is to implement a meta-analysis of GWA (MA-GWA) combining results from independent association studies. The objective of this study was to identify significant genomic regions associated with meat quality traits by performing MA-GWA for 8 different traits in 3 independent pig populations. Results from MA-GWA were used to search for genes possibly associated with the set of evaluated traits. Data from 3 pig data sets (U.S. Meat Animal Research Center, commercial, and Michigan State University Pig Resource Population) were used. A MA was implemented by combining -scores derived for each SNP in every population and then weighting them using the inverse of estimated variance of SNP effects. A search for annotated genes retrieved genes previously reported as candidates for shear force (calpain-1 catalytic subunit [] and calpastatin []), as well as for ultimate pH, purge loss, and cook loss (protein kinase, AMP-activated, γ 3 noncatalytic subunit []). In addition, novel candidate genes were identified for intramuscular fat and cook loss (acyl-CoA synthetase family member 3 mitochondrial []) and for the objective measure of muscle redness, CIE a* (glycogen synthase 1, muscle [] and ferritin, light polypeptide []). Thus, implementation of MA-GWA allowed integration of results for economically relevant traits and identified novel genes to be tested as candidates for meat quality traits in pig populations.

Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

Estimates of epistatic and pleiotropic effects of () and () genetic markers on beef heifer performance traits enhanced by selection

Genetic marker effects and type of inheritance are estimated with poor precision when minor marker allele frequencies are low. A stable composite population (MARC II) was subjected to marker assisted selection for 2 yr to equalize and genetic marker frequencies to evaluate the epistatic and pleiotropic effects of these markers on BW, reproduction, and first calf performance traits in replacement beef females ( = 171) managed under 2 postweaning development protocols. Traits evaluated on the heifers were birth BW, weaning BW, 11-mo BW, 12-mo BW, 13-mo BW, first breeding season pregnancy evaluation BW, first calving season BW, 11-mo puberty, 12-mo puberty, 13-mo puberty, first breeding season pregnancy, and first calf weaning rate. Additionally, heifer's first calf performance traits of ordinal calving date, first calf birth BW, and first calf weaning BW (with and without age adjustment) were analyzed. Selection to increase minor allele frequencies and balanced sampling across genotype classes enhanced the ability to detect all genetic effects except dominance × dominance epistasis. The × genotype effect was significant ( < 0.05) for 11-mo BW and 12-mo BW and tended to be significant ( = 0.08) for 13-mo BW. Consistently, for all 3 traits, the most significant effect among epistatic × genotype effects was the additive effect, with the G allele decreasing BW. There were no associations between × genotype and fertility related traits ( ≥ 0.46) in this study. Additionally, there were no × genotype associations with first progeny performance traits ( ≥ 0.14). The large effect of the additive × additive interaction on first calf weaning BW was imprecisely estimated, which may warrant further investigation.

Co-dependence of genotype and dietary protein intake to affect expression on amino acid/peptide transporters in porcine skeletal muscle

A total of 96 barrows (48 pure-bred Bama mini-pigs representing fatty genotype, and 48 Landrace pigs representing lean genotype) were randomly assigned to either a low- or adequate-protein treatment diet. The experimental period commenced at 5 weeks of age and extended to the finishing period. After euthanasia, blood and skeletal muscle samples were collected from pigs at the nursery, growing, and finishing phases. Our results indicate that the concentrations of free AAs in the plasma and muscle decreased as the age of the pigs increased. In addition, a strain × growth phase interaction (P < 0.05) was observed for the free AA pool in the plasma and muscle. The low-protein diet upregulated (P < 0.05) the mRNA levels for T1R1/T1R3 involved in glutamate binding, but downregulated (P < 0.05) the mRNA levels for PAT1, PAT2, and ASCT2, which transport neutral AAs into muscles. Bama mini-pigs had higher (P < 0.05) mRNA levels for LAT1, SNAT2, and EAAC1, but a lower (P < 0.05) mRNA level for PepT1, compared with Landrace pigs. Collectively, our findings indicate that adequate provision of dietary protein plays an important role in regulating profiles of free AA pools and expression of key AA/peptide transporters/transceptors in a genotype- and tissue-specific manner.