Innovative construction of the first reliable catalogue of bovine circular RNAs.

The aim of this study was to compare the circular transcriptome of divergent tissues in order to understand: i) the presence of circular RNAs (circRNAs) that are not exonic circRNAs, i.e. originated from backsplicing involving known exons and, ii) the origin of artificial circRNA (artif_circRNA), i.e. circRNA not generated in-vivo. CircRNA identification is mostly an in-silico process, and the analysis of data from the BovReg project (https://www.bovreg.eu/) provided an opportunity to explore new ways to identify reliable circRNAs. By considering 117 tissue samples, we characterized 23,926 exonic circRNAs, 337 circRNAs from 273 introns (191 ciRNAs, 146 intron circles), 108 circRNAs from small non-coding genes and nearly 36.6K circRNAs classified as other_circRNAs. Furthermore, for 63 of those samples we analysed in parallel data from total-RNAseq (ribosomal RNAs depleted prior to library preparation) with paired mRNAseq (library prepared with poly(A)-selected RNAs). The high number of circRNAs detected in mRNAseq, and the significant number of novel circRNAs, mainly other_circRNAs, led us to consider all circRNAs detected in mRNAseq as artificial. This study provided evidence of 189 false entries in the list of exonic circRNAs: 103 artif_circRNAs identified by total RNAseq/mRNAseq comparison using two circRNA tools, 26 probable artif_circRNAs, and 65 identified by deep annotation analysis. Extensive benchmarking was performed (including analyses with CIRI2 and CIRCexplorer-2) and confirmed 94% of the 23,737 reliable exonic circRNAs. Moreover, this study demonstrates the effectiveness of a panel of highly expressed exonic circRNAs (5-8%) in analysing the tissue specificity of the bovine circular transcriptome.

Transcriptomic Response of Differentiating Porcine Myotubes to Thermal Stress and Donor Piglet Age.

Climate change is a current concern that directly and indirectly affects agriculture, especially the livestock sector. Neonatal piglets have a limited thermoregulatory capacity and are particularly stressed by ambient temperatures outside their optimal physiological range, which has a major impact on their survival rate. In this study, we focused on the effects of thermal stress (35 °C, 39 °C, and 41 °C compared to 37 °C) on differentiating myotubes derived from the satellite cells of Musculus rhomboideus, isolated from two different developmental stages of thermolabile 5-day-old (p5) and thermostable 20-day-old piglets (p20). Analysis revealed statistically significant differential expression genes (DEGs) between the different cultivation temperatures, with a higher number of genes responding to cold treatment. These DEGs were involved in the macromolecule degradation and actin kinase cytoskeleton categories and were observed at lower temperatures (35 °C), whereas at higher temperatures (39 °C and 41 °C), the protein transport system, endoplasmic reticulum system, and ATP activity were more pronounced. Gene expression profiling of HSP and RBM gene families, which are commonly associated with cold and heat responses, exhibited a pattern dependent on temperature variability. Moreover, thermal stress exhibited an inhibitory effect on cell cycle, with a more pronounced downregulation during cold stress driven by ADGR genes. Additionally, our analysis revealed DEGs from donors with an undeveloped thermoregulation capacity (p5) and those with a fully developed thermoregulation capacity (p20) under various cultivation temperature. The highest number of DEGs and significant GO terms was observed under temperatures of 35 °C and 37 °C. In particular, under 35 °C, the DEGs were enriched in insulin, thyroid hormone, and calcium signaling pathways. This result suggests that the different thermoregulatory capacities of the donor piglets determined the ability of the primary muscle cell culture to differentiate into myotubes at different temperatures. This work sheds new light on the underlying molecular mechanisms that govern piglet differentiating myotube response to thermal stress and can be leveraged to develop effective thermal management strategies to enhance skeletal muscle growth.

Brain Transcriptome Responses to Dexamethasone Depending on Dose and Sex Reveal Factors Contributing to Sex-Specific Vulnerability to Stress-Induced Disorders.

BACKGROUND: Glucocorticoid (GC) receptor (GR) signaling in the hypothalamus (Hyp) and in the superordinate limbic structures, such as the hippocampus (Hip), conveys feedback regulation of the neuroendocrine stress response and acts upon other neurobiological functions that ultimately influence mental health. These responses are strongly influenced by sex, but the molecular causes are still largely unexplored. METHODS: To investigate GR targets and their GC sensitivity in the Hyp and Hip, we treated juvenile male and female piglets with 10 (D10) or 60 (D60) µg/kg dexamethasone (DEX), a selective GR agonist, and analyzed transcriptome responses compared to a saline control group using RNA sequencing. RESULTS: Both doses influenced similar biological functions, including cellular response to lipid and immune cell-related functions, but the transcriptional response to D10 was considerably weaker, particularly in the Hip. Weighted Gene Co-expression Network Analysis revealed a network of genes coordinately regulated by DEX in both structures, among which the alpha-arrestin ARRDC2 takes a central position. Distinct functional groups of genes were differentially regulated by DEX between sexes depending on the dose; at D10, these included particularly mitochondrial genes, whereas at D60 interferon signaling and lipid homeostasis genes were enriched. The general and sex-specific transcriptional responses to DEX highlight microglia as the prominent target. Several key marker genes of disease-associated microglia were regulated by DEX depending on sex, such as TREM2 and LPL. CONCLUSION: The discovered expression signatures suggest that DEX induced a dysfunctional state of microglia in males, while in females microglia were primed, which could entail predisposition for different mental disorders.

Genetic background and production periods shape the microRNA profiles of the gut in laying hens.

There is growing evidence of the importance of miRNAs for intestinal functional properties and nutritional uptake. Comparative miRNAs profiles of the jejunal mucosa were established against two genetic backgrounds (Lohmann Brown-Classic (LB) and Lohmann LSL-Classic (LSL), which are similar in egg production but differ in physiological traits including mineral utilization, along the production periods of laying hens. The target genes of miRNAs higher expressed in LB vs. LSL (miR-126-3p, miR-214, miR-24-3p, miR-726-5p, miR-29b-3p) were enriched for energy pathways at all stages. The target genes of the miRNAs higher in LSL (miR-1788-5p, miR-103-3p, miR-22-5p, miR-221-3p, miR-375) were more enriched for immune and the bone signalling pathways. The most prominent expression differences were between 16 and 24 weeks of age before and after onset of laying. Our results evidence central roles of intestinal miRNAs as regulators of gene expression, influencing intestinal homeostasis and adaptation to environment in different strains and production phases.

MicroRNA-mRNA Networks in Pregnancy Complications: A Comprehensive Downstream Analysis of Potential Biomarkers.

Pregnancy complications are a major cause of fetal and maternal morbidity and mortality in humans. The majority of pregnancy complications initiate due to abnormal placental development and function. During the last decade, the role of microRNAs (miRNAs) in regulating placental and fetal development has become evident. Dysregulation of miRNAs in the placenta not only affects placental development and function, but these miRNAs can also be exported to both maternal and fetal compartments and affect maternal physiology and fetal growth and development. Due to their differential expression in the placenta and maternal circulation during pregnancy complications, miRNAs can be used as diagnostic biomarkers. However, the differential expression of a miRNA in the placenta may not always be reflected in maternal circulation, which makes it difficult to find a reliable biomarker for placental dysfunction. In this review, we provide an overview of differentially expressed miRNAs in the placenta and/or maternal circulation during preeclampsia (PE) and intrauterine growth restriction (IUGR), which can potentially serve as biomarkers for prediction or diagnosis of pregnancy complications. Using different bioinformatics tools, we also identified potential target genes of miRNAs associated with PE and IUGR, and the role of miRNA-mRNA networks in the regulation of important signaling pathways and biological processes.

Identification of the Key Molecular Drivers of Phosphorus Utilization Based on Host miRNA-mRNA and Gut Microbiome Interactions.

Phosphorus is an essential mineral for all living organisms and a limited resource worldwide. Variation and heritability of phosphorus utilization (PU) traits were observed, indicating the general possibility of improvement. Molecular mechanisms of PU, including host and microbial effects, are still poorly understood. The most promising molecules that interact between the microbiome and host are microRNAs. Japanese quail representing extremes for PU were selected from an F2 population for miRNA profiling of the ileal tissue and subsequent association with mRNA and microbial data of the same animals. Sixty-nine differentially expressed miRNAs were found, including 21 novel and 48 known miRNAs. Combining miRNAs and mRNAs based on correlated expression and target prediction revealed enrichment of transcripts in functional pathways involved in phosphate or bone metabolism such as RAN, estrogen receptor and Wnt signaling, and immune pathways. Out of 55 genera of microbiota, seven were found to be differentially abundant between PU groups. The study reveals molecular interactions occurring in the gut of quail which represent extremes for PU including miRNA-16-5p, miR-142b-5p, miR-148a-3p, CTDSP1, SMAD3, IGSF10, Bacteroides, and Alistipes as key indicators due to their trait-dependent differential expression and occurrence as hub-members of the network of molecular drivers of PU.

Epigenome-wide skeletal muscle DNA methylation profiles at the background of distinct metabolic types and ryanodine receptor variation in pigs.

BACKGROUND: Epigenetic variation may result from selection for complex traits related to metabolic processes or appear in the course of adaptation to mediate responses to exogenous stressors. Moreover epigenetic marks, in particular the DNA methylation state, of specific loci are driven by genetic variation. In this sense, polymorphism with major gene effects on metabolic and cell signaling processes, like the variation of the ryanodine receptors in skeletal muscle, may affect DNA methylation. METHODS: DNA-Methylation profiles were generated applying Reduced Representation Bisulfite Sequencing (RRBS) on 17 Musculus longissimus dorsi samples. We examined DNA methylation in skeletal muscle of pig breeds differing in metabolic type, Duroc and Pietrain. We also included F2 crosses of these breeds to get a first clue to DNA methylation sites that may contribute to breed differences. Moreover, we compared DNA methylation in muscle tissue of Pietrain pigs differing in genotypes at the gene encoding the Ca2+ release channel (RYR1) that largely affects muscle physiology. RESULTS: More than 2000 differently methylated sites were found between breeds including changes in methylation profiles of METRNL, IDH3B, COMMD6, and SLC22A18, genes involved in lipid metabolism. Depending on RYR1 genotype there were 1060 differently methylated sites including some functionally related genes, such as CABP2 and EHD, which play a role in buffering free cytosolic Ca2+ or interact with the Na+/Ca2+ exchanger. CONCLUSIONS: The change in the level of methylation between the breeds is probably the result of the long-term selection process for quantitative traits involving an infinite number of genes, or it may be the result of a major gene mutation that plays an important role in muscle metabolism and triggers extensive compensatory processes.

Development and evaluation of a milk protein transcript depletion method for differential transcriptome analysis in mammary gland tissue.

BACKGROUND: In the mammary gland transcriptome of lactating dairy cows genes encoding milk proteins are highly abundant, which can impair the detection of lowly expressed transcripts and can bias the outcome in global transcriptome analyses. Therefore, the aim of this study was to develop and evaluate a method to deplete extremely highly expressed transcripts in mRNA from lactating mammary gland tissue. RESULTS: Selective RNA depletion was performed by hybridization of antisense oligonucleotides targeting genes encoding the caseins (CSN1S1, CSN1S2, CSN2 and CSN3) and whey proteins (LALBA and PAEP) within total RNA followed by RNase H-mediated elimination of the respective transcripts. The effect of the RNA depletion procedure was monitored by RNA sequencing analysis comparing depleted and non-depleted RNA samples from Escherichia coli (E. coli) challenged and non-challenged udder tissue of lactating cows in a proof of principle experiment. Using RNase H-mediated RNA depletion, the ratio of highly abundant milk protein gene transcripts was reduced in all depleted samples by an average of more than 50% compared to the non-depleted samples. Furthermore, the sensitivity for discovering transcripts with marginal expression levels and transcripts not yet annotated was improved. Finally, the sensitivity to detect significantly differentially expressed transcripts between non-challenged and challenged udder tissue was increased without leading to an inadvertent bias in the pathogen challenge-associated biological signaling pathway patterns. CONCLUSIONS: The implementation of selective RNase H-mediated RNA depletion of milk protein gene transcripts from the mammary gland transcriptome of lactating cows will be highly beneficial to establish comprehensive transcript catalogues of the tissue that better reflects its transcriptome complexity.

Deep sequencing of small non-coding RNA highlights brain-specific expression patterns and RNA cleavage.

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.

Unmapped reads from cattle RNAseq data: A source for missing and misassembled sequences in the reference assemblies and for detection of pathogens in the host.

Usually, reads from transcriptome sequencing data unmapped to the target species' reference genome are disregarded. A recent RNAseq project on the new fatal disease Bovine Neonatal Pancytopenia had indicated an unexplained immune response signature to a double-stranded RNA virus. To unravel its background, contigs were de novo assembled from unmapped RNAseq reads and aligned against the bovine genome assemblies and multispecies NCBI databases. Lack of genuine virus sequence contigs rejected the hypothesis of a live virus being causal for the unexplained immune response. Alignment data also demonstrated incomplete bovine reference genome assemblies. In addition, we found that several parasite and virus genome reference assemblies in NCBI were contaminated with bovine DNA and confirmed recombination of bovine DNA into BVD virus strains. Exploring unmapped reads can extract useful biological information regarding the presence of microorganisms and can highlight issues with reference genome assemblies of host and pathogen species.