Growth and reproduction genomic-polygenic and polygenic parameters and prediction trends as Brahman fraction increases in an Angus-Brahman multibreed population

Impact of Brahman genetics on skin histology characteristics with implications for heat tolerance in cattle

Cattle lose heat predominantly through cutaneous evaporation at the skin-hair coat interface when experiencing heat stress. Sweating ability, sweat gland properties, and hair coat properties are a few of the many variables determining the efficacy of evaporative cooling. Sweating is a significant heat dissipation mechanism responsible for 85% of body heat loss when temperatures rise above 86⁰F. The purpose of this study was to characterize skin morphological parameters in Angus, Brahman, and their crossbred cattle. Skin samples were taken during the summer of 2017 and 2018 from a total of 319 heifers from six breed groups ranging from 100% Angus to 100% Brahman. Epidermis thickness decreased as the percentage of Brahman genetics increased where the 100% Angus group had a significantly thicker epidermis compared to the 100% Brahman animals. A more extended epidermis layer was identified in Brahman animals due to more pronounced undulations in this skin layer. Breed groups with 75% and 100% Brahman genes were similar and had the largest sweat gland area, indicative of superior resilience to heat stress, compared to breed groups with 50% or lower Brahman genetics. There was a significant linear breed group effect on sweat gland area indicating an increase of 862.0 µm² for every 25% increase in Brahman genetics. Sweat gland length increased as the Brahman percentage increased, while the sweat gland depth showed an opposite trend, decreasing from 100% Angus to 100% Brahman. The number of sebaceous glands was highest in 100% Brahman animals which had about 1.77 more sebaceous glands (p < 0.05) per 4.6 mm²area. Conversely, the sebaceous gland area was greatest in the 100% Angus group. This study identified significant differences in skin properties related to heat exchange ability between Brahman and Angus cattle. Equally important, these differences are also accompanied by significant levels of variation within each breed, which is indicative that selection for these skin traits would improve the heat exchange ability in beef cattle. Further, selecting beef cattle for these skin traits would lead to increased resilience to heat stress without disrupting production traits.

Cattle breed type and anabolic implants impact calpastatin expression and abundance of mRNA associated with protein turnover in the longissimus thoracis of feedlot steers

Two methods that the beef cattle industry can use to improve efficiency, sustainability, and economic viability are growth promotants and crossbreeding cattle of different breed types. In the United States, over 90% of cattle receive an anabolic implant at some point during production resulting in an overall increase in skeletal muscle growth. Recent research suggests that the two main cattle breed types, Bos indicus and Bos taurus, respond differently to anabolic implants. The objective of this study was to characterize changes that occur in skeletal muscle following implanting in Bos indicus influenced steers or Bos taurus steers. Twenty steers were stratified by initial weight in a 2 × 2 factorial design examining two different breeds: Angus (AN; n = 10) or Santa Gertrudis influenced (SG; n = 10), and two implant strategies: no implant (CON; n = 10) or a combined implant containing 120 mg TBA and 24 mg E2 (IMP; n = 10; Revalor-S, Merck Animal Health). Skeletal muscle biopsies were taken from the longissimus thoracis (LT) 2 and 10 d post-implantation. The mRNA abundance of 24 genes associated with skeletal muscle growth were examined, as well as the protein expression of µ-calpain and calpastatin. Succinate dehydrogenase mRNA abundance was impacted (P = 0.05) by a breed × treatment interaction 2 d post-implanting, with SG-CON having a greater increased abundance than all other steers. A tendency for a breed × treatment interaction was observed for calpain-6 mRNA (P = 0.07), with SG-CON having greater abundance than AN-CON and SG-IMP. Additionally, calpastatin protein expression was altered (P = 0.01) by a breed × treatment interaction, with SG-CON and SG-IMP steers having increased expression (P = 0.01) compared with AN-CON steers. At 2 d post-implanting, a breed × treatment interaction was observed with SG-CON steers having greater (P = 0.05) mRNA abundance of mitogen-activated protein kinase compared with AN-CON steers. Furthermore, breed affected (P = 0.05) calpastatin abundance with AN steers having increased (P = 0.05) abundance 2 d post-implanting compared with SG steers. Meanwhile, implants tended to affect (P = 0.09) muscle RING finger protein-1 mRNA abundance, with CON steers having increased (P = 0.09) abundance compared with that of IMP steers. These findings suggest that cattle breed type and anabolic implants impact calpastatin expression and mRNA abundance associated with protein turnover in the LT of feedlot steers 2 and 10 d post-implantation.

Interactions between cattle breed type and anabolic implant strategy impact circulating serum metabolites, feedlot performance, feeding behavior, carcass characteristics, and economic return in beef steers

Introducing Bos indicus (BI) genetics into a beef herd has the potential to increase environmental sustainability. When introducing BI genetics, there are concerns regarding negative impacts on temperament, growth, and carcass characteristics. Implants are routinely used in the United States, with majority of cattle on feed receiving an anabolic implant to improve growth and efficiency, however research regarding the interaction between cattle breed type and anabolic implants is limited. This research compared the use of implants in BI influenced animals versus Bos taurus in a feedlot setting. Twenty steers were stratified by initial weight in a 2 × 2 factorial design examining two different breeds: Angus (AN; n = 10) or Santa Gertrudis influenced (SG; n = 10), and two implant strategies: no implant (CON; n = 10) or a combined implant containing 120 mg TBA and 24 mg E2 (IMP; n = 10; Revalor-S, Merck Animal Health). We hypothesized that anabolic implants would improve growth and feedlot performance of BI influenced animals. Steers were randomly placed into covered pens equipped with GrowSafe bunks and fed the same ration for 129 d. Steers were weighed every 28 d. Dry matter intake, feeding behavior, and carcass data of the steers was collected. Blood was collected and harvested as serum on d 0, 2, 10, 28 and every 28 d after that, and analyzed for serum urea nitrogen (SUN), haptoglobin, and 25HydroxyVitamin D. Angus steers tended to gain more (P = 0.06) weight than SG, while IMP tended to gain more (P = 0.10) weight than CON with no breed × treatment interaction observed (P > 0.10). A breed × treatment interaction was observed when analyzing SUN (P = 0.05) and haptoglobin (P = 0.02) concentrations. Serum 25HydroxyVitmain D concentrations tended to be increased (P = 0.09) in SG-IMP steers compared to SG-CON steers. Angus steers tended (P = 0.10) to have greater amounts of marbling compared to SG steers, while SG steers had improved (P = 0.04) yield grade. Economic return was decreased by $46 a head when introducing SG genetics, while implanting steers improved economic return by $46 a head. This research provides evidence suggesting that BI influenced animals may respond differently to anabolic implants when compared to BT animals. Economic analyses demonstrate that anabolic implants improve economic return to beef producers, while introducing SG genetics decreases economic return in animals raised in more temperate climates.

Genomic Approaches Reveal Pleiotropic Effects in Crossbred Beef Cattle

Carcass and meat quality are two important attributes for the beef industry because they drive profitability and consumer demand. These traits are of even greater importance in crossbred cattle used in subtropical and tropical regions for their superior adaptability because they tend to underperform compared to their purebred counterparts. Many of these traits are challenging and expensive to measure and unavailable until late in life or after the animal is harvested, hence unrealistic to improve through traditional phenotypic selection, but perfect candidates for genomic selection. Before genomic selection can be implemented in crossbred populations, it is important to explore if pleiotropic effects exist between carcass and meat quality traits. Therefore, the objective of this study was to identify genomic regions with pleiotropic effects on carcass and meat quality traits in a multibreed Angus-Brahman population that included purebred and crossbred animals. Data included phenotypes for 10 carcass and meat quality traits from 2,384 steers, of which 1,038 were genotyped with the GGP Bovine F-250. Single-trait genome-wide association studies were first used to investigate the relevance of direct additive genetic effects on each carcass, sensory and visual meat quality traits. A second analysis for each trait included all other phenotypes as covariates to correct for direct causal effects from identified genomic regions with pure direct effects on the trait under analysis. Five genomic windows on chromosomes BTA5, BTA7, BTA18, and BTA29 explained more than 1% of additive genetic variance of two or more traits. Moreover, three suggestive pleiotropic regions were identified on BTA10 and BTA19. The 317 genes uncovered in pleiotropic regions included anchoring and cytoskeletal proteins, key players in cell growth, muscle development, lipid metabolism and fat deposition, and important factors in muscle proteolysis. A functional analysis of these genes revealed GO terms directly related to carcass quality, meat quality, and tenderness in beef cattle, including calcium-related processes, cell signaling, and modulation of cell-cell adhesion. These results contribute with novel information about the complex genetic architecture and pleiotropic effects of carcass and meat quality traits in crossbred beef cattle.

RNA-seq analysis identifies cytoskeletal structural genes and pathways for meat quality in beef

RNA sequencing (RNA-seq) has allowed for transcriptional profiling of biological systems through the identification of differentially expressed (DE) genes and pathways. A total of 80 steers with extreme phenotypes were selected from the University of Florida multibreed Angus-Brahman herd. The average slaughter age was 12.91±8.69 months. Tenderness, juiciness and connective tissue assessed by sensory panel, along with marbling, Warner-Bratzler Shear Force (WBSF) and cooking loss, were measured in longissimus dorsi muscle. Total RNA was extracted from muscle and one RNA-seq library per sample was constructed, multiplexed, and sequenced based on protocols by Illumina HiSeq-3000 platform to generate 2×101 bp paired-end reads. The overall read mapping rate using the Btau_4.6.1 reference genome was 63%. A total of 8,799 genes were analyzed using two different methodologies, an expression association and a DE analysis. A gene and exon expression association analysis was carried out using a meat quality index on all 80 samples as a continuous response variable. The expression of 208 genes and 3,280 exons from 1,565 genes was associated with the meat quality index (p-value ≤ 0.05). A gene and isoform DE evaluation was performed analyzing two groups with extreme WBSF, tenderness and marbling. A total of 676 (adjusted p-value≤0.05), 70 (adjusted p-value≤0.1) and 198 (adjusted p-value≤0.1) genes were DE for WBSF, tenderness and marbling, respectively. A total of 106 isoforms from 98 genes for WBSF, 13 isoforms from 13 genes for tenderness and 43 isoforms from 42 genes for marbling (FDR≤0.1) were DE. Cytoskeletal and transmembrane anchoring genes and pathways were identified in the expression association, DE and the gene enrichment analyses; these proteins can have a direct effect on meat quality. Cytoskeletal proteins and transmembrane anchoring molecules can influence meat quality by allowing cytoskeletal interaction with myocyte and organelle membranes, contributing to cytoskeletal structure and architecture maintenance postmortem.

Identification of eQTLs and sQTLs associated with meat quality in beef

BACKGROUND

Transcription has a substantial genetic control and genetic dissection of gene expression could help us understand the genetic architecture of complex phenotypes such as meat quality in cattle. The objectives of the present research were: 1) to perform eQTL and sQTL mapping analyses for meat quality traits in longissimus dorsi muscle; 2) to uncover genes whose expression is influenced by local or distant genetic variation; 3) to identify expression and splicing hot spots; and 4) to uncover genomic regions affecting the expression of multiple genes.

RESULTS

Eighty steers were selected for phenotyping, genotyping and RNA-seq evaluation. A panel of traits related to meat quality was recorded in longissimus dorsi muscle. Information on 112,042 SNPs and expression data on 8588 autosomal genes and 87,770 exons from 8467 genes were included in an expression and splicing quantitative trait loci (QTL) mapping (eQTL and sQTL, respectively). A gene, exon and isoform differential expression analysis previously carried out in this population identified 1352 genes, referred to as DEG, as explaining part of the variability associated with meat quality traits. The eQTL and sQTL mapping was performed using a linear regression model in the R package Matrix eQTL. Genotype and year of birth were included as fixed effects, and population structure was accounted for by including as a covariate the first PC from a PCA analysis on genotypic data. The identified QTLs were classified as cis or trans using 1 Mb as the maximum distance between the associated SNP and the gene being analyzed. A total of 8377 eQTLs were identified, including 75.6% trans, 10.4% cis, 12.5% DEG trans and 1.5% DEG cis; while 11,929 sQTLs were uncovered: 66.1% trans, 16.9% DEG trans, 14% cis and 3% DEG cis. Twenty-seven expression master regulators and 13 splicing master regulators were identified and were classified as membrane-associated or cytoskeletal proteins, transcription factors or DNA methylases. These genes could control the expression of other genes through cell signaling or by a direct transcriptional activation/repression mechanism.

CONCLUSION

In the present analysis, we show that eQTL and sQTL mapping makes possible positional identification of gene and isoform expression regulators.

Genetic parameters for hair characteristics and core body temperature in a multibreed Brahman-Angus herd1

Thermal stress in hot humid conditions limits cattle production. The objectives for this study were to estimate genetic parameters for hair characteristics and core body temperature under low and high temperature humidity index (THI) conditions. Hair samples were collected and measured for length and diameter. Core body temperature was measured as vaginal temperature every 15 min over a 5-d period using an iButton temperature measuring device implanted in a blank CIDR in 336 heifers from the University of Florida multibreed herd (ranging from 100% Angus to 100% Brahman). Restricted maximum likelihood procedures were used to estimate heritabilities from multiple bivariate animal models using the WOMBAT program. Estimates of heritability for hair diameter, undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions were 0.50, 0.67, 0.42, 0.32, and 0.26, respectively. The genetic parameters estimated in this study indicate a large, exploitable genetic variance which can be selected upon to improve tolerance in cattle. Breed effects for differing compositions of Brahman and Angus were also estimated. As Brahman breed composition increased by 25% undercoat length, topcoat length, body temperature under low THI conditions, and body temperature under high THI conditions decreased by 1.32 mm, 2.94 mm, 0.11 °C, and 0.14 °C, respectively. Under both low and high THI conditions, cattle with 25% Brahman breed composition or greater maintained a significantly lower body temperature than the 100% Angus breed group. The incorporation of Brahman germplasm is recommended for herds that often experience heat stress conditions in order to increase resilience to heat stress.

Genome wide association and gene enrichment analysis reveal membrane anchoring and structural proteins associated with meat quality in beef

BACKGROUND

Meat quality related phenotypes are difficult and expensive to measure and predict but are ideal candidates for genomic selection if genetic markers that account for a worthwhile proportion of the phenotypic variation can be identified. The objectives of this study were: 1) to perform genome wide association analyses for Warner-Bratzler Shear Force (WBSF), marbling, cooking loss, tenderness, juiciness, connective tissue and flavor; 2) to determine enriched pathways present in each genome wide association analysis; and 3) to identify potential candidate genes with multiple quantitative trait loci (QTL) associated with meat quality.

RESULTS

The WBSF, marbling and cooking loss traits were measured in longissimus dorsi muscle from 672 steers. Out of these, 495 animals were used to measure tenderness, juiciness, connective tissue and flavor by a sensory panel. All animals were genotyped for 221,077 markers and included in a genome wide association analysis. A total number of 68 genomic regions covering 52 genes were identified using the whole genome association approach; 48% of these genes encode transmembrane proteins or membrane associated molecules. Two enrichment analysis were performed: a tissue restricted gene enrichment applying a correlation analysis between raw associated single nucleotide polymorphisms (SNPs) by trait, and a functional classification analysis performed using the DAVID Bioinformatic Resources 6.8 server. The tissue restricted gene enrichment approach identified eleven pathways including "Endoplasmic reticulum membrane" that influenced multiple traits simultaneously. The DAVID functional classification analysis uncovered eleven clusters related to transmembrane or structural proteins. A gene network was constructed where the number of raw associated uncorrelated SNPs for each gene across all traits was used as a weight. A multiple SNP association analysis was performed for the top five most connected genes in the gene-trait network. The gene network identified the EVC2, ANXA10 and PKHD1 genes as potentially harboring multiple QTLs. Polymorphisms identified in structural proteins can modulate two different processes with direct effect on meat quality: in vivo myocyte cytoskeletal organization and postmortem proteolysis.

CONCLUSION

The main result from the present analysis is the uncovering of several candidate genes associated with meat quality that have structural function in the skeletal muscle.

Fatty acid profile, mineral content, and palatability of beef from a multibreed Angus-Brahman population

Consumers demand meat that is both healthy and palatable, 2 attributes of meat that are affected by lipid content. Many cattle in the southern United States are Bos indicus influenced, as this improves the ability to survive and thrive in these subtropical regions. However, these animals tend to have leaner carcasses and less marbled meat products. Thus, the objective of this study was to examine the effect of percent Brahman genetics on carcass characteristics, palatability, fatty acids profile, and minerals content in LM of steers from a multibreed population ranging from 100% Angus to 100% Brahman. Breed effect was significant for birth weight (P = 0.0003), weaning weight (P < 0.0001), HCW (P < 0.0001), dressing percentage (P = 0.0008), ribeye area (P = 0.002), quality grade (P < 0.0001), and marbling score (P < 0.0001), and all these traits except dressing percentage decreased as the percentage of Brahman increased. Among palatability traits, breed group had a significant effect only on tenderness (TEND) and connective tissue (CT) scores (P < 0.0001). Least squares means decreased from Angus (5.75 ± 0.13 TEND score and 6.29 ± 0.14 CT score, respectively) to Brahman (4.84 ± 0.10 TEND score and 5.49 ± 0.11 CT score, respectively) as indicated by a significant linear effect. Breed group significantly affected the percentage of several individual fatty acids, saturated fatty acids (SFA), and polyunsaturated fatty acids (PUFA), but not monounsaturated fatty acids (MUFA). The 100% Angus group had the highest percentage of SFA at 49.92%, which was significantly higher (P < 0.05) than the SFA percentage in the 50%, 75%, and 100% Brahman breed groups. Brangus animals also had an increased SFA percentage compared with the 100% Brahman animals (P < 0.05). No significant effect was identified for the concentration of PUFA across the 6 breed groups (P = 0.14). Least squares means decreased from 100% Angus to 100% Brahman for concentration of total fat, SFA, and MUFA (g/mg meat). The concentration of magnesium (P < 0.0001), phosphorus (P = 0.06), and potassium (P = 0.06) increased as the percentage of Brahman increased. Our study shows that breed has a significant effect on the fatty acid profile of beef. Cattle with high Brahman percentage, which are characterized by lower marbled meat, will present a more favorable healthfulness profile with reduced content of SFA and MUFA but the same content of PUFA as purebred Angus animals.

Structural Equation Modeling and Whole-Genome Scans Uncover Chromosome Regions and Enriched Pathways for Carcass and Meat Quality in Beef

Structural equation models involving latent variables are useful tools for formulating hypothesized models defined by theoretical variables and causal links between these variables. The objectives of this study were: (1) to identify latent variables underlying carcass and meat quality traits and (2) to perform whole-genome scans for these latent variables in order to identify genomic regions and individual genes with both direct and indirect effects. A total of 726 steers from an Angus-Brahman multibreed population with records for 22 phenotypes were used. A total of 480 animals were genotyped with the GGP Bovine F-250. The single-step genomic best linear unbiased prediction method was used to estimate the amount of genetic variance explained for each latent variable by chromosome regions of 20 adjacent SNP-windows across the genome. Three types of genetic effects were considered: (1) direct effects on a single latent phenotype; (2) direct effects on two latent phenotypes simultaneously; and (3) indirect effects. The final structural model included carcass quality as an independent latent variable and meat quality as a dependent latent variable. Carcass quality was defined by quality grade, fat over the ribeye and marbling, while the meat quality was described by juiciness, tenderness and connective tissue, all of them measured through a taste panel. From 571 associated genomic regions (643 genes), each one explaining at least 0.05% of the additive variance, 159 regions (179 genes) were associated with carcass quality, 106 regions (114 genes) were associated with both carcass and meat quality, 242 regions (266 genes) were associated with meat quality, and 64 regions (84 genes) were associated with carcass quality, having an indirect effect on meat quality. Three biological mechanisms emerged from these findings: postmortem proteolysis of structural proteins and cellular compartmentalization, cellular proliferation and differentiation of adipocytes, and fat deposition.

Association of μ-Calpain and Calpastatin Polymorphisms with Meat Tenderness in a Brahman-Angus Population

Autogenous proteolytic enzymes of the calpain family are implicated in myofibrillar protein degradation. As a result, the μ-calpain gene and its specific inhibitor, calpastatin, have been repeatedly investigated for their association with meat quality traits in cattle; however, no functional mutation has been identified for these two genes. The objectives of this study were: (1) to assess breed composition effect on tenderness; (2) to perform a linkage disequilibrium (LD) analysis in μ-calpain and calpastatin genes as well as an association analyses with tenderness; and (3) to analyze putative functional SNPs inside the significant LD block for an effect on tenderness. Tenderness measurements and genotypes for 16 SNPs in μ-calpain gene and 28 SNPs in calpastatin gene from 673 steers were analyzed. A bioinformatic analysis identified “putative functional SNPs” inside the associated LD block – polymorphisms able to produce a physical and/or chemical change in the DNA, mRNA, or translated protein in silico. Breed composition had a significant (P < 0.0001) effect on tenderness where animals with more than 80% Angus composition had the most tender meat. One 11-kb LD-block and three LD-blocks of 37, 17, and 14 kb in length were identified in the μ-calpain and calpastatin genes, respectively. Out of these, the LD-block 3 in calpastatin, tagged by SNPs located at 7-98566391 and 7-98581038, had a significant effect on tenderness with the TG-CG diplotype being approximately 1 kg more tender than the toughest diplotype, TG-CG. A total of 768 SNPs in the LD-block 3 of calpastatin were included in the bioinformatic analysis, and 28 markers were selected as putative functional SNPs inside the LD-block 3 of calpastatin; however, none of them were polymorphic in this population. Out of 15 initial polymorphisms segregating inside the LD-block 3 of calpastatin in this population, markers ARSUSMARC116, Cast5, rs730723459, and rs210861835 were found to be significantly associated with tenderness.