Hematological and gene co-expression network analyses of high-risk beef cattle defines immunological mechanisms and biological complexes involved in bovine respiratory disease and weight gain

Bovine respiratory disease (BRD), the leading disease complex in beef cattle production systems, remains highly elusive regarding diagnostics and disease prediction. Previous research has employed cellular and molecular techniques to describe hematological and gene expression variation that coincides with BRD development. Here, we utilized weighted gene co-expression network analysis (WGCNA) to leverage total gene expression patterns from cattle at arrival and generate hematological and clinical trait associations to describe mechanisms that may predict BRD development. Gene expression counts of previously published RNA-Seq data from 23 cattle (2017; n = 11 Healthy, n = 12 BRD) were used to construct gene co-expression modules and correlation patterns with complete blood count (CBC) and clinical datasets. Modules were further evaluated for cross-populational preservation of expression with RNA-Seq data from 24 cattle in an independent population (2019; n = 12 Healthy, n = 12 BRD). Genes within well-preserved modules were subject to functional enrichment analysis for significant Gene Ontology terms and pathways. Genes which possessed high module membership and association with BRD development, regardless of module preservation (“hub genes”), were utilized for protein-protein physical interaction network and clustering analyses. Five well-preserved modules of co-expressed genes were identified. One module (“steelblue”), involved in alpha-beta T-cell complexes and Th2-type immunity, possessed significant correlation with increased erythrocytes, platelets, and BRD development. One module (“purple”), involved in mitochondrial metabolism and rRNA maturation, possessed significant correlation with increased eosinophils, fecal egg count per gram, and weight gain over time. Fifty-two interacting hub genes, stratified into 11 clusters, may possess transient function involved in BRD development not previously described in literature. This study identifies co-expressed genes and coordinated mechanisms associated with BRD, which necessitates further investigation in BRD-prediction research.

Genomewide association study of lung lesions in cattle using sample pooling

Bovine respiratory disease complex (BRDC) is the most expensive disease in beef cattle in the United States costing the industry at least US$1 billion annually. Bovine respiratory disease complex causes damage to lung tissue resulting in persistent lung lesions observable at slaughter. Severe lung lesions at harvest have been associated with decreased preharvest ADG and increased clinical BRDC in the feedlot. Our objective was to identify SNP that are associated with severe lung lesions observed at harvest in feedlot cattle. We conducted a genomewide association study (GWAS) using a case-control design for severe lung lesions in fed cattle at slaughter using the Illumina Bovine HD array (approximately 770,000 SNP) and sample pooling. Lung samples were collected from 11,520 young cattle, a portion of which had not been treated with antibiotics (participating in a "natural" marketing program), at a large, commercial beef processing plant in central Nebraska. Lung samples with lesions (cases) and healthy lungs (controls) were collected when both phenotypes were in close proximity on the viscera (offal) table. We constructed 60 case and 60 control pools with 96 animals per pool. Pools were constructed by sampling sequence to ensure that case and control pool pairs were matched by proximity on the processing line. The Bovine HD array (770,000 SNP) was run on all pools. Fourteen SNP on BTA 2, 3, 4, 9, 11, 14, 15, 22, 24, and 25 were significant at the genomewide experiment-wise error rate of 5% ( ≤ 1.49 × 10). Eighty-five SNP on 28 chromosomes achieved a false discovery rate of 5% ( ≤ 5.38 × 10). Significant SNP were near (±100 kb) genes involved in tissue repair and regeneration, tumor suppression, cell proliferation, apoptosis, control of organ size, and immunity. Based on 85 significantly associated SNP in or near a collection of genes with diverse function on 28 chromosomes, we conclude that the genomic footprint of lung lesions is complex. A complex genomic footprint (genes and regulatory elements that affect the trait) is consistent with what is known about the cause of the disease: complex interactions among multiple viral and bacterial pathogens along with several environmental factors including dust, commingling, transportation, and stress. Characterization of sequence variation near significant SNP will enable accurate and cost effective genome-enhanced genetic evaluations for BRDC resistance in AI bulls and seed stock populations.

Susceptibility loci revealed for bovine respiratory disease complex in pre-weaned holstein calves

BACKGROUND

Bovine respiratory disease complex (BRDC) is an infectious disease of cattle that is caused by a combination of viral and/or bacterial pathogens. Selection for cattle with reduced susceptibility to respiratory disease would provide a permanent tool for reducing the prevalence of BRDC. The objective of this study was to identify BRDC susceptibility loci in pre-weaned Holstein calves as a prerequisite to using genetic improvement as a tool for decreasing the prevalence of BRDC. High density SNP genotyping with the Illumina BovineHD BeadChip was conducted on 1257 male and 757 female Holstein calves from California (CA), and 767 calves identified as female from New Mexico (NM). Of these, 1382 were classified as BRDC cases, and 1396 were classified as controls, with all phenotypes assigned using the McGuirk health scoring system. During the acquisition of blood for DNA isolation, two deep pharyngeal and one mid-nasal diagnostic swab were obtained from each calf for the identification of bacterial and viral pathogens. Genome-wide association analyses were conducted using four analytical approaches (EIGENSTRAT, EMMAX-GRM, GBLUP and FvR). The most strongly associated SNPs from each individual analysis were ranked and evaluated for concordance. The heritability of susceptibility to BRDC in pre-weaned Holstein calves was estimated.

RESULTS

The four statistical approaches produced highly concordant results for 373 top ranked SNPs that defined 126 chromosomal regions for the CA population. Similarly, in NM, 370 SNPs defined 138 genomic regions that were identified by all four approaches. When the two populations were combined (i.e., CA + NM) and analyzed, 324 SNPs defined 116 genomic regions that were associated with BRDC across all analytical methods. Heritability estimates for BRDC were 21% for both CA and NM as individual populations, but declined to 13% when the populations were combined.

CONCLUSIONS

Four analytical approaches utilizing both single and multi-marker association methods revealed common genomic regions associated with BRDC susceptibility that can be further characterized and used for genomic selection. Moderate heritability estimates were observed for BRDC susceptibility in pre-weaned Holstein calves, thereby supporting the application of genomic selection to reduce the prevalence of BRDC in U.S. Holsteins.

Differential stress responses among newly received calves: variations in reductant capacity and Hsp gene expression

Bovine respiratory disease complex (BRD), a major economic concern to the beef cattle industry all over the world, is triggered by physical, biological and psychological stresses. It is becoming noticeable that the key to reducing BRD appears to be centered at reducing the response to stress. The aims of the present study were to detect individual variations in the stress response of newly received young calves through their leukocyte heat shock protein (Hsp) response, selected neutrophil-related gene expression and oxidative stress, and relate them to pulmonary adhesions at slaughter, an indicative sign of clinical and subclinical episodes of BRD at an early age. Differential expression patterns of Hsp60 and Hsp70A1A were revealed in newly received calves 1 h, 5 h and 1 day after arrival, distinguishing between stress-responsive and non-stress-responsive individuals. Plasma cortisol was also indicative of stress-responsive and non-stress-responsive individuals, 1 h and 5 h after arrival. At the longer term, β-glycan levels were highest 7 days after arrival and significantly correlated with an adhesion-free phenotype at slaughter. Oxidative stress responses, measured through the oxidation products of the exogenous linoleoyl tyrosine (LT) marker, revealed that hydroperoxidation and epoxidation of membranes may readily occur. Based on the LT oxidation products and levels of β-glycan, we present a discriminant analysis model, according to which vulnerable individuals may be predicted at near 100% probability 7 days after arrival. Since clinical signs of BRD may often go undetected in feedlot calves, such a model, after its examination in large-scale experiments, may be a reliable tool for an early prediction of subclinical signs of BRD.

The epidemiology of bovine respiratory disease: What is the evidence for predisposing factors?

Bovine respiratory disease (BRD) is the most costly disease of beef cattle in North America. It is multi-factorial, with a variety of physical and physiological stressors combining to predispose cattle to pneumonia. However, efforts to discern which factors are most important have frequently failed to establish definitive answers. Calves are at highest risk shortly after transport. Risk factors include purchasing from sale barns and commingling. It is unclear whether or not these practices increase susceptibility, increase exposure, or are proxies for poor management. Lighter-weight calves appear to be at greater risk, although this has not been consistent. Persistent infection (PI) with bovine virus diarrhea virus increases BRD occurrence, but it is unclear if PI calves affect other cattle in the feedlot. The complexity of BRD has made it difficult to define involvement of individual factors. Stressors may play a role as "necessary but not sufficient" components, requiring additive effects to cause disease.

Respiratory function in cattle: impact of breed, heritability and external factors

Bovine respiratory disease complex (BRDC) is the most common and costly disease of cattle. The important prevalence of this multifactorial syndrome results from interactions between different pathogens, physical constitution of the host and environmental factors. Because of the strong economical impact of BRDC, it is recommended to complete classical preventive and therapeutic measures by selection of heritable traits improving resistance against respiratory disease and by alternatives based on improved control of environmental factors. The transport of oxygen between the atmospheric air and the mitochondria essentially includes four steps, i.e. (1) respiratory function, including pulmonary ventilation, perfusion and respiratory mechanics, (2) blood circulation and oxygen transport, (3) capillary-cell diffusion and (4) oxygen combustion by mitochondria. Each step should be considered when factors improving resistance against respiratory disease are investigated. This review aims (1) to summarise briefly the anatomical, histological and physiological peculiarities of the bovine respiratory system and (2) to consider the effect of breed, heritability and external factors at each step of the oxygen transport chain, by focussing essentially on respiratory function.

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

The objective of this study was to characterize genetic and environmental factors influencing bovine respiratory disease (BRD) in beef cattle. Records from nine purebred and three composite breeds and a variety of F1 and three-way crosses, including the progeny of 12 additional different sire breeds produced over a 20-yr period (1983 to 2002), were evaluated for breed and heterozygosity effects on the observed incidence of BRD. Heterozygosity fractions for calves and dams were defined by generalized breed origins: British, Continental, and tropically adapted. Variance components were estimated for each pure and composite breed, and across all breeds and crossbreeds. The effect of incidence of observed BRD was determined by comparing groups of low and high years of incidence. Respiratory disease in this herd followed a standard epidemiological pattern of initial introduction, reaching an epidemic stage at 70 to 170 d of age, followed by a period of rapid decrease to weaning. Estimates of heritability of incidence of BRD were low, ranging from 0.00 to 0.26, with overall estimates of 0.07 and 0.19 depending on the data set analyzed. The highest incidence of BRD in preweaned calves occurred in the Braunvieh breed (18.8%). The genetic correlation between the direct and maternal genetic effects was generally large and negative, suggesting dams genetically superior for resisting BRD raise calves that are more susceptible. Perhaps maternally superior dams provide passive immunity to their calves, which delays the development of the calves' direct immune system, making them more prone to BRD during the preweaning period. Heterozygosity of calves decreased the incidence of BRD compared with purebred cattle. Calves that were Continental x British or tropically adapted x British breeds had a lower incidence of BRD than did calves of British x British breeds. As the annual incidence of BRD increased, there was an associated increase in the heritability estimate. The estimated heritability based on an underlying continuous scale was large (h2 = 0.48), inferring response to selection for BRD resistance could be large if the phenotype for BRD resistance was known.

Feedlot entry characteristics and climate: their relationship with cattle growth rate, bovine respiratory disease and mortality

OBJECTIVE

To examine the relationship between the feedlot entry characteristics of breed, gender, entry weight and the presence or absence of permanent incisors, with growth rate, the incidence of bovine respiratory disease (BRD), and mortalities of cattle destined for the Australian domestic market. The relationship between climate data and the occurrence of BRD was also investigated.

DESIGN

Two thousand four hundred and sixty eight cattle were individually identified at feedlot entry, and lot number, body weight, property identification code, breed, gender, and the presence or absence of permanent incisors were recorded for each animal. Daily climate data were collected for the duration of the experiment. Univariate correlations, multivariate analysis and survival analysis techniques were used to examine the relationship between these feedlot entry characteristics and the outcomes of BRD and growth rate.

RESULTS

The feedlot entry description of breed had a significant effect on average daily gain. The presence of permanent incisors at feedlot entry was associated with an increase in average daily gain of 0.1 kg/animal/d. Growth rate decreased by 0.702 kg/animal/d in cattle treated for BRD and by 0.302 kg/animal/d in cattle identified with anorexia. Treatment for BRD was ten times more likely with Hereford cattle. There was a moderate correlation between minimum daily temperature and the occurrence of BRD.

CONCLUSIONS

This study indicates that breed can affect growth rate and health, specifically the occurrence of BRD, in Australian domestic feedlot cattle. The moderate correlation between minimum daily temperature and the occurrence of BRD indicates potential associations between climate and BRD incidence in Australian feedlots.

Bovine respiratory disease in feedlot cattle: phenotypic, environmental, and genetic correlations with growth, carcass, and longissimus muscle palatability traits

Bovine respiratory disease (BRD) is the most costly feedlot disease in the United States. Selection for disease resistance is one of several possible interventions to prevent or reduce the economic loss associated with animal disease and to improve animal welfare. Undesirable genetic relationships, however, may exist between production and disease resistance traits. The objectives of this study were to estimate the phenotypic, environmental, and genetic correlations of BRD with growth, carcass, and LM palatability traits. Health records on 18,112 feedlot cattle over a 15-yr period and slaughter data on 1,627 steers over a 4-yr period were analyzed with bivariate animal models. Traits included ADG, adjusted carcass fat thickness at the 12th rib, marbling score, LM area, weight of retail cuts, weight of fat trim, bone weight, Warner-Bratzler shear force, tenderness score, and juiciness score. The estimated heritability of BRD incidence was 0.08 +/- 0.01. Phenotypic, environmental, and genetic correlations of the observed traits with BRD ranged from -0.35 to 0.40, -0.36 to 0.55, and -0.42 to 0.20, respectively. Most correlations were low or negligible. The percentage of carcass bone had moderate genetic, phenotypic, and environmental correlations with BRD (-0.42, -0.35, and -0.36, respectively). Hot carcass weight and weight of retail cuts had moderate, undesirable phenotypic correlations with BRD (0.37 and 0.40, respectively). Correlations of BRD with LM palatability and ADG were not detected. Low or near zero estimates of genetic correlations infer that selection to reduce BRD in feedlot cattle would have negligible correlated responses on growth, carcass, and meat palatability traits or that selection for those traits will have little effect on BRD susceptibility or resistance.