Quantitative Activity Profile and Context Dependence of All Human 5' Splice Sites.

Pre-mRNA splicing is an essential step in the expression of most human genes. Mutations at the 5' splice site (5'ss) frequently cause defective splicing and disease due to interference with the initial recognition of the exon-intron boundary by U1 small nuclear ribonucleoprotein (snRNP), a component of the spliceosome. Here, we use a massively parallel splicing assay (MPSA) in human cells to quantify the activity of all 32,768 unique 5'ss sequences (NNN/GYNNNN) in three different gene contexts. Our results reveal that although splicing efficiency is mostly governed by the 5'ss sequence, there are substantial differences in this efficiency across gene contexts. Among other uses, these MPSA measurements facilitate the prediction of 5'ss sequence variants that are likely to cause aberrant splicing. This approach provides a framework to assess potential pathogenic variants in the human genome and streamline the development of splicing-corrective therapies.

Higher-order epistasis and phenotypic prediction.

Contemporary high-throughput mutagenesis experiments are providing an increasingly detailed view of the complex patterns of genetic interaction that occur between multiple mutations within a single protein or regulatory element. By simultaneously measuring the effects of thousands of combinations of mutations, these experiments have revealed that the genotype-phenotype relationship typically reflects not only genetic interactions between pairs of sites but also higher-order interactions among larger numbers of sites. However, modeling and understanding these higher-order interactions remains challenging. Here we present a method for reconstructing sequence-to-function mappings from partially observed data that can accommodate all orders of genetic interaction. The main idea is to make predictions for unobserved genotypes that match the type and extent of epistasis found in the observed data. This information on the type and extent of epistasis can be extracted by considering how phenotypic correlations change as a function of mutational distance, which is equivalent to estimating the fraction of phenotypic variance due to each order of genetic interaction (additive, pairwise, three-way, etc.). Using these estimated variance components, we then define an empirical Bayes prior that in expectation matches the observed pattern of epistasis and reconstruct the genotype-phenotype mapping by conducting Gaussian process regression under this prior. To demonstrate the power of this approach, we present an application to the antibody-binding domain GB1 and also provide a detailed exploration of a dataset consisting of high-throughput measurements for the splicing efficiency of human pre-mRNA [Formula: see text] splice sites, for which we also validate our model predictions via additional low-throughput experiments.

Alternative splicing regulation of telomerase: a new paradigm?

Alternative splicing affects approximately 95% of eukaryotic genes, greatly expanding the coding capacity of complex genomes. Although our understanding of alternative splicing has increased rapidly, current knowledge of splicing regulation has largely been derived from studies of highly expressed mRNAs. Telomerase is a key example of a protein that is alternatively spliced, but it is expressed at very low levels and although it is known that misregulation of telomerase splicing is a hallmark of nearly all cancers, the details of this process are unclear. Here we review work showing that hTERT expression is in part regulated by atypical alternative splicing, perhaps due to its exceptionally low expression level. We propose that these differential regulatory mechanisms may be widely applicable to other genes and may provide new opportunities for the development of cancer therapeutics.

Specificity, synergy, and mechanisms of splice-modifying drugs.

Drugs that target pre-mRNA splicing hold great therapeutic potential, but the quantitative understanding of how these drugs work is limited. Here we introduce mechanistically interpretable quantitative models for the sequence-specific and concentration-dependent behavior of splice-modifying drugs. Using massively parallel splicing assays, RNA-seq experiments, and precision dose-response curves, we obtain quantitative models for two small-molecule drugs, risdiplam and branaplam, developed for treating spinal muscular atrophy. The results quantitatively characterize the specificities of risdiplam and branaplam for 5' splice site sequences, suggest that branaplam recognizes 5' splice sites via two distinct interaction modes, and contradict the prevailing two-site hypothesis for risdiplam activity at SMN2 exon 7. The results also show that anomalous single-drug cooperativity, as well as multi-drug synergy, are widespread among small-molecule drugs and antisense-oligonucleotide drugs that promote exon inclusion. Our quantitative models thus clarify the mechanisms of existing treatments and provide a basis for the rational development of new therapies.

Regulation of human telomerase splicing by RNA:RNA pairing.

Telomerase adds telomeric repeats onto chromosome ends and is almost universally upregulated in human cancers. Here we demonstrate that RNA:RNA pairing regulates splicing of the catalytic subunit of human telomerase (TERT). Human alleles contain a variable number of 38 bp repeats within TERT intron 6 (>1 kb from exon-intron junctions). At least nine repeats are required for generating the major non-functional 'minus beta' isoform, which skips exons 7 and 8. RNA:RNA pairing between the repeats and the pre-mRNA might bring exons 6 and 9 closer, thereby promoting exon skipping. To demonstrate this, we show that mutations within the repeat that abolish exon skipping are corrected by compensatory mutations in the pre-mRNA. This study thus identifies RNA:RNA pairing by repetitive sequences as a novel form of alternative splicing regulation in a gene crucial for cancer survival and sheds new light on functional roles for short repetitive sequences embedded deep within introns throughout the genome.

Regulation of telomerase alternative splicing: a target for chemotherapy.

Telomerase is present in human cancer cells but absent in most somatic tissues. The messenger RNA of human telomerase (hTERT) is alternatively spliced into mostly nonfunctional products. We sought to understand splicing so that we could decrease functional splice isoforms to reduce telomerase activity in order to complement direct enzyme inhibition. Unexpectedly, minigenes containing hTERT exons 5-10 flanked by 150-300 bp intronic sequences did not produce alternative splicing. A 1.1 kb region of 38 bp repeats ~2 kb from the exon 6/intron junction restored the exclusion of exons 7 and 8. An element within intron 8, also >1 kb from intron/exon junctions, modulated this effect. Transducing an oligonucleotide complementary to this second element increased nonfunctional hTERT messenger RNA from endogenous telomerase. These results demonstrate the potential of manipulating hTERT splicing for both chemotherapy and regenerative medicine and provide specific sequences deep within introns that regulate alternative splicing in mammalian cells by mechanisms other than the introduction of cryptic splice sites.