Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

circRIP: an accurate tool for identifying circRNA-RBP interactions.

Circular ribonucleic acids (RNAs) (circRNAs) are formed by covalently linking the downstream splice donor and the upstream splice acceptor. One of the most important functions of circRNAs is mainly exerted through binding RNA-binding proteins (RBPs). However, there is no efficient algorithm for identifying genome-wide circRNA-RBP interactions. Here, we developed a unique algorithm, circRIP, for identifying circRNA-RBP interactions from RNA immunoprecipitation sequencing (RIP-Seq) data. A simulation test demonstrated the sensitivity and specificity of circRIP. By applying circRIP, we identified 95 IGF2BP3-binding circRNAs based on the IGF2BP3 RIP-Seq dataset. We further identified 2823 and 1333 circRNAs binding to >100 RBPs in K562 and HepG2 cell lines, respectively, based on enhanced cross-linking immunoprecipitation (eCLIP) data, demonstrating the significance to survey the potential interactions between circRNAs and RBPs. In this study, we provide an accurate and sensitive tool, circRIP (https://github.com/bioinfolabwhu/circRIP), to systematically identify RBP and circRNA interactions from RIP-Seq and eCLIP data, which can significantly benefit the research community for the functional exploration of circRNAs.

CSCD2: an integrated interactional database of cancer-specific circular RNAs.

The significant function of circRNAs in cancer was recognized in recent work, so a well-organized resource is required for characterizing the interactions between circRNAs and other functional molecules (such as microRNA and RNA-binding protein) in cancer. We previously developed cancer-specific circRNA database (CSCD), a comprehensive database for cancer-specific circRNAs, which is widely used in circRNA research. Here, we updated CSCD to CSCD2 (http://geneyun.net/CSCD2 or http://gb.whu.edu.cn/CSCD2), which includes significantly more cancer-specific circRNAs identified from a large number of human cancer and normal tissues/cell lines. CSCD2 contains >1000 samples (825 tissues and 288 cell lines) and identifies a large number of circRNAs: 1 013 461 cancer-specific circRNAs, 1 533 704 circRNAs from only normal samples and 354 422 circRNAs from both cancer and normal samples. In addition, CSCD2 predicts potential miRNA-circRNA and RBP-circRNA interactions using binding motifs from >200 RBPs and 2000 microRNAs. Furthermore, the potential full-length and open reading frame sequence of these circRNAs were also predicted. Collectively, CSCD2 provides a significantly enhanced resource for exploring the function and regulation of circRNAs in cancer.

ADEIP: an integrated platform of age-dependent expression and immune profiles across human tissues.

Gene expression and immune status in human tissues are changed with aging. There is a need to develop a comprehensive platform to explore the dynamics of age-related gene expression and immune profiles across tissues in genome-wide studies. Here, we collected RNA-Seq datasets from GTEx project, containing 16 704 samples from 30 major tissues in six age groups ranging from 20 to 79 years old. Dynamic gene expression along with aging were depicted and gene set enrichment analysis was performed among those age groups. Genes from 34 known immune function categories and immune cell compositions were investigated and compared among different age groups. Finally, we integrated all the results and developed a platform named ADEIP (http://gb.whu.edu.cn/ADEIP or http://geneyun.net/ADEIP), integrating the age-dependent gene expression and immune profiles across tissues. To demonstrate the usage of ADEIP, we applied two datasets: severe acute respiratory syndrome coronavirus 2 and human mesenchymal stem cells-assoicated genes. We also included the expression and immune dynamics of these genes in the platform. Collectively, ADEIP is a powerful platform for studying age-related immune regulation in organogenesis and other infectious or genetic diseases.

Numerical and experimental investigations on the propagation property of a mid-infrared 7 × 1 multimode fiber combiner.

Mid-infrared fiber combiners have great potential in power and spectral combination. However, studies on mid-infrared transmission optical field distributions using these combiners are limited. In this study, we designed and fabricated a 7 × 1 multimode fiber combiner based on sulfur-based glass fibers and observed approximately 80% per-port transmission efficiency at 4.778 µm wavelength. We investigated the propagation properties of the prepared combiners and explored the effects of transmission wavelength, output fiber length, and fusion deviation on the transmitted optical field and beam quality factor M2. Additionally, we assessed the effect of coupling on the excitation mode and spectral combination of the mid-infrared fiber combiner for multiple light sources. Our results provide an in-depth understanding of the propagation properties of the mid-infrared multimode fiber combiners, which may find applications in high-beam-quality laser devices.

Silencing of circIgf1r plays a protective role in neuronal injury via regulating astrocyte polarization during epilepsy.

Epilepsy is a common brain disorder, repeated seizures of epilepsy may lead to a series of brain pathological changes such as neuronal or glial damage. However, whether circular RNAs are involved in neuronal injury during epilepsy is not fully understood. Here, we screened circIgf1r in the status epilepticus model through circRNA sequencing, and found that it was upregulated after the status epilepticus model through QPCR analysis. Astrocytes polarizing toward neurotoxic A1 phenotype and neurons loss were observed after status epilepticus. Through injecting circIgf1r siRNA into the lateral ventricle, it was found that knocking down circIgf1r in vivo would induce the polarization of astrocytes to phenotype A2 and reduce neuronal loss. The results in vitro further confirmed that inhibiting the expression of circIgf1r in astrocytes could protect neurons by converting reactive astrocytes from A1 to the protective A2. In addition, knocking down circIgf1r in astrocytes could functionally promote astrocyte autophagy and relieve the destruction of 4-AP-induced autophagy flux. In terms of mechanism, circIgf1r promoted the polarization of astrocytes to phenotype A1 by inhibiting autophagy. Taken together, our results reveal circIgf1r may serve as a potential target for the prevention and treatment of neuron damage after epilepsy.

CircRNA inhibits DNA damage repair by interacting with host gene.

BACKGROUND: Deregulated circular RNAs (circRNAs) are associated with the development of cancer and therapy resistance. However, functional research of circRNAs mostly focus on potential miRNA or protein binding and more potential regulation of circRNA on host gene DNA in cancers are yet to be inspected. METHOD: We performed total RNA sequencing on clinical breast cancer samples and identified the expression patterns of circRNAs and corresponding host genes in patient blood, tumor and adjacent normal tissues. qPCR, northern blot and in situ hybridization were used to validate the dysregulation of circRNA circSMARCA5. A series of procedures including R-loop dot-blotting, DNA-RNA immunoprecipitation and mass spectrum, etc. were conducted to explore the regulation of circSMARCA5 on the transcription of exon 15 of SMARCA5. Moreover, immunofluorescence and in vivo experiments were executed to investigate the overexpression of circSMARCA5 with drug sensitivities. RESULTS: We found that circRNAs has average higher expression over its host linear genes in peripheral blood. Compared to adjacent normal tissues, circSMARCA5 is decreased in breast cancer tissues, contrary to host gene SMARCA5. The enforced expression of circSMARCA5 induced drug sensitivity of breast cancer cell lines in vitro and in vivo. Furthermore, we demonstrated that circSMARCA5 can bind to its parent gene locus, forming an R-loop, which results in transcriptional pausing at exon 15 of SMARCA5. CircSMARCA5 expression resulted in the downregulation of SMARCA5 and the production of a truncated nonfunctional protein, and the overexpression of circSMARCA5 was sufficient to improve sensitivity to cytotoxic drugs. CONCLUSION: Our results revealed a new regulatory mechanism for circRNA on its host gene and provided evidence that circSMARCA5 may serve as a therapeutic target for drug-resistant breast cancer patients.

CircView: a visualization and exploration tool for circular RNAs.

Circular RNAs (circRNAs) are novel rising stars of noncoding RNAs, which are highly abundant and evolutionarily conserved across species. Number of publications related to circRNAs increased sharply in recent years, representing emerging focuses in the field. Therefore, tools, pipelines and databases have been developed to identify and store circRNAs. However, there is no existing tool to visualize and explore circRNAs. Therefore, we introduce CircView, a user-friendly visualization tool for circRNAs detected from existing tools. CircView enables users to visualize circRNAs and to quantify number of samples with detected circRNAs. CircView allows users to explore circRNAs detected by unique or multiple tools. Furthermore, CircView allows users to view the regulatory elements, such as microRNA response elements and RNA-binding protein binding sites. CircView is a unique tool to visualize and explore circRNAs, which helps users to better understand potential functions of circRNAs and design the functional experiments.

Genome-wide characterization of lncRNAs in acute myeloid leukemia.

Long noncoding RNAs (lncRNAs) are a large family of noncoding RNAs that play a critical role in various normal bioprocesses as well as tumorigenesis. However, the expression patterns and biological functions of lncRNAs in acute leukemia have not been well studied. Here, we performed transcriptome-wide lncRNA expression profiling of acute myeloid leukemia (AML) patient samples, along with non-leukemia control hematopoietic samples. We found that lncRNAs were differentially expressed in AML samples relative to control samples. Notably, we identified that lncRNAs upregulated in AML (relative to the control samples) are associated with a lower degree of DNA methylation and a higher ratio of being bound by transcription factors such as SP1, STAT4, ATF-2 and ELK-1 compared with those downregulated in AML. Moreover, an enrichment of H3K4me3 and a depletion of H3K27me3 were observed in upregulated lncRNAs in AML. Expression patterns of three types of lncRNAs (antisense, enhancer and intergenic lncRNAs) have previously been characterized. Of the identified lncRNAs, we found that high expression level lncRNA LOC285758 is associated with the poor prognosis in AML patients. Furthermore, we found that LOC285758 regulates proliferation of AML cell lines by enhancing the expression of HDAC2, a key factor in carcinogenesis. Collectively, our study depicts a landscape of important lncRNAs in AML and provides novel potential therapeutic targets and prognostic markers for AML treatment.

CSCD: a database for cancer-specific circular RNAs.

Circular RNA (circRNA) is a large group of RNA family extensively existed in cells and tissues. High-throughput sequencing provides a way to view circRNAs across different samples, especially in various diseases. However, there is still no comprehensive database for exploring the cancer-specific circRNAs. We collected 228 total RNA or polyA(-) RNA-seq samples from both cancer and normal cell lines, and identified 272 152 cancer-specific circRNAs. A total of 950 962 circRNAs were identified in normal samples only, and 170 909 circRNAs were identified in both tumor and normal samples, which could be further used as non-tumor background. We constructed a cancer-specific circRNA database (CSCD, http://gb.whu.edu.cn/CSCD). To understand the functional effects of circRNAs, we predicted the microRNA response element sites and RNA binding protein sites for each circRNA. We further predicted potential open reading frames to highlight translatable circRNAs. To understand the association between the linear splicing and the back-splicing, we also predicted the splicing events in linear transcripts of each circRNA. As the first comprehensive cancer-specific circRNA database, we believe CSCD could significantly contribute to the research for the function and regulation of cancer-associated circRNAs.