Longitudinal blood transcriptomic analysis to identify molecular regulatory patterns of bovine respiratory disease in beef cattle.

Bovine respiratory disease (BRD) is the most common disease in beef cattle and leads to considerable economic losses in both beef and dairy cattle. It is important to uncover the molecular mechanisms underlying BRD and to identify biomarkers for early identification of BRD cattle in order to address its impact on production and welfare. In this study, a longitudinal transcriptomic analysis was conducted using blood samples collected from 24 beef cattle at three production stages in the feedlot: 1) arrival (Entry group); 2) when identified as sick (diagnosed as BRD) and separated for treatment (Pulled); 3) prior to marketing (Close-out, representing healthy animals). Expressed genes were significantly different in the same animal among Entry, Pulled and Close-out stages (false discovery rate (FDR) < 0.01 & |Fold Change| > 2). Beef steers at both Entry and Pulled stages presented obvious difference in GO terms (FDR < 0.05) and affected biological functions (FDR < 0.05 & |Z-score| > 2) when compared with animals at Close-out. However, no significant functional difference was observed between Entry and Pulled animals. The interferon signaling pathway showed the most significant difference between animals at Entry/Pulled and Close-out stages (P < .001 & |Z-score| > 2), suggesting the animals initiated antiviral responses at an early stage of infection. Six key genes including IFI6, IFIT3, ISG15, MX1, and OAS2 were identified as biomarkers to predict and recognize sick cattle at Entry. A gene module with 169 co-expressed genes obtained from WGCNA analysis was most positively correlated (R = 0.59, P = 6E-08) with sickness, which was regulated by 11 transcription factors. Our findings provide an initial understanding of the BRD infection process in the field and suggests a subset of novel marker genes for identifying BRD in cattle at an early stage of infection.