key: cord-0948775-kvstm9a6
authors: Cockrum, R. R.; Speidel, S. E.; Salak-Johnson, J. L.; Chase, C. C. L.; Peel, R. K.; Weaber, R. L.; Loneagan, G. H.; Wagner, J. J.; Boddhireddy, P.; Thomas, M. G.; Prayaga, K.; DeNise, S.; Enns, R. M.
title: Genetic parameters estimated at receiving for circulating cortisol, immunoglobulin G, interleukin 8, and incidence of bovine respiratory disease in feedlot beef steers(1)
date: 2016-07-03
journal: J Anim Sci
DOI: 10.2527/jas.2015-0222
sha: d109b3d1f5ec4a0b344c6f34c1dc389a68afaccc
doc_id: 948775
cord_uid: kvstm9a6

Bovine respiratory disease complex (i.e., shipping fever and bacterial bronchopneumonia) is a multifaceted respiratory illness influenced by numerous environmental factors and microorganisms. Bovine respiratory disease (BRD) is just one component of BRD complex. Because BRD is moderately heritable, it may be possible to reduce the incidence of BRD through genetic selection. The objectives of this study were to determine the heritability and associative genetic relationships among immune system traits (i.e., cortisol, total IgG, IgG isotypes, and IL-8) in cattle monitored for BRD incidence. At an average of 83 d after weaning (219 d age and mean = 221.7 kg [SD 4.34]), crossbred Bos taurus steer calves (n = 2,869) were received at a commercial feedlot in southeastern Colorado over a 2-yr period. At receiving, jugular blood samples were collected at 212 (yr 1) and 226 d (yr 2) of age for immune trait analyses. The BRD phenotype was defined as a binomial variable (0 = no and 1 = yes) and compared with immune system traits measured at receiving (prior to illness onset). An animal identified as BRD positive exhibited ≥ 2 clinical signs (i.e., eye or nasal discharge, cough, lethargy, rapid breathing, acute interstitial pneumonia, or acute upper respiratory syndrome and/or a rectal temperature > 39.7°C). Heritability and genetic correlation estimates for categorical variable BRD, cortisol, IgG, IgG1, IgG2, and IL-8 were estimated from a sire model using ASREML. Heritability estimates were low to moderate for BRD (0.17 ± 0.08), cortisol (0.13 ± 0.05), IgG (0.15 ± 0.05), IgG1 (0.11 ± 0.05), IgG2 (0.24 ± 0.06), and IL-8 (0.30 ± 0.06). A moderate negative genetic correlation was determined between BRD and cortisol (r(g) = −0.19 ± 0.32). Moderate positive correlations were found between BRD with IgG (0.42 ± 0.28), IgG1 (0.36 ± 0.32), and IL-8 (r(g) = 0.26 ± 0.26). Variation in the BRD phenotype and immune system traits suggested herd health improvement may be achieved through genetic selection.

Prevention and treatment of bovine respiratory disease (BRD) cost the cattle industry over US$700 million annually (NASS, 2011) . Although more than 85% of feedlots in the United States use a vaccination protocol for viruses that contribute to the develop ment of BRD, the number of cases has not decreased (Gifford et al., 2012; APHIS, 2013) . Therefore, ef fective prevention and/or early detection strategies for BRD must be developed to reduce its incidence in cattle. Prevention and diagnosis of BRD is further aBSTRaCT: Bovine respiratory disease complex (i.e., shipping fever and bacterial bronchopneumonia) is a multifaceted respiratory illness influenced by numer ous environmental factors and microorganisms. Bovine respiratory disease (BRD) is just one component of BRD complex. Because BRD is moderately heritable, it may be possible to reduce the incidence of BRD through genetic selection. The objectives of this study were to determine the heritability and associative genetic rela tionships among immune system traits (i.e., cortisol, total IgG, IgG isotypes, and IL8) in cattle monitored for BRD incidence. At an average of 83 d after weaning (219 d age and mean = 221.7 kg [SD 4.34] ), crossbred Bos taurus steer calves (n = 2,869) were received at a commercial feedlot in southeastern Colorado over a 2yr period. At receiving, jugular blood samples were col lected at 212 (yr 1) and 226 d (yr 2) of age for immune trait analyses. The BRD phenotype was defined as a binomial variable (0 = no and 1 = yes) and compared with immune system traits measured at receiving (prior to illness onset). An animal identified as BRD positive exhibited ≥ 2 clinical signs (i.e., eye or nasal discharge, cough, lethargy, rapid breathing, acute interstitial pneu monia, or acute upper respiratory syndrome and/or a rectal temperature > 39.7°C). Heritability and genetic correlation estimates for categorical variable BRD, cor tisol, IgG, IgG1, IgG2, and IL8 were estimated from a sire model using ASREML. Heritability estimates were low to moderate for BRD (0.17 ± 0.08), cortisol (0.13 ± 0.05), IgG (0.15 ± 0.05), IgG1 (0.11 ± 0.05), IgG2 (0.24 ± 0.06), and IL8 (0.30 ± 0.06). A moderate negative genetic correlation was determined between BRD and cortisol (r g = −0.19 ± 0.32). Moderate posi tive correlations were found between BRD with IgG (0.42 ± 0.28), IgG1 (0.36 ± 0.32), and IL8 (r g = 0.26 ± 0.26). Variation in the BRD phenotype and immune sys tem traits suggested herd health improvement may be achieved through genetic selection. complicated by the complexity of this disease. Stress, inadequate nutrition, poor management, and/or ge netics predispose cattle to polymicrobial respiratory infections (Ellis, 2001; Duff and Galyean, 2007) , in creasing an animal's susceptibility to contracting BRD.

Because BRD is a complex trait that is heritable, economically relevant, and difficult to measure, ge nomic selection could be used to reduce its incidence. However, to improve accuracy of selection, the ge netic component of BRD must be delineated. Immune system traits such as cortisol, IgG and its subclasses (IgG1 and IgG2), and IL-8 all influence an individ ual's response to stress and immune system function (Stanley et al., 1984; Hull et al., 2001; SalakJohnson and McGlone, 2007) . Deciphering the genetic rela tionships among immune system traits will aid in pre dicting an individual's likelihood to develop BRD. We hypothesized that immune system traits measured at receiving were genetically correlated with BRD. The longterm aim of this research is to elucidate both ge notypic and phenotypic relationships of BRD with physiological, behavioral, and performance traits. Specifically, the objectives of this study were to 1) estimate the genetic variation of BRD and immune system traits (i.e., cortisol, IgG, and IL8), 2) evaluate the genetic relationships with BRD and immune traits measured prior to the onset of illness, and 3) assess the genetic relationships within immune system traits.

The Colorado State University Institutional Animal Care and Use Committee approved all animal proce dures. Crossbred Bos taurus steer calves (n = 2,869) were received at a commercial feedlot located in south eastern Colorado from a ranch in western Nebraska over 2 consecutive years. All cattle were shipped approxi mately 536 km and were received over a 3d period for each year. Descriptive statistics by year of steers studied are provided in Table 1 . To minimize the effects of ship ping stress, cattle were kept overnight in receiving pens, provided hay and water, and processed the following day.

During processing, steers received a radio frequen cy identification tag, a visual identification tag, and a lot tag. Steers were also vaccinated, dewormed, implanted with a growth promotant, and measured for BW and be havior; additionally, blood samples were collected. Prior to shipping, cattle were vaccinated in Nebraska against bovine viral diarrhea, Pasteurella multocida, and Mannheimia haemolytica (Pyramid 2 + Type II BVD and Presponse SQ; Boehringer Ingelheim Vetmedica Inc., St. Joseph, MO) before weaning and again at weaning. Year 1 cattle were not vaccinated on arrival at the feedlot. However, at the request of the commercial feedlot owners in yr 2, cattle were vaccinated on feedlot entry for common viral and bacterial agents associated with BRD: bovine viral diarrhea, Pasteurella multocida, and Mannheimia haemolytica (Pyramid 2 + Type II BVD and Presponse SQ). This difference in vaccination protocols resulted in 45% of steers developing BRD in yr 1 and 7% in yr 2 (Table 1) , which more accurately reflects the range of vaccination procedures in the beef industry (i.e., no vaccinations versus multiple vaccina tions). All steers received oral and pouron parasiticides, oxfendazole (Synanthic Suspension 9.06%; Boehringer Ingelheim Vetmedica Inc.) and ivermectin (ProMectin B Pouron; Vedco, Inc., St. Joseph, MO) on arrival to the feedlot. In yr 1, growth promotants containing 80 mg trenbolone acetate and 16 mg estradiol (RevalorIS; Intervet/ScheringPlough, Millsboro, DE) were admin istered per feedlot recommendation with a reimplant after 72 d on feed. In yr 2, cattle were implanted with a 200 d delayed releasing growth promotant of 200 mg trenbolone acetate and 40 mg estradiol (Revalor XS; Intervet/ScheringPlough). Collection of repeated phe notypic performance data (i.e., chute score, flight time, processing time, weight, and growth measurements via ultrasound) occurred on arrival and on d 73 and 130 in yr 1 and d 85 and 148 in yr 2. Approximately 10 mL of blood was collected via jugular venipuncture using va cutainers containing EDTA for wholegenome genotyp ing. Whole blood was shipped to Zoetis (Kalamazoo, MI) for DNA extraction and genotyping using the Bovine SNP50 BeadChip (Illumina Inc., San Diego, CA). Plasma was isolated from whole blood samples and frozen at −20°C until plasma immune trait analyses.

Commercial feedlot personnel monitored steer health on a daily basis. Animals identified as either sick or lame were removed and relocated to the Colorado State University South Eastern Colorado Research Center (SeCRC; Fort Collins, CO), adjacent to the com mercial feedlot, for further evaluation. A dirt lane of 6.10 m separated these 2 facilities. Steers identified as BRD positive had to exhibit ≥ 2 clinical signs and/or rectal temperature > 39.7°C. Clinical signs included eye dis charge, nasal discharge, cough, lethargy, rapid breathing, and atypical interstitial pneumonia or acute upper respi ratory syndrome. Rectal temperature was regarded to be an overall determining factor because it is less subjective than visual appraisal alone. Based on these criteria, steers were binomially classified as treated (1; n = 792) versus not treated (0; n = 2,077) for BRD.

Treatment protocols were administered in accor dance with the commercial feedlot and SECRC stan dard operating procedures. Steers with a rectal tem perature below 39.7°C and displaying no clinical signs between d 5 and 7 after treatment of tulathromycin (Draxxin; Pfizer Animal Health, Zoetis, Kalamazoo, MI) were considered recovered and returned to their respective pens at the commercial feedlot. Steers with rectal temperatures above 39.7°C or exhibiting mul tiple clinical signs of BRD within 7 d of the initial treatment received a second BRD treatment of fluoro quinolone (Baytril; Bayer Animal Health, Leverkusen, Germany). Steers were monitored for 5 d, reevaluated, and treated a third time using oxytetracycline (Bio Mycin 200; Boehringer Ingelheim Vetmedica Inc.) if exhibiting BRD signs. Animals that were previously treated and recovered from BRD but again developed clinical signs of BRD continued on a treatment proto col based on the last treatment record. Steers requiring >3 treatments were removed from the study with fur ther treatment protocols made at the discretion of the commercial feedlot. All other illnesses lacking suffi cient evidence of BRD were treated as diagnosed (e.g., bloat, lameness, and pinkeye). Similar to the feedlot, while at the SECRC, steers with clinical signs of BRD were fed a commercial diet (66 Mcal/hundredweight DM). Steers were also fed longstemmed hay to ac climate them to the bunk feeding.

Plasma isolated from jugular blood collected on receiving was used to determine circulating cortisol, IL8, IgG, IgG1, and IgG2 concentrations. Plasma samples were measured for cortisol and IL8 at the University of Illinois (UrbanaChampaign, IL). Plasma cortisol was measured using a commercially available RIA kit (CoatACount; Siemens Healthcare, Malvern, PA) following manufacturer procedures. Intra and interassay CV were 4.5 and 8.5%, respec tively, and the minimum detectable concentration was 100 ng/dL. Interleukin 8 was measured in duplicate us ing a commercially available ELISA (R&D Systems, Inc., Minneapolis, MN). Circulating IgG, IgG1, and IgG2 were measured using an ELISAbased method at South Dakota State University according to the manufacturer's instructions (Bethyl Laboratories, Inc., Montgomery, TX). Immunoglobulin levels greater than the standard curve were subsequently remeasured.

Descriptive statistics were generated using the boot and pastecs packages in R (R Core Team, 2013). Observation numbers, means, SE, SD, and minimum and maximum values are reported. Values outside the normal range were determined for IgG, IgG1, and IgG2 using Cook's distance > 4/n. Measures for IgG > 100 mg/mL and IgG1 and IgG2 > 75 mg/mL were removed, as values were outside the normal range and were biologically irrelevant according to estimates re ported in literature (Butler, 1983; Gilbert et al., 1988) . All values < 0 mg/mL were adjusted to 0 to account for nonspecific binding from the ELISA kit. After val ues outside the normal range were removed, plasma measurements for cortisol, IgG, IgG1, IgG2, and IL8 were within the expected range.

A singletrait sire model presented below was used to determine heritability for BRD, cortisol, IgG, IgG1, IgG2, and IL8 was

Specifically, y was a vector of observations, X was a known incidence matrix relating contemporary group and process order effects in b to observations in y, Z was a known incidence matrix relating random sire effects in s to observations in y, and e was a vector of random re sidual errors. Random effects in the model were assumed to have 0 means and variances represented as follows:

in which a was Wright's numerator relationship matrix and I was an identity matrix whose order was equal to the number of observations in y. For BRD, a binary trait, a probit link function was used, thereby implic itly defining a residual variance of 1. Genetic correla tions were estimated for each of the 15 trait pairs using the bivariate form of the sire model presented above with covariances fitted between traits.

The genetic parameters for all traits and their SE were estimated using ASREML (version 3.0; VSN International, Ltd., Hemel Hempstead, UK), which uses an average information REML algorithm. Additive genetic variance was calculated as 4 times the sire variance estimated through ASREML.

Contemporary group was included as a fixed effect in the above models with process order within day used as a covariate. For dependent variables measured on re ceiving (i.e., cortisol, immunoglobulins, and IL8), con temporary group consisted of ranch unit of origin (n = 3) and entry process date (n = 6). The contemporary group for BRD, a postreceiving trait, included ranch unit of or igin (n = 3), feedlot pen (n = 11), and entry process date (n = 6). Half of the steers were commercial (n = 1,395) and the other half was seed stock (n = 1,469); therefore, only the commercial steers were backgrounded before shipping. Commercial versus seed stock steers were not commingled; therefore, the fixed effect pen would have accounted for a steer's backgrounding status.

Calves in the study were the result of multiple sire breeding pastures. Sire identification was performed via DNA sampling of both sire and progeny using a commercial DNA laboratory. A sparse historical pedi gree was obtained from the ranch of origin and used to construct a 5generation pedigree starting with the sires of animals in the study. There were 308 unique sires with progeny that had data and 548 total sires in the full pedigree. In this 5generation pedigree, dam information was extremely scarce. Therefore, for the purpose of estimating heritability and genetic correla tion, a sire model was used to best describe the flow of genetic material from one generation to next.

The average weight for steers (n = 2,864) entering the feedlot was 221.75 kg (SD 24.34). In yr 1, 1,551 calves (223.90 kg [SD 23.46] ) were received at an av erage age of 212 d, and in yr 2, 1,318 calves (219.10 kg [SD 25.11] ) were received at an average age of 226 d (Table 1) . Steers were treated for BRD an average of 36 d earlier in yr 1 (SD 20.47) compared with yr 2 (SD 56.40). Additionally, a total of 698 steers were treated for BRD in yr 1 and 94 in yr 2. Descriptive statistics with corresponding SE and SD for the im mune system traits determined prior to steers develop ing BRD is provided in Table 2 . Total IgG was greater than IgG1, which was greater than IgG2. Values out side the normal range for IgG (measure > 100 mg/mL), IgG1 (measure > 75 mg/mL), and IgG2 (measure > 75 mg/mL) did not overlap among immunoglobulins. Biologically irrelevant immunoglobulin samples re moved from analyses were attributed to sample con tamination. Only 5, 6, and 6 samples for IgG, IgG1, and IgG2, respectively, were removed from analy sis. Both IgG2 (0.24 ± 0.06) and IL8 (0.30 ± 0.06; Table 2 ) had moderate heritability estimates and low heritability estimates were determined for BRD, corti sol, IgG, and IgG1 (h 2 ≤ 0.17 ± 0.08).

Moderate positive genetic correlations (r g ≥ 0.36 ± 0.32) were detected between BRD with IgG and IgG1 (Table 3) . Also, there was a nonsignificant weak posi tive correlation between IL8 and BRD. Conversely, BRD status was negatively genetically correlated with plasma cortisol (r g = −0.19 ± 0.32; Table 3 ).

A low genetic correlation between both IgG and IgG1 with cortisol (r g ≤ −0.19 ± 0.26) as well as a weak relationship with IL8 (r g = −0.08 ± 0.29) and IgG2 (r g = 0.06 ± 0.22) with cortisol was detected. There was a weak negative genetic correlation between IL8 and IgG2 (r g = −0.11 ± 0.17). As expected, a strong positive genetic correlation was observed among IgG, IgG1, and IgG2 (r g ≥ 0.83 ± 0.09). Greater than 50% of total IgG is composed of IgG1 and 20 to 40% is composed of IgG2 (Hashira et al., 2000) . 

The average steer in the United States weighs ap proximately 200 to 320 kg at weaning (6 to 10 mo age).

Calves used in the current study were slightly smaller than average at weaning (158 ± 26.48), which was at tributed to calves being only 4.5 mo of age. Because calves were younger and lighter, calves were back grounded until weights were adequate for finishing. Only half of the steers used in this study had weight and backgrounding information because half of the calves were seed stock and the other half were commercial. Upon receiving, calves were approximately 7 mo of age and had weights were similar to that of the national average at that age (Table 1 ). The majority of steers re mained healthy throughout the trial, with a little more than onethird (38%) treated for BRD. In yr 1 (not vac cinated), 45% of steers developed BRD compared with 7% in yr 2 (vaccinated). Generally, 16.2% of U.S. cattle entering a feedlot develop respiratory disease (APHIS, 2013). Rates of BRD in yr 1 were much higher than av erage; however, none of the steers were vaccinated for the BRD complex. Year 2 BRD rates were well below the average 87.5% of cattle vaccinated for BRD that still develop BRD (APHIS, 2013).

In Table 2 , more variation was observed in corti sol levels in the current study than in previous research (7.00 to 29.50 ng/mL; Bristow and Holmes, 2007) . Immunoglobulin G and isotype levels were lower in these steers than in beef calves from previous research (Waldner and Rosengren, 2009 ). Similar to our find ings, Kimura et al. (2002) determined plasma IL8 lev els in dairy cattle with retained placentas ranged from 134 to 300 pg/mL. Number of observations, calf age, breed, and environmental conditions of the stressor (i.e., shipping versus a squeeze chute) all differed with each study making it difficult to compare immune trait levels across studies. Limited information is available on im mune system traits in receiving beef steers; therefore, it was difficult to assess if plasma levels were typical of this age or if differing levels were attributed to environ mental conditions. Immune system trait plasma levels were within the expected range of the kits used, but it should be noted that levels varied greatly among steers.

Presented in Table 4 are previously reported heri tability estimates for BRD, cortisol, total IgG, IgG1, and IgG2 for comparison with those reported in this study. Previous heritability estimates for BRD confirm the estimate (h 2 = 0.17 ± 0.08) in the current research (Snowder et al., 2005 (Snowder et al., , 2006 (Snowder et al., , 2007 Schneider et al., 2010) . It is possible that heritability estimates reported in this study were inflated due to dominance; however, given the structure of the data and the unknown breed composition of the dam, these effects were unable to be accounted for. These results provide evidence to suggest that improvement on BRD is possible through genetic selection. However, accuracy of heritability estimates may be improved for BRD through a more accurate phenotype, increased incidence reporting for individu als, development and increased use of genomic mark ers, and/or incorporating genetic relationships through multiple trait evaluation into estimates of genetic merit.

Given the moderate heritability and genetic cor relation estimates between BRD and cortisol, plasma concentration of cortisol may serve as an index variable that can be used for indirect selection for BRD or used to create an index for BRD susceptibility. Limited re search was previously published describing the herita bility of cortisol; however, stress research in pigs dem onstrated cortisol is moderately heritable (h 2 = 0.40; Kadarmideen and Janss, 2007) . Although steers in the current study were not intentionally stressed, the act of shipping creates a stressed phenotype in some animals. A low to moderate heritability for cortisol suggested that Table 3 . Heritability estimates (h 2 ; diagonal) and genetic correlations (r g ; above diagonal) with SE for bovine respiratory disease (BRD), cortisol, immunoglobulins (IgG, IgG2, IgG1, and IgG2), and IL8 measured at receiving in steers varied in their innate response to environmental stressors, potentially leaving them more or less suscep tible to illness. Based on the negative genetic correla tion with BRD and the moderate heritability, cortisol levels measured after shipping may serve as a potential candidate for selection for BRD susceptibility. Low to moderate heritability estimates were ob served for IgG and isotypes in the current study. Previous reports of heritability for IgG and its subclasses ranged from 0.00 to 0.28 in dairy calves, supporting the current research estimates in beef cattle (Mallard et al., 1983; Detilleux et al., 1994) . Interleukin 8 had a moderate heritability estimate of 0.30 ± 0.06. No known research is available that described the genetic variability of IL8. However, these results provide evidence to suggest that development of a BRD susceptibility index through the incorporation of IgG, corresponding isoforms, and IL8 may decrease the incidence of BRD in cattle. Future research must be conducted to better characterize the phenotype for BRD susceptibility and identify genomic regions associated with BRD.

Results suggested that there was a genetic relation ship between circulating cortisol measured after ship ping and the development of BRD. Initially, the nega tive genetic relationship between cortisol and BRD was surprising, as prior reports indicated that cortisol levels increase after transportation, thereby increasing an ani mal's susceptibility to disease (Jericho, 1979; Slocombe et al., 1984) . However, when exposed to acute stressors, the immune system will react in an immunoprotective manner (SalakJohnson and McGlone, 2007; Carroll and Forsberg, 2007; Dhabhar, 2009 ). In our study, cor tisol was measured after shipping, which would con stitute an acute stressor. Therefore, this may, in part, explain why there was a negative genetic correlation between cortisol and BRD. Additionally, McGlone et al. (1993) determined that there was a moderate pheno typic correlation (r p = 0.35, P = 0.036) between natural killer cell cytotoxicity and plasma cortisol concentra tions. It is possible that cortisol secretion after an acute stressor versus a chronic stressor is a different trait from a genetic standpoint; therefore, phenotypic relation ships observed in previous research may not translate to a genetic observation. Therefore, more research must be conducted to elucidate these effects.

Total IgG, IgG1, and IgG2 are major components of the humoral immune system that are part of the adaptive immune response that protects cattle against bacterial and viral challenges (Stanley et al., 1984) . These data sug gested that selection for total IgG and/or IgG1 could im pact whether a steer developed BRD. Interestingly, there was no apparent genetic relationship between IgG2 and BRD. In cattle, IgG levels have been reported to decrease after shipping (Kegley et al., 1997) . Immunoglobulin G1 has been shown to increase in allergytype reactions and indicate a T helper 2 response (Kemeny et al., 1986) , whereas IgG2 increases suggest a stronger T helper 1 response and T cell memory (Estes and Brown, 2002) . Also, in the presence of encapsulated bacteria such as Mannheimia haemolytica, a key bacterium involved in (2007) IgG Holstein 1,293 0.12 ± 0.06 Mazengera et al. (1985) IgG1

Beef cattle -0.05 ± 0.17 Gilbert et al. (1988) IgG2

Holstein 190 0.28 ± 0.22 Mallard et al. (1983) Holstein 137 0.02 ± 0.31 Detilleux et al. (1994) Holstein 1,293 0.10 ± 0.07 Mazengera et al. (1985) Red Danish 758 0.12 ± 0.08 Thode Jensen and Christensen (1975) IL8 ----1 h 2 d = direct genetic effect. 2 Preweaning heritability estimate. BRD, IgG2 was increased. Immunoglobulin G2 antibod ies have also been shown to be associated with reduced incidence of cattle becoming infected with respiratory bovine coronaviruses (Lin et. al., 2000 (Lin et. al., , 2001 . As IgG1 makes up the majority of total IgG and has a differing mode of regulation compared with IgG2, it is plausible that the genetic relationship between IgG1 and IgG2 with BRD would differ. To date, genetic relationships between BRD status and total IgG or isotypes of IgG have not been reported in the literature. These results suggested that negative selection pressure for IgG and IgG1 and positive selection pressure for plasma cortisol may result in progeny less susceptible to developing BRD. This is the first description of a genetic relation ship between BRD and IL8. Interleukin 8 is a cyto kine that serves as a chemoattractant signal, recruiting neutrophils to the sight of infection (Hull et al., 2001) . Phenotypically, Caswell et al. (1998) demonstrated that IL8 expression was upregulated in lesions of pneu monic pasteurellosis in beef calves. It was concluded that IL8 is integral to the recruitment of neutrophils to establish lesions in cattle with pneumonia. In vitro, low levels of Pasteurella haemolytica, a Gramnegative bacterium associated with BRD, induced high concen trations of IL8 (Morsey et al., 1996) . Similar to previ ous phenotypic reports, there appears to be a moderate positive genetic correlation between IL8 and BRD. Therefore, genetic selection against IL8 may decrease a steer's susceptibility to developing BRD. Although a moderate genetic correlation existed between BRD and IL8, further investigation is warranted to deter mine if IL8 is useful as a predictor of BRD.

To our knowledge, this is the first study to evalu ate the genetic relationships among cortisol, IgG, IgG2, and IL8. Human maternal total IgG is com posed of 52% IgG1, 40% IgG2, 2% IgG3, and 2% IgG4 (Hashira et al., 2000) . Similar ratios would be expected for cattle, but immunoglobulins 3 and 4 were not measured; therefore, specific ratios could not be determined. However, as expected, concentrations of total IgG > IgG1 > IgG2. Because of this relationship, strong genetic correlations were expected between total IgG and isoforms, which was confirmed given the ≥0.83 ± 0.09 genetic correlation. Also, there was a strong genetic correlation between IgG1 and IgG2 (0.55 ± 0.16). Although total IgG and isoforms dif fer in phenotypic function, the strong genetic corre lation between total IgG and isoforms suggests that they are under similar genetic control. However, the genetic correlation between IgG1 and IgG2 suggested they differ genetically. Selection for total IgG will subsequently impact both IgG1 and IgG2. Selection for IgG1 or IgG2 will influence the other isoform, but not as strongly.

Limited reports were available describing the inter relationships among immunoglobulins, cortisol, and IL 8. Kimata et al. (1992) determined that IL8 had little influence on in vitro secretion of IgG1, IgG2, IgG3, or IgG4. Previous phenotypic relationships between IL8 and IgG and isoforms confirm weak genetic correlation estimates determined in the current study. Results sug gested that selection for IL8 will have little impact on total IgG or isoforms and vice versa. Interleukin 8 is a chemoattractant that impacts neutrophil production, in fluences innate immune response through angiogenesis (Murphy, 2008) , and is inhibited by cortisol (Dovio et al., 2004; Sorenson et al., 2012) . Therefore, a decrease in cortisol levels should increase IL8 production, po tentially leading to increased pain and inflammation. However, because only a low genetic correlation was observed, the genetic relationship between cortisol with IL8 appeared to be minimal at best in these animals. Overall, IgG, IgG1, IgG2, and cortisol appeared to be largely independent of IL8. Therefore, it is possible to develop a selection index for BRD susceptibility that includes cortisol, total IgG, and IL8 without confound ing genetic relationships among component traits.

Overall, reported low to moderate genetic correlations support our hypothesis that there is a genetic relationship between BRD and immune system traits. Specifically, heritability estimates and genetic relationships of cortisol and IgG with BRD suggest herd health may be improved through 1) genetic selection with breeding values and se lection indices and/or 2) improvement of the BRD phe notype through incorporation of immune traits in receiv ing management strategies. Admittedly, measurement of these parameters for use in management may be difficult given sample collection and time needed to deliver re sults from these diagnostic assays. This study described genetic relationships between treatment for BRD, stress hormones, and immune system traits. Several parameters had a favorable genetic relationship with BRD and may be useful to direct future research to reduce the incidence of BRD through genetic improvement. This is especially useful establishing a genetic selection model, because es tablishing a pedigree and relevant environmental factors to BRD is difficult; therefore, genetic improvement must be primarily driven through advancements in application and development of genomic prediction technologies.

Animal and Plant Health Inspection Service (APHIS 

Cortisol levels and anxi etyrelated behaviors in cattle

Bovine immunoglobulins: An augmented review

Influence of stress and nutri tion on cattle immunity

Expression of the neutrophil chemoattractant interleukin8 in the lesions of bovine pneumonic pasteurellosis

Immunological parameters of periparturient Holstein cattle: Genetic variation

Enhancing versus suppressive effects of stress on immune function: Implications for immunoprotection and immunopathology

Effects of physiological concentrations of steroid hormones and interleukin11 on basal and stimulated production of interleukin8 by human osteoblastlike cells with different functional profiles

Boardinvited review: Recent ad vances in management of highly stressed, newly received feedlot cattle

The immunology of the bovine respiratory disease complex

Type 1 and type 2 responses in regulation of Ig isotype expression in cattle

Growth and development symposium: Impacts of inflammation on cattle growth and carcass merit

Inbreeding and immunoglobulin G1 concentra tions in cattle

Placental transfer of IgG subclasses in a Japanese population

Unusual haplotypic structure of IL8, a susceptibility locus for a common respiratory virus

Update on pasteurellosis in young cattle

Population and sys tems genetics analyses of cortisol in pigs divergently selected for stress

Effect of shipping and chromium supplementation on performance, im mune response, and disease resistance of steers

Subclass of IgG in allergic disease. I. IgG subclass antibodies in immediate and nonimmediate food al lergy

Interleukin 8 (IL8) selectively inhibits immunoglobulin E production induced by IL4 in human B cells

Decreased neutrophil function as a cause of retained placenta in dairy cattle

Antibody responses to respiratory coronavirus infec tions of cattle during shipping fever pathogenesis

Infectivity neutralizing and hemaglutinininhibiting antibody responses to respiratory coronavirus infections of cattle in pathogenesis of shipping fever pneumonia

Variation in serum immunoglobulins in Canadian Holstein Friesians

Genetic parameters of bovine serum immuno globulins

Genetics of bovine respiratory disease in feedlot cattle

Shipping stress and social status effects on pig performance, plasma cortisol, natural killer cell activity, and leukocyte numbers

Molecular cloning and ex pression of bovine interleukin8

Genetic analysis of bovine respiratory disease in beef calves dur ing the first year of life

Chemokines

R: A language and environment for statistical computing. R Foundation for Statistical Computing

Making sense of ap parently conflicting data: Stress and immunity in swine and cattle

Evaluation of fixed sources of variation and estimation of genetic parameters for incidence of bovine respira tory disease in preweaned calves and feedlot cattle

Interactions of cold stress and Pasteurella haemolytica in the pathogenesis of pneumonic pas teurellosis in calves: Method of induction and hematologic and pathologic changes

Influence of breed, heterozygosity, and disease incidence on estimates of variance components of respiratory disease in preweaned beef calves

Bovine respiratory disease in feedlot cattle: Environmental, genetic, and economic factors

Bovine respiratory dis ease in feedlot cattle: Phenotypic, environmental, and genetic cor relations with growth, carcass, and longissimus muscle palatabil ity traits

The production of interleukin8 is increased in plasma and peripheral blood mononuclear cells of patients with fatigue

Serum IgG subclasses in chronic and recurrent respiratory infections

Genetic analysis of the serum level of IgG2 and total protein in Red Danish cattle

Factors associated with serum immunoglobulin levels in beef calves from Alberta and Saskatchewan and association between passive transfer and health outcomes