Differential Expression of Circadian Clock Genes in the Bovine Neuroendocrine Adrenal System

Abstract

Knowledge of circadian rhythm clock gene expression outside the suprachiasmatic nucleus is increasing. The purpose of this study was to determine whether expression of circadian clock genes differed within or among the bovine stress axis tissues (e.g., amygdala, hypothalamus, pituitary, adrenal cortex, and adrenal medulla). Tissues were obtained at an abattoir from eight mature nonpregnant Brahman cows that had been maintained in the same pasture and nutritional conditions. Sample tissues were stored in RNase-free sterile cryovials at −80 °C until the total RNA was extracted, quantified, assessed, and sequenced (NovaSeq 6000 system; paired-end 150 bp cycles). The trimmed reads were then mapped to a Bos taurus (B. taurus) reference genome (Umd3.1). Further analysis used the edgeR package. Raw gene count tables were read into RStudio, and low-expression genes were filtered out using the criteria of three minimum reads per gene in at least five samples. Normalization factors were then calculated using the trimmed mean of M values method to produce normalized gene counts within each sample tissue. The normalized gene counts important for a circadian rhythm were analyzed within and between each tissue of the stress axis using the GLM and CORR procedures of the Statistical Analysis System (SAS). The relative expression profiles of circadian clock genes differed (p < 0.01) within each tissue, with neuronal PAS domain protein 2 (NPAS2) having greater expression in the amygdala (p < 0.01) and period circadian regulator (PER1) having greater expression in all other tissues (p < 0.01). The expression among tissues also differed (p < 0.01) for individual circadian clock genes, with circadian locomotor output cycles protein kaput (CLOCK) expression being greater within the adrenal tissues and nuclear receptor subfamily 1 group D member 1 (NR1D1) expression being greater within the other tissues (p < 0.01). Overall, the results indicate that within each tissue, the various circadian clock genes were differentially expressed, in addition to being differentially expressed among the stress tissues of mature Brahman cows. Future use of these findings may assist in improving livestock husbandry and welfare by understanding interactions of the environment, stress responsiveness, and peripheral circadian rhythms.