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

B-cell immune deficiency in twin sisters expands the phenotype of MOPDI

Microcephalic osteodysplastic primordial dwarfism type I (MOPDI) is a very rare and severe autosomal recessive disorder characterized by marked intrauterine growth retardation, skeletal dysplasia, microcephaly and brain malformations. MOPDI is caused by biallelic mutations in RNU4ATAC, a non-coding gene involved in U12-type splicing of 1% of the introns in the genome, which are recognized by their specific splicing consensus sequences. Here, we describe a unique observation of immunodeficiency in twin sisters with mild MOPDI, who harbor a novel n.108_126del mutation, encompassing part of the U4atac snRNA 3' stem-loop and Sm protein binding site, and the previously reported n.111G>A mutation. Interestingly, both twin sisters show mild B-cell anomalies, including low naive B-cell counts and increased memory B-cell and plasmablasts counts, suggesting partial and transitory blockage of B-cell maturation and/or excessive activation of naive B-cells. Hence, the localization of a mutation in stem II of U4atac snRNA, as observed in another RNU4ATAC-opathy with immunodeficiency, that is, Roifman syndrome (RFMN), is not required for the occurrence of an immune deficiency. Finally, we emphasize the importance of considering immunodeficiency in MOPDI management to reduce the risk of serious infectious episodes.

Delineating the phenotype of RNU4ATAC-related spliceosomopathy

Biallelic pathogenic variants in RNU4ATAC cause microcephalic osteodysplastic primordial dwarfism type I (MOPD1), Roifman syndrome (RS) and Lowry-Wood syndrome (LWS). These conditions demonstrate significant phenotypic heterogeneity yet have overlapping features. Although historically described as discrete conditions they appear to represent a phenotypic spectrum with clinical features not always aligning with diagnostic categories. Clinical variability and ambiguity in diagnostic criteria exist among each disorder. Here we report an individual with a novel genotype and phenotype spanning all three disorders, expanding the phenotypic spectrum of RNU4ATAC-related spliceosomeopathies.

Mutations in the non-coding RNU4ATAC gene affect the homeostasis and function of the Integrator complex

Various genetic diseases associated with microcephaly and developmental defects are due to pathogenic variants in the U4atac small nuclear RNA (snRNA), a component of the minor spliceosome essential for the removal of U12-type introns from eukaryotic mRNAs. While it has been shown that a few RNU4ATAC mutations result in impaired binding of essential protein components, the molecular defects of the vast majority of variants are still unknown. Here, we used lymphoblastoid cells derived from RNU4ATAC compound heterozygous (g.108_126del;g.111G>A) twin patients with MOPD1 phenotypes to analyze the molecular consequences of the mutations on small nuclear ribonucleoproteins (snRNPs) formation and on splicing. We found that the U4atac108_126del mutant is unstable and that the U4atac111G>A mutant as well as the minor di- and tri-snRNPs are present at reduced levels. Our results also reveal the existence of 3'-extended snRNA transcripts in patients' cells. Moreover, we show that the mutant cells have alterations in splicing of INTS7 and INTS10 minor introns, contain lower levels of the INTS7 and INTS10 proteins and display changes in the assembly of Integrator subunits. Altogether, our results show that compound heterozygous g.108_126del;g.111G>A mutations induce splicing defects and affect the homeostasis and function of the Integrator complex.

Disruption of exon-bridging interactions between the minor and major spliceosomes results in alternative splicing around minor introns

Vertebrate genomes contain major (>99.5%) and minor (<0.5%) introns that are spliced by the major and minor spliceosomes, respectively. Major intron splicing follows the exon-definition model, whereby major spliceosome components first assemble across exons. However, since most genes with minor introns predominately consist of major introns, formation of exon-definition complexes in these genes would require interaction between the major and minor spliceosomes. Here, we report that minor spliceosome protein U11-59K binds to the major spliceosome U2AF complex, thereby supporting a model in which the minor spliceosome interacts with the major spliceosome across an exon to regulate the splicing of minor introns. Inhibition of minor spliceosome snRNAs and U11-59K disrupted exon-bridging interactions, leading to exon skipping by the major spliceosome. The resulting aberrant isoforms contained a premature stop codon, yet were not subjected to nonsense-mediated decay, but rather bound to polysomes. Importantly, we detected elevated levels of these alternatively spliced transcripts in individuals with minor spliceosome-related diseases such as Roifman syndrome, Lowry–Wood syndrome and early-onset cerebellar ataxia. In all, we report that the minor spliceosome informs splicing by the major spliceosome through exon-definition interactions and show that minor spliceosome inhibition results in aberrant alternative splicing in disease.

Clinical interpretation of variants identified in RNU4ATAC, a non-coding spliceosomal gene

Biallelic variants in RNU4ATAC, a non-coding gene transcribed into the minor spliceosome component U4atac snRNA, are responsible for three rare recessive developmental diseases, namely Taybi-Linder/MOPD1, Roifman and Lowry-Wood syndromes. Next-generation sequencing of clinically heterogeneous cohorts (children with either a suspected genetic disorder or a congenital microcephaly) recently identified mutations in this gene, illustrating how profoundly these technologies are modifying genetic testing and assessment. As RNU4ATAC has a single non-coding exon, the bioinformatic prediction algorithms assessing the effect of sequence variants on splicing or protein function are irrelevant, which makes variant interpretation challenging to molecular diagnostic laboratories. In order to facilitate and improve clinical diagnostic assessment and genetic counseling, we present i) an update of the previously reported RNU4ATAC mutations and an analysis of the genetic variations affecting this gene using the Genome Aggregation Database (gnomAD) resource; ii) the pathogenicity prediction performances of scores computed based on an RNA structure prediction tool and of those produced by the Combined Annotation Dependent Depletion tool for the 285 RNU4ATAC variants identified in patients or in large-scale sequencing projects; iii) a method, based on a cellular assay, that allows to measure the effect of RNU4ATAC variants on splicing efficiency of a minor (U12-type) reporter intron. Lastly, the concordance of bioinformatic predictions and cellular assay results was investigated.

Extending the critical regions for mutations in the non-coding gene RNU4ATAC in another patient with Roifman Syndrome

Compound heterozygosity of a previously described pathogenic variant and a second novel nucleotide substitution (NR_023343.1:n.116A>C) affecting a highly conserved nucleotide in the noncoding RNU4ATAC gene could be identified in a patient with overlapping features of Roifman Syndrome. These data extend the spectrum of pathogenic variants in RNU4ATAC.

Mutations of RNA splicing factors in hematological malignancies

Systematic large-scale cancer genomic studies have produced numerous significant findings. These studies have not only revealed new cancer-promoting genes, but they also have identified cancer-promoting functions of previously known "housekeeping" genes. These studies have identified numerous mutations in genes which play a fundamental role in nuclear precursor mRNA splicing. Somatic mutations and copy number variation in many of the splicing factors which participate in the formation of multiple spliceosomal complexes appear to play a role in many cancers and in particular in myelodysplastic syndromes (MDS). Mutated proteins seem to interfere with the recognition of the authentic splice sites (SS) leading to utilization of suboptimal alternative splicing sites generating aberrantly spliced mRNA isoforms. This short review is focusing on the function of the splice factors involved in the formation of splicing complexes and potential mechanisms which affect usage of the authentic splice site recognition.

U6atac snRNA stem-loop interacts with U12 p65 RNA binding protein and is functionally interchangeable with the U12 apical stem-loop III

Formation of catalytic core of the U12-dependent spliceosome involves U6atac and U12 interaction with the 5′ splice site and branch site regions of a U12-dependent intron, respectively. Beyond the formation of intermolecular helix I region between U6atac and U12 snRNAs, several other regions within these RNA molecules are predicted to form stem-loop structures. Our previous work demonstrated that the 3′ stem-loop region of U6atac snRNA contains a U12-dependent spliceosome-specific targeting activity. Here, we show a detailed structure-function analysis and requirement of a substructure of U6atac 3′ stem-loop in U12-dependent in vivo splicing. We show that the C-terminal RNA recognition motif of p65, a U12 snRNA binding protein, also binds to the distal 3′ stem-loop of U6atac. By using a binary splice site mutation suppressor assay we demonstrate that p65 protein-binding apical stem-loop of U12 snRNA can be replaced by this U6atac distal 3′ stem-loop. Furthermore, we tested the compatibility of the U6atac 3′ end from phylogenetically distant species in a human U6atac background, to establish the evolutionary relatedness of these structures and in vivo function. In summary, we demonstrate that RNA-RNA and RNA-protein interactions in the minor spliceosome are highly plastic as compared to the major spliceosome.

Compound heterozygous mutations in the noncoding RNU4ATAC cause Roifman Syndrome by disrupting minor intron splicing

Roifman Syndrome is a rare congenital disorder characterized by growth retardation, cognitive delay, spondyloepiphyseal dysplasia and antibody deficiency. Here we utilize whole-genome sequencing of Roifman Syndrome patients to reveal compound heterozygous rare variants that disrupt highly conserved positions of the RNU4ATAC small nuclear RNA gene, a minor spliceosome component that is essential for minor intron splicing. Targeted sequencing confirms allele segregation in six cases from four unrelated families. RNU4ATAC rare variants have been recently reported to cause microcephalic osteodysplastic primordial dwarfism, type I (MOPD1), whose phenotype is distinct from Roifman Syndrome. Strikingly, all six of the Roifman Syndrome cases have one variant that overlaps MOPD1-implicated structural elements, while the other variant overlaps a highly conserved structural element not previously implicated in disease. RNA-seq analysis confirms extensive and specific defects of minor intron splicing. Available allele frequency data suggest that recessive genetic disorders caused by RNU4ATAC rare variants may be more prevalent than previously reported.

Roifman Syndrome is a rare disorder whose disease manifestations include growth retardation, spondyloepiphyseal dysplasia and immunodeficiency. Here, the authors use whole-genome sequencing to discover that rare compound heterozygous variants disrupting the small nuclear RNA gene RNU4ATAC cause Roifman Syndrome.

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.

In vitro reconstitution of yeast splicing with U4 snRNA reveals multiple roles for the 3' stem-loop

U4 small nuclear RNA (snRNA) plays a fundamental role in the process of premessenger RNA splicing, yet many questions remain regarding the location, interactions, and roles of its functional domains. To address some of these questions, we developed the first in vitro reconstitution system for yeast U4 small nuclear ribonucleoproteins (snRNPs). We used this system to examine the functional domains of U4 by measuring reconstitution of splicing, U4/U6 base-pairing, and triple-snRNP formation. In contrast to previous work in human extracts and Xenopus oocytes, we found that the 3' stem-loop of U4 is necessary for efficient base-pairing with U6. In particular, the loop is sensitive to changes in both length and sequence. Intriguingly, a number of mutations that we tested resulted in more stable interactions with U6 than wild-type U4. Nevertheless, each of these mutants was impaired in its ability to support splicing, indicating that these regions of U4 have functions subsequent to base pair formation with U6. Our data suggest that one such function is likely to be in tri-snRNP formation, when U5 joins the U4/U6 di-snRNP. We have identified two regions, the upper stem of the 3' stem-loop and the central domain, that promote tri-snRNP formation. In addition, the loop of the 3' stem-loop promotes di-snRNP formation, while the central domain and the 3'-terminal domain appear to antagonize di-snRNP formation.

Functionally important structural elements of U12 snRNA

U12 snRNA is analogous to U2 snRNA of the U2-dependent spliceosome and is essential for the splicing of U12-dependent introns in metazoan cells. The essential region of U12 snRNA, which base pairs to the branch site of minor class introns is well characterized. However, other regions which are outside of the branch site base pairing region are not yet characterized and the requirement of these structures in U12-dependent splicing is not clear. U12 snRNA is predicted to form an intricate secondary structure containing several stem-loops and single-stranded regions. Using a previously characterized branch site genetic suppression assay, we generated second-site mutations in the suppressor U12 snRNA to investigate the in vivo requirement of structural elements in U12-dependent splicing. Our results show that stem-loop IIa is essential and required for in vivo splicing. Interestingly, an evolutionarily conserved stem-loop IIb is dispensable for splicing. We also show that stem-loop III, which binds to a p65 RNA binding protein of the U11-U12 di.snRNP complex, is essential for in vivo splicing. The data validate the existence of proposed stem-loops of U12 snRNA and provide experimental support for individual secondary structures.

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.

Association of TALS developmental disorder with defect in minor splicing component U4atac snRNA

The spliceosome, a ribonucleoprotein complex that includes proteins and small nuclear RNAs (snRNAs), catalyzes RNA splicing through intron excision and exon ligation to produce mature messenger RNAs, which, in turn serve as templates for protein translation. We identified four point mutations in the U4atac snRNA component of the minor spliceosome in patients with brain and bone malformations and unexplained postnatal death [microcephalic osteodysplastic primordial dwarfism type 1 (MOPD 1) or Taybi-Linder syndrome (TALS); Mendelian Inheritance in Man ID no. 210710]. Expression of a subgroup of genes, possibly linked to the disease phenotype, and minor intron splicing were affected in cell lines derived from TALS patients. Our findings demonstrate a crucial role of the minor spliceosome component U4atac snRNA in early human development and postnatal survival.

RNA-RNA interaction prediction based on multiple sequence alignments

MOTIVATION

Many computerized methods for RNA-RNA interaction structure prediction have been developed. Recently, O(N(6)) time and O(N(4)) space dynamic programming algorithms have become available that compute the partition function of RNA-RNA interaction complexes. However, few of these methods incorporate the knowledge concerning related sequences, thus relevant evolutionary information is often neglected from the structure determination. Therefore, it is of considerable practical interest to introduce a method taking into consideration both: thermodynamic stability as well as sequence/structure covariation.

RESULTS

We present the a priori folding algorithm ripalign, whose input consists of two (given) multiple sequence alignments (MSA). ripalign outputs (i) the partition function, (ii) base pairing probabilities, (iii) hybrid probabilities and (iv) a set of Boltzmann-sampled suboptimal structures consisting of canonical joint structures that are compatible to the alignments. Compared to the single sequence-pair folding algorithm rip, ripalign requires negligible additional memory resource but offers much better sensitivity and specificity, once alignments of suitable quality are given. ripalign additionally allows to incorporate structure constraints as input parameters.

AVAILABILITY

The algorithm described here is implemented in C as part of the rip package.

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.

Evolution of spliceosomal snRNA genes in metazoan animals

While studies of the evolutionary histories of protein families are commonplace, little is known on noncoding RNAs beyond microRNAs and some snoRNAs. Here we investigate in detail the evolutionary history of the nine spliceosomal snRNA families (U1, U2, U4, U5, U6, U11, U12, U4atac, and U6atac) across the completely or partially sequenced genomes of metazoan animals. Representatives of the five major spliceosomal snRNAs were found in all genomes. None of the minor splicesomal snRNAs were detected in nematodes or in the shotgun traces of Oikopleura dioica, while in all other animal genomes at most one of them is missing. Although snRNAs are present in multiple copies in most genomes, distinguishable paralogue groups are not stable over long evolutionary times, although they appear independently in several clades. In general, animal snRNA secondary structures are highly conserved, albeit, in particular, U11 and U12 in insects exhibit dramatic variations. An analysis of genomic context of snRNAs reveals that they behave like mobile elements, exhibiting very little syntenic conservation.

U12-dependent intron splicing in plants

U12-dependent (U12) introns have persisted in the genomes of plants since the ancestral divergence between plants and metazoans. These introns, which are rare, are found in a range of genes that include essential functions in DNA replication and RNA metabolism and are implicated in regulating the expression of their host genes. U12 introns are removed from pre-mRNAs by a U12 intron-specific spliceosome. Although this spliceosome shares many properties with the more abundant U2-dependent (U2) intron spliceosome, four of the five small nuclear RNAs (snRNAs) required for splicing are different and specific for the unique splicing of U12 introns. Evidence in plants so far indicates that splicing signals of plant U12 introns and their splicing machinery are similar to U12 intron splicing in other eukaryotes. In addition to the high conservation of splicing signals, plant U12 introns also retain unique characteristic features of plant U2 introns, such as UA-richness, which suggests a requirement for plant-specific components for both the U2 and U12 splicing reaction. This chapter compares U12 and U2 splicing and reviews what is known about plant U12 introns and their possible role in gene expression.

Free energy landscapes of RNA/RNA complexes: with applications to snRNA complexes in spliceosomes

We develop a statistical mechanical model for RNA/RNA complexes with both intramolecular and intermolecular interactions. As an application of the model, we compute the free energy landscapes, which give the full distribution for all the possible conformations, for U4/U6 and U2/U6 in major spliceosome and U4atac/U6atac and U12/U6atac in minor spliceosome. Different snRNA experiments found contrasting structures, our free energy landscape theory shows why these structures emerge and how they compete with each other. For yeast U2/U6, the model predicts that the two distinct experimental structures, the four-helix junction structure and the helix Ib-containing structure, can actually coexist and specifically compete with each other. In addition, the energy landscapes suggest possible mechanisms for the conformational switches in splicing. For instance, our calculation shows that coaxial stacking is essential for stabilizing the four-helix junction in yeast U2/U6. Therefore, inhibition of the coaxial stacking possibly by protein-binding may activate the conformational switch from the four-helix junction to the helix Ib-containing structure. Moreover, the change of the energy landscape shape gives information about the conformational changes. We find multiple (native-like and misfolded) intermediates formed through base-pairing rearrangements in snRNA complexes. For example, the unfolding of the U2/U6 undergoes a transition to a misfolded state which is functional, while in the unfolding of U12/U6atac, the functional helix Ib is found to be the last one to unfold and is thus the most stable structural component. Furthermore, the energy landscape gives the stabilities of all the possible (functional) intermediates and such information is directly related to splicing efficiency.

The ASRG database: identification and survey of Arabidopsis thaliana genes involved in pre-mRNA splicing

The database of Arabidopsis splicing related genes includes classification of genes encoding snRNAs and other splicing related proteins, together with information on gene structure, alternative splicing, gene duplications and phylogenetic relationships.

A total of 74 small nuclear RNA (snRNA) genes and 395 genes encoding splicing-related proteins were identified in the Arabidopsis genome by sequence comparison and motif searches, including the previously elusive U4atac snRNA gene. Most of the genes have not been studied experimentally. Classification of these genes and detailed information on gene structure, alternative splicing, gene duplications and phylogenetic relationships are made accessible as a comprehensive database of Arabidopsis Splicing Related Genes (ASRG) on our website.

Recycling of the U12-type spliceosome requires p110, a component of the U6atac snRNP

U12-dependent introns are spliced by the so-called minor spliceosome, requiring the U11, U12, and U4atac/U6atac snRNPs in addition to the U5 snRNP. We have recently identified U6-p110 (SART3) as a novel human recycling factor that is related to the yeast splicing factor Prp24. U6-p110 transiently associates with the U6 and U4/U6 snRNPs during the spliceosome cycle, regenerating functional U4/U6 snRNPs from singular U4 and U6 snRNPs. Here we investigated the involvement of U6-p110 in recycling of the U4atac/U6atac snRNP. In contrast to the major U6 and U4/U6 snRNPs, p110 is primarily associated with the U6atac snRNP but is almost undetectable in the U4atac/U6atac snRNP. Since p110 does not occur in U5 snRNA-containing complexes, it appears to be transiently associated with U6atac during the cycle of the minor spliceosome. The p110 binding site was mapped to U6 nucleotides 38 to 57 and U6atac nucleotides 10 to 30, which are highly conserved between these two functionally related snRNAs. With a U12-dependent in vitro splicing system, we demonstrate that p110 is required for recycling of the U4atac/U6atac snRNP.

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.