Non-coding RNAs in ovine immunity: Identification of unannotated genes and functional analyses of high throughput genomic data

Abstract

Non-coding RNAs (ncRNAs) are involved in several biological processes in mammals, including the immune system response to pathogens and vaccines. The annotation and functional characterization of two of the main classes of ncRNAs, microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), is more advanced in humans than in livestock species, and thus, there is limited knowledge about the function of these transcripts. The main objective of this work was the identification of ovine non-coding genes, concretely miRNA and lncRNA genes, that are involved in the innate and adaptive immune responses induced by vaccines, vaccine components and pathogen infections. For this purpose, high-throughput transcriptome sequencing datasets produced for this purpose and datasets publicly available were analysed with bioinformatic tools and workflows in order to identify unannotated non-coding genes, profile their expression in different tissues and perform evolutionary conservation analyses. More than 12000 unannotated ovine lncRNAs and 1000 ovine miRNAs were identified in the different analyses, with varying levels of sequence conservation. Differential expression analyses between unstimulated samples and samples stimulated with pathogen infection or vaccination resulted in hundreds of lncRNAs and miRNAs with changed expression. Gene co-expression analyses revealed immune gene-enriched clusters associated with immune system activation. These genes make up a prioritized set of potential candidates for deeper experimental analyses. Taken together, these results should help completing the sheep non-coding gene catalogue, and most importantly, they give evidence of immune state-specific ncRNA expression patterns in a livestock species.