Docker4Circ: A Framework for the Reproducible Characterization of circRNAs from RNA-Seq Data

Recent improvements in cost-effectiveness of high-throughput technologies has allowed RNA sequencing of total transcriptomes suitable for evaluating the expression and regulation of circRNAs, a relatively novel class of transcript isoforms with suggested roles in transcriptional and post-transcriptional gene expression regulation, as well as their possible use as biomarkers, due to their deregulation in various human diseases. A limited number of integrated workflows exists for prediction, characterization, and differential expression analysis of circRNAs, none of them complying with computational reproducibility requirements. We developed Docker4Circ for the complete analysis of circRNAs from RNA-Seq data. Docker4Circ runs a comprehensive analysis of circRNAs in human and model organisms, including: circRNAs prediction; classification and annotation using six public databases; back-splice sequence reconstruction; internal alternative splicing of circularizing exons; alignment-free circRNAs quantification from RNA-Seq reads; and differential expression analysis. Docker4Circ makes circRNAs analysis easier and more accessible thanks to: (i) its R interface; (ii) encapsulation of computational tasks into docker images; (iii) user-friendly Java GUI Interface availability; and (iv) no need of advanced bash scripting skills for correct use. Furthermore, Docker4Circ ensures a reproducible analysis since all its tasks are embedded into a docker image following the guidelines provided by Reproducible Bioinformatics Project

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.This work was supported by a Spanish Ministry of Economy grant AGL2013-49137-C3 to BMJ; EV-M is a predoctoral fellow from the University of the Basque Country (UPV/EHU) [PIF15/361] and received a grant for a short-term scientific mission (STSM) from the Functional Annotation of Animal Genomes—European network (FAANG-Europe) [COST Action CA15112

Patterns of Differentially Expressed circRNAs in Human Thymocytes

Circular RNAs (circRNAs) are suggested to play a discriminative role between some stages of thymocyte differentiation. However, differential aspects of the stage of mature single-positive thymocytes remain to be explored. The purpose of this study is to investigate the differential expression pattern of circRNAs in three different development stages of human thymocytes, including mature single-positive cells, and perform predictions in silico regarding the ability of specific circRNAs when controlling the expression of genes involved in thymocyte differentiation. We isolate human thymocytes at three different stages of intrathymic differentiation and determine the expression of circRNAs and mRNA by RNASeq. We show that the differential expression pattern of 50 specific circRNAs serves to discriminate between the three human thymocyte populations. Interestingly, the downregulation of RAG2, a gene involved in T-cell differentiation in the thymus, could be simultaneously controlled by the downregulation of two circRNASs (hsa_circ_0031584 and hsa_circ_0019079) through the hypothetical liberation of hsa-miR-609. Our study provides, for the first time, significant insights into the usefulness of circRNAs in discriminating between different stages of thymocyte differentiation and provides new potential circRNA-miRNA-mRNA networks capable of controlling the expression of genes involved in T-cell differentiation in the thymus.This work was financed by grants from the Spanish Ministry of Science, Innovation and Universities (MCIU)(RTI2018- 093330-B-I00; MCIU/FEDER, EU), Ramón Areces Foundation (CIVP19S7917); Autonomous Community of Madrid, Spain (B2017/BMD-3778; LINFOMAS-CM); the Spanish Association Against Cancer (AECC, 2018; PROYE18054PIRI); and the Spanish Ministry (Juan de la Cierva Grant IJCI-2016-29155). Institutional grants from the Fundación Ramón Areces and Banco de Santander to the CBMSO are also acknowledged.S

An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations.

Parietal cortex RNA-sequencing (RNA-seq) data were generated from individuals with and without Alzheimer disease (AD; ncontrol = 13; nAD = 83) from the Knight Alzheimer Disease Research Center (Knight ADRC). Using this and an independent (Mount Sinai Brain Bank (MSBB)) AD RNA-seq dataset, cortical circular RNA (circRNA) expression was quantified in the context of AD. Significant associations were identified between circRNA expression and AD diagnosis, clinical dementia severity and neuropathological severity. It was demonstrated that most circRNA-AD associations are independent of changes in cognate linear messenger RNA expression or estimated brain cell-type proportions. Evidence was provided for circRNA expression changes occurring early in presymptomatic AD and in autosomal dominant AD. It was also observed that AD-associated circRNAs co-expressed with known AD genes. Finally, potential microRNA-binding sites were identified in AD-associated circRNAs for miRNAs predicted to target AD genes. Together, these results highlight the importance of analyzing non-linear RNAs and support future studies exploring the potential roles of circRNAs in AD pathogenesis

The circular RNA landscape in multiple sclerosis: Disease-specific associated variants and exon methylation shape circular RNA expression profile

BACKGROUND: Circular RNAs (circRNAs) are a class of non-coding RNAs increasingly emerging as crucial actors in the pathogenesis of human diseases, including autoimmune and neurological disorders as multiple sclerosis (MS). Despite several efforts, the mechanisms regulating circRNAs expression are still largely unknown and the circRNA profile and regulation in MS-relevant cell models has not been completely investigated. In this work, we aimed at exploring the global landscape of circRNA expression in MS patients, also evaluating a possible correlation with their genetic and epigenetic background. METHODS: We performed RNA-seq experiments on circRNA-enriched samples, derived from peripheral blood mononuclear cells (PBMCs) of 10 MS patients and 10 matched controls and performed differential circRNA expression. The genetic background was evaluated using array genotyping, and an expression quantitative trait loci (eQTL) analysis was carried out. RESULTS: Expression analysis revealed 166 differentially expressed circRNAs in MS patients, 125 of which are downregulated. One of the top dysregulated circRNAs, hsa_circ_0007990, derives from the PGAP3 gene, encoding a protein relevant for the control of autoimmune responses. The downregulation of this circRNA was confirmed in two independent replication cohorts, suggesting its implementation as a possible RNA-based biomarker. The eQTL analysis evidenced a significant association between 89 MS-associated loci and the expression of at least one circRNA, suggesting that MS-associated variants could impact on disease pathogenesis by altering circRNA profiles. Finally, we found a significant correlation between exon methylation and circRNA expression levels, supporting the hypothesis that epigenetic features may play an important role in the definition of the cell circRNA pool. CONCLUSION: We described the circRNA expression profile of PBMCs in MS patients, suggesting that MS-associated variants may tune the expression levels of circRNAs acting as circ-QTLs , and proposing a role for exon-based DNA methylation in regulating circRNA expression

An Update on Circular RNA in Pediatric Cancers

Circular RNAs (circRNAs) are a class of single-stranded closed noncoding RNA molecules which are formed as a result of reverse splicing of mRNAs. Despite their relative abundance, only recently there appeared an increased interest in the understanding of their regulatory importance. Among their most relevant characteristics are high stability, abundance and evolutionary conservation among species. CircRNAs are implicated in several cellular functions, ranging from miRNA and protein sponges to transcriptional modulation and splicing. Additionally, circRNAs’ aberrant expression in pathological conditions is bringing to light their possible use as diagnostic and prognostic biomarkers. Their use as indicator molecules of pathological changes is also supported by their peculiar covalent closed cyclic structure which bestows resistance to RNases. Their regulatory role in cancer pathogenesis and metastasis is supported by studies involving human tumors that have investigated different expression profiles of these molecules. As endogenous competitive RNA, circRNAs can regulate tumor proliferation and invasion and they arouse great consideration as potential therapeutic biomarkers and targets for cancer. In this review, we describe the most recent findings on circRNAs in the most common pediatric solid cancers (such as brain tumors, neuroblastomas, and sarcomas) and in more rare ones (such as Wilms tumors, hepatoblastomas, and retinoblastomas)