Detection of splice isoforms and rare intermediates using multiplexed primer extension sequencing

Mechanisms and regulation of spliceosome-mediated pre-mRNA splicing in Saccharomyces cerevisiae

Pre-mRNA splicing, the removal of introns and ligation of flanking exons, is a crucial step in eukaryotic gene expression. The spliceosome, a macromolecular complex made up of five small nuclear RNAs (snRNAs) and dozens of proteins, assembles on introns via a complex pathway before catalyzing the two transesterification reactions necessary for splicing. All of these steps have the potential to be highly regulated to ensure correct mRNA isoform production for proper cellular function. While Saccharomyces cerevisiae (yeast) has a limited set of intron-containing genes, many of these genes are highly expressed, resulting in a large number of transcripts in a cell being spliced. As a result, splicing regulation is of critical importance for yeast. Just as in humans, yeast splicing can be influenced by protein components of the splicing machinery, structures and properties of the pre-mRNA itself, or by the action of trans-acting factors. It is likely that further analysis of the mechanisms and pathways of splicing regulation in yeast can reveal general principles applicable to other eukaryotes. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Timing is everything: advances in quantifying splicing kinetics

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

Specific recognition and sensitive quantification of mRNA splice variants via one-pot ligation-dependent loop-mediated isothermal amplification

Specific recognition and sensitive quantification of mRNA alternative splice variants have been a necessity for exploring the regulatory mechanism of RNA splicing and revealing the association between pre-mRNA splicing and transcriptome function, as well as disease diagnosis. However, their wide abundance range and high sequence homology pose enormous challenges for high sensitivity and selectivity quantification of splice variants. Herein, taking advantage of the excellent specificity of ligation and the powerful nucleic acid replication feature of loop-mediated isothermal amplification (LAMP), we developed a one-pot method (termed one-pot ligation-LAMP) for specific recognition and sensitive quantification of mRNA splicing variants based on two splicing junction-specific stem-loop DNA probe ligation and the subsequently initiating LAMP. The one-pot ligation-LAMP can specifically detect as low as 100 aM mRNA splice variants without any nonspecific signals and quantify them with a wide dynamics range spanning at least six orders of magnitude. We have demonstrated that the one-pot ligation-LAMP is a versatile and practical strategy for accurately quantifying different splicing variants in complex biological samples with high sensitivity all in one tube within 90 min, thereby providing an attractive tool for mRNA splice variant-related studies.

Complexities of mammalian transcriptome revealed by targeted RNA enrichment techniques

Studies using highly sensitive targeted RNA enrichment methods have shown that a large portion of the human transcriptome remains to be discovered and that most of the genome is transcribed in a complex, interleaved fashion characterized by a complex web of transcripts emanating from protein coding and noncoding loci. These results resonate with those from single-cell transcriptome profiling endeavors that reveal the existence of multiple novel, cell type-specific transcripts and clearly demonstrate that our understanding of the complexities of the human transcriptome is far from being complete. Here, we review the current status of the targeted RNA enrichment techniques, their application to the discovery of novel cell type-specific transcripts, and their impact on our understanding of the human genome and transcriptome.

Transcript-specific determinants of pre-mRNA splicing revealed through in vivo kinetic analyses of the 1st and 2nd chemical steps

We generate high-precision measurements of the in vivo rates of both chemical steps of pre-mRNA splicing across the genome-wide complement of substrates in yeast by coupling metabolic labeling, multiplexed primer-extension sequencing, and kinetic modeling. We demonstrate that the rates of intron removal vary widely, splice-site sequences are primary determinants of 1st step but have little apparent impact on 2nd step rates, and the 2nd step is generally faster than the 1st step. Ribosomal protein genes (RPGs) are spliced faster than non-RPGs at each step, and RPGs share evolutionarily conserved properties that may contribute to their faster splicing. A genetic variant defective in the 1st step of the pathway reveals a genome-wide defect in the 1st step but an unexpected, transcript-specific change in the 2nd step. Our work demonstrates that extended co-transcriptional association is an important determinant of splicing rate, a conclusion at odds with recent claims of ultra-fast splicing.

Advanced preparation of fragment libraries enabled by oligonucleotide-modified 2',3'-dideoxynucleotides

The ever-growing demand for inexpensive, rapid, and accurate exploration of genomes calls for refinement of existing sequencing techniques. The development of next-generation sequencing (NGS) was a revolutionary milestone in genome analysis. While modified nucleotides already were inherent tools in sequencing and imaging, further modification of nucleotides enabled the expansion into even more diverse applications. Herein we describe the design and synthesis of oligonucleotide-tethered 2',3'-dideoxynucleotide (dd^ONNTP) terminators bearing universal priming sites attached to the nucleobase, as well as their enzymatic incorporation and performance in read-through assays. In the context of NGS library preparation, the incorporation of dd^ONNTP fulfills two requirements at once: the fragmentation step is integrated into the workflow and the obtained fragments are readily labeled by platform-specific adapters. DNA polymerases can incorporate dd^ONNTP nucleotides, as shown by primer extension assays. More importantly, reading through the unnatural linkage during DNA synthesis was demonstrated, with 25-30% efficiency in single-cycle extension.

Baiting out a full length sequence from unmapped RNA-seq data

BACKGROUND

As a powerful tool, RNA-Seq has been widely used in various studies. Usually, unmapped RNA-seq reads have been considered as useless and been trashed or ignored.

RESULTS

We develop a strategy to mining the full length sequence by unmapped reads combining with specific reverse transcription primers design and high throughput sequencing. In this study, we salvage 36 unmapped reads from standard RNA-Seq data and randomly select one 149 bp read as a model. Specific reverse transcription primers are designed to amplify its both ends, followed by next generation sequencing. Then we design a statistical model based on power law distribution to estimate its integrality and significance. Further, we validate it by Sanger sequencing. The result shows that the full length is 1556 bp, with insertion mutations in microsatellite structure.

CONCLUSION

We believe this method would be a useful strategy to extract the sequences information from the unmapped RNA-seq data. Further, it is an alternative way to get the full length sequence of unknown cDNA.

Dynamic imaging of nascent RNA reveals general principles of transcription dynamics and stochastic splice site selection

The activities of RNA polymerase and the spliceosome are responsible for the heterogeneity in the abundance and isoform composition of mRNA in human cells. However, the dynamics of these megadalton enzymatic complexes working in concert on endogenous genes have not been described. Here, we establish a quasi-genome-scale platform for observing synthesis and processing kinetics of single nascent RNA molecules in real time. We find that all observed genes show transcriptional bursting. We also observe large kinetic variation in intron removal for single introns in single cells, which is inconsistent with deterministic splice site selection. Transcriptome-wide footprinting of the U2AF complex, nascent RNA profiling, long-read sequencing, and lariat sequencing further reveal widespread stochastic recursive splicing within introns. We propose and validate a unified theoretical model to explain the general features of transcription and pervasive stochastic splice site selection.

Targeted RNA sequencing enhances gene expression profiling of ultra-low input samples

RNA-seq is the standard method for profiling gene expression in many biological systems. Due to the wide dynamic range and complex nature of the transcriptome, RNA-seq provides an incomplete characterization, especially of lowly expressed genes and transcripts. Targeted RNA sequencing (RNA CaptureSeq) focuses sequencing on genes of interest, providing exquisite sensitivity for transcript detection and quantification. However, uses of CaptureSeq have focused on bulk samples and its performance on very small populations of cells is unknown. Here we show CaptureSeq greatly enhances transcriptomic profiling of target genes in ultra-low-input samples and provides equivalent performance to that on bulk samples. We validate the performance of CaptureSeq using multiple probe sets on samples of iPSC-derived cortical neurons. We demonstrate up to 275-fold enrichment for target genes, the detection of 10% additional genes and a greater than 5-fold increase in identified gene isoforms. Analysis of spike-in controls demonstrated CaptureSeq improved both detection sensitivity and expression quantification. Comparison to the CORTECON database of cerebral cortex development revealed CaptureSeq enhanced the identification of sample differentiation stage. CaptureSeq provides sensitive, reliable and quantitative expression measurements on hundreds-to-thousands of target genes from ultra-low-input samples and has the potential to greatly enhance transcriptomic profiling when samples are limiting.

Native molecule sequencing by nano-ID reveals synthesis and stability of RNA isoforms

Eukaryotic genes often generate a variety of RNA isoforms that can lead to functionally distinct protein variants. The synthesis and stability of RNA isoforms is poorly characterized because current methods to quantify RNA metabolism use short-read sequencing and cannot detect RNA isoforms. Here we present nanopore sequencing–based isoform dynamics (nano-ID), a method that detects newly synthesized RNA isoforms and monitors isoform metabolism. Nano-ID combines metabolic RNA labeling, long-read nanopore sequencing of native RNA molecules, and machine learning. Nano-ID derives RNA stability estimates and evaluates stability determining factors such as RNA sequence, poly(A)-tail length, secondary structure, translation efficiency, and RNA-binding proteins. Application of nano-ID to the heat shock response in human cells reveals that many RNA isoforms change their stability. Nano-ID also shows that the metabolism of individual RNA isoforms differs strongly from that estimated for the combined RNA signal at a specific gene locus. Nano-ID enables studies of RNA metabolism at the level of single RNA molecules and isoforms in different cell states and conditions.

Accessing a New Dimension in TP53 Biology: Multiplex Long Amplicon Digital PCR to Specifically Detect and Quantitate Individual TP53 Transcripts

TP53, the most commonly-mutated gene in cancer, undergoes complex alternative splicing. Different TP53 transcripts play different biological roles, both in normal function and in the progression of diseases such as cancer. The study of TP53’s alternative RNA splice forms and their use as clinical biomarkers has been hampered by limited specificity and quantitative accuracy of current methods. TP53 RNA splice variants differ at both 5’ and 3’ ends, but because they have a common central region of 618 bp, the individual TP53 transcripts are impossible to specifically detect and precisely quantitate using standard PCR-based methods or short-read RNA sequencing. Therefore, we devised multiplex probe-based long amplicon droplet digital PCR (ddPCR) assays, which for the first time allow precise end-to-end quantitation of the seven major TP53 transcripts, with amplicons ranging from 0.85 to 1.85 kb. Multiple modifications to standard ddPCR assay procedures were required to enable specific co-amplification of these long transcripts and to overcome issues with secondary structure. Using these assays, we show that several TP53 transcripts are co-expressed in breast cancers, and illustrate the potential for this method to identify novel TP53 transcripts in tumour cells. This capability will facilitate a new level of biological and clinical understanding of the alternatively-spliced TP53 isoforms.

Multiplexed primer extension sequencing: A targeted RNA-seq method that enables high-precision quantitation of mRNA splicing isoforms and rare pre-mRNA splicing intermediates

The study of pre-mRNA splicing has been greatly aided by the advent of RNA sequencing (RNA-seq), which enables the genome-wide detection of discrete splice isoforms. Quantification of these splice isoforms requires analysis of splicing informative sequencing reads, those that unambiguously map to a single splice isoform, including exon-intron spanning alignments corresponding to retained introns, as well as exon-exon junction spanning alignments corresponding to either canonically- or alternatively-spliced isoforms. Because most RNA-seq experiments are designed to produce sequencing alignments that uniformly cover the entirety of transcripts, only a comparatively small number of splicing informative alignments are generated for any given splice site, leading to a decreased ability to detect and/or robustly quantify many splice isoforms. To address this problem, we have recently described a method termed Multiplexed Primer Extension sequencing, or MPE-seq, which uses pools of reverse transcription primers to target sequencing to user selected loci. By targeting reverse transcription to pre-mRNA splice junctions, this approach enables a dramatic enrichment in the fraction of splicing informative alignments generated per splicing event, yielding an increase in both the precision with which splicing efficiency can be measured, and in the detection of splice isoforms including rare splicing intermediates. Here we provide a brief review of the shortcomings associated with RNA-seq that drove our development of MPE-seq, as well as a detailed protocol for implementation of MPE-seq.