A multimodal analysis of genomic and RNA splicing features in myeloid malignancies

RNA splicing dysfunctions are more widespread than what is believed by only estimating the effects resulting by splicing factor mutations (SF^MT) in myeloid neoplasia (MN). The genetic complexity of MN is amenable to machine learning (ML) strategies. We applied an integrative ML approach to identify co-varying features by combining genomic lesions (mutations, deletions, and copy number), exon-inclusion ratio as measure of RNA splicing (percent spliced in, PSI), and gene expression (GE) of 1,258 MN and 63 normal controls. We identified 15 clusters based on mutations, GE, and PSI. Different PSI levels were present at various extents regardless of SF^MT suggesting that changes in RNA splicing were not strictly related to SF^MT. Combination of PSI and GE further distinguished the features and identified PSI similarities and differences, common pathways, and expression signatures across clusters. Thus, multimodal features can resolve the complex architecture of MN and help identifying convergent molecular and transcriptomic pathways amenable to therapies.

Germline DDX41 mutations cause ineffective hematopoiesis and myelodysplasia

DDX41 mutations are the most common germline alterations in adult myelodysplastic syndromes (MDSs). The majority of affected individuals harbor germline monoallelic frameshift DDX41 mutations and subsequently acquire somatic mutations in their other DDX41 allele, typically missense R525H. Hematopoietic progenitor cells (HPCs) with biallelic frameshift and R525H mutations undergo cell cycle arrest and apoptosis, causing bone marrow failure in mice. Mechanistically, DDX41 is essential for small nucleolar RNA (snoRNA) processing, ribosome assembly, and protein synthesis. Although monoallelic DDX41 mutations do not affect hematopoiesis in young mice, a subset of aged mice develops features of MDS. Biallelic mutations in DDX41 are observed at a low frequency in non-dominant hematopoietic stem cell clones in bone marrow (BM) from individuals with MDS. Mice chimeric for monoallelic DDX41 mutant BM cells and a minor population of biallelic mutant BM cells develop hematopoietic defects at a younger age, suggesting that biallelic DDX41 mutant cells are disease modifying in the context of monoallelic DDX41 mutant BM.

Functional analyses of human LUC7-like proteins involved in splicing regulation and myeloid neoplasms

Vertebrates have evolved three paralogs, termed LUC7L, LUC7L2, and LUC7L3, of the essential yeast U1 small nuclear RNA (snRNA)-associated splicing factor Luc7p. We investigated the mechanistic and regulatory functions of these putative splicing factors, of which one (LUC7L2) is mutated or deleted in myeloid neoplasms. Protein interaction data show that all three proteins bind similar core but distinct regulatory splicing factors, probably mediated through their divergent arginine-serine-rich domains, which are not present in Luc7p. Knockdown of each factor reveals mostly unique sets of significantly dysregulated alternative splicing events dependent on their binding locations, which are largely non-overlapping. Notably, knockdown of LUC7L2 alone significantly upregulates the expression of multiple spliceosomal factors and downregulates glycolysis genes, possibly contributing to disease pathogenesis. RNA binding studies reveal that LUC7L2 and LUC7L3 crosslink to weak 5' splice sites and to the 5' end of U1 snRNA, establishing an evolutionarily conserved role in 5' splice site selection.

Excessive R-loops trigger an inflammatory cascade leading to increased HSPC production

Hematopoietic stem and progenitor cells (HSPCs) arise during embryonic development and are essential for sustaining the blood and immune systems throughout life. Tight regulation of HSPC numbers is critical for hematopoietic homeostasis. Here, we identified DEAD-box helicase 41 (Ddx41) as a gatekeeper of HSPC production. Using zebrafish ddx41 mutants, we unveiled a critical role for this helicase in regulating HSPC production at the endothelial-to-hematopoietic transition. We determined that Ddx41 suppresses the accumulation of R-loops, nucleic acid structures consisting of RNA:DNA hybrids and ssDNAs whose equilibrium is essential for cellular fitness. Excess R-loop levels in ddx41 mutants triggered the cGAS-STING inflammatory pathway leading to increased numbers of hemogenic endothelium and HSPCs. Elevated R-loop accumulation and inflammatory signaling were observed in human cells with decreased DDX41, suggesting possible conservation of mechanism. These findings delineate that precise regulation of R-loop levels during development is critical for limiting cGAS-STING activity and HSPC numbers.

Comprehensive database and evolutionary dynamics of U12-type introns

During nuclear maturation of most eukaryotic pre-messenger RNAs and long non-coding RNAs, introns are removed through the process of RNA splicing. Different classes of introns are excised by the U2-type or the U12-type spliceosomes, large complexes of small nuclear ribonucleoprotein particles and associated proteins. We created intronIC, a program for assigning intron class to all introns in a given genome, and used it on 24 eukaryotic genomes to create the Intron Annotation and Orthology Database (IAOD). We then used the data in the IAOD to revisit several hypotheses concerning the evolution of the two classes of spliceosomal introns, finding support for the class conversion model explaining the low abundance of U12-type introns in modern genomes.

Complex landscape of alternative splicing in myeloid neoplasms

Myeloid neoplasms are characterized by frequent mutations in at least seven components of the spliceosome that have distinct roles in the process of pre-mRNA splicing. Hotspot mutations in SF3B1, SRSF2, U2AF1 and loss of function mutations in ZRSR2 have revealed widely different aberrant splicing signatures with little overlap. However, previous studies lacked the power necessary to identify commonly mis-spliced transcripts in heterogeneous patient cohorts. By performing RNA-Seq on bone marrow samples from 1,258 myeloid neoplasm patients and 63 healthy bone marrow donors, we identified transcripts frequently mis-spliced by mutated splicing factors (SF), rare SF mutations with common alternative splicing (AS) signatures, and SF-dependent neojunctions. We characterized 17,300 dysregulated AS events using a pipeline designed to predict the impact of mis-splicing on protein function. Meta-splicing analysis revealed a pattern of reduced levels of retained introns among disease samples that was exacerbated in patients with splicing factor mutations. These introns share characteristics with “detained introns,” a class of introns that have been shown to promote differentiation by detaining pro-proliferative transcripts in the nucleus. In this study, we have functionally characterized 17,300 targets of mis-splicing by the SF mutations, identifying a common pathway by which AS may promote maintenance of a proliferative state.

The hnRNP RALY regulates transcription and cell proliferation by modulating the expression of specific factors including the proliferation marker E2F1

The heterogeneous nuclear ribonucleoproteins (hnRNP) form a large family of RNA-binding proteins that exert numerous functions in RNA metabolism. RALY is a member of the hnRNP family that binds poly-U-rich elements within several RNAs and regulates the expression of specific transcripts. RALY is up-regulated in different types of cancer, and its down-regulation impairs cell cycle progression. However, the RALY's role in regulating RNA levels remains elusive. Here, we show that numerous genes coding for factors involved in transcription and cell cycle regulation exhibit an altered expression in RALY-down-regulated HeLa cells, consequently causing impairments in transcription, cell proliferation, and cell cycle progression. Interestingly, by comparing the list of RALY targets with the list of genes affected by RALY down-regulation, we found an enrichment of RALY mRNA targets in the down-regulated genes upon RALY silencing. The affected genes include the E2F transcription factor family. Given its role as proliferation-promoting transcription factor, we focused on E2F1. We demonstrate that E2F1 mRNA stability and E2F1 protein levels are reduced in cells lacking RALY expression. Finally, we also show that RALY interacts with transcriptionally active chromatin in both an RNA-dependent and -independent manner and that this association is abolished in the absence of active transcription. Taken together, our results highlight the importance of RALY as an indirect regulator of transcription and cell cycle progression through the regulation of specific mRNA targets, thus strengthening the possibility of a direct gene expression regulation exerted by RALY.

DDX41-related myeloid neoplasia

While early presentation of familial leukemia syndromes is typical, long disease anticipation may mask cases of familial traits in seemingly spontaneous disease. Germline mutations in DDX41 gene have been discovered in several leukemia families, as well as in mostly adult patients with seemingly spontaneous disease but having strong family histories of myeloid neoplasia. As with other familial genes, DDX41 mutation carriers can develop neoplasia through acquisition of another somatic mutation, thereby affecting both DDX41 alleles. In other patients, somatic mutations of different driver genes can substitute for acquired missense DDX41 during progression. Conversely, non-familial cases with heterozygous somatic DDX41 mutations point towards other mutations that can substitute for the germ line founder DDX41 lesions. In either circumstance, total inactivation of DDX41 appears to be cell-lethal, explaining why frameshift germline lesions have not been found to be accompanied by deletions of the DDX41 locus on 5q. The precise function of the DDX41 protein is unknown; considerable evidence suggests its involvement in RNA splicing. Thus DDX41 can be included in the now large group of mutated spliceosomal genes affected in myeloid neoplasia. However, it appears that DDX4 is so far the only example of a germline spliceosomal mutation in leukemia. Clinically, recognition of DDX41 mutated cases may have implications for surveillance, assessment of prognosis, and, perhaps, for design of targeted therapies.

Inherited and Somatic Defects in DDX41 in Myeloid Neoplasms

Most cases of adult myeloid neoplasms are routinely assumed to be sporadic. Here, we describe an adult familial acute myeloid leukemia (AML) syndrome caused by germline mutations in the DEAD/H-box helicase gene DDX41. DDX41 was also found to be affected by somatic mutations in sporadic cases of myeloid neoplasms as well as in a biallelic fashion in 50% of patients with germline DDX41 mutations. Moreover, corresponding deletions on 5q35.3 present in 6% of cases led to haploinsufficient DDX41 expression. DDX41 lesions caused altered pre-mRNA splicing and RNA processing. DDX41 is exemplary of other RNA helicase genes also affected by somatic mutations, suggesting that they constitute a family of tumor suppressor genes.

Biochemical defects in minor spliceosome function in the developmental disorder MOPD I

Biallelic mutations of the human RNU4ATAC gene, which codes for the minor spliceosomal U4atac snRNA, cause the developmental disorder, MOPD I/TALS. To date, nine separate mutations in RNU4ATAC have been identified in MOPD I patients. Evidence suggests that all of these mutations lead to abrogation of U4atac snRNA function and impaired minor intron splicing. However, the molecular basis of these effects is unknown. Here, we use a variety of in vitro and in vivo assays to address this question. We find that only one mutation, 124G>A, leads to significantly reduced expression of U4atac snRNA, whereas four mutations, 30G>A, 50G>A, 50G>C and 51G>A, show impaired binding of essential protein components of the U4atac/U6atac di-snRNP in vitro and in vivo. Analysis of MOPD I patient fibroblasts and iPS cells homozygous for the most common mutation, 51G>A, shows reduced levels of the U4atac/U6atac.U5 tri-snRNP complex as determined by glycerol gradient sedimentation and immunoprecipitation. In this report, we establish a mechanistic basis for MOPD I disease and show that the inefficient splicing of genes containing U12-dependent introns in patient cells is due to defects in minor tri-snRNP formation, and the MOPD I-associated RNU4ATAC mutations can affect multiple facets of minor snRNA function.

PRPF8 defects cause missplicing in myeloid malignancies

Mutations of spliceosome components are common in myeloid neoplasms. One of the affected genes, PRPF8, encodes the most evolutionarily conserved spliceosomal protein. We identified either recurrent somatic PRPF8 mutations or hemizygous deletions in 15/447 and 24/450 cases, respectively. 50% of PRPF8 mutant and del(17p) cases were found in AML and conveyed poor prognosis. PRPF8 defects correlated with increased myeloblasts and ring sideroblasts in cases without SF3B1 mutations. Knockdown of PRPF8 in K562 and CD34+ primary bone marrow cells increased proliferative capacity. Whole RNA deep sequencing of primary cells from patients with PRPF8 abnormalities demonstrated consistent missplicing defects. In yeast models, homologous mutations introduced into Prp8 abrogated a block experimentally produced in the second step of the RNA splicing process suggesting that the mutants have defects in proof-reading functions. In sum, the exploration of clinical and functional consequences suggests that PRPF8 is a novel leukemogenic gene in myeloid neoplasms with a distinct phenotype likely manifested through aberrant splicing.

Defects in spliceosomal machinery: a new pathway of leukaemogenesis

Proper splicing of pre-mRNA is required for protein synthesis and therefore is a fundamental cellular function. The discovery of a variety of somatic spliceosomal mutations in haematological malignancies, including myeloid neoplasms and chronic lymphocytic leukaemia has pointed to a new leukaemogenic pathway involving spliceosomal dysfunction. Theoretically, spliceosomal mutations can lead to activation of incorrect splice sites, intron retention or aberrant alternative splicing occurring in patterns generated by mutations of individual spliceosomal proteins. Such events can produce a defective balance between protein isoforms leading to functional consequences including defective regulation of proliferation and differentiation. The observed pattern of occurrence of highly specific missense mutations, coupled with the lack of nonsense mutations and deletions, implies a gain-of-function or better gain-of-dysfunction mechanism. Incorrect splicing of downstream genes, such as tumour suppressor genes, may result in haploinsufficient expression through nonsense-mediated mRNA decay. Thus, spliceosomal mutations may, depending on the pattern of affected proteins, lead to similar functional effects on tumour suppressor genes as chromosomal deletions, epigenetic silencing or inactivating/hypomorphic mutations. The prognostic value of the most common mutations and their phenotypic association in the clinical setting is currently under investigation. It is likely that spliceosomal mutations may indicate sensitivity to spliceosome inhibitors applied in the form of a synthetic lethal approach. This review discusses the most current aspects of spliceosomal research in the context of haematological malignancies.

New connections between splicing and human disease

The removal by splicing of introns from the primary transcripts of most mammalian genes is an essential step in gene expression. Splicing is performed by large, complex ribonucleoprotein particles termed spliceosomes. Mammals contain two types that splice out mutually exclusive types of introns. However, the role of the minor spliceosome has been poorly studied. Recent reports have now shown that mutations in one minor spliceosomal snRNA, U4atac, are linked to a rare autosomal recessive developmental defect. In addition, very exciting recent results of exome deep-sequencing have found that recurrent, somatic, heterozygous mutations of other splicing factors occur at high frequencies in particular cancers and pre-cancerous conditions, suggesting that alterations in the core splicing machinery can contribute to tumorigenesis. Mis-splicing of crucial genes may underlie the pathologies of all of these diseases. Identifying these genes and understanding the mechanisms involved in their mis-splicing may lead to advancements in diagnosis and treatment.

Mutations in U4atac snRNA, a component of the minor spliceosome, in the developmental disorder MOPD I

Small nuclear RNAs (snRNAs) are essential factors in messenger RNA splicing. By means of homozygosity mapping and deep sequencing, we show that a gene encoding U4atac snRNA, a component of the minor U12-dependent spliceosome, is mutated in individuals with microcephalic osteodysplastic primordial dwarfism type I (MOPD I), a severe developmental disorder characterized by extreme intrauterine growth retardation and multiple organ abnormalities. Functional assays showed that mutations (30G>A, 51G>A, 55G>A, and 111G>A) associated with MOPD I cause defective U12-dependent splicing. Endogenous U12-dependent but not U2-dependent introns were found to be poorly spliced in MOPD I patient fibroblast cells. The introduction of wild-type U4atac snRNA into MOPD I cells enhanced U12-dependent splicing. These results illustrate the critical role of minor intron splicing in human development.

Biased chromatin signatures around polyadenylation sites and exons

Core RNA-processing reactions in eukaryotic cells occur cotranscriptionally in a chromatin context, but the relationship between chromatin structure and pre-mRNA processing is poorly understood. We observed strong nucleosome depletion around human polyadenylation sites (PAS) and nucleosome enrichment just downstream of PAS. In genes with multiple alternative PAS, higher downstream nucleosome affinity was associated with higher PAS usage, independently of known PAS motifs that function at the RNA level. Conversely, exons were associated with distinct peaks in nucleosome density. Exons flanked by long introns or weak splice sites exhibited stronger nucleosome enrichment, and incorporation of nucleosome density data improved splicing simulation accuracy. Certain histone modifications, including H3K36me3 and H3K27me2, were specifically enriched on exons, suggesting active marking of exon locations at the chromatin level. Together, these findings provide evidence for extensive functional connections between chromatin structure and RNA processing.

The conserved 3' end domain of U6atac snRNA can direct U6 snRNA to the minor spliceosome

U6 and U6atac snRNAs play analogous critical roles in the major U2-dependent and minor U12-dependent spliceosomes, respectively. Previous results have shown that most of the functional cores of these two snRNAs are either highly similar in sequence or functionally interchangeable. Thus, a mechanism must exist to restrict each snRNA to its own spliceosome. Here we show that a chimeric U6 snRNA containing the unique and highly conserved 3' end domain of U6atac snRNA is able to function in vivo in U12-dependent spliceosomal splicing. Function of this chimera required the coexpression of a modified U4atac snRNA; U4 snRNA could not substitute. Partial deletions of this element in vivo, as well as in vitro antisense experiments, showed that the 3' end domain of U6atac snRNA is necessary to direct the U4atac/U6atac.U5 tri-snRNP to the forming U12-dependent spliceosome. In vitro experiments also uncovered a role for U4atac snRNA in this targeting.

Mutational analysis of the U12-dependent branch site consensus sequence

Highly conserved sequences at the 5' splice site and branch site of U12-dependent introns are important determinants for splicing by U12-dependent spliceosomes. This study investigates the in vivo splicing phenotypes of mutations in the branch site consensus sequence of the U12-dependent intron F from a human NOL1 (P120) minigene. Intron F contains a fully consensus branch site sequence (UUCCUUAAC). Mutations at each position were analyzed for their effects on U12-dependent splicing in vivo. Mutations at most positions resulted in a significant reduction of correct U12-dependent splicing. Defects observed included increased unspliced RNA levels, the activation of cryptic U2-dependent 5' and 3' splice sites, and the activation of cryptic U12-dependent branch/3' splice sites. A strong correlation was observed between the predicted thermodynamic stability of the branch site: U12 snRNA interaction and correct U12-dependent splicing. The lack of a polypyrimidine tract between the branch site and 3' splice site of U12-dependent introns and the observed reliance on base-pairing interactions for correct U12-dependent splicing emphasize the importance of RNA/RNA interactions during U12-dependent intron recognition and proper splice site selection.

A glimpse into the active site of a group II intron and maybe the spliceosome, too

The X-ray crystal structure of an excised group II self-splicing intron was recently solved by the Pyle group. Here we review some of the notable features of this structure and what they may tell us about the catalytic active site of the group II ribozyme and potentially the spliceosome. The new structure validates the central role of domain V in both the structure and catalytic function of the ribozyme and resolves several outstanding puzzles raised by previous biochemical, genetic and structural studies. While lacking both exons as well as the cleavage sites and nucleophiles, the structure reveals how a network of tertiary interactions can position two divalent metal ions in a configuration that is ideal for catalysis.

Fluorescence and solution NMR study of the active site of a 160-kDa group II intron ribozyme

We have reconstructed the group II intron from Pylaiella littoralis (PL) into a hydrolytic ribozyme, comprising domains 1-3 (D123) connected in cis plus domain 5 (D5) supplied in trans that efficiently cleaves spliced exon substrates. Using a novel gel-based fluorescence assay and nuclear magnetic resonance (NMR) spectroscopy, we monitored the direct binding of D5 to D123, characterized the kinetics of the spliced exon hydrolysis reaction (which is mechanistically analogous to the reverse of the second catalytic step of splicing), and identified the binding surface of D123 on D5. This PL ribozyme acts as an RNA endonuclease even at low monovalent (100 mM KCl) and divalent ion concentrations (1-10 mM MgCl(2)). This is in contrast to other group II intron ribozyme systems that require high levels of salt, making NMR analysis problematic. D5 binds tightly to D123 with a K(d) of 650 +/- 250 nM, a K(m) of approximately 300 nM, and a K(cat) of 0.02 min(-1) under single turnover conditions. Within the approximately 160-kDa D123-D5 binary complex, site-specific binding to D123 leads to dramatic chemical shift perturbation of residues localized to the tetraloop and internal bulge within D5, suggesting a structural switch model for D5-assisted splicing. This minimal ribozyme thus recapitulates the essential features of the reverse of the second catalytic step and represents a well-behaved system for ongoing high-resolution structural work to complement folding and catalytic functional studies.

Structure of a self-splicing group II intron catalytic effector domain 5: parallels with spliceosomal U6 RNA

Domain 5 (D5) is absolutely required for all catalytic functions of group II introns. Here we describe the solution NMR structure, electrostatic calculations, and detailed magnesium ion-binding surface of D5 RNA from the Pylaiella littoralis large ribosomal RNA intron (D5-PL). The overall structure consists of a hairpin capped by a GNRA tetraloop. The stem is divided into lower and upper helices of 8 and 5 bp, respectively, separated by an internal bulge. The D5-PL internal bulge nucleotides stack into the helical junction, resulting in a coupling between the bulge A25 and the closing base pair (G8-C27) of the lower helix. Comparison of the D5-PL structure to previously reported related structures indicates that our structure is most similar, in the helical regions, to the crystal structure of D5 from yeast Ai5gamma (D5-Ai5gamma) and the NMR structure of the U6 snRNA stem-loop region. Our structure differs in many respects from both the NMR and X-ray structures of D5-Ai5gamma in the bulge region. Electrostatic calculations and NMR chemical shift perturbation analyses reveal magnesium ion-binding sites in the tetraloop, internal bulge, and the AGC triad in the lower stem. Our results suggest that the structure, electrostatic environment, and the magnesium ion-binding sites within the tetraloop, bulge, and triad regions are conserved features of the splicing machinery of both the group II introns and the spliceosome that are likely key for catalytic function.

A mutational analysis of U12-dependent splice site dinucleotides

Introns spliced by the U12-dependent minor spliceosome are divided into two classes based on their splice site dinucleotides. The /AU-AC/ class accounts for about one-third of U12-dependent introns in humans, while the /GU-AG/ class accounts for the other two-thirds. We have investigated the in vivo and in vitro splicing phenotypes of mutations in these dinucleotide sequences. A 5' A residue can splice to any 3' residue, although C is preferred. A 5' G residue can splice to 3' G or U residues with a preference for G. Little or no splicing was observed to 3' A or C residues. A 5' U or C residue is highly deleterious for U12-dependent splicing, although some combinations, notably 5' U to 3' U produced detectable spliced products. The dependence of 3' splice site activity on the identity of the 5' residue provides evidence for communication between the first and last nucleotides of the intron. Most mutants in the second position of the 5' splice site and the next to last position of the 3' splice site were defective for splicing. Double mutants of these residues showed no evidence of communication between these nucleotides. Varying the distance between the branch site and the 3' splice site dinucleotide in the /GU-AG/ class showed that a somewhat larger range of distances was functional than for the /AU-AC/ class. The optimum branch site to 3' splice site distance of 11-12 nucleotides appears to be the same for both classes.

U4 small nuclear RNA can function in both the major and minor spliceosomes

U4 small nuclear RNA (snRNA) and U6 snRNA form a base-paired di-snRNP complex that is essential for pre-mRNA splicing of the major class of metazoan nuclear introns. The functionally analogous but highly diverged U4atac and U6atac snRNAs form a similar complex that is involved in splicing of the minor class of introns. Previous results with mutants of U6atac in which a substructure was replaced by the analogous structure from U6 snRNA suggested that wild-type U4 snRNA might be able to interact productively with the mutant U6atac snRNA. Here we show that a mutant U4 snRNA designed to base pair with a mutant U6atac snRNA can activate U12-dependent splicing when coexpressed in an in vivo genetic suppression assay. This genetic interaction could also be demonstrated in an in vitro crosslinking assay. These results show that a U4/U6atac di-snRNP can correctly splice a U12-dependent intron and suggest that the specificity for spliceosome type resides in the U6 and U6atac small nuclear ribonucleoproteins. Further experiments suggest that expression of a mutant U4 snRNA that can bind to wild-type U6atac snRNA alters the specificity of some splice sites, providing an evolutionary rationale for maintaining two U4-like snRNAs.

Domains of human U4atac snRNA required for U12-dependent splicing in vivo

U4atac snRNA forms a base-paired complex with U6atac snRNA. Both snRNAs are required for the splicing of the minor U12-dependent class of eukaryotic nuclear introns. We have developed a new genetic suppression assay to investigate the in vivo roles of several regions of U4atac snRNA in U12-dependent splicing. We show that both the stem I and stem II regions, which have been proposed to pair with U6atac snRNA, are required for in vivo splicing. Splicing activity also requires U4atac sequences in the 5' stem-loop element that bind a 15.5 kDa protein that also binds to a similar region of U4 snRNA. In contrast, mutations in the region immediately following the stem I interaction region, as well as a deletion of the distal portion of the 3' stem-loop element, were active for splicing. Complete deletion of the 3' stem-loop element abolished in vivo splicing function as did a mutation of the Sm protein binding site. These results show that the in vivo sequence requirements of U4atac snRNA are similar to those described previously for U4 snRNA using in vitro assays and provide experimental support for models of the U4atac/U6atac snRNA interaction.

A catalytically active group II intron domain 5 can function in the U12-dependent spliceosome

Both spliceosomal and self-splicing group II introns require the function of similar small, metal binding RNA stem-loop elements located in U6 or U6atac snRNAs of the spliceosome or domain 5 (D5) of group II introns. Here we report that two different D5 elements can functionally replace the U6atac snRNA stem-loop in an in vivo splicing assay. For efficient function in vivo, a single base pair from the upper helical section of the D5 sequence had to be removed. Introducing the equivalent base pair deletion into the D5 element of a group II intron reduced but did not eliminate self-splicing activity. Our results strengthen the case that these RNA elements play similar roles in the catalytic centers of both the spliceosome and a self-splicing ribozyme.

Alternative splicing of U12-dependent introns in vivo responds to purine-rich enhancers

Alternative splicing increases the coding capacity of genes through the production of multiple protein isoforms by the conditional use of splice sites and exons. Many alternative splice sites are regulated by the presence of purine-rich splicing enhancer elements (ESEs) located in the downstream exon. Although the role of ESEs in alternative splicing of the major class U2-dependent introns is well established, no alternatively spliced minor class U12-dependent introns have so far been described. Although in vitro studies have shown that ESEs can stimulate splicing of individual U12-dependent introns, there is no direct evidence that the U12-dependent splicing system can respond to ESEs in vivo. To investigate the ability of U12-dependent introns to use alternative splice sites and to respond to ESEs in an in vivo context, we have constructed two sets of artificial minigenes with alternative splicing pathways and evaluated the effects of ESEs on their alternative splicing patterns. In minigenes with alternative U12-dependent 3' splice sites, a purine-rich ESE promotes splicing to the immediately upstream 3' splice site. As a control, a mutant ESE has no stimulatory effect. In minigene constructs with two adjacent U12-dependent introns, the predominant in vivo splicing pattern results in the skipping of the internal exon. Insertion of a purine-rich ESE into the internal exon promotes the inclusion of the internal exon. These results show that U12-dependent introns can participate in alternative splicing pathways and that U12-dependent splice sites can respond to enhancer elements in vivo.

Role of the 3' splice site in U12-dependent intron splicing

U12-dependent introns containing alterations of the 3' splice site AC dinucleotide or alterations in the spacing between the branch site and the 3' splice site were examined for their effects on splice site selection in vivo and in vitro. Using an intron with a 5' splice site AU dinucleotide, any nucleotide could serve as the 3'-terminal nucleotide, although a C residue was most active, while a U residue was least active. The penultimate A residue, by contrast, was essential for 3' splice site function. A branch site-to-3' splice site spacing of less than 10 or more than 20 nucleotides strongly activated alternative 3' splice sites. A strong preference for a spacing of about 12 nucleotides was observed. The combined in vivo and in vitro results suggest that the branch site is recognized in the absence of an active 3' splice site but that formation of the prespliceosomal complex A requires an active 3' splice site. Furthermore, the U12-type spliceosome appears to be unable to scan for a distal 3' splice site.

The intramolecular stem-loop structure of U6 snRNA can functionally replace the U6atac snRNA stem-loop

The U6 spliceosomal snRNA forms an intramolecular stem-loop structure during spliceosome assembly that is required for splicing and is proposed to be at or near the catalytic center of the spliceosome. U6atac snRNA, the analog of U6 snRNA used in the U12-dependent splicing of the minor class of spliceosomal introns, contains a similar stem-loop whose structure but not sequence is conserved between humans and plants. To determine if the U6 and U6atac stem-loops are functionally analogous, the stem-loops from human and budding yeast U6 snRNAs were substituted for the U6atac snRNA structure and tested in an in vivo genetic suppression assay. Both chimeric U6/U6atac snRNA constructs were active for splicing in vivo. In contrast, several mutations of the native U6atac stem-loop that either delete putatively unpaired residues or disrupt the putative stem regions were inactive for splicing. Compensatory mutations that are expected to restore base pairing within the stem regions restored splicing activity. However, other mutants that retained base pairing potential were inactive, suggesting that functional groups within the stem regions may contribute to function. These results show that the U6atac snRNA stem-loop structure is required for in vivo splicing within the U12-dependent spliceosome and that its role is likely to be similar to that of the U6 snRNA intramolecular stem-loop.

Conservation of functional features of U6atac and U12 snRNAs between vertebrates and higher plants

Splicing of U12-dependent introns requires the function of U11, U12, U6atac, U4atac, and U5 snRNAs. Recent studies have suggested that U6atac and U12 snRNAs interact extensively with each other, as well as with the pre-mRNA by Watson-Crick base pairing. The overall structure and many of the sequences are very similar to the highly conserved analogous regions of U6 and U2 snRNAs. We have identified the homologs of U6atac and U12 snRNAs in the plant Arabidopsis thaliana. These snRNAs are significantly diverged from human, showing overall identities of 65% for U6atac and 55% for U12 snRNA. However, there is almost complete conservation of the sequences and structures that are implicated in splicing. The sequence of plant U6atac snRNA shows complete conservation of the nucleotides that base pair to the 5' splice site sequences of U12-dependent introns in human. The immediately adjacent AGAGA sequence, which is found in human U6atac and all U6 snRNAs, is also conserved. High conservation is also observed in the sequences of U6atac and U12 that are believed to base pair with each other. The intramolecular U6atac stem-loop structure immediately adjacent to the U12 interaction region differs from the human sequence in 9 out of 21 positions. Most of these differences are in base pairing regions with compensatory changes occurring across the stem. To show that this stem-loop was functional, it was transplanted into a human suppressor U6atac snRNA expression construct. This chimeric snRNA was inactive in vivo but could be rescued by coexpression of a U4atac snRNA expression construct containing compensatory mutations that restored base pairing to the chimeric U6atac snRNA. These data show that base pairing of U4atac snRNA to U6atac snRNA has a required role in vivo and that the plant U6atac intramolecular stem-loop is the functional analog of the human sequence.

Evolutionary fates and origins of U12-type introns

U2-type and U12-type introns are spliced by distinct spliceosomes in eukaryotic nuclei. A classification method was devised to distinguish these two types of introns based on splice site sequence properties and was used to identify 56 different genes containing U12-type introns in available genomic sequences. U12-type introns occur with consistently low frequency in diverse eukaryotic taxa but have almost certainly been lost from C. elegans. Comparisons with available homologous sequences demonstrate subtype switching of U12 introns between termini of AT-AC and GT-AG as well as conversion of introns from U12-type to U2-type and provide evidence for a fission/fusion model in which the two splicing systems evolved in separate lineages that were fused in a eukaryotic progenitor.

The two steps of group II intron self-splicing are mechanistically distinguishable

The two transesterification reactions catalyzed by self-splicing group II introns take place in either two active sites or two conformations of a single active site involving rearrangements of the positions of the reacting groups. We have investigated the effects on the rates of the chemical steps of the two reactions due to sulfur substitution of nonbridging oxygens at both the 5' and 3' splice sites as well as the deoxyribose substitution of the ribose 2' hydroxyl group at the 5' splice site. The data suggest that the two active sites differ in their interactions with several of these groups. Specifically, sulfur substitution of the pro-Sp nonbridging oxygen at the 5' splice site reduces the chemical rate of the step one branching reaction by at least 250-fold, whereas substitution of the pro-Sp oxygen at the 3' splice site has only a 4.5-fold effect on the chemical rate of step two. Previous work demonstrated that the Rp phosphorothioate substitutions at both the 5' and 3' splice sites reduced the rate of both steps of splicing to an undetectable level. These results suggest that either two distinct active sites catalyze the two steps or that more significant alterations must be made in a single bifunctional active site to accommodate the two different reactions.

Base pairing with U6atac snRNA is required for 5' splice site activation of U12-dependent introns in vivo

The minor U12-dependent class of eukaryotic nuclear pre-mRNA introns is spliced by a distinct spliceosomal mechanism that requires the function of U11, U12, U5, U4atac, and U6atac snRNAs. Previous work has shown that U11 snRNA plays a role similar to U1 snRNA in the major class spliceosome by base pairing to the conserved 5' splice site sequence. Here we show that U6atac snRNA also base pairs to the 5' splice site in a manner analogous to that of U6 snRNA in the major class spliceosome. We show that splicing defective mutants of the 5' splice site can be activated for splicing in vivo by the coexpression of compensatory U6atac snRNA mutants. In some cases, maximal restoration of splicing required the coexpression of compensatory U11 snRNA mutants. The allelic specificity of mutant phenotype suppression is consistent with Watson-Crick base pairing between the pre-mRNA and the snRNAs. These results provide support for a model of the RNA-RNA interactions at the core of the U12-dependent spliceosome that is strikingly similar to that of the major class U2-dependent spliceosome.

Terminal intron dinucleotide sequences do not distinguish between U2- and U12-dependent introns

Two types of eukaryotic nuclear introns are known: the common U2-dependent class with /GU and AG/ terminal intron dinucleotides, and the rare U12-dependent class with /AU and AC/ termini. Here we show that the U12-dependent splicing system can splice introns with /GU and AG/ termini and that such introns occur naturally. Further, U2-dependent introns with /AU and AC/termini also occur naturally and are evolutionarily conserved. Thus, the sequence of the terminal dinucleotides does not determine which spliceosomal system removes an intron. Rather, the four classes of introns described here can be sorted into two mechanistic classes (U2- or U12-dependent) by inspection of the complete set of conserved splice site sequences.

U11 snRNA interacts in vivo with the 5' splice site of U12-dependent (AU-AC) pre-mRNA introns

A notable feature of the newly described U12 snRNA-dependent class of eukaryotic nuclear pre-mRNA introns is the highly conserved 8-nt 5' splice site sequence. This sequence is virtually invariant in all known members of this class from plants to mammals. Based on sequence complementarity between this sequence and the 5' end of the U11 snRNA, we proposed that U11 snRNP may play a role in identifying and/or activating the 5' splice site for splicing. Here we show that mutations of the conserved 5' splice site sequence of a U12-dependent intron severely reduce correct splicing in vivo and that compensatory mutations in U11 snRNA can suppress the effects of the 5' splice site mutations to varying extents. This provides evidence for a required interaction between U11 snRNA and the 5' splice site sequence involving Watson-Crick base pairing. This data, in addition to a report that U11 snRNP is bound transiently to the U12-dependent spliceosome, suggests that U11 snRNP is the analogue of U1 snRNP in splicing this rare class of introns.

Requirement of U12 snRNA for in vivo splicing of a minor class of eukaryotic nuclear pre-mRNA introns

A conserved sequence element in a minor class of eukaryotic pre-messenger RNA (pre-mRNA) introns was previously proposed to base pair with a complementary sequence in the U12 small nuclear RNA (snRNA) in a manner analogous to the pairing of US snRNA with the branch site sequence of the major class of introns. Here, mutations generated in this conserved sequence element block the splicing of a member of this minor intron class in vivo. The block was relieved by coexpression of a U12 snRNA containing compensatory mutations that restore the proposed base pairing interaction. These results show that this minor class of pre-mRNA introns is a distinct class existing alongside the major class of introns in animal genomes, and these results also establish an in vivo function for U12 snRNA.

Stereochemical selectivity of group II intron splicing, reverse splicing, and hydrolysis reactions

We have previously shown, using phosphorothioate substitutions at splice site, that both transesterification steps of group II intron self-splicing proceed, by stereochemical inversion, with an Sp but not an Rp phosphorothioate. Under alternative reaction conditions or with various intron fragments, group II introns can splice following hydrolysis at the 5' splice site and can also hydrolyze the bond between spliced exons (the spliced-exon reopening reaction). In this study, we have determined the stereochemical specificities of all of the major model hydrolytic reactions carried out by the aI5 gamma intron from Saccharomyces cerevisiae mitochondria. For all substrates containing exon 1 and most of the intron, the stereospecificity of hydrolysis is the same as for the step 1 transesterification reaction. In contrast, the spliced-exon reopening reaction proceeds with an Rp but not an Sp phosphorothioate at the scissile bond, as does true reverse splicing. Thus, by stereochemistry, this reaction appears to be related to the reverse of step 2 of self-splicing. Finally, a substrate RNA that contains the first exon and nine nucleotides of the intron, when reacted with the intron ribozyme, releases the first exon regardless of the configuration of the phosphorothioate at the 5' splice site, suggesting that this substrate can be cleaved by either the step 1 or the step 2 reaction site. Our findings clarify the relationships of these model reactions to the transesterification reactions of the intact self-splicing system and permit new studies to be interpreted more rigorously.

The stereochemical course of group II intron self-splicing

The stereochemical specificities and reaction courses for both self-splicing steps of a group II intron have been determined by phosphorothioate substitution at the 5' and 3' splice site phosphodiester bonds. Both steps of the splicing reaction proceeded with a phosphorothioate in the Sp configuration but were blocked by the Rp diastereomer. Both steps also proceeded with inversion of stereochemical configuration around phosphorus, consistent with a concerted transesterification reaction. These results are identical to those found for nuclear precursor mRNA (pre-mRNA) splicing and provide support for the hypothesis that group II introns and nuclear pre-mRNA introns share a common evolutionary history.

Conserved sequences in a class of rare eukaryotic nuclear introns with non-consensus splice sites

Eukaryotic nuclear genomes contain a rare class of pre-mRNA introns with consensus sequence features that differ markedly from most pre-mRNA introns. Four genes have so far been identified that contain one copy each of this rare intron class in addition to several standard introns. These introns and homologous introns from several species were compared to identify conserved sequence elements and to establish consensus sequences for these elements. The only well-conserved elements are found at the 5' and 3' ends of the introns. The 5' splice site sequence is ATATCCTT beginning with the first nucleotide of the intron and is invariant in the introns examined to date. The 3' splice site consensus sequence is YCCAC ending at the last nucleotide of the intron. An almost invariant sequence of TCCTTAAC is also found near the 3' end of the intron (the 3' upstream element). The length of the introns varies between 95 and 2940 nucleotides. The sequence organization of these introns suggests that they represent a variant class of pre-mRNA introns that might be spliced via a spliceosome mechanism employing factors distinct from those used by other pre-mRNA introns. A search of small nuclear RNA (snRNA) sequences for regions complementary to the conserved elements of this rare class of introns found a strong match between U12 snRNA and the 3' upstream element and a weaker match between U11 snRNA and the 5' splice site sequence.

The stereochemical course of the first step of pre-mRNA splicing

We have determined the effects on splicing of sulfur substitution of the non-bridging oxygens in the phosphodiester bond at the 5' splice site of a pre-mRNA intron. Pre-mRNAs containing stereochemically pure Rp and Sp phosphorothioate isomers were produced by ligation of a chemically synthesized modified RNA oligonucleotide to enzymatically synthesized RAs. When these modified pre-mRNA substrates were tested for in vitro splicing activity in a HeLa cell nuclear extract system, the RNA with the Rp diastereomeric phosphorothioate was not spliced while the Sp diastereomeric RNA spliced readily. The sulfur-containing branched trinucleotide was purified from the splicing reaction of the Sp RNA and analyzed by cleavage with a stereospecific nuclease. The results showed that the Sp phosphorothioate was inverted during the splicing reaction to the Rp configuration; a finding previously obtained for a Group I self-splicing RNA. This inversion of configuration is consistent with a transesterification mechanism for pre-mRNA splicing. The lack of splicing of the Rp modified RNA also suggests that the pro-Rp oxygen at the 5' splice site is involved in a critical chemical contact in the splicing mechanism. Additionally, we have found that the HeLa cell RNA debranching enzyme is inactive on branches containing an Rp phosphorothioate.

Phosphorothioate substitution identifies phosphate groups important for pre-mRNA splicing

Substitution of pre-mRNA in vitro splicing substrates with alpha-phosphorothioate ribonucleotide analogs has multiple effects on the processes of spliceosome formation and splicing. A major effect of substitution is on the splicing cleavage/ligation reactions. Substitution at the 5' splice junction blocks the first cleavage/ligation reaction while substitution at the 3' splice junction blocks the second cleavage/ligation reaction. A second effect of phosphorothioate substitution is the inhibition of spliceosome formation. A substitution/interference assay was used to determine positions where substitution inhibits spliceosome formation or splicing. Substitution in the 3' splice site polypyrimidine tract was found to inhibit spliceosome formation and splicing. This effect was enhanced with multiple substitutions in the region. No sites of substitution within the exons were found which affected spliceosome formation or splicing.

Hydroxyl radical "footprinting" of RNA: application to pre-mRNA splicing complexes

We present an adaptation of the hydroxyl radical DNA "footprinting" technique that permits high-resolution mapping of protected regions of RNA. Hydroxyl radical cleaves RNA independently of base sequence and secondary structure of the RNAs examined and allows resolution of protected regions at the single nucleotide level. By using this technique, we show that several regions of the 3' splice site of mRNA precursors are protected during the formation of splicing-specific ribonucleoprotein complexes in an in vitro splicing system. These regions include the 3' intron/exon junction and a portion of the adjacent exon, the polypyrimidine tract, and the site of branch formation. These protections appear to be due to splicing specific complexes since their formation is sensitive to point mutations at crucial residues and requires ATP and incubation. The formation of these protected regions is independent of the presence of a 5' splice site.

Sequence requirements for splicing of higher eukaryotic nuclear pre-mRNA

We determined the effect on splicing of 24 point mutations in the 5' and 3' splice region of the large rabbit beta-globin intron. In vitro, 3' AG mutations drastically reduce 5' cleavage and abolish splicing. In vivo, the same mutations elicit efficient splicing at a cryptic, rather than the correct, 3' splice site. In vitro, mutations at all but 2 positions of the consensus 5' splice region impair correct splicing and promote joining of exon 1 to exon 3. In vivo, the same mutations show no effect, except for those converting 5' GT to AT or GA, which cause accumulation of lariat intermediate in vitro and in vivo. We conclude that the 5' GT need not be conserved for 5' cleavage and that it plays an important role in cleavage and exon joining at the 3' splice site.

Nonconsensus branch-site sequences in the in vitro splicing of transcripts of mutant rabbit beta-globin genes

Mutants of the rabbit beta-globin gene lacking the natural site of branch formation in the second intervening sequence have been analyzed for in vitro splicing activity. RNAs transcribed from these mutants were spliced, via lariat formation, at a reduced rate compared to wild-type RNA. The sites of branch formation were mapped by direct RNA analysis and primer-extension analysis. The sequences at the branch sites in the three mutants examined did not conform to the previously determined consensus sequence, nor were the 5' splice sites and branch sites complementary.

Trans splicing of mRNA precursors in vitro

Two exon segments from two separate RNA molecules can be joined in a trans splicing process. In trans splicing reactions, an RNA molecule containing an exon, a 5' splice site, and adjacent intron sequences was mixed with an RNA molecule containing an exon, a 3' splice site, and adjacent intron sequences. The efficiency of trans splicing of these two RNAs increased if the two termini of the intervening sequences were paired in a short RNA duplex. However, trans splicing of two RNA molecules with no significant complementarity was also observed. These results strongly suggest that significant secondary structures within intervening sequences could affect the splicing of flanking exons. Similarly, RNAs that are complementary to segments within the intervening sequences could potentially regulate the selection of splice sites. Finally, some organisms might use trans splicing to distribute a single exon to many different mRNAs.

Bioluminescence screening for bacteriuria

A rapid screening test (45 min) for bacteriuria was evaluated in 1,000 clinical urine specimens. The test procedure is based upon firefly luciferase analysis of bacterial ATP and uses the Lumac kit and Lumac M2010 Biocounter (3M Co., St. Paul, Minn.). The procedure allows for removal and destruction of nonbacterial ATP and subsequent analysis of bacterial ATP by firefly luciferase with a single photon counter. Results, expressed in relative light units, were compared with actual CFU by the calibrated loop technique. Sensitivities and specificities were calculated separately for clean-catch midstream specimens and for urines obtained by catheterization. The sensitivity for 719 clean-catch specimens with a Lumac cutoff of greater than or equal to 500 relative light units, representing greater than or equal to 10(5) CFU/ml, was 93%. The sensitivity for 281 catheterized specimens with a Lumac cutoff of greater than or equal to 200 relative light units, representing greater than or equal to 10(4) CFU/ml, was 95%. There were 19 false-negative results in the 1,000 specimens tested; more than 50% of these were contaminated cultures and were not considered significant in determining bacteriuria. In conclusion, the Lumac bioluminescence assay is a reliable, rapid bacteriuria screening technique with the potential of reducing the laboratory cost and for reducing the turnaround time in processing negative urine cultures.

Characterization of the branch site in lariat RNAs produced by splicing of mRNA precursors

The branch site of lariat RNAs produced during the splicing of the first two late leader exons of adenovirus-2 has a structure of (formula; see text) There is a distinct complementarity between sequences preceding the adenosine at the branch site and the 5' terminus of the intervening sequence.