Assessing long-distance RNA sequence connectivity via RNA-templated DNA-DNA ligation

Computational approaches and challenges in the analysis of circRNA data

Circular RNAs (circRNA) are a class of non-coding RNA, forming a single-stranded covalently closed loop structure generated via back-splicing. Advancements in sequencing methods and technologies in conjunction with algorithmic developments of bioinformatics tools have enabled researchers to characterise the origin and function of circRNAs, with practical applications as a biomarker of diseases becoming increasingly relevant. Computational methods developed for circRNA analysis are predicated on detecting the chimeric back-splice junction of circRNAs whilst mitigating false-positive sequencing artefacts. In this review, we discuss in detail the computational strategies developed for circRNA identification, highlighting a selection of tool strengths, weaknesses and assumptions. In addition to circRNA identification tools, we describe methods for characterising the role of circRNAs within the competing endogenous RNA (ceRNA) network, their interactions with RNA-binding proteins, and publicly available databases for rich circRNA annotation.

Droplet Digital PCR for the Detection and Quantification of Bona Fide CircRNAs

CircRNAs are covalently closed RNA molecules gaining increasing attention over the years. Initially considered mere splicing errors, circRNAs are now recognized as a novel class of endogenous, conserved RNAs, expressed in many different species. The unique structure, the low levels of expression, and the almost complete sequence overlap with the cognate linear RNA make their detection and quantification challenging. Moreover, it has become crucial to prove the circular nature of the targeted transcript and unequivocally distinguish the circRNA from its linear counterpart. Nowadays, the most widely used technique to quantify circRNA expression is real-time quantitative PCR (qPCR). However, in the particular case of quantification of circles, it shows several technical shortcomings which affect the accuracy of the quantification. To precisely assess circRNA expression level, droplet digital PCR (ddPCR) is rapidly taking over for the more popular qPCR. In this chapter, we describe the detailed procedure based on droplets partitioning to quantify both linear and circRNA abundancy and demonstrate the circularity of the transcript under study with high precision, in a single experiment.

Recent advances in cancer fusion transcript detection

Extensive investigation of gene fusions in cancer has led to the discovery of novel biomarkers and therapeutic targets. To date, most studies have neglected chromosomal rearrangement-independent fusion transcripts and complex fusion structures such as double or triple-hop fusions, and fusion-circRNAs. In this review, we untangle fusion-related terminology and propose a classification system involving both gene and transcript fusions. We highlight the importance of RNA-level fusions and how long-read sequencing approaches can improve detection and characterization. Moreover, we discuss novel bioinformatic tools to identify fusions in long-read sequencing data and strategies to experimentally validate and functionally characterize fusion transcripts.

The Use of circRNAs as Biomarkers of Cancer

CircRNAs are a subclass of lncRNAs that have been found to be abundantly present in a wide range of species, including humans. CircRNAs are generally produced by a noncanonical splicing event called backsplicing that is dependent on the canonical splicing machinery, giving rise to circRNAs classified into three main categories: exonic circRNA, circular intronic RNA, and exon-intron circular RNA. Notably, circRNAs possess functional importance and display their functions through different mechanisms of action including sponging miRNAs, or even being translated into functional proteins. In addition, circRNAs also have great potential as biomarkers, particularly in cancer, thanks to their high stability, tissue type and developmental stage specificity, and their presence in biological fluids, which make them promising candidates as noninvasive biomarkers. In this chapter, we describe the most commonly used techniques for the study of circRNAs as cancer biomarkers, including high-throughput techniques such as RNA-Seq and microarrays, and other methods to analyze the presence of specific circRNAs in patient samples.

Intron-targeted mutagenesis reveals roles for Dscam1 RNA pairing architecture-driven splicing bias in neuronal wiring

Drosophila melanogaster Down syndrome cell adhesion molecule (Dscam1) can generate 38,016 different isoforms through largely stochastic, yet highly biased, alternative splicing. These isoforms are required for nervous functions. However, the functional significance of splicing bias remains unknown. Here, we provide evidence that Dscam1 splicing bias is required for mushroom body (MB) axonal wiring. We generate mutant flies with normal overall protein levels and an identical number but global changes in exon 4 and 9 isoform bias (DscamΔ4D^-/- and DscamΔ9D^-/-), respectively. In contrast to DscamΔ4D^-/-, DscamΔ9D^-/- exhibits remarkable MB defects, suggesting a variable domain-specific requirement for isoform bias. Importantly, changes in isoform bias cause axonal defects but do not influence the self-avoidance of axonal branches. We conclude that, in contrast to the isoform number that provides the molecular basis for neurite self-avoidance, isoform bias may play a role in MB axonal wiring by influencing non-repulsive signaling.

A spatially resolved brain region- and cell type-specific isoform atlas of the postnatal mouse brain

Splicing varies across brain regions, but the single-cell resolution of regional variation is unclear. We present a single-cell investigation of differential isoform expression (DIE) between brain regions using single-cell long-read sequencing in mouse hippocampus and prefrontal cortex in 45 cell types at postnatal day 7 ( www.isoformAtlas.com ). Isoform tests for DIE show better performance than exon tests. We detect hundreds of DIE events traceable to cell types, often corresponding to functionally distinct protein isoforms. Mostly, one cell type is responsible for brain-region specific DIE. However, for fewer genes, multiple cell types influence DIE. Thus, regional identity can, although rarely, override cell-type specificity. Cell types indigenous to one anatomic structure display distinctive DIE, e.g. the choroid plexus epithelium manifests distinct transcription-start-site usage. Spatial transcriptomics and long-read sequencing yield a spatially resolved splicing map. Our methods quantify isoform expression with cell-type and spatial resolution and it contributes to further our understanding of how the brain integrates molecular and cellular complexity.

Circular RNAs in Sepsis: Biogenesis, Function, and Clinical Significance

Sepsis is a life-threatening condition that occurs when the body responds to an infection that damages it is own tissues. The major problem in sepsis is rapid, vital status deterioration in patients, which can progress to septic shock with multiple organ failure if not properly treated. As there are no specific treatments, early diagnosis is mandatory to reduce high mortality. Despite more than 170 different biomarkers being postulated, early sepsis diagnosis and prognosis remain a challenge for clinicians. Recent findings propose that circular RNAs (circRNAs) may play a prominent role in regulating the patients’ immune system against different pathogens, including bacteria and viruses. Mounting evidence also suggests that the misregulation of circRNAs is an early event in a wide range of diseases, including sepsis. Despite circRNA levels being altered in sepsis, the specific mechanisms controlling the dysregulation of these noncoding RNAs are not completely elucidated, although many factors are known to affect circRNA biogenesis. Therefore, there is a need to explore the molecular pathways that lead to this disorder. This review describes the role of this new class of regulatory RNAs in sepsis and the feasibility of using circRNAs as diagnostic biomarkers for sepsis, opening up new avenues for circRNA-based medicine.

Plasmid-based gap-repair recombineered transgenes reveal a central role for introns in mutually exclusive alternative splicing in Down Syndrome Cell Adhesion Molecule exon 4

Alternative splicing is a key feature of human genes, yet studying its regulation is often complicated by large introns. The Down Syndrome Cell Adhesion Molecule (Dscam) gene from Drosophila is one of the most complex genes generating vast molecular diversity by mutually exclusive alternative splicing. To resolve how alternative splicing in Dscam is regulated, we first developed plasmid-based UAS reporter genes for the Dscam variable exon 4 cluster and show that its alternative splicing is recapitulated by GAL4-mediated expression in neurons. We then developed gap-repair recombineering to very efficiently manipulate these large reporter plasmids in Escherichia coli using restriction enzymes or sgRNA/Cas9 DNA scission to capitalize on the many benefits of plasmids in phiC31 integrase-mediated transgenesis. Using these novel tools, we show that inclusion of Dscam exon 4 variables differs little in development and individual flies, and is robustly determined by sequences harbored in variable exons. We further show that introns drive selection of both proximal and distal variable exons. Since exon 4 cluster introns lack conserved sequences that could mediate robust long-range base-pairing to bring exons into proximity for splicing, our data argue for a central role of introns in mutually exclusive alternative splicing of Dscam exon 4 cluster.

Getting the Entire Message: Progress in Isoform Sequencing

The advent of second-generation sequencing and its application to RNA sequencing have revolutionized the field of genomics by allowing quantification of gene expression, as well as the definition of transcription start/end sites, exons, splice sites and RNA editing sites. However, due to the sequencing of fragments of cDNAs, these methods have not given a reliable picture of complete RNA isoforms. Third-generation sequencing has filled this gap and allows end-to-end sequencing of entire RNA/cDNA molecules. This approach to transcriptomics has been a "niche" technology for a couple of years but now is becoming mainstream with many different applications. Here, we review the background and progress made to date in this rapidly growing field. We start by reviewing the progressive realization that alternative splicing is omnipresent. We then focus on long-noncoding RNA isoforms and the distinct combination patterns of exons in noncoding and coding genes. We consider the implications of the recent technologies of direct RNA sequencing and single-cell isoform RNA sequencing. Finally, we discuss the parameters that define the success of long-read RNA sequencing experiments and strategies commonly used to make the most of such data.

CircRNAs and cancer: Biomarkers and master regulators

Circular RNAs (circRNAs) are a novel class of regulatory RNAs that despite being relatively abundant have only recently begun to be explored. There are many thousands of genes that appear capable of producing circRNAs, however the function of all but a handful remain to be determined. What is emerging about these highly conserved molecules is that they play important roles in biology and cancer biology in particular. The most explored function of circRNAs is as master regulators of gene expression that act to sequester or ´sponge´ other gene expression regulators, in particular miRNAs. They have also been demonstrated to function via direct modulation of transcription, and by interfering with splicing mechanisms. Although generally expressed in low abundance when compared to their linear counterparts, they are often expressed in a tissue- and developmental stage- specific manner. Coupled with their remarkable resistance to RNAse activity due to a covalent closed cyclic structure, circRNAs show great promise as novel biomarkers of cancer and other diseases. In this review we consider the current state of knowledge regarding these molecules, their synthesis, function, and association with cancer. We will also review some of the challenges that remain to be resolved if this emerging class of RNAs are really to become useful in the clinic.

Single-cell isoform RNA sequencing characterizes isoforms in thousands of cerebellar cells

Full-length RNA sequencing (RNA-Seq) has been applied to bulk tissue, cell lines and sorted cells to characterize transcriptomes, but applying this technology to single cells has proven to be difficult, with less than ten single-cell transcriptomes having been analyzed thus far. Although single splicing events have been described for ≤200 single cells with statistical confidence, full-length mRNA analyses for hundreds of cells have not been reported. Single-cell short-read 3' sequencing enables the identification of cellular subtypes, but full-length mRNA isoforms for these cell types cannot be profiled. We developed a method that starts with bulk tissue and identifies single-cell types and their full-length RNA isoforms without fluorescence-activated cell sorting. Using single-cell isoform RNA-Seq (ScISOr-Seq), we identified RNA isoforms in neurons, astrocytes, microglia, and cell subtypes such as Purkinje and Granule cells, and cell-type-specific combination patterns of distant splice sites. We used ScISOr-Seq to improve genome annotation in mouse Gencode version 10 by determining the cell-type-specific expression of 18,173 known and 16,872 novel isoforms.

Molecular functions and specific roles of circRNAs in the cardiovascular system

As part of the superfamily of long noncoding RNAs, circular RNAs (circRNAs) are emerging as a new type of regulatory molecules that partake in gene expression control. Here, we review the current knowledge about circRNAs in cardiovascular disease. CircRNAs are not only associated with different types of cardiovascular disease, but they have also been identified as intracellular effector molecules for pathophysiological changes in cardiovascular tissues, and as cardiovascular biomarkers. This evidence is put in the context of the current understanding of general circRNA biogenesis and of known interactions of circRNAs with DNA, RNA, and proteins.

Molecular roles and function of circular RNAs in eukaryotic cells

Protein-coding and noncoding genes in eukaryotes are typically expressed as linear messenger RNAs, with exons arranged colinearly to their genomic order. Recent advances in sequencing and in mapping RNA reads to reference genomes have revealed that thousands of genes express also covalently closed circular RNAs. Many of these circRNAs are stable and contain exons, but are not translated into proteins. Here, we review the emerging understanding that both, circRNAs produced by co- and posttranscriptional head-to-tail "backsplicing" of a downstream splice donor to a more upstream splice acceptor, as well as circRNAs generated from intronic lariats during colinear splicing, may exhibit physiologically relevant regulatory functions in eukaryotes. We describe how circRNAs impact gene expression of their host gene locus by affecting transcriptional initiation and elongation or splicing, and how they partake in controlling the function of other molecules, for example by interacting with microRNAs and proteins. We conclude with an outlook how circRNA dysregulation affects disease, and how the stability of circRNAs might be exploited in biomedical applications.

Microfluidic isoform sequencing shows widespread splicing coordination in the human transcriptome

Understanding transcriptome complexity is crucial for understanding human biology and disease. Technologies such as Synthetic long-read RNA sequencing (SLR-RNA-seq) delivered 5 million isoforms and allowed assessing splicing coordination. Pacific Biosciences and Oxford Nanopore increase throughput also but require high input amounts or amplification. Our new droplet-based method, sparse isoform sequencing (spISO-seq), sequences 100k–200k partitions of 10–200 molecules at a time, enabling analysis of 10–100 million RNA molecules. SpISO-seq requires less than 1 ng of input cDNA, limiting or removing the need for prior amplification with its associated biases. Adjusting the number of reads devoted to each molecule reduces sequencing lanes and cost, with little loss in detection power. The increased number of molecules expands our understanding of isoform complexity. In addition to confirming our previously published cases of splicing coordination (e.g., BIN1), the greater depth reveals many new cases, such as MAPT. Coordination of internal exons is found to be extensive among protein coding genes: 23.5%–59.3% (95% confidence interval) of highly expressed genes with distant alternative exons exhibit coordination, showcasing the need for long-read transcriptomics. However, coordination is less frequent for noncoding sequences, suggesting a larger role of splicing coordination in shaping proteins. Groups of genes with coordination are involved in protein–protein interactions with each other, raising the possibility that coordination facilitates complex formation and/or function. We also find new splicing coordination types, involving initial and terminal exons. Our results provide a more comprehensive understanding of the human transcriptome and a general, cost-effective method to analyze it.

Circular RNAs: analysis, expression and potential functions

Just a few years ago, it had been assumed that the dominant RNA isoforms produced from eukaryotic genes were variants of messenger RNA, functioning as intermediates in gene expression. In early 2012, however, a surprising discovery was made: circular RNA (circRNA) was shown to be a transcriptional product in thousands of human and mouse genes and in hundreds of cases constituted the dominant RNA isoform. Subsequent studies revealed that the expression of circRNAs is developmentally regulated, tissue and cell-type specific, and shared across the eukaryotic tree of life. These features suggest important functions for these molecules. Here, we describe major advances in the field of circRNA biology, focusing on the regulation of and functional roles played by these molecules.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Detecting circular RNAs: bioinformatic and experimental challenges

The pervasive expression of circular RNAs (circRNAs) is a recently discovered feature of gene expression in highly diverged eukaryotes. Numerous algorithms that are used to detect genome-wide circRNA expression from RNA sequencing (RNA-seq) data have been developed in the past few years, but there is little overlap in their predictions and no clear gold-standard method to assess the accuracy of these algorithms. We review sources of experimental and bioinformatic biases that complicate the accurate discovery of circRNAs and discuss statistical approaches to address these biases. We conclude with a discussion of the current experimental progress on the topic.

Sensitive and specific miRNA detection method using SplintR Ligase

We describe a simple, specific and sensitive microRNA (miRNA) detection method that utilizes Chlorella virus DNA ligase (SplintR^® Ligase). This two-step method involves ligation of adjacent DNA oligonucleotides hybridized to a miRNA followed by real-time quantitative PCR (qPCR). SplintR Ligase is 100X faster than either T4 DNA Ligase or T4 RNA Ligase 2 for RNA splinted DNA ligation. Only a 4–6 bp overlap between a DNA probe and miRNA was required for efficient ligation by SplintR Ligase. This property allows more flexibility in designing miRNA-specific ligation probes than methods that use reverse transcriptase for cDNA synthesis of miRNA. The qPCR SplintR ligation assay is sensitive; it can detect a few thousand molecules of miR-122. For miR-122 detection the SplintR qPCR assay, using a FAM labeled double quenched DNA probe, was at least 40× more sensitive than the TaqMan assay. The SplintR method, when coupled with NextGen sequencing, allowed multiplex detection of miRNAs from brain, kidney, testis and liver. The SplintR qPCR assay is specific; individual let-7 miRNAs that differ by one nucleotide are detected. The rapid kinetics and ability to ligate DNA probes hybridized to RNA with short complementary sequences makes SplintR Ligase a useful enzyme for miRNA detection.

New Analytical Tool for the Detection of Ractopamine Abuse in Goat Skeletal Muscle by Potential Gene Expression Biomarkers

In this study, quantification of mRNA gene expression was examined as biomarkers to detect ractopamaine abuse and ractopamaine residues in cashmere goats. It was focused on the identification of potential gene expression biomarkers and describing the coreletionship between gene expression and residue level by 58 animals for 49 days. The results showed that administration periods and residue levels significantly influenced mRNA expressions of the β2-adrenergic receptor (β2AR), the enzymes PRKACB, ADCY3, ATP1A3, ATP2A3, PTH, and MYLK, and the immune factors IL-1β and TNF-α. Statistical analysis like principal components analysis (PCA), hierarchical cluster analysis (HCA), and discriminant analysis (DA) showed that these genes can serve as potential biomarkers for ractopamine in skeletal muscle and that they are also suitable for describing different residue levels separately.

Determining exon connectivity in complex mRNAs by nanopore sequencing

Short-read high-throughput RNA sequencing, though powerful, is limited in its ability to directly measure exon connectivity in mRNAs that contain multiple alternative exons located farther apart than the maximum read length. Here, we use the Oxford Nanopore MinION sequencer to identify 7,899 ‘full-length’ isoforms expressed from four Drosophila genes, Dscam1, MRP, Mhc, and Rdl. These results demonstrate that nanopore sequencing can be used to deconvolute individual isoforms and that it has the potential to be a powerful method for comprehensive transcriptome characterization.