Dbr1 functions in mRNA processing, intron turnover and human diseases

Rapid exome sequencing for children with severe acute encephalopathy - A case series

Increasingly, next-generation sequencing (NGS) is becoming an invaluable tool in the diagnosis of unexplained acute neurological disorders, such as acute encephalopathy/encephalitis. Here, we describe a brief series of pediatric patients who presented at the pediatric intensive care unit with severe acute encephalopathy, initially suspected as infectious or inflammatory but subsequently diagnosed with a monogenic disorder. Rapid exome sequencing was performed during the initial hospitalization of three unrelated patients, and results were delivered within 7-21 days. All patients were previously healthy, 1.5-3 years old, of Muslim Arab descent, with consanguineous parents. One patient presenting with acute necrotizing encephalopathy (ANEC). Her sister presented with ANEC one year prior. Exome sequencing was diagnostic in all three patients. All were homozygous for pathogenic and likely-pathogenic variants associated with recessive disorders; MOCS2, NDUFS8 and DBR1. Surprisingly, the initial workup was not suggestive of the final diagnosis. This case series demonstrates that the use of rapid exome sequencing is shifting the paradigm of diagnostics even in critical care situations and should be considered early on in children with acute encephalopathy. A timely diagnosis can direct initial treatment as well as inform decisions regarding long-term care.

Trichothiodystrophy-associated MPLKIP maintains DBR1 levels for proper lariat debranching and ectodermal differentiation

The brittle hair syndrome Trichothiodystrophy (TTD) is characterized by variable clinical features, including photosensitivity, ichthyosis, growth retardation, microcephaly, intellectual disability, hypogonadism, and anaemia. TTD-associated mutations typically cause unstable mutant proteins involved in various steps of gene expression, severely reducing steady-state mutant protein levels. However, to date, no such link to instability of gene-expression factors for TTD-associated mutations in MPLKIP/TTDN1 has been established. Here, we present seven additional TTD individuals with MPLKIP mutations from five consanguineous families, with a newly identified MPLKIP variant in one family. By mass spectrometry-based interaction proteomics, we demonstrate that MPLKIP interacts with core splicing factors and the lariat debranching protein DBR1. MPLKIP-deficient primary fibroblasts have reduced steady-state DBR1 protein levels. Using Human Skin Equivalents (HSEs), we observed impaired keratinocyte differentiation associated with compromised splicing and eventually, an imbalanced proteome affecting skin development and, interestingly, also the immune system. Our data show that MPLKIP, through its DBR1 stabilizing role, is implicated in mRNA splicing, which is of particular importance in highly differentiated tissue.

Activation of human RNA lariat debranching enzyme Dbr1 by binding protein TTDN1 occurs though an intrinsically disordered C-terminal domain

In eukaryotic cells, the introns are excised from pre-mRNA by the spliceosome. These introns typically have a lariat configuration due to the 2'-5' phosphodiester bond between an internal branched residue and the 5' terminus of the RNA. The only enzyme known to selectively hydrolyze the 2'-5' linkage of these lariats is the RNA lariat debranching enzyme Dbr1. In humans, Dbr1 is involved in processes such as class-switch recombination of immunoglobulin genes, and its dysfunction is implicated in viral encephalitis, HIV, ALS, and cancer. However, mechanistic details of precisely how Dbr1 affects these processes are missing. Here we show that human Dbr1 contains a disordered C-terminal domain through sequence analysis and nuclear magnetic resonance. This domain stabilizes Dbr1 in vitro by reducing aggregation but is dispensable for debranching activity. We establish that Dbr1 requires Fe2+ for efficient catalysis and demonstrate that the noncatalytic protein Drn1 and the uncharacterized protein trichothiodystrophy nonphotosensitive 1 directly bind to Dbr1. We demonstrate addition of trichothiodystrophy nonphotosensitive 1 to in vitro debranching reactions increases the catalytic efficiency of human Dbr1 19-fold but has no effect on the activity of Dbr1 from the amoeba Entamoeba histolytica, which lacks a disordered C-terminal domain. Finally, we systematically examine how the identity of the branchpoint nucleotide affects debranching rates. These findings describe new aspects of Dbr1 function in humans and further clarify how Dbr1 contributes to human health and disease.

Circular stable intronic RNAs possess distinct biological features and are deregulated in bladder cancer

Until recently, intronic lariats were regarded as short-lasting splicing byproducts with no apparent function; however, increasing evidence of stable derivatives suggests regulatory roles. Yet little is known about their characteristics, functions, distribution, and expression in healthy and tumor tissue. Here, we profiled and characterized circular stable intronic sequence RNAs (sisRNAs) using total RNA-Seq data from bladder cancer (BC; n = 457, UROMOL cohort), healthy tissue (n = 46), and fractionated cell lines (n = 5). We found that the recently-discovered full-length intronic circles and the stable lariats formed distinct subclasses, with a surprisingly high intronic circle fraction in BC (∼45%) compared to healthy tissues (0-20%). The stable lariats and their host introns were characterized by small transcript sizes, highly conserved BP regions, enriched BP motifs, and localization in multiple cell fractions. Additionally, circular sisRNAs showed tissue-specific expression patterns. We found nine circular sisRNAs as differentially expressed across early-stage BC patients with different prognoses, and sisHNRNPK expression correlated with progression-free survival. In conclusion, we identify distinguishing biological features of circular sisRNAs and point to specific candidates (incl. sisHNRNPK, sisWDR13 and sisMBNL1) that were highly expressed, had evolutionary conserved sequences, or had clinical correlations, which may facilitate future studies and further insights into their functional roles.

Recursive splicing discovery using lariats in total RNA sequencing

Recursive splicing is a non-canonical splicing mechanism in which an intron is removed in segments via multiple splicing reactions. Relatively few recursive splice sites have been identified with high confidence in human introns, and more comprehensive analyses are needed to better characterize where recursive splicing happens and whether or not it has a regulatory function. In this study, we use an unbiased approach using intron lariats to search for recursive splice sites in constitutive introns and alternative exons in the human transcriptome. We find evidence for recursive splicing in a broader range of intron sizes than previously reported and detail a new location for recursive splicing at the distal ends of cassette exons. In addition, we identify evidence for the conservation of these recursive splice sites among higher vertebrates and the use of these sites to influence alternative exon exclusion. Together, our data demonstrate the prevalence of recursive splicing and its potential influence on gene expression through alternatively spliced isoforms.

Transcriptome sequencing suggests that pre-mRNA splicing counteracts widespread intronic cleavage and polyadenylation

Alternative splicing (AS) and alternative polyadenylation (APA) are two crucial steps in the post-transcriptional regulation of eukaryotic gene expression. Protocols capturing and sequencing RNA 3'-ends have uncovered widespread intronic polyadenylation (IPA) in normal and disease conditions, where it is currently attributed to stochastic variations in the pre-mRNA processing. Here, we took advantage of the massive amount of RNA-seq data generated by the Genotype Tissue Expression project (GTEx) to simultaneously identify and match tissue-specific expression of intronic polyadenylation sites with tissue-specific splicing. A combination of computational methods including the analysis of short reads with non-templated adenines revealed that APA events are more abundant in introns than in exons. While the rate of IPA in composite terminal exons and skipped terminal exons expectedly correlates with splicing, we observed a considerable fraction of IPA events that lack AS support and attributed them to spliced polyadenylated introns (SPI). We hypothesize that SPIs represent transient byproducts of a dynamic coupling between APA and AS, in which the spliceosome removes the intron while it is being cleaved and polyadenylated. These findings indicate that cotranscriptional pre-mRNA splicing could serve as a rescue mechanism to suppress premature transcription termination at intronic polyadenylation sites.

An unknown essential function of tRNA splicing endonuclease is linked to the integrated stress response and intron debranching

tRNA splicing endonuclease (TSEN) has a well-characterized role in transfer RNA (tRNA) splicing but also other functions. For yeast TSEN, these other functions include degradation of a subset of mRNAs that encode mitochondrial proteins and an unknown essential function. In this study, we use yeast genetics to characterize the unknown tRNA-independent function(s) of TSEN. Using a high-copy suppressor screen, we found that sen2 mutants can be suppressed by overexpression of SEN54. This effect was seen both for tRNA-dependent and tRNA-independent functions indicating that SEN54 is a general suppressor of sen2, likely through structural stabilization. A spontaneous suppressor screen identified mutations in the intron-debranching enzyme, Dbr1, as tRNA splicing-independent suppressors. Transcriptome analysis showed that sen2 mutation activates the Gcn4 stress response. These Gcn4 target transcripts decreased considerably in the sen2 dbr1 double mutant. We propose that Dbr1 and TSEN may compete for a shared substrate, which TSEN normally processes into an essential RNA, while Dbr1 initiates its degradation. These data provide further insight into the essential function(s) of TSEN. Importantly, single amino acid mutations in TSEN cause the generally fatal neuronal disease pontocerebellar hypoplasia (PCH). The mechanism by which defects in TSEN cause this disease is unknown, and our results reveal new possible mechanisms.

Maturation of small nucleolar RNAs: from production to function

Small Nucleolar RNAs (snoRNAs) are an abundant group of non-coding RNAs with well-defined roles in ribosomal RNA processing, folding and chemical modification. Besides their classic roles in ribosome biogenesis, snoRNAs are also implicated in several other cellular activities including regulation of splicing, transcription, RNA editing, cellular trafficking, and miRNA-like functions. Mature snoRNAs must undergo a series of processing steps tightly regulated by transiently associating factors and coordinated with other cellular processes including transcription and splicing. In addition to their mature forms, snoRNAs can contribute to gene expression regulation through their derivatives and degradation products. Here, we review the current knowledge on mechanisms of snoRNA maturation, including the different pathways of processing, and the regulatory mechanisms that control snoRNA levels and complex assembly. We also discuss the significance of studying snoRNA maturation, highlight the gaps in the current knowledge and suggest directions for future research in this area.

The microRNA Lifecycle in Health and Cancer

MicroRNAs (miRNAs) are small non-coding RNAs of ~22 nucleotides that regulate gene expression at the post-transcriptional level. They can bind to around 60% of all protein-coding genes with an average of 200 targets per miRNA, indicating their important function within physiological and pathological cellular processes. miRNAs can be quickly produced in high amounts through canonical and non-canonical pathways that involve a multitude of steps and proteins. In cancer, miRNA biogenesis, availability and regulation of target expression can be altered to promote tumour progression. This can be due to genetic causes, such as single nucleotide polymorphisms, epigenetic changes, differences in host gene expression, or chromosomal remodelling. Alternatively, post-transcriptional changes in miRNA stability, and defective or absent components and mediators of the miRNA-induced silencing complex can lead to altered miRNA function. This review provides an overview of the current knowledge on the lifecycle of miRNAs in health and cancer. Understanding miRNA function and regulation is fundamental prior to potential future application of miRNAs as cancer biomarkers.

From the comparative study of a circRNA originating from an mammalian ATXN2L intron to understanding the genesis of intron lariat-derived circRNAs

Circular intronic RNAs (ciRNAs) are still unexplored regarding mechanisms for their emergence. We considered the ATXN2L intron lariat-derived circular RNA (ciRNA-ATXN2L) as an opportunity to conduct a cross-species examination of ciRNA genesis. To this end, we investigated 207 datasets from 4 tissues and from 13 mammalian species. While in eight species, ciRNA-ATXN2L was never detected, in pigs and rabbits, ciRNA-ATXN2L was expressed in all tissues and sometimes at very high levels. Bovine tissues were an intermediate case and in macaques and cats, only ciRNA-ATXN2L traces were detected. The pattern of ciRNA-ATXN2L restricted to only five species is not related to a particular evolution of intronic sequences. To empower our analysis, we considered 221 additional introns including 80 introns where a lariat-derived ciRNA was previously described. The primary driver of micro-ciRNA genesis (< 155 nt as ciRNA-ATXN2L) appears to be the absence of a canonical "A" (i.e. a "tnA" located in the usual branching region) to build the lariat around this adenosine. The balance between available "non canonical-A" (no ciRNA genesis) and "non-A" (ciRNA genesis) for use as a branch point to build the lariat could modify the expression level of ciRNA-ATXN2L. In addition, the rare localization of the 2'-5' bond in an open RNA secondary structure could also negatively affect the lifetime of ciRNAs (macaque ciRNA-ATXN2L). Our analyses suggest that ciRNA-ATXN2L is likely a functionless splice remnant. This study provides a better understanding of the ciRNAs origin, especially drivers for micro ciRNA genesis.

The Role of DBR1 as a Candidate Prognosis Biomarker in Esophageal Squamous Cell Carcinoma

Aims: Esophageal squamous cell carcinoma (ESCC) is one of the most prevalent malignancies with unfavorable clinical outcomes and limited therapeutic methods. As a key enzyme in RNA metabolism, debranching RNA Lariats 1 (DBR1) is involved in intron turnover and biogenesis of noncoding RNA. Although cancer cells often show disorder of nucleic acid metabolism, it is unclear whether DBR1 has any effect on the carcinogenesis and progression of ESCC. Methods: Here we detected DBR1 expression in 112 ESCC samples by immunohistochemistry and analyzed its correlation with clinical parameters and survival. Results: DBR1 is mainly located in the nucleus of ESCC tissue. And DBR1 was associated with several malignant clinical features in patients, including tumor location (χ² = 9.687, P = .021), pathologic T stage (χ² = 5.771, P = .016), lymph node metastasis (χ² = 8.215, P = .004) and N classification (χ² = 10.066, P = .018). Moreover, Kaplan-Meier analysis showed that ESCC patients harboring lower DBR1 expression had a worse prognosis in comparison with those with higher DBR1 expression (P = .005). Univariate and multivariate Cox proportional hazards regression analyses indicated that decreased DBR1 might act as an independent predictor of poor prognosis for ESCC patients. Conclusion: Abnormal RNA metabolism might play a critical role in promoting the progression of ESCC, and DBR1 may be a promising potential biomarker for predicting the prognosis of ESCC patients.

Internally Symmetrical Stwintrons and Related Canonical Introns in Hypoxylaceae Species

Spliceosomal introns are pervasive in eukaryotes. Intron gains and losses have occurred throughout evolution, but the origin of new introns is unclear. Stwintrons are complex intervening sequences where one of the sequence elements (5′-donor, lariat branch point element or 3′-acceptor) necessary for excision of a U2 intron (external intron) is itself interrupted by a second (internal) U2 intron. In Hypoxylaceae, a family of endophytic fungi, we uncovered scores of donor-disrupted stwintrons with striking sequence similarity among themselves and also with canonical introns. Intron–exon structure comparisons suggest that these stwintrons have proliferated within diverging taxa but also give rise to proliferating canonical introns in some genomes. The proliferated (stw)introns have integrated seamlessly at novel gene positions. The recently proliferated (stw)introns appear to originate from a conserved ancestral stwintron characterised by terminal inverted repeats (45–55 nucleotides), a highly symmetrical structure that may allow the formation of a double-stranded intron RNA molecule. No short tandem duplications flank the putatively inserted intervening sequences, which excludes a DNA transposition-based mechanism of proliferation. It is tempting to suggest that this highly symmetrical structure may have a role in intron proliferation by (an)other mechanism(s).

Linking circular intronic RNA degradation and function in transcription by RNase H1

Circular intronic RNAs (ciRNAs) escaping from DBR1 debranching of intron lariats are co-transcriptionally produced from pre-mRNA splicing, but their turnover and mechanism of action have remained elusive. We report that RNase H1 degrades a subgroup of ciRNAs in human cells. Many ciRNAs contain high GC% and tend to form DNA:RNA hybrids (R-loops) for RNase H1 cleavage, a process that appears to promote Pol II transcriptional elongation at ciRNA-producing loci. One ciRNA, ciankrd52, shows a stronger ability of R-loop formation than that of its cognate pre-mRNA by maintaining a locally open RNA structure in vitro. This allows the release of pre-mRNA from R-loops by ci-ankrd52 replacement and subsequent ciRNA removal via RNase H1 for efficient transcriptional elongation. We propose that such an R-loop dependent ciRNA degradation likely represents a mechanism that on one hand limits ciRNA accumulation by recruiting RNase H1 and on the other hand resolves R-loops for transcriptional elongation at some GC-rich ciRNA-producing loci.

Mechanistic Insights of Aberrant Splicing with Splicing Factor Mutations Found in Myelodysplastic Syndromes

Pre-mRNA splicing is an essential process for gene expression in higher eukaryotes, which requires a high order of accuracy. Mutations in splicing factors or regulatory elements in pre-mRNAs often result in many human diseases. Myelodysplastic syndrome (MDS) is a heterogeneous group of chronic myeloid neoplasms characterized by many symptoms and a high risk of progression to acute myeloid leukemia. Recent findings indicate that mutations in splicing factors represent a novel class of driver mutations in human cancers and affect about 50% of Myelodysplastic syndrome (MDS) patients. Somatic mutations in MDS patients are frequently found in genes SF3B1, SRSF2, U2AF1, and ZRSR2. Interestingly, they are involved in the recognition of 3' splice sites and exons. It has been reported that mutations in these splicing regulators result in aberrant splicing of many genes. In this review article, we first describe molecular mechanism of pre-mRNA splicing as an introduction and mainly focus on those four splicing factors to describe their mutations and their associated aberrant splicing patterns.