Genome-wide association study for response to vaccination in Angus calves1.

BACKGROUND: Bovine respiratory disease complex (BRDC) is one of the most important sources of loss within the beef cattle industry in the USA. Steps have been taken to reduce the incidence of BRDC through vaccination. Despite the effectiveness of vaccines, large proportions of cattle still experience morbidity and mortality. Identification of genomic regions that are associated with variation in response to vaccination would allow for the selection of individuals genetically predisposed to respond to vaccination based on specific markers, while heritability and accuracy estimates would help facilitate genomic selection. This in turn may lead to selection for beef cattle herds that may have lower incidence rate of BRDC after vaccination. This study utilizes an Angus herd of more than 2000 head of cattle to identify these regions of association. RESULTS: Genome wide association studies were performed for viral neutralization antibody level and response to vaccination traits against four different viruses associated with BRDC: bovine viral diarrhea virus 1 and 2 (BVDV1 and BVDV2), bovine respiratory syncytial virus (BRSV), and bovine herpesvirus (BHV1). A total of six 1-Mb windows were associated with greater than 1% of the genetic variance for the analyzed vaccination response traits. Heritabilities ranged from 0.08 to 0.21 and prediction accuracy ranged from 0.01 to 0.33 across 7 different vaccination traits. CONCLUSIONS: Although six 1-Mb windows were identified as associated with 1% or greater genetic variance for viral neutralization antibody level and response to vaccination traits, few genes around these windows could readily be considered candidates. This indicates the need for further functional genomic annotation, as these regions appear to be gene deserts. Traits ranged from lowly to moderately heritable, which indicated the potential for selection of individuals that are genetically pre-disposed to respond to vaccination. The relatively low amount of genetic variance accounted for by any 1-Mb window indicated that viral neutralization antibody level and response to vaccination traits are polygenic in nature. Selection for these traits is possible, but likely to be slow due to the low heritabilities and absence of markers with high genetic variation associated with them.

Evaluation of responses to vaccination of Angus cattle for four viruses that contribute to bovine respiratory disease complex.

Although vaccination is an effective measure in reducing the risk of bovine respiratory disease complex (BRDC) in cattle, BRDC losses remain significant. Increasing the efficacy of vaccination depends on elucidating the protective immune response to different antigens included in vaccines, determining the best timing for vaccination, and understanding the impact of the age of the calf on vaccination. This study measured the serum antibodies present in calves following vaccination against 4 viruses commonly associated with BRDC: bovine viral diarrhea virus type 1 and 2 (BVDV1 and BVDV2), bovine respiratory syncytial virus (BRSV), and bovine herpesvirus 1 (BHV1). Serum antibody titers were measured in more than 1,600 calves at 3-wk intervals starting at the time of the first vaccination. This first vaccination occurred at weaning for approximately half of the individuals and 3 wk before weaning for the other half. Dam age (years), time of weaning (initial vaccination or booster vaccination), and age of calf within year-season (days within year-season) classification all were found to have a significant effect on measured traits such as the initial titer and overall response. An increased initial titer was negatively correlated with each response trait (initial, booster, and overall response). Calves that were weaned at initial vaccination had greater overall antibody response to BVDV1 and BVDV2 compared with calves weaned 3 wk before initial vaccination. In contrast, calves given their initial vaccination 3 wk before weaning had greater overall antibody response to BRSV and BHV1 compared with calves that were vaccinated at weaning. Furthermore, the circulating antibody titer at which each virus needed to be below for an individual calf to positively respond to vaccination was determined (log titer of 0.38 for BVDV1, 1.5 for BVDV2, 3.88 for BRSV, and 1.5 for BHV1). This information can be used to improve vaccination protocols to allow for a greater response rate of individuals to vaccination and, hopefully, improved protection.

Epistatic interactions associated with fatty acid concentrations of beef from angus sired beef cattle.

BACKGROUND: Consumers are becoming increasingly conscientious about the nutritional value of their food. Consumption of some fatty acids has been associated with human health traits such as blood pressure and cardiovascular disease. Therefore, it is important to investigate genetic variation in content of fatty acids present in meat. Previously publications reported regions of the cattle genome that are additively associated with variation in fatty acid content. This study evaluated epistatic interactions, which could account for additional genetic variation in fatty acid content. RESULTS: Epistatic interactions for 44 fatty acid traits in a population of Angus beef cattle were evaluated with EpiSNPmpi. False discovery rate (FDR) was controlled at 5 % and was limited to well-represented genotypic combinations. Epistatic interactions were detected for 37 triacylglyceride (TAG), 36 phospholipid (PL) fatty acid traits, and three weight traits. A total of 6,181, 7,168, and 0 significant epistatic interactions (FDR < 0.05, 50-animals per genotype combination) were associated with Triacylglyceride fatty acids, Phospholipid fatty acids, and weight traits respectively and most were additive-by-additive interactions. A large number of interactions occurred in potential regions of regulatory control along the chromosomes where genes related to fatty acid metabolism reside. CONCLUSIONS: Many fatty acids were associated with epistatic interactions. Despite a large number of significant interactions, there are a limited number of genomic locations that harbored these interactions. While larger population sizes are needed to accurately validate and quantify these epistatic interactions, the current findings point towards additional genetic variance that can be accounted for within these fatty acid traits.

Evaluation of response to bovine viral diarrhea virus type 2 vaccination and timing of weaning on yearling ultrasound body composition, performance, and carcass quality traits in Angus calves.

There are concerns about antagonisms between immunity and animal productivity in livestock production. The objective of this study was to evaluate the effect of antibody levels through a response to vaccination protocol, weaning timing, and their interaction on performance and carcass quality traits in Angus beef cattle. Final antibody level and response to vaccination were based on neutralizing serum antibodies against bovine viral diarrhea virus type 2 (BVDV2). Calves were followed through development and the feedlot phase, with collection of yearling ultrasound (n=957), preharvest (n=762), and carcass (n=673) data. In this study, 48% of the animals were observed to have positively responded to the vaccine, as evidenced by higher final antibody levels compared to prevaccination antibody levels. Increased final antibody levels were significantly (P<0.05) associated with increased yearling weight and increased subcutaneous fat over the rump. An interaction between final antibody level and weaning time also was associated (P<0.05) with Warner-Bratzler shear force (WBSF) and meat pH, with a favorable, negative relationship between final antibody and WBSF in calves weaned at initial vaccination. Overall antibody response by wean time interaction had a significant (P<0.05) association with ADG and meat pH, with calves weaned at initial vaccination having a favorable, positive relationship between overall antibody response and ADG. Under both the final antibody and overall antibody response models, animals weaned at initial vaccination had significantly (P<0.05) lower intramuscular fat at yearling time and conversely higher harvest weight than animals weaned at the booster vaccination. When antibody response was grouped (none, low, high), a significant interaction (P<0.05) between antibody response group and weaning time was identified for ADG, harvest weight, and HCW. Animals weaned at the initial vaccination in the high antibody response group had the advantage for ADG, harvest weight, and HCW compared to animals in the high-response group that were weaned at booster vaccination. Linear increases in antibody response generally did not have negative effects on performance or carcass quality traits in finished cattle (P>0.05). Therefore, producers should not be concerned about decreased production or quality attributes as a result of developing a robust antibody response to vaccination for BVDV2 in beef cattle.

An evaluation of circulating bovine viral diarrhea virus type 2 maternal antibody level and response to vaccination in Angus calves.

Vaccination against viruses has been shown to help prevent bovine respiratory disease in cattle. However, both passively acquired maternal antibody concentration and calf age have been shown to impact the ability of the immune system of a calf to respond to vaccination. The objectives of this study were to identify and evaluate environmental and management factors that affect 1) passively acquired bovine viral diarrhea virus (BVDV) type 2 antibody level, 2) decay rate of passively acquired BVDV type 2 antibody level, and 3) responses to BVDV type 2 vaccinations. A 2-shot modified live vaccine was administered to 1,004 Angus calves that were weaned at either the initial vaccination (n = 508) or the booster vaccination (n = 496). Calves weaned at the initial vaccination averaged 139 d whereas calves weaned at booster vaccination averaged 128 d of age. Bovine viral diarrhea virus type 2 antibodies were measured in 3 approximately 21-d intervals, serially collected serum samples to quantify antibody levels at initiation and end of vaccination protocol in addition to responses to initial, booster, and overall vaccination protocol. Amount of passively transferred antibody in the calf increased as dam age increased from 2 to 6 yr (P < 0.05) with no differences after dams reached 6 yr (P > 0.05). Calf age nested within birth year-season and dam age affected both initial and final antibody level, initial response, booster response, and overall antibody response to vaccination. The level of circulating, passively acquired maternal antibodies present at the time of vaccination had a significant (P < 0.05) negative effect on antibody responses to vaccination (initial response, booster response, and overall response). Calves that were weaned at the time of initial vaccination had significantly (P < 0.05) greater final antibody level, initial response, and overall response to vaccination than animals weaned at booster vaccination. In order for a calf to mount an overall antibody response to vaccination, maternal antibodies in circulation need to be less than 3.12 titers. However, the age at which a calf reached this antibody threshold was dependent on dam age. This information will help cattle managers and consultants design vaccination protocols to successfully mount an antibody response to vaccination.

Genome-wide association study of concentrations of iron and other minerals in longissimus muscle of Angus cattle.

The objective of this study was to characterize variation and identify SNP and chromosomal regions associated with mineral concentrations in LM of Angus beef cattle. Samples of LM from 2,285 Angus cattle were obtained, and concentrations of 7 minerals, including iron, magnesium, manganese, phosphorus, potassium, sodium, and zinc, were quantified. Genomic DNA extracted from the ground beef sample used for mineral composition was genotyped with the Bovine SNP50 Infinium II BeadChip, and effects of SNP on each trait were estimated using the Bayes-Cπ module of GenSel software. Pedigree-based estimates of heritabilities and corresponding genetic variances indicate iron was the only mineral concentration that could be considered a good candidate for manipulation by genomic selection. The amount of variation that could be accounted for by SNP genotypes was concordant with pedigree-based heritabilities and varied from very low for potassium and sodium (<0.09) to medium high (0.37) for iron. Although significant chromosomal regions were identified for all minerals analyzed in this study, further study focused on iron. Seven regions on 6 chromosomes (1, 2, 7, 10, 15, and 28) were identified to have a major effect on the iron content of LM in Angus cattle. The accuracy of direct genomic values (DGV) for iron concentration was estimated using a 5-fold cross-validation strategy. The accuracy of DGV estimated as the genetic correlation between DGV and the phenotype (iron concentration) adjusted for contemporary groups was 0.59. A bivariate animal model was used to estimate genetic correlations between iron concentrations and a reduced set of economically important carcass traits: HCW, rib eye area, calculated USDA yield grade, percent KPH, and marbling score. The genetic correlations between iron concentration and HCW, percentage KPH, marbling score, and rib eye area were small (-0.19 to 0.15) and nonsignificant. Although still weak (0.22), a positive significant genetic correlation was identified between iron content and USDA calculated yield grade. Beef is a major contributor of iron and zinc in the human diet, and this study found that iron content might be effectively manipulated through marker-assisted selection programs, without compromising other carcass and palatability traits.

Genetic parameters for concentrations of minerals in longissimus muscle and their associations with palatability traits in Angus cattle.

The objective of this study was to estimate genetic parameters for concentrations of minerals in LM and to evaluate their associations with beef palatability traits. Samples of LM from 2,285 Angus cattle were obtained and fabricated into steaks for analysis of mineral concentrations and for trained sensory panel assessments. Nine minerals, including calcium, copper, iron, magnesium, manganese, phosphorus, potassium, sodium, and zinc, were quantified. Restricted maximum likelihood procedures were used to obtain estimates of variance and covariance components under a multiple-trait animal model. Estimates of heritability for mineral concentrations in LM varied from 0.01 to 0.54. Iron and sodium were highly and moderately heritable, respectively, whereas the other minerals were lowly heritable except for calcium, copper, and manganese, which exhibited no genetic variation. Strong positive genetic correlations existed between iron and zinc (0.49, P < 0.05), between magnesium and phosphorus (0.88, P < 0.05), between magnesium and sodium (0.68, P < 0.05), and between phosphorus and potassium (0.69, P < 0.05). Overall tenderness assessed by trained sensory panelists was positively associated with manganese, potassium, and sodium and negatively associated with phosphorus and zinc concentrations (P < 0.05). Juiciness assessed by trained sensory panelists was negatively associated with magnesium and positively associated with manganese and sodium concentrations (P < 0.05). Livery or metallic flavor was not associated with any of the minerals (P > 0.05). Beefy flavor was positively associated with calcium, iron, and zinc and negatively associated with sodium concentration, whereas a painty or fishy flavor was positively associated with sodium and negatively associated with calcium and potassium concentrations (P < 0.05). Beef is a major contributor of iron and zinc in the human diet, and these results demonstrate sufficient genetic variation for these traits to be improved through marker-assisted selection programs without compromising beef palatability.

Genetic parameters for carnitine, creatine, creatinine, carnosine, and anserine concentration in longissimus muscle and their association with palatability traits in Angus cattle.

The objective of this study was to estimate genetic parameters for carnitine, creatine, creatinine, carnosine, and anserine concentration in LM and to evaluate their associations with Warner-Bratzler shear force (WBSF) and beef palatability traits. Longissimus muscle samples from 2,285 Angus cattle were obtained and fabricated into steaks for analysis of carnitine, creatine, creatinine, carnosine, anserine, and other nutrients, and for trained sensory panel and WBSF assessments. Restricted maximum likelihood procedures were used to obtain estimates of variance and covariance components under a multiple-trait animal model. Estimates of heritability for carnitine, creatine, creatinine, carnosine, and anserine concentrations in LM from Angus cattle were 0.015, 0.434, 0.070, 0.383, and 0.531, respectively. Creatine, carnosine, and anserine were found to be moderately heritable, whereas almost no genetic variation was observed in carnitine and creatinine. Moderate positive genetic (0.25, P < 0.05) and phenotypic correlations (0.25, P < 0.05) were identified between carnosine and anserine. Medium negative genetic correlations were identified between creatine and both carnosine (-0.53, P < 0.05) and anserine (-0.46, P < 0.05). Beef and livery/metallic flavor were not associated with any of the 5 compounds analyzed (P > 0.10), and carnitine concentrations were not associated (P > 0.10) with any of the meat palatability traits analyzed. Carnosine was negatively associated with overall tenderness as assessed by trained sensory panelists. Similar negative associations with overall tenderness were identified for creatinine and anserine. Painty/fishy was the only flavor significantly and negatively associated with creatinine and carnosine. These results provide information regarding the concentration of these compounds, the amount of genetic variation, and evidence for negligible associations with beef palatability traits in LM of beef cattle.

Genetic polymorphisms in bovine transferrin receptor 2 (TFR2) and solute carrier family 40 (iron-regulated transporter), member 1 (SLC40A1) genes and their association with beef iron content.

Beef is considered to be an excellent source of dietary iron. However, little is known about the genetic control of beef iron content. We hypothesized that genetic polymorphisms in transferrin receptor 2 (TFR2) and solute carrier family 40 (iron-regulated transporter), member 1 (SLC40A1) could influence skeletal muscle iron content. The objective of this study was to use Angus cattle to identify single-nucleotide polymorphisms (SNPs) in the exons and flanking regions of the bovine TFR2 and SLC40A1 genes and to evaluate the extent to which genetic variation in them was associated with bovine longissimus dorsi muscle iron content. Ten novel SNPs were identified in TFR2, of which one SNP tended to be associated (P < 0.013) with skeletal muscle iron content. Nine novel SNPs in SLC40A1, NC007300: rs133108154, rs137140497, rs135205621, rs136600836, rs134388440, rs136347850, rs134186279, rs134621419 and rs137555693, were identified, of which SNPs rs134388440, rs136347850 and rs137555693 were significantly associated (P < 0.007) with skeletal muscle iron content. High linkage disequilibrium was observed among SLC40A1 SNPs rs134388440, rs136347850 and rs137555693 (R(2) > 0.99), from which two haplotypes, TGC and CAT, were defined. Beef from individuals that were homozygous for the TGC haplotype had significantly (P < 0.001) higher iron content than did beef from CAT homozygous or heterozygous individuals. The estimated size of effect of the identified haplotypes was 0.3% of the phenotypic variance. In conclusion, our study provides evidence for genetic control of beef iron concentration. Moreover, SNPs identified in SLC40A1, rs134388440, rs136347850 and rs137555693 might be useful markers for the selection of Angus cattle for altered iron content.

Mapping genetic loci that interact with myostatin to affect growth traits.

Myostatin, or GDF8, is an inhibitor of skeletal muscle growth. A non-functional myostatin mutation leads to a double muscling phenotype in some species, for example, mice, cattle and humans. Previous studies have indicated that there are loci in the genome that interact with myostatin to control backfat depth and other complex traits. We now report a quantitative trait loci (QTL) mapping study designed to identify loci that interact with myostatin to impact growth traits in mice. Body weight and average daily gain traits were collected on F2 progeny derived from a myostatin-null C57BL/6 strain by M16i cross. In all, 44 main effect QTL were detected above a 5% genome-wide significance threshold when an interval mapping method was used. An additional 37 QTL were identified to significantly interact with myostatin, sex or reciprocal cross. A total of 12 of these QTL interacted with myostatin genotype. These results provide a foundation for the further fine mapping of genome regions that harbor loci that interact with myostatin.

Estimation of relationships between mineral concentration and fatty acid composition of longissimus muscle and beef palatability traits.

The objective of this study was to determine the influence of beef LM nutrient components on beef palatability traits and evaluate the impact of USDA quality grade on beef palatability. Longissimus muscle samples from related Angus cattle (n = 1,737) were obtained and fabricated into steaks for trained sensory panel, Warner-Bratzler shear force (WBSF), lipid oxidation measured by thiobarbituric acid reactive substances (TBARS), fatty acid, and mineral composition analysis. Pearson phenotypic correlations were obtained by the correlation procedure of SAS. Beef palatability data were analyzed by the GLM procedure of SAS with USDA quality grade as the main effect. Specific mineral concentrations did not demonstrate strong correlations with WBSF or sensory traits (r = -0.14 to 0.16). However, minerals appeared to have a stronger relationship with flavor; all minerals evaluated except Ca and Mn were positively correlated (P < 0.05) with beef flavor. Stearic acid (C18:0), C18:2, C20:4, and PUFA were negatively correlated (P < 0.05) with all 3 panelist tenderness traits (r = -0.09 to -0.22) and were positively correlated (P < 0.05) with WBSF (r = 0.09 to 0.15). The MUFA were positively correlated (P < 0.05) with panelist tenderness ratings (r = 0.07 to 0.10) and negatively associated (P < 0.05) with WBSF (r = -0.11). The strongest correlations with juiciness were negative relationships (P < 0.05) with C18:2, C18:3, C20:4, and PUFA (r = -0.08 to -0.20). Correlations with beef flavor were weak, but the strongest was a positive relationship with MUFA (r = 0.13). Quality grade affected (P < 0.05) WBSF, TBARS, and all trained sensory panel traits, except livery/metallic flavor. As quality grade increased, steaks were more tender (P < 0.05), as evidenced by both WBSF and sensory panel tenderness ratings. Prime steaks were rated juiciest (P < 0.05) by panelists, whereas Select and Low Choice were similarly rated below Top Choice for sustained juiciness. Quality grade influenced (P < 0.05) beef flavor, but not in a linear fashion. Although there were significant correlations, these results indicate tenderness, juiciness, and flavor are not strongly influenced by individual nutrient components in beef LM. Furthermore, the positive linear relationships between USDA quality grade and beef palatability traits suggest quality grade is still one of the most valuable tools available to predict beef tenderness.

Early postmortem biochemical factors influence tenderness and water-holding capacity of three porcine muscles.

The objective of this study was to determine whether differences in pork tenderness and water-holding capacity could be explained by factors influencing calpain activity and proteolysis. Halothane-negative (HAL-1843 normal) Duroc pigs (n = 16) were slaughtered, and temperature and pH of the longissimus dorsi (LD), semimembranosus (SM), and psoas major (PM) were measured at 30 and 45 min and 1, 6, 12, and 24 h postmortem. Calpastatin activity; mu-calpain activity; and autolysis and proteolysis of titin, nebulin, desmin, and troponin-T were determined on muscle samples from the LD, SM, and PM at early times postmortem. Myofibrils from each muscle were purified to assess myofibril-bound (mu-calpain. Percentage drip loss was determined, and Warner-Bratzler shear (WBS) force was analyzed. Myosin heavy-chain (MHC) isoforms were examined using SDS-PAGE. The pH of PM was lower (P < 0.01) than the pH of LD and SM at 30 and 45 min and 1 h postmortem. The PM had a higher (P < 0.01) percentage of the MHC type IIa/IIx isoforms than the LD. The-LD had the greatest proportion of (P < 0.01) MHC IIb isoforms of any of the muscles. The PM had the lowest (P < 0.01) percentage of MHC IIb isoforms and a greater (P < 0.05) percentage of type I MHC isoforms than the LD and SM. The PM had less (P < 0.01) drip loss after 96 h of storage than the SM and LD. The PM had more desmin degradation (P < 0.01) than the LD and SM at 45 min and 6 h postmortem. Degradation of titin occurred earlier in the PM than the LD and SM. At 45 min postmortem, the PM consistently had some autolysis of mu-calpain, whereas the LD and SM did not. At 6 h postmortem, some autolysis of mu-calpain (80-kDa subunit) was observed in all three muscles. The rapid pH decline and increased rate of autolysis in the PM paralleled an earlier appearance of myofibril-bound mu-calpain. The SM had higher calpastatin activity (P < 0.05) at 45 min, 6 h, and 24 h and had higher WBS values at 48 h (P < 0.01) and 120 h (P < 0.05) postmortem than the LD. At 48 and 120 h postmortem, more degradation of desmin, titin, and nebulin were observed in the LD than in the SM. These results show that mu-calpain activity, mu-calpain autolysis, and protein degradation are associated with differences in pork tenderness and water-holding capacity observed in different muscles.