Retrospective analysis of cohort risk factors and feeding phase timing associated with noninfectious heart disease deaths in U.S. feedlot cattle

Heart disease, specifically congestive heart failure, has become of increased interest to geneticists and cattle feeders. Data on cohort associations of risk factors related to heart disease and when heart disease deaths occur in U.S. feedlot cattle are limited. The study objectives were to 1) determine potential associations between feedlot cohort demographics and the risk of at least one noninfectious heart disease (NIHD) death occurrence and 2) determine potential association between feedlot cohort demographics and the timing of NIHD deaths during the feeding phase. Data were downloaded from commercial feedyard software and analyzed by constructing a generalized linear mixed model for both analyses. A binomial and Gaussian distributions for risk of NIHD death and timing of NIHD were utilized as link functions for their respective models. Our study population consisted of 28,950 cohorts (representing 4,596,205 cattle) that were placed in 22 U.S. commercial feedlots from January 01, 2016, to January 01, 2019. There were 3,282 cases of NIHD deaths from a population of 75,963 cattle that died during the 3-yr study period. Average cohort arrival weight’s effect on NIHD probability was influenced by arrival quarter and arrival year of placement (P < 0.01). Cohorts with steers were associated with a greater probability of at least one NIHD death (2.38%) compared with heifers (1.95%; P < 0.01). Increasing cohort size was associated with an increased probability of a cohort having at least one NIHD death (P < 0.01). The probability of a cohort having at least one NIHD death increased with increasing DOF categories from 1.51% in cattle fed 100 to 175 d, to 2.12% in cattle fed 176 to 250 d, and 2.87% for cattle fed 251 to 326 d. Cattle > 326 DOF were no different in the probability of a NIHD death compared with the other feeding categories. Timing of a NIHD death had a mean and median occurrence of 110 DOF with an interquartile range of 64 to 153 DOF. The effect of arrival weight on days at death was influenced by year placed with heavier cattle generally decreasing the model adjusted means of DOF at NIHD death. Arrival quarter was influenced by year placed on model adjusted means on the timing of a NIHD death. Steers with NIHD died later compared with heifers (P < 0.01) diagnosed with NIHD. In conclusion, multiple factors are associated with the probability and timing of a NIHD death. Probability of having at least one NIHD death within a cohort was low, and half of the deaths occurred before 110 DOF.

Brisket Disease Is Associated with Lower Volatile Fatty Acid Production and Altered Rumen Microbiome in Holstein Heifers

Simple Summary

Development of the dairy industry in the high-altitude plateau environment through incorporation of Holstein cows is complicated by the risk of brisket disease. While the physiological effects of brisket disease are well-studied, its effects on rumen function and microbial community composition are not. There are clear shifts in volatile fatty acids production and rumen microbial community composition in Holstein heifers suffering from brisket disease. Observed shifts reveal key genera associated with healthy and disease states and suggest that bovine brisket disease is associated with impaired rumen functioning. This work supports further understanding of the roles of key rumen taxa in bovine brisket disease, with particular focus on candidate rumen biomarkers in healthy animals that may be able to reduce economic losses for farmers.

Abstract

Brisket disease is heritable but is also associated with non-genetic risk factors and effects of the disease on the rumen microbiome are unknown. Ten Holstein heifers were exposed to the plateau environment for three months and divided into two groups according to the index of brisket disease, the mean pulmonary arterial pressure (mPAP): brisket disease group (BD, n = 5, mPAP > 63 mmHg) and healthy heifer group (HH, n = 5, mPAP < 41 mmHg). Rumen fluid was collected for analysis of the concentrations of volatile fatty acids (VFAs). Extracted DNA from rumen contents was analyzed using Illumina MiSeq 16S rRNA sequencing technology. The concentration of total VFA and alpha-diversity metrics were significantly lower in BD group (p < 0.05). Ruminococcus and Treponema were significantly decreased in BD heifers (p < 0.05). Correlation analysis indicated that 10 genera were related to the mPAP (p < 0.05). Genera of Anaerofustis, Campylobacter, and Catonella were negatively correlated with total VFA and acetic acid (R < −0.7, p < 0.05), while genera of Blautia, YRC22, Ruminococcus, and Treponema were positively related to total VFA and acetic acid (R > 0.7; p < 0.05). Our findings may be a useful biomarker in future brisket disease work.

Evaluation of EPAS1 variants for association with bovine congestive heart failure

Background:  Bovine congestive heart failure (BCHF) has become increasingly prevalent in feedlot cattle in the Western Great Plains of North America. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. A protein variant of hypoxia-inducible factor 2 alpha (HIF2α, encoded by the endothelial PAS domain-containing protein 1 gene, EPAS1) was previously reported to be associated with pulmonary hypertension at altitudes exceeding 2,000 m. Our aim was to evaluate EPAS1 haplotypes for association with BCHF in feedlot cattle raised at moderate altitudes (1,200 m).

Methods: Paired samples of clinical cases and unaffected controls were collected at four feedlots in Nebraska and Wyoming. Each pair (n =102) was matched for source, pen, breed type, sex, arrival date, and management conditions. Cases were identified by animal caretakers, euthanized, and diagnosis was confirmed at necropsy. Cases were derived from 30 different ranch operations, with the largest source contributing 32. Animals were tested for eight EPAS1 haplotypes encoding 36 possible different diploid combinations.

Results: The common, ancestral EPAS1 haplotype encoding HIF2α with alanine (A) at position 606 and glycine (G) at position 610 was equally frequent in cases and controls (0.67). The EPAS1 variant haplotype reported to be associated with disease (encoding threonine (T) at position 606 and serine (S) at position 610) was not enriched in cases compared with controls (0.21 and 0.25, respectively). Frequencies of other EPAS1 haplotypes (e.g., encoding Q270, L362, or G671) were each less than 0.05 overall. McNemar’s test with 45 discordant pairs showed the linked T606/S610 variant was not associated with BCHF (OR = 0.73, CI ₉₅ 0.38 -1.4, p-value = 0.37).

Conclusions: HIF2α polypeptide variants were not significantly associated with BCHF in feedlot cattle at moderate altitudes. Thus, a wider search is needed to identify genetic risk factors underlying this disease.

High-altitude pulmonary hypertension in cattle (brisket disease): Candidate genes and gene expression profiling of peripheral blood mononuclear cells

High-altitude pulmonary hypertension (HAPH) is a consequence of chronic alveolar hypoxia, leading to hypoxic vasoconstriction and remodeling of the pulmonary circulation. Brisket disease in cattle is a naturally occurring animal model of hypoxic pulmonary hypertension. Genetically susceptible cattle develop severe pulmonary hypertension and right heart failure at altitudes >7,000 ft. No information currently exists regarding the identity of the pathways and gene(s) responsible for HAPH or influencing severity. We hypothesized that initial insights into the pathogenesis of the disease could be discovered by a strategy of (1) sequencing of functional candidates revealed by single nucleotide polymorphism (SNP) analysis and (2) gene expression profiling of affected cattle compared with altitude-matched normal controls, with gene set enrichment analysis (GSEA) and Ingenuity pathway analysis (IPA). We isolated blood from a single herd of Black Angus cattle of both genders, aged 12-18 months, by jugular vein puncture. Mean pulmonary arterial pressures were 85.6±13 mmHg STD in the 10 affected and 35.3±1.2 mmHg STD in the 10 resistant cattle, P<0.001. From peripheral blood mononuclear cells, DNA was hybridized to an Affymetrix 10K Gene Chip SNP, and RNA was used to probe an Affymetrix Bovine genome array. SNP loci were remapped using the Btau 4.0 bovine genome assembly. mRNA data was analyzed by the Partek software package to identify sets of genes with an expression that was statistically different between the two groups. GSEA and IPA were conducted on the refined expression data to identify key cellular pathways and to generate networks and conduct functional analyses of the pathways and networks. Ten SNPs were identified by allelelic association and four candidate genes were sequenced in the cohort. Neither endothelial nitric oxide synthetase, NADH dehydrogenase, TG-interacting factor-2 nor BMPR2 were different among affected and resistant cattle. A 60-gene mRNA signature was identified that differentiated affected from unaffected cattle. Forty-six genes were overexpressed in the affected and 14 genes were downregulated in the affected cattle by at least 20%. GSEA and Ingenuity analysis identified respiratory diseases, inflammatory diseases and pathways as the top diseases and disorders (P<5.14×10(-14)), cell development and cell signaling as the top cellular functions (P<1.20×10(-08)), and IL6, TREM, PPAR, NFkB cell signaling (P<8.69×10(-09)) as the top canonical pathways associated with this gene signature. This study provides insights into differences in RNA expression in HAPH at a molecular level, and eliminates four functional gene candidates. Further studies are needed to validate and refine these preliminary findings and to determine the role of transcribed genes in the development of HAPH.