With the advance of high-throughput sequencing technology numerous new regulatory small RNAs have been identified, that broaden the variety of processing mechanisms and functions of non-coding RNA. Here we explore small non-coding RNA (sncRNA) expression in central parts of the physiological stress and anxiety response system. Therefore, we characterize the sncRNA profile of tissue samples from Amygdala, Hippocampus, Hypothalamus and Adrenal Gland, obtained from 20 pigs. Our analysis reveals that all tissues but Amygdala and Hippocampus possess distinct, tissue-specific expression pattern of miRNA that are associated with Hypoxia, stress responses as well as memory and fear conditioning. In particular, we observe marked differences in the expression profile of limbic tissues compared to those associated to the HPA/stress axis, with a surprisingly high aggregation of 3´-tRNA halves in Amygdala and Hippocampus. Since regulation of sncRNA and RNA cleavage plays a pivotal role in the central nervous system, our work provides seminal insights in the role/involvement of sncRNA in the transcriptional and post-transcriptional regulation of negative emotion, stress and coping behaviour in pigs, and mammals in general.

The hypothalamic-pituitary-adrenal (HPA) axis is an important component of neuroendocrine stress regulation and coping behavior. Transcriptome profiles of the hypothalamus and adrenal gland were assessed to identify molecular pathways and candidate genes for coping behavior in pigs. Ten each of high- (HR) and low- (LR) reactive pigs (n = 20) were selected for expression profiling based haplotype information of a prominent QTL-region on SSC12 discovered in our previous genome-wide association study (GWAS) on coping behavior. Comparing the HR and LR pigs showed 692 differentially expressed genes (DEGs) in the adrenal gland and 853 DEGs in the hypothalamus, respectively. Interestingly, 47% (17 out of 36) of DEGs found in both tissues were located in GWAS regions identified on SSC12, indicating that there are significant functional positional candidate genes for coping behaviour. Pathway analysis assigned DEGs to glucocorticoid receptor signaling in the adrenal gland. Furthermore, oxidative phosphorylation, mitochondrial dysfunction, and NGF signaling as well as cholecystokinin/Gastrin-mediated were identified in the hypothalamus. We narrowed the list of candidate genes in GWAS regions by analyzing their DEGs in the HPA axis. The top identified transcripts, including ATP1B2, AURKB, MPDU1 and NDEL1 provide evidence for molecular correlates of coping behavior in GWAS regions.

No abstract

No abstract

The molecular basis of porcine coping behavior (CB) relies on a sophisticated interplay of genetic and epigenetic features. Deep sequencing technologies allowed the identification of a plethora of new regulatory small non-coding RNA (sncRNA). We characterized mRNA and sncRNA profiles of central parts of the physiological stress response system including amygdala, hippocampus, hypothalamus and adrenal gland using systems biology for integration. Therefore, ten each of high- (HR) and low- (LR) reactive pigs (n = 20) carrying a CB associated haplotype in a prominent QTL-region on SSC12 were selected for mRNA and sncRNA expression profiling. The molecular markers related to the LR group included ATP1B2, MPDU1, miR-19b-5p, let-7g-5p, and 5′-tiRNALeu in the adrenal gland, miR-194a-5p, miR-125a-5p, miR-7-1-5p, and miR-107-5p in the hippocampus and CBL and PVRL1 in the hypothalamus. Interestingly, amygdalae of the LR group showed 5′-tiRNA and 5′-tRF (5′-tRFLys, 5′-tiRNALys, 5′-tiRNACys, and 5′-tiRNAGln) enrichment. Contrarily, molecular markers associated with the HR group encompassed miR-26b-5p, tRNAArg, tRNAGlyiF in the adrenal gland, IGF1 and APOD in the amygdala and PBX1, TOB1, and C18orf1 in the hippocampus and miR-24 in the hypothalamus. In addition, hypothalami of the HR group were characterized by 3′-tiRNA enrichment (3′-tiRNAGln, 3′-tiRNAAsn, 3′-tiRNAVal, 3′-tRFPro, 3′-tiRNACys, and 3′-tiRNAAla) and 3′-tRFs enrichment (3′-tRFAsn, 3′-tRFGlu, and 3′-tRFVal). These evidence suggest that tRNA-derived fragments and their cleavage activity are a specific marker for coping behavior. Data integration revealed new bio-signatures of important molecular interactions on a multi-transcript level in HPA axis and limbic system of pigs carrying a CB-associated haplotype.