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.

Use of nCounter mRNA profiling to identify at-arrival gene expression patterns for predicting bovine respiratory disease in beef cattle.

BACKGROUND: Transcriptomics has identified at-arrival differentially expressed genes associated with bovine respiratory disease (BRD) development; however, their use as prediction molecules necessitates further evaluation. Therefore, we aimed to selectively analyze and corroborate at-arrival mRNA expression from multiple independent populations of beef cattle. In a nested case-control study, we evaluated the expression of 56 mRNA molecules from at-arrival blood samples of 234 cattle across seven populations via NanoString nCounter gene expression profiling. Analysis of mRNA was performed with nSolver Advanced Analysis software (p < 0.05), comparing cattle groups based on the diagnosis of clinical BRD within 28 days of facility arrival (n = 115 Healthy; n = 119 BRD); BRD was further stratified for severity based on frequency of treatment and/or mortality (Treated_1, n = 89; Treated_2+, n = 30). Gene expression homogeneity of variance, receiver operator characteristic (ROC) curve, and decision tree analyses were performed between severity cohorts. RESULTS: Increased expression of mRNAs involved in specialized pro-resolving mediator synthesis (ALOX15, HPGD), leukocyte differentiation (LOC100297044, GCSAML, KLF17), and antimicrobial peptide production (CATHL3, GZMB, LTF) were identified in Healthy cattle. BRD cattle possessed increased expression of CFB, and mRNA related to granulocytic processes (DSG1, LRG1, MCF2L) and type-I interferon activity (HERC6, IFI6, ISG15, MX1). Healthy and Treated_1 cattle were similar in terms of gene expression, while Treated_2+ cattle were the most distinct. ROC cutoffs were used to generate an at-arrival treatment decision tree, which classified 90% of Treated_2+ individuals. CONCLUSIONS: Increased expression of complement factor B, pro-inflammatory, and type I interferon-associated mRNA hallmark the at-arrival expression patterns of cattle that develop severe clinical BRD. Here, we corroborate at-arrival mRNA markers identified in previous transcriptome studies and generate a prediction model to be evaluated in future studies. Further research is necessary to evaluate these expression patterns in a prospective manner.

Multipopulational transcriptome analysis of post-weaned beef cattle at arrival further validates candidate biomarkers for predicting clinical bovine respiratory disease.

Bovine respiratory disease (BRD) remains the leading infectious disease in post-weaned beef cattle. The objective of this investigation was to contrast the at-arrival blood transcriptomes from cattle derived from two distinct populations that developed BRD in the 28 days following arrival versus cattle that did not. Forty-eight blood samples from two populations were selected for mRNA sequencing based on even distribution of development (n = 24) or lack of (n = 24) clinical BRD within 28 days following arrival; cattle which developed BRD were further stratified into BRD severity cohorts based on frequency of antimicrobial treatment: treated once (treated_1) or treated twice or more and/or died (treated_2+). Sequenced reads (~ 50 M/sample, 150 bp paired-end) were aligned to the ARS-UCD1.2 bovine genome assembly. One hundred and thirty-two unique differentially expressed genes (DEGs) were identified between groups stratified by disease severity (healthy, n = 24; treated_1, n = 13; treated_2+, n = 11) with edgeR (FDR ≤ 0.05). Differentially expressed genes in treated_1 relative to both healthy and treated_2+ were predicted to increase neutrophil activation, cellular cornification/keratinization, and antimicrobial peptide production. Differentially expressed genes in treated_2+ relative to both healthy and treated_1 were predicted to increase alternative complement activation, decrease leukocyte activity, and increase nitric oxide production. Receiver operating characteristic (ROC) curves generated from expression data for six DEGs identified in our current and previous studies (MARCO, CFB, MCF2L, ALOX15, LOC100335828 (aka CD200R1), and SLC18A2) demonstrated good-to-excellent (AUC: 0.800-0.899; ≥ 0.900) predictability for classifying disease occurrence and severity. This investigation identifies candidate biomarkers and functional mechanisms in at arrival blood that predicted development and severity of BRD.

Genes and regulatory mechanisms associated with experimentally-induced bovine respiratory disease identified using supervised machine learning methodology.

Bovine respiratory disease (BRD) is a multifactorial disease involving complex host immune interactions shaped by pathogenic agents and environmental factors. Advancements in RNA sequencing and associated analytical methods are improving our understanding of host response related to BRD pathophysiology. Supervised machine learning (ML) approaches present one such method for analyzing new and previously published transcriptome data to identify novel disease-associated genes and mechanisms. Our objective was to apply ML models to lung and immunological tissue datasets acquired from previous clinical BRD experiments to identify genes that classify disease with high accuracy. Raw mRNA sequencing reads from 151 bovine datasets (n = 123 BRD, n = 28 control) were downloaded from NCBI-GEO. Quality filtered reads were assembled in a HISAT2/Stringtie2 pipeline. Raw gene counts for ML analysis were normalized, transformed, and analyzed with MLSeq, utilizing six ML models. Cross-validation parameters (fivefold, repeated 10 times) were applied to 70% of the compiled datasets for ML model training and parameter tuning; optimized ML models were tested with the remaining 30%. Downstream analysis of significant genes identified by the top ML models, based on classification accuracy for each etiological association, was performed within WebGestalt and Reactome (FDR ≤ 0.05). Nearest shrunken centroid and Poisson linear discriminant analysis with power transformation models identified 154 and 195 significant genes for IBR and BRSV, respectively; from these genes, the two ML models discriminated IBR and BRSV with 100% accuracy compared to sham controls. Significant genes classified by the top ML models in IBR (154) and BRSV (195), but not BVDV (74), were related to type I interferon production and IL-8 secretion, specifically in lymphoid tissue and not homogenized lung tissue. Genes identified in Mannheimia haemolytica infections (97) were involved in activating classical and alternative pathways of complement. Novel findings, including expression of genes related to reduced mitochondrial oxygenation and ATP synthesis in consolidated lung tissue, were discovered. Genes identified in each analysis represent distinct genomic events relevant to understanding and predicting clinical BRD. Our analysis demonstrates the utility of ML with published datasets for discovering functional information to support the prediction and understanding of clinical BRD.

Comprehensive at-arrival transcriptomic analysis of post-weaned beef cattle uncovers type I interferon and antiviral mechanisms associated with bovine respiratory disease mortality.

BACKGROUND: Despite decades of extensive research, bovine respiratory disease (BRD) remains the most devastating disease in beef cattle production. Establishing a clinical diagnosis often relies upon visual detection of non-specific signs, leading to low diagnostic accuracy. Thus, post-weaned beef cattle are often metaphylactically administered antimicrobials at facility arrival, which poses concerns regarding antimicrobial stewardship and resistance. Additionally, there is a lack of high-quality research that addresses the gene-by-environment interactions that underlie why some cattle that develop BRD die while others survive. Therefore, it is necessary to decipher the underlying host genomic factors associated with BRD mortality versus survival to help determine BRD risk and severity. Using transcriptomic analysis of at-arrival whole blood samples from cattle that died of BRD, as compared to those that developed signs of BRD but lived (n = 3 DEAD, n = 3 ALIVE), we identified differentially expressed genes (DEGs) and associated pathways in cattle that died of BRD. Additionally, we evaluated unmapped reads, which are often overlooked within transcriptomic experiments. RESULTS: 69 DEGs (FDR<0.10) were identified between ALIVE and DEAD cohorts. Several DEGs possess immunological and proinflammatory function and associations with TLR4 and IL6. Biological processes, pathways, and disease phenotype associations related to type-I interferon production and antiviral defense were enriched in DEAD cattle at arrival. Unmapped reads aligned primarily to various ungulate assemblies, but failed to align to viral assemblies. CONCLUSION: This study further revealed increased proinflammatory immunological mechanisms in cattle that develop BRD. DEGs upregulated in DEAD cattle were predominantly involved in innate immune pathways typically associated with antiviral defense, although no viral genes were identified within unmapped reads. Our findings provide genomic targets for further analysis in cattle at highest risk of BRD, suggesting that mechanisms related to type I interferons and antiviral defense may be indicative of viral respiratory disease at arrival and contribute to eventual BRD mortality.

Whole blood transcriptomic analysis of beef cattle at arrival identifies potential predictive molecules and mechanisms that indicate animals that naturally resist bovine respiratory disease.

Bovine respiratory disease (BRD) is a multifactorial disease complex and the leading infectious disease in post-weaned beef cattle. Clinical manifestations of BRD are recognized in beef calves within a high-risk setting, commonly associated with weaning, shipping, and novel feeding and housing environments. However, the understanding of complex host immune interactions and genomic mechanisms involved in BRD susceptibility remain elusive. Utilizing high-throughput RNA-sequencing, we contrasted the at-arrival blood transcriptomes of 6 beef cattle that ultimately developed BRD against 5 beef cattle that remained healthy within the same herd, differentiating BRD diagnosis from production metadata and treatment records. We identified 135 differentially expressed genes (DEGs) using the differential gene expression tools edgeR and DESeq2. Thirty-six of the DEGs shared between these two analysis platforms were prioritized for investigation of their relevance to infectious disease resistance using WebGestalt, STRING, and Reactome. Biological processes related to inflammatory response, immunological defense, lipoxin metabolism, and macrophage function were identified. Production of specialized pro-resolvin mediators (SPMs) and endogenous metabolism of angiotensinogen were increased in animals that resisted BRD. Protein-protein interaction modeling of gene products with significantly higher expression in cattle that naturally acquire BRD identified molecular processes involving microbial killing. Accordingly, identification of DEGs in whole blood at arrival revealed a clear distinction between calves that went on to develop BRD and those that resisted BRD. These results provide novel insight into host immune factors that are present at the time of arrival that confer protection from BRD.

Modeling the pasture-associated severe equine asthma bronchoalveolar lavage fluid proteome identifies molecular events mediating neutrophilic airway inflammation.

Background: Pasture-associated severe equine asthma is a warm season, environmentally-induced respiratory disease characterized by reversible airway obstruction, persistent and non-specific airway hyper-responsiveness, and chronic neutrophilic airway inflammation. During seasonal exacerbation, signs vary from mild to life-threatening episodes of wheezing, coughing, and chronic debilitating labored breathing. Purpose: In human asthma, neutrophilic airway inflammation is associated with more severe and steroid-refractory asthma phenotypes, highlighting a need to decipher the mechanistic basis of this disease characteristic. We hypothesize that the collective biological activities of proteins in bronchoalveolar lavage fluid (BALF) of horses with pasture-associated severe asthma predict changes in neutrophil functions that contribute to airway neutrophilic inflammation. Methods: Using shotgun proteomics, we identified 1,003 unique proteins in cell-free BALF from six horses experiencing asthma exacerbation and six control herdmates. Contributions of each protein to ten neutrophil functions were modeled using manual biocuration to determine each protein's net effect on the respective neutrophil functions. Results: A total of 417 proteins were unique to asthmatic horses, 472 proteins were unique to control horses (p<0.05), and 114 proteins were common in both groups. Proteins whose biological activities are responsible for increasing neutrophil migration, chemotaxis, cell spreading, transmigration, and infiltration, which would collectively bring neutrophils to airways, were over-represented in the BALF of asthmatic relative to control horses. By contrast, proteins whose biological activities support neutrophil activation, adhesion, phagocytosis, respiratory burst, and apoptosis, which would collectively shorten neutrophil lifespan, were under-represented in BALF of asthmatic relative to control horses. Interaction networks generated using Ingenuity® Pathways Analysis further support the results of our biocuration. Conclusion: Congruent with our hypothesis, the collective biological functions represented in differentially expressed proteins of BALF from horses with pasture-associated severe asthma support neutrophilic airway inflammation. This illustrates the utility of systems modeling to organize functional genomics data in a manner that characterizes complex molecular events associated with clinically relevant disease.

Effect of a syringe aspiration technique versus a mechanical suction technique and use of N-butylscopolammonium bromide on the quantity and quality of bronchoalveolar lavage fluid samples obtained from horses with the summer pasture endophenotype of equine asthma.

OBJECTIVE To evaluate the effect of 2 bronchoalveolar lavage (BAL) sampling techniques and the use of N-butylscopolammonium bromide (NBB) on the quantity and quality of BAL fluid (BALF) samples obtained from horses with the summer pasture endophenotype of equine asthma. ANIMALS 8 horses with the summer pasture endophenotype of equine asthma. PROCEDURES BAL was performed bilaterally (right and left lung sites) with a flexible videoendoscope passed through the left or right nasal passage. During lavage of the first lung site, a BALF sample was collected by means of either gentle syringe aspiration or mechanical suction with a pressure-regulated wall-mounted suction pump. The endoscope was then maneuvered into the contralateral lung site, and lavage was performed with the alternate fluid retrieval technique. For each horse, BAL was performed bilaterally once with and once without premedication with NBB (21-day interval). The BALF samples retrieved were evaluated for volume, total cell count, differential cell count, RBC count, and total protein concentration. RESULTS Use of syringe aspiration significantly increased total BALF volume (mean volume increase, 40 mL [approx 7.5% yield]) and decreased total RBC count (mean decrease, 142 cells/µL), compared with use of mechanical suction. The BALF nucleated cell count and differential cell count did not differ between BAL procedures. Use of NBB had no effect on BALF retrieval. CONCLUSIONS AND CLINICAL RELEVANCE Results indicated that retrieval of BALF by syringe aspiration may increase yield and reduce barotrauma in horses at increased risk of bronchoconstriction and bronchiolar collapse. Further studies to determine the usefulness of NBB and other bronchodilators during BAL procedures in horses are warranted.

Horses With Pasture Asthma Have Airway Remodeling That Is Characteristic of Human Asthma.

Severe equine asthma, formerly recurrent airway obstruction (RAO), is the horse counterpart of human asthma, affecting horses maintained indoors in continental climates. Equine pasture asthma, formerly summer pasture RAO, is clinically similar but affects grazing horses during hot, humid conditions in the southeastern United States and United Kingdom. To advance translational relevance of equine pasture asthma to human asthma, histologic features of airway remodeling in human asthma were scored in lung lobes from 15 pasture asthma-affected and 9 control horses of mixed breeds. All noncartilaginous airways were scored using a standardized grading rubric (0-3) in hematoxylin and eosin (HE) and Movat's pentachrome-stained sections; 15 airways were chosen randomly from each lobe for analysis. Logistic regression identified disease, age, and lobe effects on probability of histologic outcomes. Airway smooth muscle (odds ratio [OR] = 2.5, P < .001), goblet cell hyperplasia/metaplasia (OR = 37.6, P < .0001), peribronchiolar elastic system fibers (OR = 4.2, P < .001), peribronchiolar fibrosis (OR = 3.8, P = .01), airway occlusion by mucus/inflammation (OR = 4.2, P = .04), and airway adventitial inflammation (OR = 3.0, P = .01) were significantly greater in diseased airways. A novel complex tissue disorganization, designated terminal bronchiolar remodeling, was overrepresented in diseased airways (OR = 3.7, P < .0001). Distribution of terminal bronchiolar remodeling corresponded to putative sites of air trapping in human asthma, at secondary pulmonary lobules. Age (>15 years) was an independent risk factor for increased peribronchiolar fibrosis, elastic system fibers, and terminal bronchiolar remodeling. Remodeling differed significantly between lung lobes, congruent with nonhomogeneous remodeling in human asthma. Equine pasture asthma recapitulates airway remodeling in human asthma in a manner not achieved in induced animal asthma models, endorsing its translational relevance for human asthma investigation.

Structural and functional-annotation of an equine whole genome oligoarray.

BACKGROUND: The horse genome is sequenced, allowing equine researchers to use high-throughput functional genomics platforms such as microarrays; next-generation sequencing for gene expression and proteomics. However, for researchers to derive value from these functional genomics datasets, they must be able to model this data in biologically relevant ways; to do so requires that the equine genome be more fully annotated. There are two interrelated types of genomic annotation: structural and functional. Structural annotation is delineating and demarcating the genomic elements (such as genes, promoters, and regulatory elements). Functional annotation is assigning function to structural elements. The Gene Ontology (GO) is the de facto standard for functional annotation, and is routinely used as a basis for modelling and hypothesis testing, large functional genomics datasets. RESULTS: An Equine Whole Genome Oligonucleotide (EWGO) array with 21,351 elements was developed at Texas A&M University. This 70-mer oligoarray was designed using the approximately 7 x assembled and annotated sequence of the equine genome to be one of the most comprehensive arrays available for expressed equine sequences. To assist researchers in determining the biological meaning of data derived from this array, we have structurally annotated it by mapping the elements to multiple database accessions, including UniProtKB, Entrez Gene, NRPD (Non-Redundant Protein Database) and UniGene. We next provided GO functional annotations for the gene transcripts represented on this array. Overall, we GO annotated 14,531 gene products (68.1% of the gene products represented on the EWGO array) with 57,912 annotations. GAQ (GO Annotation Quality) scores were calculated for this array both before and after we added GO annotation. The additional annotations improved the meanGAQ score 16-fold. This data is publicly available at AgBase http://www.agbase.msstate.edu/. CONCLUSION: Providing additional information about the public databases which link to the gene products represented on the array allows users more flexibility when using gene expression modelling and hypothesis-testing computational tools. Moreover, since different databases provide different types of information, users have access to multiple data sources. In addition, our GO annotation underpins functional modelling for most gene expression analysis tools and enables equine researchers to model large lists of differentially expressed transcripts in biologically relevant ways.

Biomechanical and molecular characteristics of hereditary equine regional dermal asthenia in Quarter Horses.

Hereditary equine regional dermal asthenia (HERDA) is an autosomal recessive skin disorder that has yet to be fully characterized. HERDA is predominately expressed in Quarter Horses, with the majority of these disseminating from elite cutting horse bloodlines, leading to the increased incidence of HERDA in recent years. Affected horses have loose, hyper-extensible, fragile skin and are frequently euthanized due to poor wound healing and disfiguring scars. This study sought to better characterize HERDA by analysis of the biomechanical parameters of tensile strength, modulus of elasticity, energy to failure and thickness of skin from 10 affected and 6 unaffected horses using an Instron Universal Testing Instrument. In addition, total soluble collagen and glycosaminoglycan concentrations of skin were analysed from 13 affected and 12 unaffected horses using Sircol Soluble Collagen and Blyscan Sulfated Glycosaminoglycan assays respectively. Affected horses exhibited a two to threefold reduction in tensile strength versus unaffected horses with statistically significant differences at six of seven sample locations (P < or = 0.05). The modulus of elasticity proved to be significantly different at six of seven sample locations, energy to failure at six of seven sample locations, and skin thickness at one of seven sample locations (P < or = 0.05). Affected horses exhibited significantly higher amounts of total soluble collagen than unaffected horses (P < or = 0.05). No significant difference was demonstrated between groups for glycosaminoglycan concentration. Affected horses demonstrated uniformly weaker skin across sample locations, indicating the biomechanical properties of HERDA are not regionally confined to specific areas of the horses' skin.