Splice site m6A methylation prevents binding of U2AF35 to inhibit RNA splicing.

The N6-methyladenosine (m6A) RNA modification is used widely to alter the fate of mRNAs. Here we demonstrate that the C. elegans writer METT-10 (the ortholog of mouse METTL16) deposits an m6A mark on the 3' splice site (AG) of the S-adenosylmethionine (SAM) synthetase pre-mRNA, which inhibits its proper splicing and protein production. The mechanism is triggered by a rich diet and acts as an m6A-mediated switch to stop SAM production and regulate its homeostasis. Although the mammalian SAM synthetase pre-mRNA is not regulated via this mechanism, we show that splicing inhibition by 3' splice site m6A is conserved in mammals. The modification functions by physically preventing the essential splicing factor U2AF35 from recognizing the 3' splice site. We propose that use of splice-site m6A is an ancient mechanism for splicing regulation.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing.

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing.

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution.

High-throughput sequencing-based methods and their applications in the study of transcriptomes have revolutionized our understanding of alternative splicing. Networks of functionally coordinated and biologically important alternative splicing events continue to be discovered in an ever-increasing diversity of cell types in the context of physiologically normal and disease states. These studies have been complemented by efforts directed at defining sequence codes governing splicing and their cognate trans-acting factors, which have illuminated important combinatorial principles of regulation. Additional studies have revealed critical roles of position-dependent, multivalent protein-RNA interactions that direct splicing outcomes. Investigations of evolutionary changes in RNA binding proteins, splice variants, and associated cis elements have further shed light on the emergence, mechanisms, and functions of splicing networks. Progress in these areas has emphasized the need for a coordinated, community-based effort to systematically address the functions of individual splice variants associated with normal and disease biology.

Reevaluating the splice-altering variant in TECTA as a cause of nonsyndromic hearing loss DFNA8/12 by functional analysis of RNA.

PURPOSE: The aim of this study was to determine the genetic cause of early onset autosomal dominant hearing loss segregating in five-generation kindred of Chinese descent and provide preimplantation genetic testing (PGT)for them. METHODS: Clinical examination, pedigree analysis and exome sequencing were carried out on the family. Minigene-based splicing analysis, in vivo RNA analysis and protein structure prediction by molecular modeling were conducted on the candidate variant. PGT for the causative variation and chromosome aneuploidis based on SNP analysis has been used for avoidance of hearing loss in this family. RESULTS: All the affected individuals presented with moderate down-sloping hearing loss and whole-exome sequencing identified a novel splice-site variant c.5383+6T>A in the tested subjects within the TECTA locus. Genotyping of all the 32 family members confirmed segregation of this variant and the hearing loss phenotype in the extended family. Functional analysis of RNA and molecular modeling indicates that c.5383+6T>A is a pathogenic splice-site variant and should be considered as genetic cause of the hearing loss. Furthermore, a successful singleton pregnancy with no variation in TECTA c.5383+6 was established and a healthy male child was born by PGT. CONCLUSION: We have identified a novel variant c.5383+6T>A in TECTA ZA-ZP inter-domain, which could be attributable to the early-onset autosomal dominant hearing loss. The implications of our study are valuable in elucidating the disrupted RNA splicing and uncovering the genetic cause of hearing loss with TECTA pathogenic variants, as well as providing reproductive approaches to healthy offspring.

Examining the capacity of human U1 snRNA variants to facilitate pre-mRNA splicing.

The human U1 snRNA is encoded by a multigene family consisting of transcribed variants and defective pseudogenes. Many variant U1 (vU1) snRNAs have been demonstrated to not only be transcribed but also processed by the addition of a trimethylated guanosine cap, packaged into snRNPs, and assembled into spliceosomes; however, their capacity to facilitate pre-mRNA splicing has, so far, not been tested. A recent systematic analysis of the human snRNA genes identified 178 U1 snRNA genes that are present in the genome as either tandem arrays or single genes on multiple chromosomes. Of these, 15 were found to be expressed in human tissues and cell lines, although at significantly low levels from their endogenous loci, <0.001% of the canonical U1 snRNA. In this study, we found that placing the variants in the context of the regulatory elements of the RNU1-1 gene improves the expression of many variants to levels comparable to the canonical U1 snRNA. Application of a previously established HeLa cell-based minigene reporter assay to examine the capacity of the vU1 snRNAs to support pre-mRNA splicing revealed that even though the exogenously expressed variant snRNAs were enriched in the nucleus, only a few had a measurable effect on splicing.

Selected humanization of yeast U1 snRNP leads to global suppression of pre-mRNA splicing and mitochondrial dysfunction in the budding yeast.

The recognition of 5' splice site (5' ss) is one of the earliest steps of pre-mRNA splicing. To better understand the mechanism and regulation of 5' ss recognition, we selectively humanized components of the yeast U1 snRNP to reveal the function of these components in 5' ss recognition and splicing. We targeted U1C and Luc7, two proteins that interact with and stabilize the yeast U1 (yU1) snRNA and the 5' ss RNA duplex. We replaced the Zinc-Finger (ZnF) domain of yU1C with its human counterpart, which resulted in a cold-sensitive growth phenotype and moderate splicing defects. We next added an auxin-inducible degron to yLuc7 protein (to mimic the lack of Luc7Ls in human U1 snRNP) and found that Luc7-depleted yU1 snRNP resulted in the concomitant loss of PRP40 and Snu71 (two other essential yeast U1 snRNP proteins), and further biochemical analyses suggest a model of how these three proteins interact with each other in the U1 snRNP. The loss of these proteins resulted in a significant growth retardation accompanied by a global suppression of pre-mRNA splicing. The splicing suppression led to mitochondrial dysfunction as revealed by a release of Fe2+ into the growth medium and an induction of mitochondrial reactive oxygen species. Together, these observations indicate that the human U1C ZnF can substitute that of yeast, Luc7 is essential for the incorporation of the Luc7-Prp40-Snu71 trimer into yeast U1 snRNP, and splicing plays a major role in the regulation of mitochondrial function in yeast.

Genetic Modulation of RNA Splicing with a CRISPR-Guided Cytidine Deaminase.

RNA splicing is a critical mechanism by which to modify transcriptome, and its dysregulation is the underlying cause of many human diseases. It remains challenging, however, to genetically modulate a splicing event in its native context. Here, we demonstrate that a CRISPR-guided cytidine deaminase (i.e., targeted-AID mediated mutagenesis [TAM]) can efficiently modulate various forms of mRNA splicing. By converting invariant guanines to adenines at either 5' or 3' splice sites (SS), TAM induces exon skipping, activation of alternative SS, switching between mutually exclusive exons, or targeted intron retention. Conversely, TAM promotes downstream exon inclusion by mutating cytidines into thymines at the polypyrimidine tract. Applying this approach, we genetically restored the open reading frame and dystrophin function of a mutant DMD gene in patient-derived induced pluripotent stem cells (iPSCs). Thus, the CRISPR-guided cytidine deaminase provides a versatile genetic platform to modulate RNA splicing and to correct mutations associated with aberrant splicing in human diseases.

Cryptic splice mutation in the fumarate hydratase gene in patients with clinical manifestations of Hereditary Leiomyomatosis and Renal Cell Cancer.

Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is an autosomal dominant condition characterized by the development of cutaneous and uterine leiomyomas and risk for development of an aggressive form of papillary renal cell cancer. HLRCC is caused by germline inactivating pathogenic variants in the fumarate hydratase (FH) gene, which encodes the enzyme that catalyzes the interconversion of fumarate and L-malate. We utilized enzyme and protein mobility assays to evaluate the FH enzyme in a cohort of patients who showed clinical manifestations of HLRCC but were negative for known pathogenic FH gene variants. FH enzyme activity and protein levels were decreased by 50% or greater in three family members, despite normal FH mRNA expression levels as measured by quantitative PCR. Direct Nanopore RNA sequencing demonstrated 57 base pairs of retained intron sequence between exons 9 and 10 of polyadenylated FH mRNA in these patients, resulting in a truncated FH protein. Genomic sequencing revealed a heterozygous intronic alteration of the FH gene (chr1: 241498239 T/C) resulting in formation of a splice acceptor site near a polypyrimidine tract, and a uterine fibroid obtained from a patient showed loss of heterozygosity at this site. The same intronic FH variant was identified in an unrelated patient who also showed a clinical phenotype of HLRCC. These data demonstrate that careful clinical assessment as well as biochemical characterization of FH enzyme activity, protein expression, direct RNA sequencing, and genomic DNA sequencing of patient-derived cells can identify pathogenic variants outside of the protein coding regions of the FH gene.

Mouse nuclear RNAi-defective 2 promotes splicing of weak 5' splice sites.

Removal of introns during pre-mRNA splicing, which is central to gene expression, initiates by base pairing of U1 snRNA with a 5' splice site (5'SS). In mammals, many introns contain weak 5'SSs that are not efficiently recognized by the canonical U1 snRNP, suggesting alternative mechanisms exist. Here, we develop a cross-linking immunoprecipitation coupled to a high-throughput sequencing method, BCLIP-seq, to identify NRDE2 (nuclear RNAi-defective 2), and CCDC174 (coiled-coil domain-containing 174) as novel RNA-binding proteins in mouse ES cells that associate with U1 snRNA and 5'SSs. Both proteins bind directly to U1 snRNA independently of canonical U1 snRNP-specific proteins, and they are required for the selection and effective processing of weak 5'SSs. Our results reveal that mammalian cells use noncanonical splicing factors bound directly to U1 snRNA to effectively select suboptimal 5'SS sequences in hundreds of genes, promoting proper splice site choice, and accurate pre-mRNA splicing.

Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants.

BACKGROUND: Single-nucleotide variants (SNVs) within gene coding sequences can significantly impact pre-mRNA splicing, bearing profound implications for pathogenic mechanisms and precision medicine. In this study, we aim to harness the well-established full-length gene splicing assay (FLGSA) in conjunction with SpliceAI to prospectively interpret the splicing effects of all potential coding SNVs within the four-exon SPINK1 gene, a gene associated with chronic pancreatitis. RESULTS: Our study began with a retrospective analysis of 27 SPINK1 coding SNVs previously assessed using FLGSA, proceeded with a prospective analysis of 35 new FLGSA-tested SPINK1 coding SNVs, followed by data extrapolation, and ended with further validation. In total, we analyzed 67 SPINK1 coding SNVs, which account for 9.3% of the 720 possible coding SNVs. Among these 67 FLGSA-analyzed SNVs, 12 were found to impact splicing. Through detailed comparison of FLGSA results and SpliceAI predictions, we inferred that the remaining 653 untested coding SNVs in the SPINK1 gene are unlikely to significantly affect splicing. Of the 12 splice-altering events, nine produced both normally spliced and aberrantly spliced transcripts, while the remaining three only generated aberrantly spliced transcripts. These splice-impacting SNVs were found solely in exons 1 and 2, notably at the first and/or last coding nucleotides of these exons. Among the 12 splice-altering events, 11 were missense variants (2.17% of 506 potential missense variants), and one was synonymous (0.61% of 164 potential synonymous variants). Notably, adjusting the SpliceAI cut-off to 0.30 instead of the conventional 0.20 would improve specificity without reducing sensitivity. CONCLUSIONS: By integrating FLGSA with SpliceAI, we have determined that less than 2% (1.67%) of all possible coding SNVs in SPINK1 significantly influence splicing outcomes. Our findings emphasize the critical importance of conducting splicing analysis within the broader genomic sequence context of the study gene and highlight the inherent uncertainties associated with intermediate SpliceAI scores (0.20 to 0.80). This study contributes to the field by being the first to prospectively interpret all potential coding SNVs in a disease-associated gene with a high degree of accuracy, representing a meaningful attempt at shifting from retrospective to prospective variant analysis in the era of exome and genome sequencing.

The polyA tail facilitates splicing of last introns with weak 3' splice sites via PABPN1.

The polyA tail of mRNAs is important for many aspects of RNA metabolism. However, whether and how it regulates pre-mRNA splicing is still unknown. Here, we report that the polyA tail acts as a splicing enhancer for the last intron via the nuclear polyA binding protein PABPN1 in HeLa cells. PABPN1-depletion induces the retention of a group of introns with a weaker 3' splice site, and they show a strong 3'-end bias and mainly locate in nuclear speckles. The polyA tail is essential for PABPN1-enhanced last intron splicing and functions in a length-dependent manner. Tethering PABPN1 to nonpolyadenylated transcripts also promotes splicing, suggesting a direct role for PABPN1 in splicing regulation. Using TurboID-MS, we construct the PABPN1 interactome, including many spliceosomal and RNA-binding proteins. Specifically, PABPN1 can recruit RBM26&27 to promote splicing by interacting with the coiled-coil and RRM domain of RBM27. PABPN1-regulated terminal intron splicing is conserved in mice. Together, our study establishes a novel mode of post-transcriptional splicing regulation via the polyA tail and PABPN1.

Unique features of conventional and nonconventional introns in Euglena gracilis.

BACKGROUND: Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns. RESULTS: A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans. CONCLUSIONS: Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Spectrum of splicing variants in disease genes and the ability of RNA analysis to reduce uncertainty in clinical interpretation.

The complexities of gene expression pose challenges for the clinical interpretation of splicing variants. To better understand splicing variants and their contribution to hereditary disease, we evaluated their prevalence, clinical classifications, and associations with diseases, inheritance, and functional characteristics in a 689,321-person clinical cohort and two large public datasets. In the clinical cohort, splicing variants represented 13% of all variants classified as pathogenic (P), likely pathogenic (LP), or variants of uncertain significance (VUSs). Most splicing variants were outside essential splice sites and were classified as VUSs. Among all individuals tested, 5.4% had a splicing VUS. If RNA analysis were to contribute supporting evidence to variant interpretation, we estimated that splicing VUSs would be reclassified in 1.7% of individuals in our cohort. This would result in a clinically significant result (i.e., P/LP) in 0.1% of individuals overall because most reclassifications would change VUSs to likely benign. In ClinVar, splicing VUSs were 4.8% of reported variants and could benefit from RNA analysis. In the Genome Aggregation Database (gnomAD), splicing variants comprised 9.4% of variants in protein-coding genes; most were rare, precluding unambiguous classification as benign. Splicing variants were depleted in genes associated with dominant inheritance and haploinsufficiency, although some genes had rare variants at essential splice sites or had common splicing variants that were most likely compatible with normal gene function. Overall, we describe the contribution of splicing variants to hereditary disease, the potential utility of RNA analysis for reclassifying splicing VUSs, and how natural variation may confound clinical interpretation of splicing variants.

Characterization of NEB pathogenic variants in patients reveals novel nemaline myopathy disease mechanisms and omecamtiv mecarbil force effects.

Nebulin, a critical protein of the skeletal muscle thin filament, plays important roles in physiological processes such as regulating thin filament length (TFL), cross-bridge cycling, and myofibril alignment. Pathogenic variants in the nebulin gene (NEB) cause NEB-based nemaline myopathy (NEM2), a genetically heterogeneous disorder characterized by hypotonia and muscle weakness, currently lacking curative therapies. In this study, we examined a cohort of ten NEM2 patients, each with unique pathogenic variants, aiming to understand their impact on mRNA, protein, and functional levels. Results show that pathogenic truncation variants affect NEB mRNA stability and lead to nonsense-mediated decay of the mutated transcript. Moreover, a high incidence of cryptic splice site activation was found in patients with pathogenic splicing variants that are expected to disrupt the actin-binding sites of nebulin. Determination of protein levels revealed patients with either relatively normal or markedly reduced nebulin. We observed a positive relation between the reduction in nebulin and a reduction in TFL, or reduction in tension (both maximal and submaximal tension). Interestingly, our study revealed a pathogenic duplication variant in nebulin that resulted in a four-copy gain in the triplicate region of NEB and a much larger nebulin protein and longer TFL. Additionally, we investigated the effect of Omecamtiv mecarbil (OM), a small-molecule activator of cardiac myosin, on force production of type 1 muscle fibers of NEM2 patients. OM treatment substantially increased submaximal tension across all NEM2 patients ranging from 87 to 318%, with the largest effects in patients with the lowest level of nebulin. In summary, this study indicates that post-transcriptional or post-translational mechanisms regulate nebulin expression. Moreover, we propose that the pathomechanism of NEM2 involves not only shortened but also elongated thin filaments, along with the disruption of actin-binding sites resulting from pathogenic splicing variants. Significantly, our findings highlight the potential of OM treatment to improve skeletal muscle function in NEM2 patients, especially those with large reductions in nebulin levels.

Loss-of-function mutations in MYO15A and OTOF cause non-syndromic hearing loss in two Yemeni families.

BACKGROUND: Hearing loss is a rare hereditary deficit that is rather common among consanguineous populations. Autosomal recessive non-syndromic hearing loss is the predominant form of hearing loss worldwide. Although prevalent, hearing loss is extremely heterogeneous and poses a pitfall in terms of diagnosis and screening. Using next-generation sequencing has enabled a rapid increase in the identification rate of genes and variants in heterogeneous conditions, including hearing loss. We aimed to identify the causative variants in two consanguineous Yemeni families affected with hearing loss using targeted next-generation sequencing (clinical exome sequencing). The proband of each family was presented with sensorineural hearing loss as indicated by pure-tone audiometry results. RESULTS: We explored variants obtained from both families, and our analyses collectively revealed the presence and segregation of two novel loss-of-function variants: a frameshift variant, c.6347delA in MYO15A in Family I, and a splice site variant, c.5292-2A > C, in OTOF in Family II. Sanger sequencing and PCR-RFLP of DNA samples from 130 deaf and 50 control individuals confirmed that neither variant was present in our in-house database. In silico analyses predicted that each variant has a pathogenic effect on the corresponding protein. CONCLUSIONS: We describe two novel loss-of-function variants in MYO15A and OTOF that cause autosomal recessive non-syndromic hearing loss in Yemeni families. Our findings are consistent with previously reported pathogenic variants in the MYO15A and OTOF genes in Middle Eastern individuals and suggest their implication in hearing loss.

Two patients with KDM3B variants and new presentations of Diets-Jongmans syndrome.

KDM3B is located on chromosome 5q31 and encodes KDM3B, which is involved in histone demethylation and epigenetic regulation. Pathogenic KDM3B variants cause a dominantly inherited disorder presenting with intellectual disability (ID), short stature, and facial dysmorphism, named Diets-Jongmans syndrome. We describe two patients with KDM3B variants presenting with Diets-Jongmans syndrome. Genetic testing was performed because of the clinical data and a lack of a clear diagnosis in both patients. Candidate variants were verified by Sanger sequencing. After KDM3B variants were detected, in silico tools were used to predict the pathogenicity of the missense variants. A minigene assay was performed to evaluate the splicing effects of the c.5070 + 1G > A variant on KDM3B. Patient 1 mainly presented with repetitive upper respiratory tract infection and patient 2 presented with palpitation, shortness of breath, and pitting edema; both had ID. Whole exome sequencing identified variants of KDM3B. Patient 1 had the de novo KDM3B c.5070 + 1G > A variant, whereas patient 2 had the c.2828G > A (p.R943Q) variant. Transcriptional experiments of the splicing variant c.5070 + 1G > A revealed aberrant transcripts leading to truncated protein products. We found two pathogenic variants in KDM3B, one of which is novel. Both patients had additional clinical presentations, and patient 1 had transient neutropenia. KDM3B c.5070 + 1G > A is the first KDM3B splice-site variant and was identified as a germline variant. Neutropenia and cardiomyopathy are newly found presentations of Diets-Jongmans syndrome. Our report enriches our knowledge of the genotypic spectrum of the KDM3B variants and phenotypic diversity of Diets-Jongmans syndrome.

A rare non-canonical splice site in Trema orientalis SYMRK does not affect its dual symbiotic functioning in endomycorrhiza and rhizobium nodulation.

BACKGROUND: Nitrogen-fixing nodules occur in ten related taxonomic lineages interspersed with lineages of non-nodulating plant species. Nodules result from an endosymbiosis between plants and diazotrophic bacteria; rhizobia in the case of legumes and Parasponia and Frankia in the case of actinorhizal species. Nodulating plants share a conserved set of symbiosis genes, whereas related non-nodulating sister species show pseudogenization of several key nodulation-specific genes. Signalling and cellular mechanisms critical for nodulation have been co-opted from the more ancient plant-fungal arbuscular endomycorrhizal symbiosis. Studies in legumes and actinorhizal plants uncovered a key component in symbiotic signalling, the LRR-type SYMBIOSIS RECEPTOR KINASE (SYMRK). SYMRK is essential for nodulation and arbuscular endomycorrhizal symbiosis. To our surprise, however, despite its arbuscular endomycorrhizal symbiosis capacities, we observed a seemingly critical mutation in a donor splice site in the SYMRK gene of Trema orientalis, the non-nodulating sister species of Parasponia. This led us to investigate the symbiotic functioning of SYMRK in the Trema-Parasponia lineage and to address the question of to what extent a single nucleotide polymorphism in a donor splice site affects the symbiotic functioning of SYMRK. RESULTS: We show that SYMRK is essential for nodulation and endomycorrhization in Parasponia andersonii. Subsequently, it is revealed that the 5'-intron donor splice site of SYMRK intron 12 is variable and, in most dicotyledon species, doesn't contain the canonical dinucleotide 'GT' signature but the much less common motif 'GC'. Strikingly, in T. orientalis, this motif is converted into a rare non-canonical 5'-intron donor splice site 'GA'. This SYMRK allele, however, is fully functional and spreads in the T. orientalis population of Malaysian Borneo. A further investigation into the occurrence of the non-canonical GA-AG splice sites confirmed that these are extremely rare. CONCLUSION: SYMRK functioning is highly conserved in legumes, actinorhizal plants, and Parasponia. The gene possesses a non-common 5'-intron GC donor splice site in intron 12, which is converted into a GA in T. orientalis accessions of Malaysian Borneo. The discovery of this functional GA-AG splice site in SYMRK highlights a gap in our understanding of splice donor sites.

Inducible nuclear import by TetR aptamer-controlled 3' splice site selection.

Correct cellular localization is essential for the function of many eukaryotic proteins and hence cell physiology. Here, we present a synthetic genetic device that allows the control of nuclear and cytosolic localization based on controlled alternative splicing in human cells. The device is based on the fact that an alternative 3' splice site is located within a TetR aptamer that in turn is positioned between the branch point and the canonical splice site. The novel splice site is only recognized when the TetR repressor is bound. Addition of doxycycline prevents TetR aptamer binding and leads to recognition of the canonical 3' splice site. It is thus possible to produce two independent splice isoforms. Since the terminal loop of the aptamer may be replaced with any sequence of choice, one of the two isoforms may be extended by the respective sequence of choice depending on the presence of doxycycline. In a proof-of-concept study, we fused a nuclear localization sequence to a cytosolic target protein, thus directing the protein into the nucleus. However, the system is not limited to the control of nuclear localization. In principle, any target sequence can be integrated into the aptamer, allowing not only the production of a variety of different isoforms on demand, but also to study the function of mislocalized proteins. Moreover, it also provides a valuable tool for investigating the mechanism of alternative splicing in human cells.