Species-specific signals for the splicing of a short Drosophila intron in vitro

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Pathogenic Abnormal Splicing Due to Intronic Deletions that Induce Biophysical Space Constraint for Spliceosome Assembly

A precise genetic diagnosis is the single most important step for families with genetic disorders to enable personalized and preventative medicine. In addition to genetic variants in coding regions (exons) that can change a protein sequence, abnormal pre-mRNA splicing can be devastating for the encoded protein, inducing a frameshift or in-frame deletion/insertion of multiple residues. Non-coding variants that disrupt splicing are extremely challenging to identify. Stemming from an initial clinical discovery in two index Australian families, we define 25 families with genetic disorders caused by a class of pathogenic non-coding splice variant due to intronic deletions. These pathogenic intronic deletions spare all consensus splice motifs, though they critically shorten the minimal distance between the 5' splice-site (5'SS) and branchpoint. The mechanistic basis for abnormal splicing is due to biophysical constraint precluding U1/U2 spliceosome assembly, which stalls in A-complexes (that bridge the 5'SS and branchpoint). Substitution of deleted nucleotides with non-specific sequences restores spliceosome assembly and normal splicing, arguing against loss of an intronic element as the primary causal basis. Incremental lengthening of 5'SS-branchpoint length in our index EMD case subject defines 45-47 nt as the critical elongation enabling (inefficient) spliceosome assembly for EMD intron 5. The 5'SS-branchpoint space constraint mechanism, not currently factored by genomic informatics pipelines, is relevant to diagnosis and precision medicine across the breadth of Mendelian disorders and cancer genomics.

Mathematical modeling identifies potential gene structure determinants of co-transcriptional control of alternative pre-mRNA splicing

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons' inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon-intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. 'intron definition') or across short exons (i.e. 'exon definition'), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture

Production of most eukaryotic mRNAs requires splicing of introns from pre-mRNA. The splicing reaction requires definition of splice sites, which are initially recognized in either intron-spanning (‘intron definition’) or exon-spanning (‘exon definition’) pairs. To understand how exon and intron length and splice site recognition mode impact splicing, we measured splicing rates genome-wide in Drosophila, using metabolic labeling/RNA sequencing and new mathematical models to estimate rates. We found that the modal intron length range of 60–70 nt represents a local maximum of splicing rates, but that much longer exon-defined introns are spliced even faster and more accurately. We observed unexpectedly low variation in splicing rates across introns in the same gene, suggesting the presence of gene-level influences, and we identified multiple gene level variables associated with splicing rate. Together our data suggest that developmental and stress response genes may have preferentially evolved exon definition in order to enhance the rate or accuracy of splicing.

Cardinal Epigenetic Role of non-coding Regulatory RNAs in Circadian Rhythm

Circadian rhythm which governs basic physiological activities like sleeping, feeding and energy consumption is regulated by light-controlled central clock genes in the pacemaker neuron. The timekeeping machinery with unique transcriptional and post-transcriptional feedback loops is controlled by different small regulatory RNAs in the brain. Roles of the multiple neuronal genes, especially post-transcriptional regulation, splicing, polyadenylation, mature mRNA editing, and stability of translation products, are controlled by epigenetic activities orchestrated via small RNAs. Collectively, these mechanisms regulate clock and light-controlled genes for effecting pacemaker activity and entrainment. Regulatory small RNAs of the circadian circuit, timekeeping mechanism, synchronization of regular entrainment, oscillation, and rhythmicity are regulated by diversified RNA molecules. Regulatory small RNAs operate critical roles in brain activities including the neuronal clock activity. In this report, we propose the emergence of the earlier unexpected small RNAs for a historic perspective of epigenetic regulation of the brain clock system.

Prp40 and early events in splice site definition

The alternative splicing (AS) of precursor messenger RNA (pre-mRNA) is a tightly regulated process through which introns are removed to leave the resulting exons in the mRNA appropriately aligned and ligated. The AS of pre-mRNA is a key mechanism for increasing the complexity of proteins encoded in the genome. In humans, more than 90% of genes undergo AS, underscoring the importance of this process in RNA biogenesis. As such, AS misregulation underlies multiple human diseases. The splicing reaction is catalyzed by the spliceosome, a highly dynamic complex that assembles at or near the intron/exon boundaries and undergoes sequential conformational and compositional changes during splicing. The initial recognition of splice sites defines the exons that are going to be removed, which is a critical step in the highly regulated splicing process. Although the available lines of evidence are increasing, the molecular mechanisms governing AS, including the initial interactions occurring at intron/exon boundaries, and the factors that modulate these critical connections by functioning as a scaffold for active-site RNAs or proteins, remain poorly understood. In this review, we summarize the major hallmarks of the initial steps in the splicing process and the role of auxiliary factors that contribute to the assembly of the spliceosomal complex. We also discuss the role of the essential yeast Prp40 protein and its mammalian homologs in the specificity of this pre-mRNA processing event. In addition, we provide the first exhaustive phylogenetic analysis of the molecular evolution of Prp40 family members. WIREs RNA 2016, 7:17-32. doi: 10.1002/wrna.1312 For further resources related to this article, please visit the WIREs website.

Differential GC content between exons and introns establishes distinct strategies of splice-site recognition

During evolution segments of homeothermic genomes underwent a GC content increase. Our analyses reveal that two exon-intron architectures have evolved from an ancestral state of low GC content exons flanked by short introns with a lower GC content. One group underwent a GC content elevation that abolished the differential exon-intron GC content, with introns remaining short. The other group retained the overall low GC content as well as the differential exon-intron GC content, and is associated with longer introns. We show that differential exon-intron GC content regulates exon inclusion level in this group, in which disease-associated mutations often lead to exon skipping. This group's exons also display higher nucleosome occupancy compared to flanking introns and exons of the other group, thus "marking" them for spliceosomal recognition. Collectively, our results reveal that differential exon-intron GC content is a previously unidentified determinant of exon selection and argue that the two GC content architectures reflect the two mechanisms by which splicing signals are recognized: exon definition and intron definition.

A 31 bp VNTR in the cystathionine beta-synthase (CBS) gene is associated with reduced CBS activity and elevated post-load homocysteine levels

Molecular defects in genes encoding enzymes involved in homocysteine metabolism may account for mild hyperhomocysteinaemia, an independent and graded risk factor for cardiovascular disease (CVD). Although heterozygosity for cystathionine beta-synthase (CBS) deficiency has been excluded as a major genetic cause of mild hyperhomocysteinaemia in vascular disease, mutations in (non-)coding DNA sequences may lead to a mildly decreased CBS expression and, consequently, to elevated plasma homocysteine levels. We assessed the association between a 31 bp VNTR, that spans the exon 13-intron 13 boundary of the CBS gene, and fasting, post-methionine load and increase upon methionine load plasma homocysteine levels in 190 patients with arterial occlusive disease, and in 381 controls. The 31 bp VNTR consists of 16, 17, 18, 19 or 21 repeat units and shows a significant increase in plasma homocysteine concentrations with an increasing number of repeat elements, in particular after methionine loading. In 26 vascular disease patients the relationship between this 31 bp VNTR and CBS enzyme activity in cultured fibroblasts was studied. The CBS enzyme activity decreased with increasing number of repeat units of the 31 bp VNTR. RT-PCR experiments showed evidence of alternative splicing at the exon 13-intron 13 splice junction site. The 31 bp VNTR in the CBS gene is associated with post-methionine load hyperhomocysteinaemia that may predispose individuals to an increased risk of cardiovascular diseases.

Vertebrate exon trapping methods: implications for transcript mapping with mosquito DNA

Exon trapping methods have played an important role in the development of transcript maps. In one in vivo vertebrate method, exons in a genomic DNA clone are transcribed, and they are recovered without any a priori information on the nature of the expressed transcript. The only requirement is that the genomic DNA clone contains exons separated by intervening introns that are removed by splicing during mRNA transcription and that the splice donor and acceptor site sequences follow those used by vertebrates. It is not known whether invertebrate splice donor and acceptor sites from genes that contain short introns will be processed correctly using an in vivo vertebrate exon trapping method. In this report, an analysis of mosquito splice sites using software designed to identify exons in genomic DNA sequence suggested that the vertebrate exon trapping method could recognize mosquito introns and exons. When a mosquito genomic DNA clone containing the D7 gene was tested experimentally, this method failed to recognize and process small introns (< 63 bp) faithfully. In spite of this failure, exons and exon fragments were recovered. The implications of these findings and their application to map-based positional cloning in mosquito genomics is discussed.

Human and fungal 3' splice sites are used by Trypanosoma brucei for trans splicing

In Trypanosoma brucei, pre-mRNAs are joined to a 5' 39 nt spliced leader sequence by trans splicing, a process that has not been well characterized. We have asked whether the 3' splice site regions of human and yeast introns are able to substitute in vivo for the 3' spliced leader acceptor regions of trypanosome pre-mRNA sequences. The ability of heterologous sequences to participate in trans splicing in trypanosomes was assayed by chloramphenicol acetyltransferase (CAT) enzyme activity and/or the detection of spliced CAT mRNA. Four out of the six heterologous 3' splice site regions (human beta-globin intervening sequence (IVS)2, human c-myc IVS2, human factor-VIII IVS1, and yeast actin IVS) functioned as 3' spliced leader acceptor regions in T. brucei, while two did not show significant or detectable levels of CAT activity (human beta-globin IVS1 and human c-myc IVS1). In the case of the human beta-globin IVS1 however, lengthening of the polypyrimidine tract as a result of single purine to pyrimidine transversions produced an active acceptor in which the spliced leader addition site coincides with the 3' splice site of the beta-globin exon 2. These studies indicate that some, but not all 3' acceptor regions in humans can function as spliced leader addition sites in trypansomes.

Distribution and evolution of introns in Drosophila amylase genes

While the two amylase genes of Drosophila melanogaster are intronless, the three genes of D. pseudoobscura harbor a short intron. This raises the question of the common structure of the Amy gene in Drosophila species. We have investigated the presence or absence of an intron in the amylase genes of 150 species of Drosophilids. Using polymerase chain reaction (PCR), we have amplified a region that surrounds the intron site reported in D. pseudoobscura and a few other species. The results revealed that most species contain an intron, with a variable size ranging from 50 to 750 bp, although the very majoritary size was around 60-80 bp. Several species belonging to different lineages were found to lack an intron. This loss of intervening sequence was likely due to evolutionarily independent and rather frequent events. Some other species had both types of genes: In the obscura group, and to a lesser extent in the ananassae subgroup, intronless copies had much diverged from intron-containing genes. Base composition of short introns was found to be variable and correlated with that of the surrounding exons, whereas long introns were all A-T rich. We have extended our study to non-Drosophilid insects. In species from other orders of Holometaboles, Lepidoptera and Hymenoptera, an intron was found at an identical position in the Amy gene, suggesting that the intron was ancestral.

Sequence divergence associated with species-specific splicing of the nonmuscle beta-tropomyosin alternative exon

Alternative splicing of vertebrate beta-tropomyosin transcripts ensures mutually exclusive expression of internal exons 6A and 6B in nonmuscle and skeletal muscle cells, respectively. Recently, we reported that this splicing regulation requires species-specific elements, since the splicing profile for the chicken, rat, and Xenopus beta-tropomyosin alternative exons is not reproduced in transfection experiments when heterologous myogenic cells are used. By analyzing the splicing pattern of hybrid chicken/rat beta-TM constructions transfected into both quail and mouse cell lines, we demonstrate that chicken beta-tropomyosin exon 6A is flanked by stronger splicing signals than rat exon 6A, thus leading to the misregulation of splicing in heterologous cells. We have characterized three splicing signals that contribute to this difference: 1) nonconsensus nucleotide differences at positions +4 and +6 in the donor site downstream of exon 6A, 2) differences in the pyrimidine composition between the branch site and acceptor site upstream of exon 6A, and 3) a pyrimidine-rich intronic exon 6A splicing enhancer present upstream of exon 6A only in the chicken beta-TM gene. The functional divergence between splicing signals in two homologous vertebrate genes reveals species-specific strategies for proper modulation of splicing of alternative exons.

Suppressor U1 snRNAs in Drosophila

Although the role of U1 small nuclear RNAs (snRNAs) in 5' splice site recognition is well established, suppressor U1 snRNAs active in intact multicellular animals have been lacking. Here we describe suppression of a 5' splice site mutation in the Drosophila melanogaster white gene (wDR18) by compensatory changes in U1 snRNA. Mutation of positions -1 and +6 of the 5' splice site of the second intron (ACG[GTGAGT to ACC]GTGAGC) results in the accumulation of RNA retaining this 74-nucleotide intron in both transfected cells and transgenic flies. U1-3G, a suppressor U1 snRNA which restores base-pairing at position +6 of the mutant intron, increases the ratio of spliced to unspliced wDR18 RNA up to fivefold in transfected Schneider cells and increases eye pigmentation in wDR18 flies. U1-9G, which targets position -1, suppresses wDR18 in transfected cells less well. U1-3G,9G has the same effect as U1-3G although it accumulates to lower levels. Suppression of wDR18 has revealed that the U1b embryonic variant (G134 to U) is active in Schneider cells and pupal eye discs. However, the combination of 9G with 134U leads to reduced accumulation of both U1b-9G and U1b-3G,9G, possibly because nucleotides 9 and 134 both participate in a potential long-range intramolecular base-pairing interaction. High levels of functional U1-3G suppressor reduce both viability and fertility in transformed flies. These results show that, despite the difficulties inherent in stably altering splice site selection in multicellular organisms, it is possible to obtain suppressor U1 snRNAs in flies.