Upstream Splice Research Articles

State-of-the-Art Circular RNA Analytics Using the Circtools Software Suite.

Circular RNAs (circRNAs) are types of RNA molecules that have been discovered relatively recently and have been found to be widely expressed in eukaryotic cells. Unlike canonical linear RNA molecules, circRNAs form a covalently closed continuous loop structure without a 5' or 3' end. They are generated by a process called back-splicing, in which a downstream splice donor site is joined to an upstream splice acceptor site. CircRNAs have been found to play important roles in various biological processes, including gene regulation, alternative splicing, and protein translation. They can act as sponges for microRNAs or RNA-binding proteins and can also encode peptides or proteins. Additionally, circRNAs have been implicated in several diseases, including cancer, neurological disorders, and cardiovascular diseases.This protocol provides all necessary steps to detect and analyze circRNAs in silico from RNA sequencing data using the circtools circRNA analytics software suite. The protocol starts from raw sequencing data with circRNA detection via back-splice events and includes statistical testing of circRNAs as well as primer design for follow-up wet lab experiments.

A recurrent de novo splice site variant involving DNM1 exon 10a causes developmental and epileptic encephalopathy through a dominant-negative mechanism

Heterozygous pathogenic variants in DNM1 cause developmental and epileptic encephalopathy (DEE) as a result of a dominant-negative mechanism impeding vesicular fission. Thus far, pathogenic variants in DNM1 have been studied with a canonical transcript that includes the alternatively spliced exon 10b. However, after performing RNA sequencing in 39 pediatric brain samples, we find the primary transcript expressed in the brain includes the downstream exon 10a instead. Using this information, we evaluated genotype-phenotype correlations of variants affecting exon 10a and identified a cohort of eleven previously unreported individuals. Eight individuals harbor a recurrent de novo splice site variant, c.1197-8G>A (GenBank: NM_001288739.1), which affects exon 10a and leads to DEE consistent with the classical DNM1 phenotype. We find this splice site variant leads to disease through an unexpected dominant-negative mechanism. Functional testing reveals an in-frame upstream splice acceptor causing insertion of two amino acids predicted to impair oligomerization-dependent activity. This is supported by neuropathological samples showing accumulation of enlarged synaptic vesicles adherent to the plasma membrane consistent with impaired vesicular fission. Two additional individuals with missense variants affecting exon 10a, p.Arg399Trp and p.Gly401Asp, had a similar DEE phenotype. In contrast, one individual with a missense variant affecting exon 10b, p.Pro405Leu, which is less expressed in the brain, had a correspondingly less severe presentation. Thus, we implicate variants affecting exon 10a as causing the severe DEE typically associated with DNM1-related disorders. We highlight the importance of considering relevant isoforms for disease-causing variants as well as the possibility of splice site variants acting through a dominant-negative mechanism.

Zinc supplementation prevents arsenic-induced dysregulation of ZRANB2 splice function

Environmentally relevant (100 nM) inorganic arsenic (iAs) exposure displaces zinc from zinc fingers of upstream splice regulator ZRANB2 disrupting the splicing of its target TRA2B. Excess zinc displaced iAs from ZRANB2 zinc fingers in cell free system. Thus, the hypothesis that zinc supplementation could prevent iAs-mediated disruption of ZRANB2 splice function in human keratinocytes was tested. The data show that zinc supplementation prevented iAs-induced dysregulation of TRA2B splicing by ZRANB2 as well as the induction of ZRANB2 protein expression. These results provide additional support for the hypothesis that zinc supplementation could prevent iAs-mediated disease in iAs-exposed populations.

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.

Comprehensive comparison of two types of algorithm for circRNA detection from short-read RNA-Seq.

Circular RNA is generally formed by the 'back-splicing' process between the upstream splice acceptor and the downstream donor in/not in the regulation of the corresponding RNA-binding proteins or cis-elements. Therefore, more and more software packages have been developed and they are mostly based on the identification of the back-spliced junction reads. However, recent studies developed two software tools that can detect circRNA candidates by constructing k-mer table or/and de Bruijn graph rather than reads mapping. Here, we compared the precision, sensitivity and detection efficiency between software tools based on different algorithms. Eleven representative detection tools with two types of algorithm were selected for the overall pipeline analysis of RNA-seq datasets with/without RNase R treatment in two cell lines. Precision, sensitivity, AUC, F1 score and detection efficiency metrics were assessed to compare the prediction tools. Meanwhile, the sensitivity and distribution of highly expressed circRNAs before and after RNase R treatment were also revealed by their enrichment, unaffected and depleted candidate frequencies. Eventually, we found that compared to the k-mer based tools, CIRI2 and KNIFE based on reads mapping had relatively superior and more balanced detection performance regardless of the cell line or RNase R (-/+) datasets. All predicted results and source codes can be retrieved from https://github.com/luffy563/circRNA_tools_comparison. Supplementary data are available at Bioinformatics online.

Nucleotides in both donor and acceptor splice sites are responsible for choice in NAGNAG tandem splice sites.

Among alternative splicing events in the human transcriptome, tandem NAGNAG acceptor splice sites represent an appreciable proportion. Both proximal and distal NAG can be used to produce two splicing isoforms differing by three nucleotides. In some cases, the upstream exon can be alternatively spliced as well, which further increases the number of possible transcripts. In this study, we showed that NAG choice in tandem splice site depends considerably not only on the concerned acceptor, but also on the upstream donor splice site sequence. Using an extensive set of experiments with systematically modified two-exonic minigene systems of AFAP1L2 or CSTD gene, we recognized the third and fifth intronic upstream donor splice site position and the tandem acceptor splice site region spanning from -10 to +2, including NAGNAG itself, as the main drivers. In addition, competition between different branch points and their composition were also shown to play a significant role in NAG choice. All these nucleotide effects appeared almost additive, which explained the high variability in proximal versus distal NAG usage.

The genome position of a therapeutic transgene strongly influences the level of expression in an armed oncolytic human adenovirus vector

Efficacy of oncolytic, conditionally-replicating adenovirus (CRAd) vectors can be enhanced by “arming” the vector with therapeutic transgenes. We examined whether inclusion of an intact early region 3 (E3) and the reptilian reovirus fusogenic p14 fusion-associated small transmembrane (FAST) protein enhanced vector efficacy. The p14 FAST transgene was cloned between the fiber gene and E4 region, with an upstream splice acceptor for replication-dependent expression from the major late promoter. In A549 cells, this vector expressed p14 FAST protein at very low levels, and showed a poor ability to mediate cell-cell fusion, relative to a similar vector encoding p14 FAST within the E3 deletion. Although expression of E3 proteins from the CRAd increased plaque size, poor expression of p14 FAST protein compromised the fusogenic capacity of the vector. Thus, location of a therapeutic transgene within a CRAd can significantly impact expression of the transgene and is an important consideration in vector design.

An intronic splice site alteration in combination with a large deletion affecting VPS13B (COH1) causes Cohen syndrome

Cohen syndrome (CS) is a rare, autosomal recessive disorder characterized by intellectual disability, postnatal microcephaly, facial abnormalities, abnormal truncal fat distribution, myopia, and pigmentary retinopathy. It is often considered an underdiagnosed condition, especially in children with developmental delay and intellectual disability. Here we report on four individuals from a large Jordanian family clinically diagnosed with CS. Using Trio Exome Sequencing (Trio-WES) and MLPA analyses we identified a maternally inherited novel intronic nucleotide substitution c.3446-23T>G leading to the activation of a cryptic splice site and a paternally inherited multi-exon deletion in VPS13B (previously termed COH1) in the index patient. Expression analysis showed a strong decrease of VPS13B mRNA levels and direct sequencing of cDNA confirmed splicing at a cryptic upstream splice acceptor site, resulting in the inclusion of 22 intronic bases. This extension results in a frameshift and a premature stop of translation (p.Gly1149Valfs*9). Segregation analysis revealed that three affected maternal cousins were homozygous for the intronic splice site variant. Our data show causality of both alterations and strongly suggest the expansion of the diagnostic strategy to search for intronic splice variants in molecularly unconfirmed patients affected by CS.

Oxaliplatin-Induced DHX9 Phosphorylation Promotes Oncogenic Circular RNA CCDC66 Expression and Development of Chemoresistance.

Circular RNA (circRNA), generated through backsplicing in which the downstream splice donor joins the upstream splice acceptor, is a novel class of RNA molecules. Our previous study found that a novel oncogenic circRNA—consisting exon 8–10 of CCDC66—is aberrantly expressed in colorectal cancer (CRC) tissues and cells. The failure of treatment for colorectal cancer is typically associated with recurrent and chemoresistant cancerous tissues. In this study, we aimed to investigate the role(s) of circCCDC66 during the development of chemoresistance. We discovered that the expression level of circCCDC66 is elevated in colorectal cancer cells with resistance to oxaliplatin. Knockdown of circCCDC66 caused the downregulation of a subset of genes which are regulated by circCCDC66-associated miRNAs and related to the modulation of apoptosis and the cell cycle, suppressing cell survival, promoting oxaliplatin-induced apoptosis and, thus, hindering the development of oxaliplatin-resistance (OxR). The induction of circCCDC66 was dependent on the time-course and dose of oxaliplatin treatment. Our analyses revealed that DHX9 harbors two phosphorylation sites of phosphatidylinositol 3-kinase-related kinases (PI3KKs) close to substrate-binding domains. Blockage of phosphorylation by either PI3KK inhibitors or nonphosphorable mutants of DHX9 decreased the oxaliplatin-induced circCCDC66 expression and the ability to develop chemoresistant cells. Taken together, we demonstrated and linked the functional role of DHX9 phosphorylation to oncogenic circCCDC66 expression during the development of resistance to oxaliplatin, providing a mechanistic insight for the development of therapeutic strategies to recurring/chemoresistant colorectal cancer.

Kaposi's Sarcoma-Associated Herpesvirus-Encoded circRNAs Are Expressed in Infected Tumor Tissues and Are Incorporated into Virions.

Kaposi's sarcoma-associated herpesvirus (KSHV) has recently been found to generate circular RNAs (circRNAs) from several KSHV genes, most abundantly from K10 (viral interferon regulatory factor 4 [vIRF4]), K7.3, and polyadenylated nuclear (PAN) RNA. To define expression of these circRNAs, KSHV-infected cell lines, patient tissues, and purified virions were examined. KSHV circRNA expression was universally detected in tests of six primary effusion lymphoma (PEL) cell lines but ranged from low-level expression in BC-1 cells dually infected with tightly latent KSHV and Epstein-Barr virus to abundant expression in KSHV-only BCBL-1 cells with spontaneous virus production. Generally, the PAN/K7.3 locus broadly and bidirectionally generated circRNA levels that paralleled the corresponding linear RNA levels. However, RNA corresponding to a particular KSHV circularization site (circ-vIRF4) was minimally induced, despite linear vIRF4 RNA being activated by virus induction. In situ hybridization showed abundant circ-vIRF4 in noninduced PEL cells. All three KSHV circRNAs were isolated as nuclease-protected forms from gradient-purified virions collected from BrK.219 cells infected with a KSHV molecular clone. For circ-vIRF4, the fully processed form that is exported to the cytoplasm was incorporated into virus particles but the nuclear, intron-retaining form was not. The half-life of circ-vIRF4 was twice as long as that of its linear counterpart. The KSHV circRNAs could be detected at a higher rate than their corresponding linear counterparts by in situ hybridization in archival tissues and by reverse transcription-PCR (RT-PCR) in sera stored for over 25 years. In summary, KSHV circRNAs are expressed in infection-associated diseases, can be regulated depending on virus life cycle, and are incorporated into viral particles for preformed delivery, suggesting a potential function in early infection.IMPORTANCE KSHV has recently been found to encode circRNAs. circRNAs result from back-splicing of an upstream pre-mRNA splice donor exon-intron junction to an acceptor site, generating a covalently closed circle. This study revealed that for one KSHV region, the PAN/K7.3 locus, broadly and bidirectionally generated circRNA levels parallel corresponding linear RNA levels. Another KSHV circularization site (circ-vIRF4), however, showed expression that differed from that of the corresponding linear RNA. All KSHV circRNAs are incorporated into KSHV virions and are potentially expressed as immediate early products in newly infected cells.

SpliceV: analysis and publication quality printing of linear and circular RNA splicing, expression and regulation

BackgroundIn eukaryotes, most genes code for multiple transcript isoforms that are generated through the complex and tightly regulated process of RNA splicing. Despite arising from identical precursor transcripts, alternatively spliced RNAs can have dramatically different functions. Transcriptome complexity is elevated further by the production of circular RNAs (circRNAs), another class of mature RNA that results from the splicing of a downstream splice donor to an upstream splice acceptor. While there has been a rapid expansion of circRNA catalogs in the last few years through the utilization of next generation sequencing approaches, our understanding of the mechanisms and regulation of circular RNA biogenesis, the impact that circRNA generation has on parental transcript processing, and the functions carried out by circular RNAs remains limited.ResultsHere, we present a visualization and analysis tool, SpliceV, that rapidly plots all relevant forward- and back-splice data, with exon and single nucleotide level coverage information from RNA-seq experiments in a publication quality format. SpliceV also integrates analysis features that assist investigations into splicing regulation and transcript functions through the display of predicted RNA binding protein sites and the configuration of repetitive elements along the primary transcript.ConclusionsSpliceV is an easy-to-use splicing visualization tool, compatible with both Python 2.7 and 3+, and distributed under the GNU Public License. The source code is freely available for download at https://github.com/flemingtonlab/SpliceV and can be installed from PyPI using pip.

Efficient system for upstream mRNA trans-splicing to generate covalent, head-to-tail, protein multimers

We present a plasmid-based system in which upstream trans-splicing efficiently generates mRNAs that encode head-to-tail protein multimers. In this system, trans-splicing occurs between one of two downstream splice donors in the sequence encoding a C-terminal V5 epitope tag and an upstream splice acceptor in the 5′ region of the pCS2(+) host plasmid. Using deletion and fusion constructs of the DUX4 protein as an example, we found that this system produced trans-spliced mRNAs in which coding regions from independent transcripts were fused in phase such that covalent head-to-tail protein multimers were translated. For a cDNA of ~450 bp, about half of the expressed proteins were multimeric, with the efficiency of trans-splicing and extent of multimer expression decreasing as cDNA length increased. This system generated covalent heterodimeric proteins upon co-transfections of plasmids encoding separate proteins and did not require a long complementary binding domain to position mRNAs for trans-splicing. This plasmid-based trans-splicing system is adaptable to multiple gene delivery systems, and it presents new opportunities for investigating molecular mechanisms of trans-splicing, generating covalent protein multimers with novel functions within cells, and producing mRNAs encoding large proteins from split precursors.

Dual Mechanisms of LYN Kinase Dysregulation Drive Aggressive Behavior in Breast Cancer Cells

SummaryThe SRC-family kinase LYN is highly expressed in triple-negative/basal-like breast cancer (TNBC) and in the cell of origin of these tumors, c-KIT-positive luminal progenitors. Here, we demonstrate LYN is a downstream effector of c-KIT in normal mammary cells and protective of apoptosis upon genotoxic stress. LYN activity is modulated by PIN1, a prolyl isomerase, and in BRCA1 mutant TNBC PIN1 upregulation activates LYN independently of c-KIT. Furthermore, the full-length LYN splice isoform (as opposed to the Δaa25–45 variant) drives migration and invasion of aggressive TNBC cells, while the ratio of splice variants is informative for breast cancer-specific survival across all breast cancers. Thus, dual mechanisms—uncoupling from upstream signals and splice isoform ratios—drive the activity of LYN in aggressive breast cancers.

A comprehensive overview and evaluation of circular RNA detection tools.

Circular RNA (circRNA) is mainly generated by the splice donor of a downstream exon joining to an upstream splice acceptor, a phenomenon known as backsplicing. It has been reported that circRNA can function as microRNA (miRNA) sponges, transcriptional regulators, or potential biomarkers. The availability of massive non-polyadenylated transcriptomes data has facilitated the genome-wide identification of thousands of circRNAs. Several circRNA detection tools or pipelines have recently been developed, and it is essential to provide useful guidelines on these pipelines for users, including a comprehensive and unbiased comparison. Here, we provide an improved and easy-to-use circRNA read simulator that can produce mimicking backsplicing reads supporting circRNAs deposited in CircBase. Moreover, we compared the performance of 11 circRNA detection tools on both simulated and real datasets. We assessed their performance regarding metrics such as precision, sensitivity, F1 score, and Area under Curve. It is concluded that no single method dominated on all of these metrics. Among all of the state-of-the-art tools, CIRI, CIRCexplorer, and KNIFE, which achieved better balanced performance between their precision and sensitivity, compared favorably to the other methods.

Monoallelic Large Intragenic KRT5 Deletions Account for Genetically Unsolved Cases of Epidermolysis Bullosa Simplex

Monoallelic Large Intragenic KRT5 Deletions Account for Genetically Unsolved Cases of Epidermolysis Bullosa Simplex

Circular RNAs: An emerging type of RNA in cancer.

Circular RNAs (circRNAs), a novel type of widespread and diverse endogenous non-coding RNAs (ncRNAs), which are different from the linear RNAs, form a covalently closed continuous loop without 5’ or 3’ polarities. The majority of circRNAs are abundant, conserved and stable across different species, and exhibit tissue/developmental-stage-specific characteristics. They are generated primarily through a type of alternative RNA splicing called “back-splicing,” in which a downstream splice donor is joined to an upstream splice acceptor through splice skipping or direct splice. Recent studies have discovered circRNAs function as microRNA sponges, binding with RNA-associated proteins to form RNA–protein complexes and then regulating gene transcription and translation into polypeptides. Emerging evidence indicates that circRNAs play important roles in the regulation of the development and progression of multiple cancers by serving as potential diagnostic and predictive biomarkers involved in tumor growth and invasion and providing new strategies for cancer diagnosis and targeted therapy. In this review, we briefly delineate the diversity and characteristics of circRNAs and discuss the highlights of the biogenesis of circRNAs and their potential functions in tumor.

The Research on Identification of Gene Splice Sites by Support Vector Machine

The recognition of splicing sites is a very important step in the eukaryotic DNA se-quence analysis. Many scholars are working hard to improve the accuracy of identifi-cation. Our team carried out research on this issue based on support vector machine, which is one famous algorithm in data mining. The training and testing data is from the HS3D dataset, and excellent accuracy rate is achieved by nucleic acid sequence orthogonal coding and RBF core function, and the cross validation experiment hints that base pattern information is mainly located within 20 nucleotides upstream and downstream splice sites.

Evolutionarily conserved exon definition interactions with U11 snRNP mediate alternative splicing regulation on U11–48K and U11/U12–65K genes

Many splicing regulators bind to their own pre-mRNAs to induce alternative splicing that leads to formation of unstable mRNA isoforms. This provides an autoregulatory feedback mechanism that regulates the cellular homeostasis of these factors. We have described such an autoregulatory mechanism for two core protein components, U11–48K and U11/U12–65K, of the U12-dependent spliceosome. This regulatory system uses an atypical splicing enhancer element termed USSE (U11 snRNP-binding splicing enhancer), which contains two U12-type consensus 5′ splice sites (5′ss). Evolutionary analysis of the USSE element from a large number of animal and plant species indicate that USSE sequence must be located 25–50 nt downstream from the target 3′ splice site (3′ss). Together with functional evidence showing a loss of USSE activity when this distance is reduced and a requirement for RS-domain of U11–35K protein for 3′ss activation, our data suggests that U11 snRNP bound to USSE uses exon definition interactions for regulating alternative splicing. However, unlike standard exon definition where the 5′ss bound by U1 or U11 will be subsequently activated for splicing, the USSE element functions similarly as an exonic splicing enhancer and is involved only in upstream splice site activation but does not function as a splicing donor. Additionally, our evolutionary and functional data suggests that the function of the 5′ss duplication within the USSE elements is to allow binding of two U11/U12 di-snRNPs that stabilize each others' binding through putative mutual interactions.

Coupling of Transcription with mRNA Processing in time and Space

The cellular transcription and mRNA processing machineries appear to have co-evolved to allow coupling of the reactions they perform in space and time. Coupling in space is achieved by “recruitment” of processing factors to the nascent transcript, often by binding to the pol II CTD. Much less is known about how coupling in time is achieved, that is to say the coordination of rates of elongation and processing. Slow elongation is predicted to lengthen the window of opportunity for an event to occur at an upstream site on the nascent transcript before it must compete with an RNA sequence element further downstream. Consistent with this prediction, in a few cases slow elongation is associated with enhanced inclusion of alternative exons (ref. 1) through processing at upstream splice sites and processing at upstream alternative poly(A) sites (ref. 2). We investigated kinetic coupling in human cell lines that express pol II large subunits with mutations in funnel and trigger loop domains that slow down and speed up elongation rates genome-wide. Slow and fast elongation affected constitutive and alternative splicing, frequently altering exon inclusion and intron retention in ways not predicted by the “window of opportunity” model. Unexpectedly, slow and fast elongation often both increased or both decreased inclusion of a particular exon or retained intron. Together our results suggest that an optimal rate of transcriptional elongation is required for normal co-transcriptional pre-mRNA splicing.

Pre-mRNA splicing is facilitated by an optimal RNA polymerase II elongation rate.

Alternative splicing modulates expression of most human genes. The kinetic model of cotranscriptional splicing suggests that slow elongation expands and that fast elongation compresses the "window of opportunity" for recognition of upstream splice sites, thereby increasing or decreasing inclusion of alternative exons. We tested the model using RNA polymerase II mutants that change average elongation rates genome-wide. Slow and fast elongation affected constitutive and alternative splicing, frequently altering exon inclusion and intron retention in ways not predicted by the model. Cassette exons included by slow and excluded by fast elongation (type I) have weaker splice sites, shorter flanking introns, and distinct sequence motifs relative to "slow-excluded" and "fast-included" exons (type II). Many rate-sensitive exons are misspliced in tumors. Unexpectedly, slow and fast elongation often both increased or both decreased inclusion of a particular exon or retained intron. These results suggest that an optimal rate of transcriptional elongation is required for normal cotranscriptional pre-mRNA splicing.