The landscape of long noncoding RNA expression in the goat brain.

The brain regulates multiple metabolic processes, such as food intake, energy expenditure, insulin secretion, hepatic glucose production, and glucose and fatty acid metabolism in adipose tissue, which are fundamental for the maintenance of energy and glucose homeostasis during lactation and pregnancy. In addition, brain expression has a fundamental impact on the development of maternal behavior. Although brain functions are partly regulated by long noncoding RNAs (lncRNAs), their expression profiles have not been characterized in depth in any ruminant species. We have sequenced the transcriptome of 12 brain tissues from 3 goats that were 1 mo pregnant and 4 nonpregnant goats to investigate their lncRNA expression patterns. Between 4,363 (adenohypophysis) and 4,604 (olfactory bulb) lncRNAs were expressed in brain tissues, leading us to establish a set of 794 already annotated lncRNAs and 5,098 novel lncRNA candidates. The detected lncRNAs shared features with those of other mammals, and tissue-specific lncRNAs were enriched in brain development-related terms. Differential expression analyses between goats that were 1 mo pregnant and nonpregnant goats showed that the lncRNA expression profiles of certain brain regions experience substantial changes associated with early pregnancy (238 lncRNAs are differentially expressed in the olfactory bulb), but others do not. Enrichment analysis showed that differentially expressed lncRNAs from the olfactory bulb are co-expressed with genes previously linked to behavioral changes related to pregnancy. These findings provide a first characterization of the landscape of lncRNA expression in the goat brain and provides valuable clues to understand the molecular events triggered by early pregnancy in the central nervous system.

Identification and characterization of miRNAs in spleens of sheep subjected to repetitive vaccination.

Accumulative evidence has shown that short non-coding RNAs such as miRNAs can regulate the innate and adaptive immune responses. Aluminium hydroxide is a commonly used adjuvant in human and veterinary vaccines. Despite its extended use, its mechanism of action is not fully understood and very few in vivo studies have been done to enhance understanding at the molecular level. In this work, we took advantage of a previous long-term experiment in which lambs were exposed to three different treatments by parallel subcutaneous inoculations with aluminium-containing commercial vaccines, an equivalent dose of aluminium or mock injections. Spleen samples were used for miRNA-seq. A total of 46 and 16 miRNAs were found differentially expressed when animals inoculated with commercial vaccines or the adjuvant alone were compared with control animals, respectively. Some miRNAs previously related to macrophage polarization were found dysregulated exclusively by the commercial vaccine treatment but not in the aluminium inoculated animals. The dysregulated miRNAs in vaccine group let-7b-5p, miR-29a-3p, miR-27a and miR-101-3p are candidates for further research, since they may play key roles in the immune response induced by aluminium adjuvants added to vaccines. Finally, protein-protein interaction network analysis points towards leucocyte transendothelial migration as a specific mechanism in animals receiving adjuvant only.

The sheep miRNAome: Characterization and distribution of miRNAs in 21 tissues.

MicroRNAs (miRNAs) are evolutionarily conserved small non-coding RNAs that regulate gene expression. In livestock species, miRNAs also show great potential as biomarkers for animal health and product quality. In sheep, few miRNAs have been described in comparison with other livestock species or model organisms. We uniformly analyzed 172 public ovine small RNA sequencing datasets from 21 different tissues in order to predict conserved and novel miRNA precursors and profile their expression patterns. In addition to the 106 annotated sheep miRNAs, 1047 previously unannotated miRNA precursor sequences were detected and 41 % of them were assigned an orthologue from other close species. Considering expression levels, a set of miRNAs with high sequence conservation were detected in all tissues, while 733 mature miRNAs were robustly expressed in at least one tissue. 270 miRNAs showed high tissue specificity index values. Brain, male reproductive tissues and PBMCs showed the most distinct expression patterns. Strikingly, over one hundred precursors from the ruminant specific family of mir-2284/mir-2285 miRNAs were found, which were enriched in immune related tissues. This work supports the known high conservation of many miRNAs, but also highlights the potential of clade-specific innovations in ruminant evolution.

Identification of sheep lncRNAs related to the immune response to vaccines and aluminium adjuvants.

BACKGROUND: Long non-coding RNAs (lncRNAs) are involved in several immune processes, including the immune response to vaccination, but most of them remain uncharacterised in livestock species. The mechanism of action of aluminium adjuvants as vaccine components is neither not fully understood. RESULTS: We built a transcriptome from sheep PBMCs RNA-seq data in order to identify unannotated lncRNAs and analysed their expression patterns along protein coding genes. We found 2284 novel lncRNAs and assessed their conservation in terms of sequence and synteny. Differential expression analysis performed between animals inoculated with commercial vaccines or aluminium adjuvant alone and the co-expression analysis revealed lncRNAs related to the immune response to vaccines and adjuvants. A group of co-expressed genes enriched in cytokine signalling and production highlighted the differences between different treatments. A number of differentially expressed lncRNAs were correlated with a divergently located protein-coding gene, such as the OSM cytokine. Other lncRNAs were predicted to act as sponges of miRNAs involved in immune response regulation. CONCLUSIONS: This work enlarges the lncRNA catalogue in sheep and puts an accent on their involvement in the immune response to repetitive vaccination, providing a basis for further characterisation of the non-coding sheep transcriptome within different immune cells.

The dystrophia myotonica WD repeat-containing protein DMWD and WDR20 differentially regulate USP12 deubiquitinase.

Despite its potential clinical relevance, the product of the DMWD (dystrophia myotonica, WD repeat containing) gene is a largely uncharacterized protein. The DMWD amino acid sequence is similar to that of WDR20, a known regulator of the USP12 and USP46 deubiquitinases (DUBs). Here, we apply a combination of in silico and experimental methods to investigate several aspects of DMWD biology. Molecular evolution and phylogenetic analyses reveal that WDR20 and DMWD, similar to USP12 and USP46, arose by duplication of a common ancestor during the whole genome duplication event in the vertebrate ancestor lineage. The analysis of public human gene expression datasets suggests that DMWD expression is positively correlated with USP12 expression in normal tissues and negatively correlated with WDR20 expression in tumors. Strikingly, a survey of the annotated interactome for DMWD and WDR20 reveals a largely nonoverlapping set of interactors for these proteins. Experimentally, we first confirmed that DMWD binds both USP12 and USP46 through direct coimmunoprecipitation of epitope-tagged proteins. We found that DMWD and WDR20 share the same binding interface in USP12, suggesting that their interaction with the DUB may be mutually exclusive. Finally, we show that both DMWD and WDR20 promote USP12 enzymatic activity, but they differentially modulate the subcellular localization of the DUB. Altogether, our findings suggest a model whereby mutually exclusive binding of DMWD and WDR20 to USP12 may lead to formation of deubiquitinase complexes with distinct subcellular localization, potentially targeting different substrate repertoires.

Novel circRNA discovery in sheep shows evidence of high backsplice junction conservation.

Circular RNAs (circRNAs) are covalently closed circular non-coding RNAs. Due to their structure, circRNAs are more stable and have longer half-lives than linear RNAs making them good candidates for disease biomarkers. Despite the scientific relevance of these molecules, the study of circRNAs in non-model organisms is still in its infancy. Here, we analyse total RNA-seq data to identify circRNAs in sheep from peripheral blood mononuclear cells (PBMCs) and parietal lobe cortex. Out of 2510 and 3403 circRNAs detected in parietal lobe cortex and in PBMCs, a total of 1379 novel circRNAs were discovered. Remarkably, around 63% of all detected circRNAs were found to be completely homologous to a circRNA annotated in human. Functional enrichment analysis was conducted for both tissues based on GO terms and KEGG pathways. The enriched terms suggest an important role of circRNAs from encephalon in synaptic functions and the involvement of circRNAs from PBMCs in basic immune system functions. In addition to this, we investigated the role of circRNAs in repetitive vaccination experiments via differential expression analysis and did not detect any significant relationship. At last, our results support both the miRNA sponge and the miRNA shuttle functions of CDR1-AS in sheep brain. To our knowledge, this is the first study on circRNA annotation in sheep PBMCs or parietal lobe cortex samples.

Whole transcriptome approach to evaluate the effect of aluminium hydroxide in ovine encephalon.

Aluminium hydroxide adjuvants are crucial for livestock and human vaccines. Few studies have analysed their effect on the central nervous system in vivo. In this work, lambs received three different treatments of parallel subcutaneous inoculations during 16 months with aluminium-containing commercial vaccines, an equivalent dose of aluminium hydroxide or mock injections. Brain samples were sequenced by RNA-seq and miRNA-seq for the expression analysis of mRNAs, long non-coding RNAs and microRNAs and three expression comparisons were made. Although few differentially expressed genes were identified, some dysregulated genes by aluminium hydroxide alone were linked to neurological functions, the lncRNA TUNA among them, or were enriched in mitochondrial energy metabolism related functions. In the same way, the miRNA expression was mainly disrupted by the adjuvant alone treatment. Some differentially expressed miRNAs had been previously linked to neurological diseases, oxidative stress and apoptosis. In brief, in this study aluminium hydroxide alone altered the transcriptome of the encephalon to a higher degree than commercial vaccines that present a milder effect. The expression changes in the animals inoculated with aluminium hydroxide suggest mitochondrial disfunction. Further research is needed to elucidate to which extent these changes could have pathological consequences.

RNA-Seq analysis of ileocecal valve and peripheral blood from Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis revealed dysregulation of the CXCL8/IL8 signaling pathway.

Paratuberculosis is chronic granulomatous enteritis of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). Whole RNA-sequencing (RNA-Seq) is a promising source of novel biomarkers for early MAP infection and disease progression in cattle. Since the blood transcriptome is widely used as a source of biomarkers, we analyzed whether it recapitulates, at least in part, the transcriptome of the ileocecal valve (ICV), the primary site of MAP colonization. Total RNA was prepared from peripheral blood (PB) and ICV samples, and RNA-Seq was used to compare gene expression between animals with focal or diffuse histopathological lesions in gut tissues versus control animals with no detectable signs of infection. Our results demonstrated both shared, and PB and ICV-specific gene expression in response to a natural MAP infection. As expected, the number of differentially expressed (DE) genes was larger in the ICV than in the PB samples. Among the DE genes in the PB and ICV samples, there were some common genes irrespective of the type of lesion including the C-X-C motif chemokine ligand 8 (CXCL8/IL8), apolipoprotein L (APOLD1), and the interferon inducible protein 27 (IFI27). The biological processes (BP) enriched in the PB gene expression profiles from the cows with diffuse lesions included the killing of cells of other organism, defense response, immune response and the regulation of neutrophil chemotaxis. Two of these BP, the defense and immune response, were also enriched in the ICV from the cows with diffuse lesions. Metabolic analysis of the DE genes revealed that the N-glycan biosynthesis, bile secretion, one-carbon pool by folate and purine metabolism were significantly enriched in the ICV from the cows with focal lesions. In the ICV from cows with diffuse lesions; the valine, leucine and isoleucine degradation route, purine metabolism, vitamin digestion and absorption and the cholesterol routes were enriched. Some of the identified DE genes, BP and metabolic pathways will be studied further to develop novel diagnostic tools, vaccines and immunotherapeutics.

Expression analysis of lung miRNAs responding to ovine VM virus infection by RNA-seq.

BACKGROUND: MicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease. RESULTS: Using RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways. CONCLUSIONS: To the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep's lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

Molecular Signature of Aluminum Hydroxide Adjuvant in Ovine PBMCs by Integrated mRNA and microRNA Transcriptome Sequencing.

There have been few in vivo studies on the effect of aluminum hydroxide adjuvant and its influence on the immune response to vaccination. In this study, lambs received a parallel subcutaneous treatment with either commercial vaccines containing aluminum hydroxide or an equivalent dose of this compound only with the aim of identifying the activated molecular signature. Blood samples were taken from each animal at the beginning and at the end of the experiment and PBMCs isolated. Total RNA and miRNA libraries were prepared and sequenced. After alignment to the Oar3.1 reference genome and differential expression with 3 programs, gene enrichment modeling was performed. For miRNAs, miRBase and RNAcentral databases were used for detection and characterization. Three expression comparisons were made: vaccinated animals at the beginning and at the end of the treatment, adjuvanted animals at the same times, and animals of both treatments at the end of the experiment. After exposure to both treatments, a total of 2,473; 2,980 and 429 differentially expressed genes were identified in vaccinated animals, adjuvanted animals and animals at the end of both treatments, respectively. In both adjuvant and vaccine treated animals the NF-κB signaling pathway was enriched. On the other hand, it can be observed a downregulation of cytokines and cytokine receptors in the adjuvanted group compared to the vaccinated group at the final time, suggesting a milder induction of the immune response when the adjuvant is alone. As for the miRNA analysis, 95 miRNAs were detected: 64 previously annotated in Ovis aries, 11 annotated in Bos taurus and 20 newly described. Interestingly, 6 miRNAs were differentially expressed in adjuvant treated animals, and 3 and 1 in the other two comparisons. Lastly, an integrated miRNA-mRNA expression profile was developed, in which a miRNA-mediated regulation of genes related to DNA damage stimulus was observed. In brief, it seems that aluminum-containing adjuvants are not simple delivery vehicles for antigens, but also induce endogenous danger signals that can stimulate the immune system. Whether this contributes to long-lasting immune activation or to the overstimulation of the immune system remains to be elucidated.