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.

A unique circular RNA expression pattern in the peripheral blood of myalgic encephalomyelitis/chronic fatigue syndrome patients.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a debilitating disease with obscure aetiology. The underdiagnosis rate of ME/CFS is high due to the lack of diagnostic criteria based on objective markers. In recent years, circRNAs have emerged as potential genetic biomarkers for neurological diseases, including Parkinson's disease and Alzheimer's disease, making them likely to have the same prospect of being biomarkers in ME/CFS. However, despite the extensive amount of research that has been performed on the transcriptomes of ME/CFS patients, all of them are solely focused on linear RNAs, and the profiling of circRNAs in ME/CFS has been completely omitted. In this study, we investigated the expression profiles of circRNAs, comparing ME/CFS patients and controls before and after two sessions of cardiopulmonary exercise longitudinally. In patients with ME/CFS, the number of detected circRNAs was higher compared to healthy controls, indicating potential differences in circRNA expression associated with the disease. Additionally, healthy controls showed an increase in the number of circRNAs following exercise testing, while no similar pattern was evident in ME/CFS patients, further highlighting physiological differences between the two groups. A lack of correlation was observed between differentially expressed circRNAs and their corresponding coding genes in terms of expression and function, suggesting the potential of circRNAs as independent biomarkers in ME/CFS. Specifically, 14 circRNAs were highly expressed in ME/CFS patients but absent in controls throughout the exercise study, indicating a unique molecular signature specific to ME/CFS patients and providing potential diagnostic biomarkers for the disease. Significant enrichment of protein and gene regulative pathways were detected in relation to five of these 14 circRNAs based on their predicted miRNA target genes. Overall, this is the first study to describe the circRNA expression profile in peripheral blood of ME/CFS patients, providing valuable insights into the molecular mechanisms underlying the disease.

Comparative Analysis of the Circular Transcriptome in Muscle, Liver, and Testis in Three Livestock Species.

Circular RNAs have been observed in a large number of species and tissues and are now recognized as a clear component of the transcriptome. Our study takes advantage of functional datasets produced within the FAANG consortium to investigate the pervasiveness of circular RNA transcription in farm animals. We describe here the circular transcriptional landscape in pig, sheep and bovine testicular, muscular and liver tissues using total 66 RNA-seq datasets. After an exhaustive detection of circular RNAs, we propose an annotation of exonic, intronic and sub-exonic circRNAs and comparative analyses of circRNA content to evaluate the variability between individuals, tissues and species. Despite technical bias due to the various origins of the datasets, we were able to characterize some features (i) (ruminant) liver contains more exonic circRNAs than muscle (ii) in testis, the number of exonic circRNAs seems associated with the sexual maturity of the animal. (iii) a particular class of circRNAs, sub-exonic circRNAs, are produced by a large variety of multi-exonic genes (protein-coding genes, long non-coding RNAs and pseudogenes) and mono-exonic genes (protein-coding genes from mitochondrial genome and small non-coding genes). Moreover, for multi-exonic genes there seems to be a relationship between the sub-exonic circRNAs transcription level and the linear transcription level. Finally, sub-exonic circRNAs produced by mono-exonic genes (mitochondrial protein-coding genes, ribozyme, and sno) exhibit a particular behavior. Caution has to be taken regarding the interpretation of the unannotated circRNA proportion in a given tissue/species: clusters of circRNAs without annotation were characterized in genomic regions with annotation and/or assembly problems of the respective animal genomes. This study highlights the importance of improving genome annotation to better consider candidate circRNAs and to better understand the circular transcriptome. Furthermore, it emphasizes the need for considering the relative "weight" of circRNAs/parent genes for comparative analyses of several circular transcriptomes. Although there are points of agreement in the circular transcriptome of the same tissue in two species, it will be not possible to do without the characterization of it in both species.

Cell type-specific circular RNA expression in human glial cells.

The circular transcriptome of human glial cells is an area of neuroscience that has not been thoroughly elucidated. Circular RNAs (circRNAs) have the potential to facilitate the understanding of vast, complex and unknown mechanisms derived from the human transcriptome, including elements of the human brain that are not known and the evolution of the human brain, the complexities of which are not well understood. Moreover, the glial cells have been determined to contribute to human brain evolution. This study presents the first comprehensive analysis of the human brain glia circRNA transcriptome, that is, astrocytes, microglia and oligodendrocytes. After stringent criteria applied to the detection of circRNAs, it was found that the circular transcriptomes of these glia are unique from one another, and hence might be indicative of distinct roles for circRNAs within the brain. This study found 265, 239 and 442 circRNAs comprising the unique circular transcriptome of astrocytes, microglia and oligodendrocytes, respectively. The most abundant circRNAs in these glial cell types are expressed by parent genes co-expressing linear RNAs in low abundance, suggesting spliceosome activity favorable to the back-splicing mechanism instead of canonical splicing activity.

Neural circular transcriptomes across mammalian species.

Circular RNAs (circRNAs) have recently attracted significant interest in the realm of science and the evolution of species. Given the lack of information available on circRNAs due to various barriers related to sequencing techniques and bioinformatics tools, little regarding their function is known. It has been predicted that circRNAs contribute to gene expression regulation, but aside from a few specific cases, this contention has yet to be proven. Although the role of circRNAs in evolution remains elusive, from the few studies that have shown circRNA conservation in mammalian species, tissue specificity in brain regions, and the abundance of circRNAs in the brains of various species, the concept is becoming more likely with much gravitas. The proposed functional role of circRNAs being gene regulators is of great interest and would provide a basis to further understand not only the functional capabilities of organisms, but also the evolution of mammalian species.

Sequence and structural properties of circular RNAs in the brain of nurse and forager honeybees (Apis mellifera).

BACKGROUND: The honeybee (Apis mellifera) represents a model organism for social insects displaying behavioral plasticity. This is reflected by an age-dependent task allocation. The most protruding tasks are performed by young nurse bees and older forager bees that take care of the brood inside the hive and collect food from outside the hive, respectively. The molecular mechanism leading to the transition from nurse bees to foragers is currently under intense research. Circular RNAs, however, were not considered in this context so far. As of today, this group of non-coding RNAs was only known to exist in two other insects, Drosophila melanogaster and Bombyx mori. Here we complement the state of circular RNA research with the first characterization in a social insect. RESULTS: We identified numerous circular RNAs in the brain of A. mellifera nurse bees and forager bees using RNA-Seq with exonuclease enrichment. Presence and circularity were verified for the most abundant representatives. Back-splicing in honeybee occurs further towards the end of transcripts and in transcripts with a high number of exons. The occurrence of circularized exons is correlated with length and CpG-content of their flanking introns. The latter coincides with increased DNA-methylation in the respective loci. For two prominent circular RNAs the abundance in worker bee brains was quantified in TaqMan assays. In line with previous findings of circular RNAs in Drosophila, circAmrsmep2 accumulates with increasing age of the insect. In contrast, the levels of circAmrad appear age-independent and correlate with the bee's task. Its parental gene is related to amnesia-resistant memory. CONCLUSIONS: We provide the first characterization of circRNAs in a social insect. Many of the RNAs identified here show homologies to circular RNAs found in Drosophila and Bombyx, indicating that circular RNAs are a common feature among insects. We find that exon circularization is correlated to DNA-methylation at the flanking introns. The levels of circAmrad suggest a task-dependent abundance that is decoupled from age. Moreover, a GO term analysis shows an enrichment of task-related functions. We conclude that circular RNAs could be relevant for task allocation in honeybee and should be investigated further in this context.

Identification of circular RNA in the Bombyx mori silk gland.

Bombyx mori is an economically important holometabolous lepidopteran insect. In B. mori endogenous noncoding RNAs such as microRNAs (miRNAs) and Piwi-interacting RNAs play crucial biological functions in metamorphosis and sex determination. In addition, circular RNAs (circRNAs) have been recently identified as noncoding RNAs in most common model organisms and show potential as gene regulators. However, to date, there have been few studies on the circRNAs present in the B. mori genome conducted to date. Here, we identified 3916 circRNAs by deep circular transcriptome sequencing using the silk gland of B. mori. 3155 circRNAs were found to be derived from 1727 parental genes. The circRNAs displayed tissue-specific expression between the middle silk gland (MSG) and posterior silk gland (PSG), with 2532 and 880 being upregulated circRNAs in the MSG and PSG, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the parental genes from the MSG and PSG were generally annotated to similar categories and pathways. The interaction network of circRNAs and miRNAs showed that circRNAs might act as miRNA sponges or interact with miRNAs in some other way. Overall, the results revealed the complicated patterns of circRNAs in the B. mori silk gland providing a new angle from which to explore the mechanisms of complex gene regulation and efficient silk protein synthesis.