Intron lariat spliceosomes convert lariats to true circles: implications for intron transposition

Rare, full-length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envisioned and tested a hypothesis for their formation using Saccharomyces cerevisiae, documenting full-length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full-length and processed circles. Postsplicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Intron-lariat spliceosomes convert lariats to true circles: implications for intron transposition

Rare, full length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envision and test a hypothesis for their formation using Saccharomyces cerevisiae, documenting full length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron-lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full length and processed circles. Post-splicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Broad variation in response of individual introns to splicing inhibitors in a humanized yeast strain

Intron branchpoint (BP) recognition by the U2 snRNP is a critical step of splicing, vulnerable to recurrent cancer mutations and bacterial natural product inhibitors. The BP binds a conserved pocket in the SF3B1 (human) or Hsh155 (yeast) U2 snRNP protein. Amino acids that line this pocket affect the binding of splicing inhibitors like Pladienolide-B (Plad-B), such that organisms differ in their sensitivity. To study the mechanism of splicing inhibitor action in a simplified system, we modified the naturally Plad-B resistant yeast Saccharomyces cerevisiae by changing 14 amino acids in the Hsh155 BP pocket to those from human. This humanized yeast grows normally, and splicing is largely unaffected by the mutation. Splicing is inhibited within minutes after the addition of Plad-B, and different introns appear inhibited to different extents. Intron-specific inhibition differences are also observed during cotranscriptional splicing in Plad-B using single-molecule intron tracking to minimize gene-specific transcription and decay rates that cloud estimates of inhibition by standard RNA-seq. Comparison of Plad-B intron sensitivities to those of the structurally distinct inhibitor Thailanstatin-A reveals intron-specific differences in sensitivity to different compounds. This work exposes a complex relationship between the binding of different members of this class of inhibitors to the spliceosome and intron-specific rates of BP recognition and catalysis. Introns with variant BP sequences seem particularly sensitive, echoing observations from mammalian cells, where monitoring individual introns is complicated by multi-intron gene architecture and alternative splicing. The compact yeast system may hasten the characterization of splicing inhibitors, accelerating improvements in selectivity and therapeutic efficacy.

Broad variation in response of individual introns to splicing inhibitors in a humanized yeast strain

Intron branch point (BP) recognition by the U2 snRNP is a critical step of splicing, vulnerable to recurrent cancer mutations and bacterial natural product inhibitors. The BP binds a conserved pocket in the SF3B1 (human) or Hsh155 (yeast) U2 snRNP protein. Amino acids that line this pocket affect binding of splicing inhibitors like Pladienolide-B (Plad-B), such that organisms differ in their sensitivity. To study the mechanism of splicing inhibitor action in a simplified system, we modified the naturally Plad-B resistant yeast Saccharomyces cerevisiae by changing 14 amino acids in the Hsh155 BP pocket to those from human. This humanized yeast grows normally, and splicing is largely unaffected by the mutation. Splicing is inhibited within minutes after addition of Plad-B, and different introns appear inhibited to different extents. Intron-specific inhibition differences are also observed during co-transcriptional splicing in Plad-B using single-molecule intron tracking (SMIT) to minimize gene-specific transcription and decay rates that cloud estimates of inhibition by standard RNA-seq. Comparison of Plad-B intron sensitivities to those of the structurally distinct inhibitor Thailanstatin-A reveals intron-specific differences in sensitivity to different compounds. This work exposes a complex relationship between binding of different members of this class of inhibitors to the spliceosome and intron-specific rates of BP recognition and catalysis. Introns with variant BP sequences seem particularly sensitive, echoing observations from mammalian cells, where monitoring individual introns is complicated by multi-intron gene architecture and alternative splicing. The compact yeast system may hasten characterization of splicing inhibitors, accelerating improvements in selectivity and therapeutic efficacy.

Transposable elements drive intron gain in diverse eukaryotes

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Concerted modification of nucleotides at functional centers of the ribosome revealed by single-molecule RNA modification profiling

Nucleotides in RNA and DNA are chemically modified by numerous enzymes that alter their function. Eukaryotic ribosomal RNA (rRNA) is modified at more than 100 locations, particularly at highly conserved and functionally important nucleotides. During ribosome biogenesis, modifications are added at various stages of assembly. The existence of differently modified classes of ribosomes in normal cells is unknown because no method exists to simultaneously evaluate the modification status at all sites within a single rRNA molecule. Using a combination of yeast genetics and nanopore direct RNA sequencing, we developed a reliable method to track the modification status of single rRNA molecules at 37 sites in 18 S rRNA and 73 sites in 25 S rRNA. We use our method to characterize patterns of modification heterogeneity and identify concerted modification of nucleotides found near functional centers of the ribosome. Distinct, undermodified subpopulations of rRNAs accumulate upon loss of Dbp3 or Prp43 RNA helicases, suggesting overlapping roles in ribosome biogenesis. Modification profiles are surprisingly resistant to change in response to many genetic and acute environmental conditions that affect translation, ribosome biogenesis, and pre-mRNA splicing. The ability to capture single-molecule RNA modification profiles provides new insights into the roles of nucleotide modifications in RNA function.

Synthesis of modified nucleotide polymers by the poly(U) polymerase Cid1: application to direct RNA sequencing on nanopores

Understanding transcriptomes requires documenting the structures, modifications, and abundances of RNAs as well as their proximity to other molecules. The methods that make this possible depend critically on enzymes (including mutant derivatives) that act on nucleic acids for capturing and sequencing RNA. We tested two 3' nucleotidyl transferases, Saccharomyces cerevisiae poly(A) polymerase and Schizosaccharomyces pombe Cid1, for the ability to add base and sugar modified rNTPs to free RNA 3' ends, eventually focusing on Cid1. Although unable to polymerize ΨTP or 1meΨTP, Cid1 can use 5meUTP and 4thioUTP. Surprisingly, Cid1 can use inosine triphosphate to add poly(I) to the 3' ends of a wide variety of RNA molecules. Most poly(A) mRNAs efficiently acquire a uniform tract of about 50 inosine residues from Cid1, whereas non-poly(A) RNAs acquire longer, more heterogeneous tails. Here we test these activities for use in direct RNA sequencing on nanopores, and find that Cid1-mediated poly(I)-tailing permits detection and quantification of both mRNAs and non-poly(A) RNAs simultaneously, as well as enabling the analysis of nascent RNAs associated with RNA polymerase II. Poly(I) produces a different current trace than poly(A), enabling recognition of native RNA 3' end sequence lost by in vitro poly(A) addition. Addition of poly(I) by Cid1 offers a broadly useful alternative to poly(A) capture for direct RNA sequencing on nanopores.

The UCSC SARS-CoV-2 Genome Browser

Background:

Researchers are generating molecular data pertaining to the SARS-CoV-2 RNA genome and its proteins at an unprecedented rate during the COVID-19 pandemic. As a result, there is a critical need for rapid and continuously updated access to the latest molecular data in a format in which all data can be quickly cross-referenced and compared. We adapted our genome browser visualization tool to the viral genome for this purpose. Molecular data, curated from published studies or from database submissions, are mapped to the viral genome and grouped together into “annotation tracks” where they can be visualized along the linear map of the viral genome sequence and programmatically downloaded in standard format for analysis.

Results:

The UCSC Genome Browser for SARS-CoV-2 (https://genome.ucsc.edu/covid19.html ) provides continuously updated access to the mutations in the many thousands of SARS-CoV-2 genomes deposited in GISAID and the international nucleotide sequencing databases, displayed alongside phylogenetic trees. These data are augmented with alignments of bat, pangolin, and other animal and human coronavirus genomes, including per-base evolutionary rate analysis. All available annotations are cross-referenced on the virus genome, including those from major databases (PDB, RFAM, IEDB, UniProt) as well as up-to-date individual results from preprints. Annotated data include predicted and validated immune epitopes, promising antibodies, RT-PCR and sequencing primers, CRISPR guides (from research, diagnostics, vaccines, and therapies), and points of interaction between human and viral genes. As a community resource, any user can add manual annotations which are quality checked and shared publicly on the browser the next day.

Conclusions:

We invite all investigators to contribute additional data and annotations to this resource to accelerate research and development activities globally. Contact us at genome-www@soe.ucsc.edu with data suggestions or requests for support for adding data. Rapid sharing of data will accelerate SARS-CoV-2 research, especially when researchers take time to integrate their data with those from other labs on a widely-used community browser platform with standardized machine-readable data formats, such as the SARS-CoV-2 Genome Browser.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these “hungry spliceosome” conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of “intronization”, whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

The protein coding information in eukaryotic genes is broken by intervening sequences called introns that are removed from RNA during transcription by a large protein-RNA complex called the spliceosome. Where introns come from and how the spliceosome contributes to genome evolution are open questions. In this study, we find more than 150 new places in the yeast genome that are recognized by the spliceosome and spliced out as introns. Since they appear to have arisen very recently in evolution by sequence drift and do not appear to contribute to gene expression or its regulation, we call these protointrons. Protointrons are found in both protein-coding and non-coding RNAs and are not efficiently removed by the splicing machinery. Although most protointrons are not conserved and will likely disappear as evolution proceeds, a few are spliced more efficiently, and are located where they might begin to play functional roles in gene expression, as predicted by the proposed process of intronization. The challenge now is to understand how spontaneously appearing splicing events like protointrons might contribute to the creation of new genes, new genetic controls, and new protein isoforms as genomes evolve.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these "hungry spliceosome" conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of "intronization", whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

Cus2 enforces the first ATP-dependent step of splicing by binding to yeast SF3b1 through a UHM-ULM interaction

Stable recognition of the intron branchpoint (BP) by the U2 snRNP to form the pre-spliceosome is the first ATP-dependent step of splicing. Genetic and biochemical data from yeast indicate that Cus2 aids U2 snRNA folding into the stem IIa conformation prior to pre-spliceosome formation. Cus2 must then be removed by an ATP-dependent function of Prp5 before assembly can progress. However, the location from which Cus2 is displaced and the nature of its binding to the U2 snRNP are unknown. Here, we show that Cus2 contains a conserved UHM (U2AF homology motif) that binds Hsh155, the yeast homolog of human SF3b1, through a conserved ULM (U2AF ligand motif). Mutations in either motif block binding and allow pre-spliceosome formation without ATP. A 2.0 Å resolution structure of the Hsh155 ULM in complex with the UHM of Tat-SF1, the human homolog of Cus2, and complementary binding assays show that the interaction is highly similar between yeast and humans. Furthermore, we show that Tat-SF1 can replace Cus2 function by enforcing ATP dependence of pre-spliceosome formation in yeast extracts. Cus2 is removed before pre-spliceosome formation, and both Cus2 and its Hsh155 ULM binding site are absent from available cryo-EM structure models. However, our data are consistent with the apparent location of the disordered Hsh155 ULM between the U2 stem-loop IIa and the HEAT repeats of Hsh155 that interact with Prp5. We propose a model in which Prp5 uses ATP to remove Cus2 from Hsh155 such that extended base-pairing between U2 snRNA and the intron BP can occur.

Autogenous cross-regulation of Quaking mRNA processing and translation balances Quaking functions in splicing and translation

Quaking protein isoforms arise from a single Quaking gene and bind the same RNA motif to regulate splicing, translation, decay, and localization of a large set of RNAs. However, the mechanisms by which Quaking expression is controlled to ensure that appropriate amounts of each isoform are available for such disparate gene expression processes are unknown. Here we explore how levels of two isoforms, nuclear Quaking-5 (Qk5) and cytoplasmic Qk6, are regulated in mouse myoblasts. We found that Qk5 and Qk6 proteins have distinct functions in splicing and translation, respectively, enforced through differential subcellular localization. We show that Qk5 and Qk6 regulate distinct target mRNAs in the cell and act in distinct ways on their own and each other's transcripts to create a network of autoregulatory and cross-regulatory feedback controls. Morpholino-mediated inhibition of Qk translation confirms that Qk5 controls Qk RNA levels by promoting accumulation and alternative splicing of Qk RNA, whereas Qk6 promotes its own translation while repressing Qk5. This Qk isoform cross-regulatory network responds to additional cell type and developmental controls to generate a spectrum of Qk5/Qk6 ratios, where they likely contribute to the wide range of functions of Quaking in development and cancer.

Protein-RNA Networks Regulated by Normal and ALS-Associated Mutant HNRNPA2B1 in the Nervous System

HnRNPA2B1 encodes an RNA binding protein associated with neurodegeneration. However, its function in the nervous system is unclear. Transcriptome-wide crosslinking and immunoprecipitation in mouse spinal cord discover UAGG motifs enriched within ∼2,500 hnRNP A2/B1 binding sites and an unexpected role for hnRNP A2/B1 in alternative polyadenylation. HnRNP A2/B1 loss results in alternative splicing (AS), including skipping of an exon in amyotrophic lateral sclerosis (ALS)-associated D-amino acid oxidase (DAO) that reduces D-serine metabolism. ALS-associated hnRNP A2/B1 D290V mutant patient fibroblasts and motor neurons differentiated from induced pluripotent stem cells (iPSC-MNs) demonstrate abnormal splicing changes, likely due to increased nuclear-insoluble hnRNP A2/B1. Mutant iPSC-MNs display decreased survival in long-term culture and exhibit hnRNP A2/B1 localization to cytoplasmic granules as well as exacerbated changes in gene expression and splicing upon cellular stress. Our findings provide a cellular resource and reveal RNA networks relevant to neurodegeneration, regulated by normal and mutant hnRNP A2/B1. VIDEO ABSTRACT.

Distinct and shared functions of ALS-associated proteins TDP-43, FUS and TAF15 revealed by multisystem analyses

The RNA-binding protein (RBP) TAF15 is implicated in amyotrophic lateral sclerosis (ALS). To compare TAF15 function to that of two ALS-associated RBPs, FUS and TDP-43, we integrate CLIP-seq and RNA Bind-N-Seq technologies, and show that TAF15 binds to ∼4,900 RNAs enriched for GGUA motifs in adult mouse brains. TAF15 and FUS exhibit similar binding patterns in introns, are enriched in 3' untranslated regions and alter genes distinct from TDP-43. However, unlike FUS and TDP-43, TAF15 has a minimal role in alternative splicing. In human neural progenitors, TAF15 and FUS affect turnover of their RNA targets. In human stem cell-derived motor neurons, the RNA profile associated with concomitant loss of both TAF15 and FUS resembles that observed in the presence of the ALS-associated mutation FUS R521G, but contrasts with late-stage sporadic ALS patients. Taken together, our findings reveal convergent and divergent roles for FUS, TAF15 and TDP-43 in RNA metabolism.

Abstract 47: Quaking Post-Transcriptionally Promotes Differentiation of Monocytes Into Pro-Atherogenic Macrophages by Controling Pre-mRNA Splicing and Gene Expression

Aim: Atherosclerosis is accelerated by excessive monocyte recruitment, influx and differentiation into pro-inflammatory macrophages. While the genome-wide mRNA expression profiles of human monocytes and pro-inflammatory macrophages are well-established, their transcriptomes are ultimately defined by factors, such as RNA-binding proteins, that modulate pre-mRNA splicing patterns and mRNA transcript abundance. This prompted us to investigate the role of the RNA-binding protein Quaking in regulating global changes in pre-mRNA splicing and mature mRNA expression as human monocytes acquire the pro-inflammatory macrophage identity.

Methods: We employed RNA-sequencing and splicing-sensitive microarrays to determine genome-wide changes in pre-mRNA splicing and mRNA expression upon conversion of human monocytes into pro-inflammatory macrophages, including those derived from a unique Quaking haploinsufficient patient.

Results: Using laser-capture micro-dissection and immunohistochemistry, we discovered that expression levels of Quaking mRNA and protein are low in monocytes of early human atherosclerotic lesions, but abundant in macrophages of advanced plaques. Depletion of Quaking protein using both siRNA and GapmeR technology significantly impaired monocyte adhesion and migration; delayed differentiation into pro-inflammatory macrophages while maintaining the capacity to adopt the anti-inflammatory phenotype; and diminished foam cell formation in vitro and in vivo. RNA-sequencing and microarray analysis of human monocyte and macrophage transcriptomes revealed striking changes in Quaking-dependent pre-mRNA splicing and mRNA transcript levels, with gene ontology analyses identifying an enrichment in transcripts involved in cellular migration and lipid metabolism. Furthermore, these studies uncovered common as well as novel alternatively spliced transcripts with unknown biological functions in monocytes and macrophages.

Conclusions: Our studies illustrate a central role for Quaking in post-transcriptionally guiding pro-inflammatory macrophage identity and function.

Intron Invasions Trace Algal Speciation and Reveal Nearly Identical Arctic and Antarctic Micromonas Populations

Spliceosomal introns are a hallmark of eukaryotic genes that are hypothesized to play important roles in genome evolution but have poorly understood origins. Although most introns lack sequence homology to each other, new families of spliceosomal introns that are repeated hundreds of times in individual genomes have recently been discovered in a few organisms. The prevalence and conservation of these introner elements (IEs) or introner-like elements in other taxa, as well as their evolutionary relationships to regular spliceosomal introns, are still unknown. Here, we systematically investigate introns in the widespread marine green alga Micromonas and report new families of IEs, numerous intron presence–absence polymorphisms, and potential intron insertion hot-spots. The new families enabled identification of conserved IE secondary structure features and establishment of a novel general model for repetitive intron proliferation across genomes. Despite shared secondary structure, the IE families from each Micromonas lineage bear no obvious sequence similarity to those in the other lineages, suggesting that their appearance is intimately linked with the process of speciation. Two of the new IE families come from an Arctic culture (Micromonas Clade E2) isolated from a polar region where abundance of this alga is increasing due to climate induced changes. The same two families were detected in metagenomic data from Antarctica—a system where Micromonas has never before been reported. Strikingly high identity between the Arctic isolate and Antarctic coding sequences that flank the IEs suggests connectivity between populations in the two polar systems that we postulate occurs through deep-sea currents. Recovery of Clade E2 sequences in North Atlantic Deep Waters beneath the Gulf Stream supports this hypothesis. Our research illuminates the dynamic relationships between an unusual class of repetitive introns, genome evolution, speciation, and global distribution of this sentinel marine alga.

Abstract 217: Quaking Post-Transcriptionally Guides Monocyte Adhesion and Differentiation into the Pro-Inflammatory Macrophage

A hallmark of inflammatory diseases is the excessive recruitment and influx of monocytes to sites of tissue damage and their ensuing differentiation into macrophages. Numerous stimuli are known to induce new transcription necessary for macrophage identity, but post-transcriptional control of human macrophage differentiation is less well understood. Here, we detail our discovery that levels of the RNA-binding protein Quaking (QKI) are low in monocytes of early atherosclerotic lesions, but abundant in macrophages of advanced plaques. Specific depletion of QKI protein impaired monocyte adhesion, migration and differentiation into macrophages, and lesion formation. RNA-seq and microarray analysis of human monocyte and macrophage transcriptomes, including those of a unique QKI haploinsufficient patient, reveal developmental changes in RNA levels and alternative splicing of RNA transcripts enriched in QKI-bound sequence elements. The importance of these transcripts and requirement for QKI during differentiation illustrates a central role for QKI in post-transcriptionally guiding macrophage identity and function. These studies implicate QKI as a novel target for therapeutic intervention in inflammatory diseases.

Methods for processing high-throughput RNA sequencing data

High-throughput sequencing (HTS) methods for analyzing RNA populations (RNA-Seq) are gaining rapid application to many experimental situations. The steps in an RNA-Seq experiment require thought and planning, especially because the expense in time and materials is currently higher and the protocols are far less routine than those used for other high-throughput methods, such as microarrays. As always, good experimental design will make analysis and interpretation easier. Having a clear biological question, an idea about the best way to do the experiment, and an understanding of the number of replicates needed will make the entire process more satisfying. Whether the goal is capturing transcriptome complexity from a tissue or identifying small fragments of RNA cross-linked to a protein of interest, conversion of the RNA to cDNA followed by direct sequencing using the latest methods is a developing practice, with new technical modifications and applications appearing every day. Even more rapid are the development and improvement of methods for analysis of the very large amounts of data that arrive at the end of an RNA-Seq experiment, making considerations regarding reproducibility, validation, visualization, and interpretation increasingly important. This introduction is designed to review and emphasize a pathway of analysis from experimental design through data presentation that is likely to be successful, with the recognition that better methods are right around the corner.

Basic quantitative polymerase chain reaction using real-time fluorescence measurements

This protocol uses quantitative polymerase chain reaction (qPCR) to measure the number of DNA molecules containing a specific contiguous sequence in a sample of interest (e.g., genomic DNA or cDNA generated by reverse transcription). The sample is subjected to fluorescence-based PCR amplification and, theoretically, during each cycle, two new duplex DNA molecules are produced for each duplex DNA molecule present in the sample. The progress of the reaction during PCR is evaluated by measuring the fluorescence of dsDNA-dye complexes in real time. In the early cycles, DNA duplication is not detected because inadequate amounts of DNA are made. At a certain threshold cycle, DNA-dye complexes double each cycle for 8-10 cycles, until the DNA concentration becomes so high and the primer concentration so low that the reassociation of the product strands blocks efficient synthesis of new DNA and the reaction plateaus. There are two types of measurements: (1) the relative change of the target sequence compared to a reference sequence and (2) the determination of molecule number in the starting sample. The first requires a reference sequence, and the second requires a sample of the target sequence with known numbers of the molecules of sequence to generate a standard curve. By identifying the threshold cycle at which a sample first begins to accumulate DNA-dye complexes exponentially, an estimation of the numbers of starting molecules in the sample can be extrapolated.

RBPmap: a web server for mapping binding sites of RNA-binding proteins

Regulation of gene expression is executed in many cases by RNA-binding proteins (RBPs) that bind to mRNAs as well as to non-coding RNAs. RBPs recognize their RNA target via specific binding sites on the RNA. Predicting the binding sites of RBPs is known to be a major challenge. We present a new webserver, RBPmap, freely accessible through the website http://rbpmap.technion.ac.il/ for accurate prediction and mapping of RBP binding sites. RBPmap has been developed specifically for mapping RBPs in human, mouse and Drosophila melanogaster genomes, though it supports other organisms too. RBPmap enables the users to select motifs from a large database of experimentally defined motifs. In addition, users can provide any motif of interest, given as either a consensus or a PSSM. The algorithm for mapping the motifs is based on a Weighted-Rank approach, which considers the clustering propensity of the binding sites and the overall tendency of regulatory regions to be conserved. In addition, RBPmap incorporates a position-specific background model, designed uniquely for different genomic regions, such as splice sites, 5’ and 3’ UTRs, non-coding RNA and intergenic regions. RBPmap was tested on high-throughput RNA-binding experiments and was proved to be highly accurate.

Tips on hybridizing, washing, and scanning affymetrix microarrays

Starting in the late 1990s, Affymetrix, Inc. produced a commercial system for hybridizing, washing, and scanning microarrays that was designed to be easy to operate and reproducible. The system used arrays packaged in a plastic cassette or chamber in which the prefabricated array was mounted and could be filled with fluid through resealable membrane ports either by hand or by an automated "fluidics station" specially designed to handle the arrays. A special rotating hybridization oven and a specially designed scanner were also required. Primarily because of automation and standardization the Affymetrix system was and still remains popular. Here, we provide a skeleton protocol with the potential pitfalls identified. It is designed to augment the protocols provided by Affymetrix.

Methods for processing microarray data

Quality control must be maintained at every step of a microarray experiment, from RNA isolation through statistical evaluation. Here we provide suggestions for analyzing microarray data. Because the utility of the results depends directly on the design of the experiment, the first critical step is to ensure that the experiment can be properly analyzed and interpreted. What is the biological question? What is the best way to perform the experiment? How many replicates will be required to obtain the desired statistical resolution? Next, the samples must be prepared, pass quality controls for integrity and representation, and be hybridized and scanned. Also, slides with defects, missing data, high background, or weak signal must be rejected. Data from individual slides must be normalized and combined so that the data are as free of systematic bias as possible. The third phase is to apply statistical filters and tests to the data to determine genes (1) expressed above background, (2) whose expression level changes in different samples, and (3) whose RNA-processing patterns or protein associations change. Next, a subset of the data should be validated by an alternative method, such as reverse transcription-polymerase chain reaction (RT-PCR). Provided that this endorses the general conclusions of the array analysis, gene sets whose expression, splicing, polyadenylation, protein binding, etc. change in different samples can be classified with respect to function, sequence motif properties, as well as other categories to extract hypotheses for their biological roles and regulatory logic.

Microarray slide hybridization using fluorescently labeled cDNA

Microarray hybridization is used to determine the amount and genomic origins of RNA molecules in an experimental sample. Unlabeled probe sequences for each gene or gene region are printed in an array on the surface of a slide, and fluorescently labeled cDNA derived from the RNA target is hybridized to it. This protocol describes a blocking and hybridization protocol for microarray slides. The blocking step is particular to the chemistry of "CodeLink" slides, but it serves to remind us that almost every kind of microarray has a treatment step that occurs after printing but before hybridization. We recommend making sure of the precise treatment necessary for the particular chemistry used in the slides to be hybridized because the attachment chemistries differ significantly. Hybridization is similar to northern or Southern blots, but on a much smaller scale.

Preparation of fluorescent-dye-labeled cDNA from RNA for microarray hybridization

This protocol describes how to prepare fluorescently labeled cDNA for hybridization to microarrays. It consists of two steps: first, a mixture of anchored oligo(dT) and random hexamers is used to prime amine-modified cDNA synthesis by reverse transcriptase using a modified deoxynucleotide with a reactive amine group (aminoallyl-dUTP) and an RNA sample as a template. Second, the cDNA is purified and exchanged into bicarbonate buffer so that the amine groups in the cDNA react with the dye N-hydroxysuccinimide (NHS) esters, covalently joining the dye to the cDNA. The dye-coupled cDNA is purified again, and the amount of dye incorporated per microgram of cDNA is determined.

Safer one-pot synthesis of the 'SHAPE' reagent 1-methyl-7-nitroisatoic anhydride (1m7)

Estimating the reactivity of 2'-hydroxyl groups along an RNA chain of interest aids in the modeling of the folded RNA structure; flexible loops tend to be reactive, whereas duplex regions are generally not. Among the most useful reagents for probing 2'-hydroxyl reactivity is 1-methyl-7-nitroisatoic anhydride (1m7), but the absence of a reliable, inexpensive source has prevented widespread adoption. An existing protocol for the conversion of an inexpensive precursor 4-nitroisatoic anhydride (4NIA) recommends the use of NaH in dimethylformamide (DMF), a reagent combination that most molecular biology labs are not equipped to handle, and that does not scale safely in any case. Here we describe a safer, one-pot method for bulk conversion of 4NIA to 1m7 that reduces costs and bypasses the use of NaH. We show that 1m7 produced by this method is free of side products and can be used to probe RNA structure in vitro.

High-yield synthesis of RNA using T7 RNA polymerase and plasmid DNA or oligonucleotide templates

When a large amount (milligram quantities) of RNA is desired (e.g., for making RNAs for RNA interference studies or for structural analyses), it is advantageous to use optimized transcription conditions. The method described here uses high millimolar concentrations of rNTPs, pyrophosphatase, and a specially formulated high-yield transcription buffer. These conditions prevent product (pyrophosphate) inhibition of RNA polymerase, which permits the accumulation of large quantities of RNA.

Preparation of cell-free splicing extracts from Saccharomyces cerevisiae

Much of our understanding of the mechanism of splicing comes from the analysis of cell extracts able to carry out splicing complex formation and splicing reactions in vitro using exogenously added synthetic model pre-mRNA transcripts. This protocol describes the preparation of whole-cell extracts from the budding yeast Saccharomyces cerevisiae. These extracts can be used to dissect the biochemical steps of the splicing reaction and to determine the macromolecules, cofactors, and substrate features necessary for successful splicing.

Analysis of splicing complexes on native gels

Splicing requires a complex set of ATP-dependent macromolecular associations that lead to the rearrangement of just a few covalent bonds in the pre-mRNA substrate. Seeing only the covalent bonds breaking and forming is seeing only a very small part of the process. Analysis of native splicing complexes into which the radiolabeled substrate has been assembled, but not necessarily completely reacted, has provided a good understanding of the process. This protocol describes a gel method for detecting and analyzing yeast splicing complexes formed in vitro on a radiolabeled pre-mRNA substrate. Complexes formed during the splicing reaction are quenched by dilution and addition of an excess of RNA, which is thought to strip nonspecifically bound proteins from the labeled substrate RNA. After loading on a low-percentage, low-cross-linking ratio composite agarose-acrylamide gel (in 10% glycerol), labeled bands are detected. These can be extracted and shown to contain small nuclear RNAs (snRNAs) and partly reacted pre-mRNA.

Analysis of splicing in vitro using extracts of Saccharomyces cerevisiae

In vitro splicing studies are a powerful means of investigating the requirements and mechanisms of action of the many components of the splicing apparatus. The ability to add and subtract components, purify activities, and reconstitute activity, as well as to expose the apparatus to chemical probes of various types, allows a far more mechanistically detailed view of the process to emerge than is available from genetic or in vivo studies alone. Two kinds of activities are assayed during in vitro splicing. The first concerns the chemical conversion of the substrate pre-mRNA into splicing intermediates and products and is usually visualized using a labeled substrate followed by separation on a denaturing gel. The second concerns the assembly of noncovalent complexes between the substrate and the myriad components of the splicing apparatus. This is also visualized using a labeled substrate, but the separation of complexes is achieved using native gel electrophoresis or gradient sedimentation. In this protocol, we describe the splicing reaction and its preparation for analysis by denaturing gels and native splicing complex gels. We also provide conditions for depletion of ATP, a critical cofactor that is hydrolyzed during numerous key steps in spliceosome assembly and splicing progression.

Competition between pre-mRNAs for the splicing machinery drives global regulation of splicing

During meiosis in yeast, global splicing efficiency increases and then decreases. Here we provide evidence that splicing improves due to reduced competition for the splicing machinery. The timing of this regulation corresponds to repression and reactivation of ribosomal protein genes (RPGs) during meiosis. In vegetative cells, RPG repression by rapamycin treatment also increases splicing efficiency. Downregulation of the RPG-dedicated transcription factor gene IFH1 genetically suppresses two spliceosome mutations, prp11-1 and prp4-1, and globally restores splicing efficiency in prp4-1 cells. We conclude that the splicing apparatus is limiting and that pre-messenger RNAs compete. Splicing efficiency of a pre-mRNA therefore depends not just on its own concentration and affinity for limiting splicing factor(s), but also on those of competing pre-mRNAs. Competition between RNAs for limiting processing factors appears to be a general condition in eukaryotes for a variety of posttranscriptional control mechanisms including microRNA (miRNA) repression, polyadenylation, and splicing.

A high-throughput splicing assay identifies new classes of inhibitors of human and yeast spliceosomes

The spliceosome is the macromolecular machine responsible for pre-mRNA splicing, an essential step in eukaryotic gene expression. During splicing, myriad subunits join and leave the spliceosome as it works on the pre-mRNA substrate. Strikingly, there are very few small molecules known to interact with the spliceosome. Splicing inhibitors are needed to capture transient spliceosome conformations and probe important functional components. Such compounds may also have chemotherapeutic applications, as links between splicing and cancer are increasingly uncovered. To identify new splicing inhibitors, we developed a high-throughput assay for in vitro splicing using a reverse transcription followed by quantitative PCR readout. In a pilot screen of 3080 compounds, we identified three small molecules that inhibit splicing in HeLa extract by interfering with different stages of human spliceosome assembly. Two of the compounds similarly affect spliceosomes in yeast extracts, suggesting selective targeting of conserved components. By examining related molecules, we identified chemical features required for the activity of two of the splicing inhibitors. In addition to verifying our assay procedure and paving the way to larger screens, these studies establish new compounds as chemical probes for investigating the splicing machinery.

Quaking and PTB control overlapping splicing regulatory networks during muscle cell differentiation

Alternative splicing contributes to muscle development, but a complete set of muscle-splicing factors and their combinatorial interactions are unknown. Previous work identified ACUAA ("STAR" motif) as an enriched intron sequence near muscle-specific alternative exons such as Capzb exon 9. Mass spectrometry of myoblast proteins selected by the Capzb exon 9 intron via RNA affinity chromatography identifies Quaking (QK), a protein known to regulate mRNA function through ACUAA motifs in 3' UTRs. We find that QK promotes inclusion of Capzb exon 9 in opposition to repression by polypyrimidine tract-binding protein (PTB). QK depletion alters inclusion of 406 cassette exons whose adjacent intron sequences are also enriched in ACUAA motifs. During differentiation of myoblasts to myotubes, QK levels increase two- to threefold, suggesting a mechanism for QK-responsive exon regulation. Combined analysis of the PTB- and QK-splicing regulatory networks during myogenesis suggests that 39% of regulated exons are under the control of one or both of these splicing factors. This work provides the first evidence that QK is a global regulator of splicing during muscle development in vertebrates and shows how overlapping splicing regulatory networks contribute to gene expression programs during differentiation.

Fragmentation of whole-transcriptome RNA using E. coli RNase III

High-throughput sequencing (HTS) methods can provide short sequence reads for many millions of individual molecules in a sample, allowing the use of sequencing to measure the abundance of RNA molecules. To quantify the amount of a particular sequence in a sample of large RNAs (e.g., mRNAs), it is important to fragment the RNA into short pieces that can be ligated to oligonucleotides that allow polymerase chain reaction (PCR) amplification and sequencing. The most desired end structure of RNA for such ligation steps is a 5' phosphate and a 3' OH. Thus, enzymes that leave these groups after cleavage are of particular utility, avoiding the need to dephosphorylate the 3' end with phosphatases or phosphorylate the 5' end with kinase before proceeding. One such enzyme, RNase III, is widely available. Although it primarily cuts duplex RNA, this specificity is salt- and concentration-dependent, and many RNAs that lack strong extended duplexes are nonetheless susceptible to cleavage at many spots. RNA fragmentation by RNase III does not seem to grossly affect the distribution of RNA sequencing reads. Thus, it has become a standard method for creating nominally representative pools of transcriptome sequences with 5' phosphates and 3' OH for library construction. Three steps in preparing fragmented transcriptome RNA for sequencing library construction are described here: (1) fragmenting the RNA with RNase III to the extent that ~60-100-nucleotide fragments are created, (2) purifying the RNA from the RNase III reaction, and (3) analyzing the digestion products for their suitability in library production.

Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing

SR proteins are well-characterized RNA binding proteins that promote exon inclusion by binding to exonic splicing enhancers (ESEs). However, it has been unclear whether regulatory rules deduced on model genes apply generally to activities of SR proteins in the cell. Here, we report global analyses of two prototypical SR proteins, SRSF1 (SF2/ASF) and SRSF2 (SC35), using splicing-sensitive arrays and CLIP-seq on mouse embryo fibroblasts (MEFs). Unexpectedly, we find that these SR proteins promote both inclusion and skipping of exons in vivo, but their binding patterns do not explain such opposite responses. Further analyses reveal that loss of one SR protein is accompanied by coordinated loss or compensatory gain in the interaction of other SR proteins at the affected exons. Therefore, specific effects on regulated splicing by one SR protein actually depend on a complex set of relationships with multiple other SR proteins in mammalian genomes.

Rbfox1 downregulation and altered calpain 3 splicing by FRG1 in a mouse model of Facioscapulohumeral muscular dystrophy (FSHD)

Facioscapulohumeral muscular dystrophy (FSHD) is a common muscle disease whose molecular pathogenesis remains largely unknown. Over-expression of FSHD region gene 1 (FRG1) in mice, frogs, and worms perturbs muscle development and causes FSHD–like phenotypes. FRG1 has been implicated in splicing, and we asked how splicing might be involved in FSHD by conducting a genome-wide analysis in FRG1 mice. We find that splicing perturbations parallel the responses of different muscles to FRG1 over-expression and disease progression. Interestingly, binding sites for the Rbfox family of splicing factors are over-represented in a subset of FRG1-affected splicing events. Rbfox1 knockdown, over-expression, and RNA-IP confirm that these are direct Rbfox1 targets. We find that FRG1 is associated to the Rbfox1 RNA and decreases its stability. Consistent with this, Rbfox1 expression is down-regulated in mice and cells over-expressing FRG1 as well as in FSHD patients. Among the genes affected is Calpain 3, which is mutated in limb girdle muscular dystrophy, a disease phenotypically similar to FSHD. In FRG1 mice and FSHD patients, the Calpain 3 isoform lacking exon 6 (Capn3 E6–) is increased. Finally, Rbfox1 knockdown and over-expression of Capn3 E6- inhibit muscle differentiation. Collectively, our results suggest that a component of FSHD pathogenesis may arise by over-expression of FRG1, reducing Rbfox1 levels and leading to aberrant expression of an altered Calpain 3 protein through dysregulated splicing.

Alternative splicing is a major contributor to the complexity of human cells, and its disruption can lead to a wide range of human disorders. FSHD is one of the most important muscle diseases. While muscle differentiation defects have been widely reported in the disease, the molecular mechanisms responsible are largely unknown. We found that expression of the alternative splicing factor Rbfox1 is a direct FRG1 target, and its expression decreased in the muscles of a mouse model of FSHD and FSHD patients. Moreover, alternative splicing of Calpain 3, encoding for a protease involved in muscle differentiation, is regulated by Rbfox1 and is altered in the muscles of the mouse model of FSHD and FSHD patients. Interestingly, we found that Rbfox1 is required for muscle differentiation and that this activity is likely mediated by Calpain 3 alternative splicing. Hence, our results suggest that decreased expression of Rbfox1 and aberrant Calpain 3 splicing contribute to the muscle differentiation defects of FSHD patients.

Structural analysis of the quaking homodimerization interface

Quaking (QkI) is a prototypical member of the STAR (signal transducer and activator of RNA) protein family, which plays key roles in posttranscriptional gene regulation by controlling mRNA translation, stability and splicing. QkI-5 has been shown to regulate mRNA expression in the central nervous system, but little is known about its roles in other tissues. STAR proteins function as dimers and bind to bipartite RNA sequences; however, the structural and functional roles of homodimerization and heterodimerization are still unclear. Here, we present the crystal structure of the QkI dimerization domain, which adopts a similar stacked helix-turn-helix arrangement as its homologs GLD-1 (germ line development defective-1) and Sam68 (Src-associated protein during mitosis, 68kDa) but differs by an additional helix inserted in the dimer interface. Variability of the dimer interface residues likely ensures selective homodimerization by preventing association with non-cognate STAR family proteins in the cell. Mutations that inhibit dimerization also significantly impair RNA binding in vitro, alter QkI-5 protein levels and impair QkI function in a splicing assay in vivo. Together, our results indicate that a functional Qua1 homodimerization domain is required for QkI-5 function in mammalian cells.

Muscleblind-like 2-mediated alternative splicing in the developing brain and dysregulation in myotonic dystrophy

The RNA-mediated disease model for myotonic dystrophy (DM) proposes that microsatellite C(C)TG expansions express toxic RNAs that disrupt splicing regulation by altering MBNL1 and CELF1 activities. While this model explains DM manifestations in muscle, less is known about the effects of C(C)UG expression on the brain. Here, we report that Mbnl2 knockout mice develop several DM-associated central nervous system (CNS) features including abnormal REM sleep propensity and deficits in spatial memory. Mbnl2 is prominently expressed in the hippocampus and Mbnl2 knockouts show a decrease in NMDA receptor (NMDAR) synaptic transmission and impaired hippocampal synaptic plasticity. While Mbnl2 loss did not significantly alter target transcript levels in the hippocampus, misregulated splicing of hundreds of exons was detected using splicing microarrays, RNA-seq, and HITS-CLIP. Importantly, the majority of the Mbnl2-regulated exons examined were similarly misregulated in DM. We propose that major pathological features of the DM brain result from disruption of the MBNL2-mediated developmental splicing program.

Isolation of total RNA from yeast cell cultures

This article describes two procedures for isolating total RNA from yeast cell cultures. The first allows the convenient isolation of total RNA from early log-phase cultures (vegetative cells). RNA isolated in this way is intact and sufficiently pure for use in microarray experiments, primer extension, and RNase protection mapping. With additional treatment to remove contaminating genomic DNA, the preparation is suitable for reverse transcription-polymerase chain reaction (RT-PCR), quantitative PCR (qPCR), cDNA library construction, high-throughput sequencing of RNA, or other manipulations. However, compared to vegetative cells, the isolation of RNA from cells late in meiosis (asci and ascospores) requires additional effort. This is because a tough cell wall composed of heavily cross-linked polysaccharides and proteins is built around the four spores during meiosis and ascospore development. Therefore, an alternative protocol is presented for extracting RNA from cells late in meiosis. This alternative may also be preferable for cells from stationary cultures or from yeast strains and other fungal species isolated from the environment.

Bacterial RNA isolation

In this bacterial RNA isolation protocol, an "RNA-protective" treatment is followed by lysozyme digestion of the peptidoglycan component of the cell wall. EDTA promotes the loss of the outer membrane of Gram-negative bacteria and allows the lysozyme better access to the peptidoglycan. Cells begin to lyse during digestion in hypotonic lysozyme buffer and lysis is completed by sodium dodecyl sulfate (SDS) and hot phenol:chloroform:isoamyl alcohol (PCA) extraction. SDS and hot phenol disrupt membranes, denature protein (including RNase), and strip proteins from RNA. The separation of the organic phase from the aqueous phase is achieved using Phase Lock Gel, an inert material with a density intermediate between the organic and aqueous samples. The sample is split into three phases: from bottom to top, these are phenol and chloroform (organic phase), the inert gel with the interface material, and the aqueous phase with the RNA. The gel acts as a physical barrier between the sample and the organic phase plus interface. Following organic extraction, the RNA is concentrated by ethanol precipitation.

Muscleblind-like 1 knockout mice reveal novel splicing defects in the myotonic dystrophy brain

Myotonic dystrophy type 1 (DM1) is a multi-systemic disorder caused by a CTG trinucleotide repeat expansion (CTG(exp)) in the DMPK gene. In skeletal muscle, nuclear sequestration of the alternative splicing factor muscleblind-like 1 (MBNL1) explains the majority of the alternative splicing defects observed in the HSA(LR) transgenic mouse model which expresses a pathogenic range CTG(exp). In the present study, we addressed the possibility that MBNL1 sequestration by CUG(exp) RNA also contributes to splicing defects in the mammalian brain. We examined RNA from the brains of homozygous Mbnl1(ΔE3/ΔE3) knockout mice using splicing-sensitive microarrays. We used RT-PCR to validate a subset of alternative cassette exons identified by microarray analysis with brain tissues from Mbnl1(ΔE3/ΔE3) knockout mice and post-mortem DM1 patients. Surprisingly, splicing-sensitive microarray analysis of Mbnl1(ΔE3/ΔE3) brains yielded only 14 candidates for mis-spliced exons. While we confirmed that several of these splicing events are perturbed in both Mbnl1 knockout and DM1 brains, the extent of splicing mis-regulation in the mouse model was significantly less than observed in DM1. Additionally, several alternative exons, including Grin1 exon 4, App exon 7 and Mapt exons 3 and 9, which have previously been reported to be aberrantly spliced in human DM1 brain, were spliced normally in the Mbnl1 knockout brain. The sequestration of MBNL1 by CUG(exp) RNA results in some of the aberrant splicing events in the DM1 brain. However, we conclude that other factors, possibly other MBNL proteins, likely contribute to splicing mis-regulation in the DM1 brain.