Transcriptomic analysis of the stress response to weaning at housing in bovine leukocytes using RNA-seq technology

Background: Weaning of beef calves is a necessary husbandry practice and involves separating the calf from its mother, resulting in numerous stressful events including dietary change, social reorganisation and the cessation of the maternal-offspring bond and is often accompanied by housing. While much recent research has focused on the physiological response of the bovine immune system to stress in recent years, little is known about the molecular mechanisms modulating the immune response. Therefore, the objective of this study was to provide new insights into the molecular mechanisms underlying the physiological response to weaning at housing in beef calves using Illumina RNA-seq. Results: The leukocyte transcriptome was significantly altered for at least 7 days following either housing or weaning at housing. Analysis of differentially expressed genes revealed that four main pathways, cytokine signalling, transmembrane transport, haemostasis and G-protein-coupled receptor (GPRC) signalling, were differentially regulated between control and weaned calves and underwent significant transcriptomic alterations in response to weaning stress on day 1, 2 and 7. Of particular note, chemokines, cytokines and integrins were consistently found to be up-regulated on each day following weaning. Evidence for alternative splicing of genes was also detected, indicating that a number of genes involved in the innate and adaptive immune response may be alternatively transcribed, including those responsible for toll receptor cascades and T cell receptor signalling. Conclusions: This study represents the first application of RNA-Seq technology for genomic studies in bovine leukocytes in response to weaning stress. Weaning stress induces the activation of a number of cytokine, chemokine and integrin transcripts and may alter the immune system whereby the ability of a number of cells of the innate and adaptive immune system to locate and destroy pathogens is transcriptionally enhanced. Stress alters the homeostasis of the transcriptomic environment of leukocytes for at least 7 days following weaning, indicating long-term effects of stress exposure in the bovine. The identification of gene signature networks that are stress activated provides a mechanistic framework to characterise the multifaceted nature of weaning stress adaptation in beef calves. Thus, capturing subtle transcriptomic changes provides insight into the molecular mechanisms that underlie the physiological response to weaning stress. Overall design: Examination of a time course (day 0, 1, 2 and 7) for 2 treatments, calves either housed with their dam (control) or housed and simultaneously weaned, using RNA-seq. The supplementary processed data file 'read_counts.txt' contains unnormalized read counts for each Ensembl bovine gene in each of the 48 samples. Unnormalized counts are required for input to EdgeR. Genome build: Btau4.0

48

O'Loughlin A, Lynn DJ, McGee M, Doyle S, McCabe M, Earley B