Specific SR protein-dependent splicing substrates identified through genomic SELEX

Therapeutic manipulation of IKBKAP mis-splicing with a small molecule to cure familial dysautonomia

Approximately half of genetic disease-associated mutations cause aberrant splicing. However, a widely applicable therapeutic strategy to splicing diseases is yet to be developed. Here, we analyze the mechanism whereby IKBKAP-familial dysautonomia (FD) exon 20 inclusion is specifically promoted by a small molecule splice modulator, RECTAS, even though IKBKAP-FD exon 20 has a suboptimal 5' splice site due to the IVS20 + 6 T > C mutation. Knockdown experiments reveal that exon 20 inclusion is suppressed in the absence of serine/arginine-rich splicing factor 6 (SRSF6) binding to an intronic splicing enhancer in intron 20. We show that RECTAS directly interacts with CDC-like kinases (CLKs) and enhances SRSF6 phosphorylation. Consistently, exon 20 splicing is bidirectionally manipulated by targeting cellular CLK activity with RECTAS versus CLK inhibitors. The therapeutic potential of RECTAS is validated in multiple FD disease models. Our study indicates that small synthetic molecules affecting phosphorylation state of SRSFs is available as a new therapeutic modality for mechanism-oriented precision medicine of splicing diseases.

Key Aspects of Nucleic Acid Library Design for in Vitro Selection

Nucleic acid aptamers capable of selectively recognizing their target molecules have nowadays been established as powerful and tunable tools for biospecific applications, be it therapeutics, drug delivery systems or biosensors. It is now generally acknowledged that in vitro selection enables one to generate aptamers to almost any target of interest. However, the success of selection and the affinity of the resulting aptamers depend to a large extent on the nature and design of an initial random nucleic acid library. In this review, we summarize and discuss the most important features of the design of nucleic acid libraries for in vitro selection such as the nature of the library (DNA, RNA or modified nucleotides), the length of a randomized region and the presence of fixed sequences. We also compare and contrast different randomization strategies and consider computer methods of library design and some other aspects.

Single neuron transcriptomics identify SRSF/SR protein B52 as a regulator of axon growth and Choline acetyltransferase splicing

We removed single identified neurons from living Drosophila embryos to gain insight into the transcriptional control of developing neuronal networks. The microarray analysis of the transcriptome of two sibling neurons revealed seven differentially expressed transcripts between both neurons (threshold: log₂1.4). One transcript encodes the RNA splicing factor B52. Loss of B52 increases growth of axon branches. B52 function is also required for Choline acetyltransferase (ChAT ) splicing. At the end of embryogenesis, loss of B52 function impedes splicing of ChAT, reduces acetylcholine synthesis, and extends the period of uncoordinated muscle twitches during larval hatching. ChAT regulation by SRSF proteins may be a conserved feature since changes in SRSF5 expression and increased acetylcholine levels in brains of bipolar disease patients have been reported recently.

Aptamers: A Feasible Technology in Cancer Immunotherapy

Aptamers are single-chained RNA or DNA oligonucleotides (ODNs) with three-dimensional folding structures which allow them to bind to their targets with high specificity. Aptamers normally show affinities comparable to or higher than that of antibodies. They are chemically synthesized and therefore less expensive to manufacture and produce. A variety of aptamers described to date have been shown to be reliable in modulating immune responses against cancer by either blocking or activating immune receptors. Some of them have been conjugated to other molecules to target the immune system and reduce off-target side effects. Despite the success of first-line treatments against cancer, the elevated number of relapsing cases and the tremendous side effects shown by the commonly used agents hinder conventional treatments against cancer. The advantages provided by aptamers could enhance the therapeutic index of a given strategy and therefore enhance the antitumor effect. Here we recapitulate the provided benefits of aptamers with immunomodulatory activity described to date in cancer therapy and the benefits that aptamer-based immunotherapy could provide either alone or combined with first-line treatments in cancer therapy.

The Wright stuff: reimagining path analysis reveals novel components of the sex determination hierarchy in Drosophila melanogaster

BACKGROUND

The Drosophila sex determination hierarchy is a classic example of a transcriptional regulatory hierarchy, with sex-specific isoforms regulating morphology and behavior. We use a structural equation modeling approach, leveraging natural genetic variation from two studies on Drosophila female head tissues--DSPR collection (596 F1-hybrids from crosses between DSPR sub-populations) and CEGS population (75 F1-hybrids from crosses between DGRP/Winters lines to a reference strain w1118)--to expand understanding of the sex hierarchy gene regulatory network (GRN). This approach is completely generalizable to any natural population, including humans.

RESULTS

We expanded the sex hierarchy GRN adding novel links among genes, including a link from fruitless (fru) to Sex-lethal (Sxl) identified in both populations. This link is further supported by the presence of fru binding sites in the Sxl locus. 754 candidate genes were added to the pathway, including the splicing factors male-specific lethal 2 and Rm62 as downstream targets of Sxl which are well-supported links in males. Independent studies of doublesex and transformer mutants support many additions, including evidence for a link between the sex hierarchy and metabolism, via Insulin-like receptor.

CONCLUSIONS

The genes added in the CEGS population were enriched for genes with sex-biased splicing and components of the spliceosome. A common goal of molecular biologists is to expand understanding about regulatory interactions among genes. Using natural alleles we can not only identify novel relationships, but using supervised approaches can order genes into a regulatory hierarchy. Combining these results with independent large effect mutation studies, allows clear candidates for detailed molecular follow-up to emerge.

Identification of Nucleic Acid High Affinity Binding Sequences of Proteins by SELEX

A technique is described for the identification of nucleic acid sequences bound with high affinity by proteins or by other molecules suitable for a partitioning assay. Here, a histidine-tagged protein is allowed to interact with a pool of nucleic acids and the protein-nucleic acid complexes formed are retained on a Ni-NTA matrix. Nucleic acids with a low level of recognition by the protein are washed away. The pool of recovered nucleic acids is amplified by the polymerase chain reaction and is submitted to further rounds of selection. Each round of selection increases the proportion of sequences that are avidly bound by the protein of interest. The cloning and sequencing of these sequences finally completes their identification.

Aptamer and its potential applications for food safety

Accompanied by industrial globalization, rapid urbanization, and population increment, mass production and staple trading for food consumption are upsoaring continuously, foodborne disease resulted from various food safety issues is currently a crucial public health concern worldwide, which has not only created a great burden on both economy and society, but also greatly threatened the sustainability of mankind's livelihood and human reproduction. In order to better ensure food safety and thus effectively curb the occurrence of foodborne diseases, the development and evolving of inspection strategies are indispensable measures for quality assurance and conformity assessment. Nowadays, as complementary measures to and with advantageous merits over classic analytical methods, highly specific and selective aptamer-based assays have found their increasingly important roles in various domains of food analysis. This critical review summarizes the advantages of aptamer as compared with antibody, introduces important evolving variants of systematic evolution of ligands by exponential enrichment (SELEX), and presents an overview of potential aptamer applications for food safety.

Shrimp Dscam and its cytoplasmic tail splicing activator serine/arginine (SR)-rich protein B52 were both induced after white spot syndrome virus challenge

The serine/arginine (SR)-rich protein family is phylogenetically conserved and plays significant roles in mRNA maturation, including alternative splicing (AS). In Drosophila, SR protein B52 functions as a splicing activator to regulate AS events in several genes, including the Down syndrome cell adhesion molecule (Dscam). In this study, the B52 gene from Litopenaeus vannamei (LvB52) was isolated and characterized. The open reading frame of LvB52 contains 1149 bp encoding 382 amino acids. The deduced LvB52 protein includes two RNA recognition motifs (RRM) at the N terminus and an arginine/serine rich domain (RS rich domain) at the C terminus, and thus shows the expected RRM1-RRM2-RS domain architecture. Tissue tropism analysis revealed that LvB52 is expressed in most tissues and at high levels in stomach and muscle. After white spot syndrome virus (WSSV) infection, a parallel increase in the expression of total LvDscam, tail-less LvDscam, membrane-bound LvDscam and LvB52 was observed after 24 hpi. Conversely, there was no obvious change in the expression of the AS repressor Lvhrp36. In vivo dsRNA silencing of LvB52 induced element 3 exclusion in the LvDscam cytoplasmic tail, but no abnormal exclusions in the Ig2-Ig3 region or the transmembrane region. We also found that the exon of the Ig7 region was quite often excluded, even in normal shrimp, and that LvB52 silencing was associated with a decrease in the variability of this region. Taken together, our data suggest that LvB52 acts as a splicing activator that regulates AS events in LvDscam.

Identification of RNA targets for the nuclear multidomain cyclophilin atCyp59 and their effect on PPIase activity

AtCyp59 is a multidomain cyclophilin containing a peptidyl-prolyl cis/trans isomerase (PPIase) domain and an evolutionarily highly conserved RRM domain. Deregulation of this class of cyclophilins has been shown to affect transcription and to influence phosphorylation of the C-terminal repeat domain of the largest subunit of the RNA polymerase II. We used a genomic SELEX method for identifying RNA targets of AtCyp59. Analysis of the selected RNAs revealed an RNA-binding motif (G[U/C]N[G/A]CC[A/G]) and we show that it is evolutionarily conserved. Binding to this motif was verified by gel shift assays in vitro and by RNA immunopreciptation assays of AtCyp59 in vivo. Most importantly, we show that binding also occurs on unprocessed transcripts in vivo and that binding of specific RNAs inhibits the PPIase activity of AtCyp59 in vitro. Surprisingly, genome-wide analysis showed that the RNA motif is present in about 70% of the annotated transcripts preferentially in exons. Taken together, the available data suggest that these cyclophilins might have an important function in transcription regulation.

VERSE: a varying effect regression for splicing elements discovery

Identification of splicing regulatory elements (SREs) deserves special attention because these cis-acting short sequences are vital parts of splicing code. The fact that a variety of other biological signals cooperatively govern the splicing pattern indicates the necessity of developing novel tools to incorporate information from multiple sources to improve splicing factor binding sites prediction. Under this context, we proposed a Varying Effect Regression for Splicing Elements (VERSE) to discover intronic SREs in the proximity of exon junctions by integrating other biological features. As a result, 1562 intronic SREs were identified in 16 human tissues, many of which overlapped with experimentally verified binding motifs for several well-known splicing factors, including FOX-1, PTB, hnRNP A/B, hnRNP F/H, and so on. The discovered tissue, region, and conservation preferences of the putative motifs demonstrate that splice site selection is a complicated process that needs subtle and delicate regulation. VERSE may serve as a powerful tool to not only discover SREs by incorporating additional informative signals but also precisely quantify their varying contribution under different biological contexts.

Regulated functional alternative splicing in Drosophila

Alternative splicing expands the coding capacity of metazoan genes, and it was largely genetic studies in the fruit-fly Drosophila melanogaster that established the principle that regulated alternative splicing results in tissue- and stage-specific protein isoforms with different functions in development. Alternative splicing is particularly prominent in germ cells, muscle and the central nervous system where it modulates the expression of various proteins including cell-surface molecules and transcription factors. Studies in flies have given us numerous insights into alternative splicing in terms of upstream regulation, the exquisite diversity of their forms and the key differential cellular functions of alternatively spliced gene products. The current inundation of transcriptome sequencing data from Drosophila provides an unprecedented opportunity to gain a comprehensive view of alternative splicing.

Isolation of putative FOXO1 genomic elements using an improved in vitro technique to isolate genomic regulatory sequences

The regulation of gene expression drives many biological processes and alterations in normal regulation are integral in the development of the diseased state. Therefore, the ability to screen genomic DNA for direct targets of DNA binding proteins (DNA-BP) would provide valuable information about the mechanisms underlying these processes. At present chromatin immunoprecipitation (ChIP) and its variants are the primary methods for identifying regulatory elements. However, some DNA-BPs, such as the winged-helix transcription factor FOXO1, are difficult to ChIP thereby detracting from the use of this technique as a nonbiased screen to isolate regulatory sequences. In this report we use an improved in vitro method to Pull Out Regulatory Elements (PORE), which uses purified protein with a stable genomic library to isolate regulatory elements directly bound by a DNA-BP, to identify putative FOXO1 genomic regulatory sequences. We first validate this technique using two known DNA-BP (FOXO1 and Pax3) by demonstrating their ability to bind and amplify identified promoter elements when present in a genomic DNA context or when present in the context of our stable genomic library. Subsequent use of this technique with FOXO1 isolated regulatory elements associated with several genes known to be regulated in a FOXO1-dependent manner and multiple genes whose biological functions are consistent with the known biological functions of FOXO1 proving that the in vitro PORE is a valuable and easy to use alternative to ChIP for the isolation of genomic regulatory elements.

Alternative splicing: global insights

Following the original reports of pre-mRNA splicing in 1977, it was quickly realized that splicing together of different combinations of splice sites--alternative splicing--allows individual genes to generate more than one mRNA isoform. The full extent of alternative splicing only began to be revealed once large-scale genome and transcriptome sequencing projects began, rapidly revealing that alternative splicing is the rule rather than the exception. Recent technical innovations have facilitated the investigation of alternative splicing at a global scale. Splice-sensitive microarray platforms and deep sequencing allow quantitative profiling of very large numbers of alternative splicing events, whereas global analysis of the targets of RNA binding proteins reveals the regulatory networks involved in post-transcriptional gene control. Combined with sophisticated computational analysis, these new approaches are beginning to reveal the so-called 'RNA code' that underlies tissue and developmentally regulated alternative splicing, and that can be disrupted by disease-causing mutations.

Analysis of RNA folding by native polyacrylamide gel electrophoresis

Polyacrylamide gel electrophoresis under native conditions (native PAGE) is a well-established and versatile method for probing nucleic acid conformation and nucleic acid-protein interactions. Native PAGE has been used to measure RNA folding equilibria and kinetics under a wide variety of conditions. Advantages of this method are its adaptability, absolute determination of reaction endpoints, and direct analysis of conformational hetereogeneity within a sample. Native PAGE is also useful for resolving ligand-induced structural changes.

Selection and elution of aptamers using nanoporous sol-gel arrays with integrated microheaters

RNA and DNA aptamers that bind to target molecules with high specificity and affinity have been a focus of diagnostics and therapeutic research. These aptamers are obtained by SELEX (Systematic Evolution of Ligands by EXponential enrichment) often requiring more than 10 successive cycles of selection and amplification, where each cycle normally takes 2 days per cycle of SELEX. Here, we have demonstrated the use of sol-gel arrays of proteins in a microfluidic system for efficient selection of RNA aptamers against multiple target molecules. The microfluidic chip incorporates five sol-gel binding droplets, within which specific target proteins are imbedded. The droplets are patterned on top of individually addressable electrical microheaters used for selective elution of aptamers bound to target proteins in the sol-gel droplets. We demonstrate that specific aptamers bind their respective protein targets and can be selectively eluted by micro-heating. Finally, our microfluidic SELEX system greatly improved selection efficiency, reducing the number of selection cycles needed to produce high affinity aptamers. The process is readily scalable to larger arrays of sol-gel-embedded proteins. To our knowledge, this is the first demonstration of a chip-based selection of aptamers using microfluidics, thereby allowing development of a high throughput and efficient SELEX procedures.

The SR protein family of splicing factors: master regulators of gene expression

The SR protein family comprises a number of phylogenetically conserved and structurally related proteins with a characteristic domain rich in arginine and serine residues, known as the RS domain. They play significant roles in constitutive pre-mRNA splicing and are also important regulators of alternative splicing. In addition they participate in post-splicing activities, such as mRNA nuclear export, nonsense-mediated mRNA decay and mRNA translation. These wide-ranging roles of SR proteins highlight their importance as pivotal regulators of mRNA metabolism, and if these functions are disrupted, developmental defects or disease may result. Furthermore, animal models have shown a highly specific, non-redundant role for individual SR proteins in the regulation of developmental processes. Here, we will review the current literature to demonstrate how SR proteins are emerging as one of the master regulators of gene expression.

Genome-wide analysis of RNA-protein interactions in plants

RNA-protein interactions profoundly impact organismal development and function through their contributions to the basal gene expression machineries and their regulation of post-transcriptional processes. The repertoire of predicted RNA binding proteins (RBPs) in plants is particularly large, suggesting that the RNA-protein interactome in plants may be more complex and dynamic even than that in metazoa. To dissect RNA-protein interaction networks, it is necessary to identify the RNAs with which each RBP interacts and to determine how those interactions influence RNA fate and downstream processes. Identification of the native RNA ligands of RBPs remains a challenge, but several high-throughput methods for the analysis of RNAs that copurify with specific RBPs from cell extract have been reported recently. This chapter reviews approaches for defining the native RNA ligands of RBPs on a genome-wide scale and provides a protocol for a method that has been used to this end for RBPs that localize to the chloroplast.

Identification of nucleic acid high-affinity binding sequences of proteins by SELEX

A technique is described for the identification of nucleic acid sequences bound with high affinity by proteins or by other molecules suitable for a partitioning assay. Here, a histidine-tagged protein is allowed to interact with a pool of nucleic acids and the protein-nucleic acid complexes formed are retained on a Ni-NTA matrix. Nucleic acids with a low level of recognition by the protein are washed away. The pool of recovered nucleic acids is amplified by the polymerase chain reaction and is submitted to further rounds of selection. Each round of selection increases the proportion of sequences that are avidly bound by the protein of interest. The cloning and sequencing of these sequences finally completes their identification.

A whole-genome approach to identifying protein binding sites: promoters in Methanocaldococcus (Methanococcus) jannaschii

We have adapted an electrophoretic mobility shift assay (EMSA) to isolate genomic DNA fragments that bind the archaeal transcription initiation factors TATA-binding protein (TBP) and transcription factor B (TFB) to perform a genome-wide search for promoters. Mobility-shifted fragments were cloned, tested for their ability to compete with known promoter-containing fragments for a limited concentration of transcription factors, and sequenced. We applied the method to search for promoters in the genome of Methanocaldococcus jannaschii. Selection was most efficient for promoters of tRNA genes and genes for several presumed small non-coding RNAs (ncRNA). Protein-coding gene promoters were dramatically underrepresented relative to their frequency in the genome. The repeated isolation of these genomic regions was partially rectified by including a hybridization-based screening. Sequence alignment of the affinity-selected promoters revealed previously identified TATA box, BRE, and the putative initiator element. In addition, the conserved bases immediately upstream and downstream of the BRE and TATA box suggest that the composition and structure of archaeal natural promoters are more complicated.

Searching for splicing motifs

Intron removal during pre-mRNA splicing in higher eukaryotes requires the accurate identification of the two splice sites at the ends of the exons, or exon definition. The sequences constituting the splice sites provide insufficient information to distinguish true splice sites from the greater number of false splice sites that populate transcripts. Additional information used for exon recognition resides in a large number of positively or negatively acting elements that lie both within exons and in the adjacent introns. The identification of such sequence motifs has progressed rapidly in recent years, such that extensive lists are now available for exonic splicing enhancers and exonic splicing silencers. These motifs have been identified both by empirical experiments and by computational predictions, the validity of the latter being confirmed by experimental verification. Molecular searches have been carried out either by the selection of sequences that bind to splicing factors, or enhance or silence splicing in vitro or in vivo. Computational methods have focused on sequences of 6 or 8 nucleotides that are over- or under-represented in exons, compared to introns or transcripts that do not undergo splicing. These various methods have sought to provide global definitions of motifs, yet the motifs are distinctive to the method used for identification and display little overlap. Astonishingly, at least three-quarters of a typical mRNA would be comprised of these motifs. A present challenge lies in understanding how the cell integrates this surfeit of information to generate what is usually a binary splicing decision.

Gene regulations in HBV-related liver cirrhosis closely correlate with disease severity

Liver cirrhosis (LC) is defined as comprising diffuse fibrosis and regenerating nodules of the liver. The biochemical and anatomical dysfunction in LC results from both reduced liver cell number and portal vascular derangement. Although several studies have investigated dysregulated genes in cirrhotic nodules, little is known about the genes implicated in the pathophysiologic change of LC or about their relationship with the degree of decompensation. Here, we applied cDNA microarray analysis using 38 HBsAg-positive LC specimens to identify the genes dysregulated in HBV-associated LC and to evaluate their relation to disease severity. Among 1063 known cancer- and apoptosis-related genes, we identified 104 genes that were significantly up- (44) or down- (60) regulated in LC. Interestingly, this subset of 104 genes was characteristically correlated with the degree of decompensation, called the Pugh-Child classification (20 Pugh-Child A, 10 Pugh-Child B, and 8 Pugh-Child C). Patient samples from Pugh-Child C exhibited a distinct pattern of gene expression relative to those of Pugh-Child A and B. Especially in Pugh-Child C, genes encoding hepatic proteins and metabolizing enzymes were significantly down-regulated, while genes encoding various molecules related to cell replication were up-regulated. Our results suggest that subsets of genes in liver cells correspond to the pathophysiologic change of LC according to disease severity and possibly to hepatocarcinogenesis.

Plant SR proteins and their functions

SR proteins are a family of splicing factors important for splice site recognition and spliceosome assembly. Their ability to bind to RNA and to interact with proteins as well identifies them as important players in splice site choice and alternative splicing. Plants possess twice as many SR proteins as animals, and some of the subfamilies are plant specific. Arabidopsis SR proteins are involved in different aspects of plant growth and development as well as in responses to environmental cues. The plant-specific subfamilies have been shown to be regulated by alternative splicing events, which are highly conserved in evolution. The tight regulation of splicing factors by alternative splicing might allow coordinated responses of their target genes.

Caffeine regulates alternative splicing in a subset of cancer-associated genes: a role for SC35

Alternative splicing of pre-mRNA contributes significantly to human proteomic complexity, playing a key role in development, gene expression and, when aberrant, human disease onset. Many of the factors involved in alternative splicing have been identified, but little is known about their regulation. Here we report that caffeine regulates alternative splicing of a subset of cancer-associated genes, including the tumor suppressor KLF6. This regulation is at the level of splice site selection, occurs rapidly and reversibly, and is concentration dependent. We have recapitulated caffeine-induced alternative splicing of KLF6 using a cell-based minigene assay and identified a "caffeine response element" within the KLF6 intronic sequence. Significantly, a chimeric minigene splicing assay demonstrated that this caffeine response element is functional in a heterologous context; similar elements exist within close proximity to caffeine-regulated exons of other genes in the subset. Furthermore, the SR splicing factor, SC35, was shown to be required for induction of the alternatively spliced KLF6 transcript. Importantly, SC35 is markedly induced by caffeine, and overexpression of SC35 is sufficient to mimic the effect of caffeine on KLF6 alternative splicing. Taken together, our data implicate SC35 as a key player in caffeine-mediated splicing regulation. This novel effect of caffeine provides a valuable tool for dissecting the regulation of alternative splicing of a large gene subset and may have implications with respect to splice variants associated with disease states.

Genomic systematic evolution of ligands by exponential enrichment (Genomic SELEX) for the identification of protein-binding RNAs independent of their expression levels

Genomic systematic evolution of ligands by exponential enrichment (Genomic SELEX) is an experimental procedure for the expression condition-independent identification of protein-binding RNAs. RNA libraries derived from genomic DNA are generated via random priming, PCR amplification and in vitro transcription. Libraries consist of genomic sequences of selected size, and fragments are flanked by constant sequences required for amplification and transcription. This RNA pool is then subjected to several rounds of selection and amplification to enrich for RNAs meeting the selection criteria. Various selection criteria are possible. Here we describe selection by affinity to a protein of interest. High-affinity ligands can then be cloned and sequenced to allow their identification. With this method, protein-binding RNAs can be discovered, nucleic acid-protein interactions can be identified, and whole protein-nucleic acid networks can be defined. This method is also suitable for discovering novel genes, including non-protein-coding RNAs, and it complements in silico approaches. It is better suited to detect protein-binding RNAs that are differentially expressed (and therefore absent from many tissues) and low-abundance RNAs than experimental procedures that start from the isolation of expressed RNAs. The protocol takes approximately 3 months to complete.

Evidence for a trade-off between translational efficiency and splicing regulation in determining synonymous codon usage in Drosophila melanogaster

In Drosophila melanogaster, synonymous codons corresponding to the most abundant cognate tRNAs are used more frequently, especially in highly expressed genes. Increased use of such "optimal" codons is considered an adaptation for translational efficiency. Need it always be the case that selection should favor the use of a translationally optimal codon? Here, we investigate one possible confounding factor, namely, the need to specify information in exons necessary to enable correct splicing. As expected from such a model, in Drosophila many codons show different usage near intron-exon boundaries versus exon core regions. However, this finding is in principle also consistent with Hill-Robertson effects modulating usage of translationally optimal codons. However, several results support the splice model over the translational selection model: 1) the trends in codon usage are strikingly similar to those in mammals in which codon usage near boundaries correlates with abundance in exonic splice enhancers (ESEs), 2) codons preferred near boundaries tend to be enriched for A and avoid C (conversely those avoided near boundaries prefer C rather than A), as expected were ESEs involved, and 3) codons preferred near boundaries are typically not translationally optimal. We conclude that usage of translationally optimal codons usage is compromised in the vicinity of splice junctions in intron-containing genes, to the effect that we observe higher levels of usage of translationally optimal codons at the center of exons. On the gene level, however, controlling for known correlates of codon bias, the impact on codon usage patterns is quantitatively small. These results have implications for inferring aspects of the mechanism of splicing given nothing more than a well-annotated genome.

The impact of alternative splicing in vivo: mouse models show the way

Alternative splicing is widely believed to have a major impact on almost all biological processes since it increases proteome complexity and thereby controls protein function. Recently, gene targeting in mice has been used to create in vivo models to study the regulation and consequences of alternative splicing. The evidence accumulated so far argues for a nonredundant, highly specific role of individual splicing factors in mammalian development, and furthermore, demonstrates the importance of distinct protein isoforms in vivo. In this review, we will compare phenotypes of mouse models for alternative splicing to crystallize common themes and to put them into perspective with the available in vitro data.

Identification of novel universal housekeeping genes by statistical analysis of microarray data

Housekeeping genes are widely used as internal controls in a variety of study types, including real time RT-PCR, microarrays, Northern analysis and RNase protection assays. However, even commonly used housekeeping genes may vary in stability depending on the cell type or disease being studied. Thus, it is necessary to identify additional housekeeping-type genes that show sample-independent stability. Here, we used statistical analysis to examine a large human microarray database, seeking genes that were stably expressed in various tissues, disease states and cell lines. We further selected genes that were expressed at different levels, because reference and target genes should be present in similar copy numbers to achieve reliable quantitative results. Real time RT-PCR amplification of three newly identified reference genes, CGI-119, CTBP1 and GOLGAl, alongside three well-known housekeeping genes, B2M, GAPD, and TUBB, confirmed that the newly identified genes were more stably expressed in individual samples with similar ranges. These results collectively suggest that statistical analysis of microarray data can be used to identify new candidate housekeeping genes showing consistent expression across tissues and diseases. Our analysis identified three novel candidate housekeeping genes (CGI-119, GOLGA1, and CTBP1) that could prove useful for normalization across a variety of RNA-based techniques.

Chip-based detection of hepatitis C virus using RNA aptamers that specifically bind to HCV core antigen

The development of reagents with high affinity and specificity to the antigens of hepatitis C virus (HCV) is important for the early stage diagnosis of its infection. Aptamers are short, single-stranded oligonucleotides with the ability to specifically recognize target molecules with high affinity. Herein, we report the selection of RNA aptamers that bind to the core antigen of HCV. High affinity aptamers were isolated from a 10(15) random library of 60 mer RNAs using the SELEX procedure. Importantly, the selected aptamers specifically bound to the core antigen, but not to another HCV antigen, NS5, in a protein chip-based assay. Using these aptamers, we developed an aptamer-based biosensor for HCV diagnosis and detected the core antigen from HCV infected patients' sera with good specificity. This novel aptamer-based antigen detection sensor could be applied to the early diagnosis of HCV infection.

SELEX experiments: new prospects, applications and data analysis in inferring regulatory pathways

Systematic Evolution of Ligands by EXponential enrichment (SELEX) is an experimental procedure that allows extraction, from an initially random pool of oligonucleotides, of the oligomers with a desired binding affinity for a given molecular target. The procedure can be used to infer the strongest binders for a given DNA or RNA binding protein, and the highest affinity binding sequences isolated through SELEX can have numerous research, diagnostic and therapeutic applications. Recently, important new modifications of the SELEX protocol have been proposed. In particular, a modification of the standard SELEX procedure allows generating a dataset from which protein-DNA interaction parameters can be determined with unprecedented accuracy. Another variant of SELEX allows investigating interactions of a protein with nucleic-acid fragments derived from the entire genome of an organism. We review here different SELEX-based methods, with particular emphasis on the experimental design and on the applications aimed at inferring protein-DNA interactions. In addition to the experimental issues, we also review relevant methods of data analysis, as well as theoretical modeling of SELEX.

The gene encoding the splicing factor SF2/ASF is a proto-oncogene

Alternative splicing modulates the expression of many oncogene and tumor-suppressor isoforms. We have tested whether some alternative splicing factors are involved in cancer. We found that the splicing factor SF2/ASF is upregulated in various human tumors, in part due to amplification of its gene, SFRS1. Moreover, slight overexpression of SF2/ASF is sufficient to transform immortal rodent fibroblasts, which form sarcomas in nude mice. We further show that SF2/ASF controls alternative splicing of the tumor suppressor BIN1 and the kinases MNK2 and S6K1. The resulting BIN1 isoforms lack tumor-suppressor activity; an isoform of MNK2 promotes MAP kinase-independent eIF4E phosphorylation; and an unusual oncogenic isoform of S6K1 recapitulates the transforming activity of SF2/ASF. Knockdown of either SF2/ASF or isoform-2 of S6K1 is sufficient to reverse transformation caused by the overexpression of SF2/ASF in vitro and in vivo. Thus, SF2/ASF can act as an oncoprotein and is a potential target for cancer therapy.

The SR family proteins B52 and dASF/SF2 modulate development of the Drosophila visual system by regulating specific RNA targets

Deciphering the role of alternative splicing in developmental processes relies on the identification of key genes whose expression is controlled by splicing regulators throughout the growth of a whole organism. Modulating the expression levels of five SR proteins in the developing eye of Drosophila melanogaster revealed that these splicing factors induce various phenotypic alterations in eye organogenesis and also affect viability. Although the SR proteins dASF/SF2 and B52 caused defects in ommatidia structure, only B52 impaired normal axonal projections of photoreceptors and neurogenesis in visual ganglia. Microarray analyses revealed that many transcripts involved in brain organogenesis have altered splicing profiles upon both loss and gain of B52 function. Conversely, a large proportion of transcripts regulated by dASF/SF2 are involved in eye development. These differential and specific effects of SR proteins indicate that they function to confer accuracy to developmental gene expression programs by facilitating the cell lineage decisions that underline the generation of tissue identities.

Alternative splicing of pre-messenger RNAs in plants in the genomic era

Primary transcripts (precursor-mRNAs) with introns can undergo alternative splicing to produce multiple transcripts from a single gene by differential use of splice sites, thereby increasing the transcriptome and proteome complexity within and between cells and tissues. Alternative splicing in plants is largely an unexplored area of gene expression, as this phenomenon used to be considered rare. However, recent genome-wide computational analyses have revealed that alternative splicing in flowering plants is far more prevalent than previously thought. Interestingly, pre-mRNAs of many spliceosomal proteins, especially serine/arginine-rich (SR) proteins, are extensively alternatively spliced. Furthermore, stresses have a dramatic effect on alternative splicing of pre-mRNAs including those that encode many spliceosomal proteins. Although the mechanisms that regulate alternative splicing in plants are largely unknown, several reports strongly suggest a key role for SR proteins in spliceosome assembly and regulated splicing. Recent studies suggest that alternative splicing in plants is an important posttranscriptional regulatory mechanism in modulating gene expression and eventually plant form and function.

Specific role of the SR protein splicing factor B52 in cell cycle control in Drosophila

E2F and retinoblastoma tumor suppressor protein pRB are important regulators of cell proliferation; however, the regulation of these proteins in vivo is not well understood. In Drosophila there are two E2F genes, an activator, de2f1, and a repressor, de2f2. The loss of de2f1 gives rise to the G(1)/S block accompanied by the repression of E2F-dependent transcription. These defects can be suppressed by mutation of de2f2. In this work, we show that the de2f1 mutant phenotype is rescued by the loss of the pre-mRNA splicing factor SR protein B52. Mutations in B52 restore S phase in clones of de2f1 mutant cells and phenocopy the loss of the de2f2 function. B52 acts upstream of de2f2 and plays a specific role in regulation of de2f2 pre-mRNA splicing. In B52-deficient cells, the level of dE2F2 protein is severely reduced and the expression of dE2F2-dependent genes is deregulated. Reexpression of the intronless copy of dE2F2 in B52-deficient cells restores the dE2F2-mediated repression. These results uncover a previously unrecognized role of the splicing factor in maintaining the G(1)/S block in vivo by specific regulation of the dE2F2 repressor function.

Experimental approaches to identify non-coding RNAs

Cellular RNAs that do not function as messenger RNAs (mRNAs), transfer RNAs (tRNAs) or ribosomal RNAs (rRNAs) comprise a diverse class of molecules that are commonly referred to as non-protein-coding RNAs (ncRNAs). These molecules have been known for quite a while, but their importance was not fully appreciated until recent genome-wide searches discovered thousands of these molecules and their genes in a variety of model organisms. Some of these screens were based on biocomputational prediction of ncRNA candidates within entire genomes of model organisms. Alternatively, direct biochemical isolation of expressed ncRNAs from cells, tissues or entire organisms has been shown to be a powerful approach to identify ncRNAs both at the level of individual molecules and at a global scale. In this review, we will survey several such wet-lab strategies, i.e. direct sequencing of ncRNAs, shotgun cloning of small-sized ncRNAs (cDNA libraries), microarray analysis and genomic SELEX to identify novel ncRNAs, and discuss the advantages and limits of these approaches.

RNomics: identification and function of small non-protein-coding RNAs in model organisms

In the recent past, our knowledge on small non-protein-coding RNAs (ncRNAs) has exponentially grown. Different approaches to identify novel ncRNAs that include computational and experimental RNomics have led to a plethora of novel ncRNAs. A picture emerges, in which ncRNAs have a variety of roles during regulation of gene expression. Thereby, many of these ncRNAs appear to function in guiding specific protein complexes to target nucleic acids. The concept of RNA guiding seems to be a widespread and very effective regulatory mechanism. In addition to guide RNAs, numerous RNAs were identified by RNomics screens, lacking known sequence and structure motifs; hence no function could be assigned to them as yet. Future challenges in the field of RNomics will include elucidation of their biological roles in the cell.

Quantitative modeling and data analysis of SELEX experiments

SELEX (systematic evolution of ligands by exponential enrichment) is an experimental procedure that allows the extraction, from an initially random pool of DNA, of those oligomers with high affinity for a given DNA-binding protein. We address what is a suitable experimental and computational procedure to infer parameters of transcription factor-DNA interaction from SELEX experiments. To answer this, we use a biophysical model of transcription factor-DNA interactions to quantitatively model SELEX. We show that a standard procedure is unsuitable for obtaining accurate interaction parameters. However, we theoretically show that a modified experiment in which chemical potential is fixed through different rounds of the experiment allows robust generation of an appropriate dataset. Based on our quantitative model, we propose a novel bioinformatic method of data analysis for such a modified experiment and apply it to extract the interaction parameters for a mammalian transcription factor CTF/NFI. From a practical point of view, our method results in a significantly improved false positive/false negative trade-off, as compared to both the standard information theory based method and a widely used empirically formulated procedure.

Identification of a splicing enhancer in MLH1 using COMPARE, a new assay for determination of relative RNA splicing efficiencies

Exonic splicing enhancers (ESEs) are sequences that facilitate recognition of splice sites and prevent exon-skipping. Because ESEs are often embedded within protein-coding sequences, alterations in them can also often be interpreted as nonsense, missense or silent mutations. To correctly interpret exonic mutations and their roles in diseases, it is important to develop strategies that identify ESE mutations. Potential ESEs can be found computationally in many exons but it has proven difficult to predict whether a given mutation will have effects on splicing based on sequence alone. Here, we describe a flexible in vitro method that can be used to functionally compare the effects of multiple sequence variants on ESE activity in a single in vitro splicing reaction. We have applied this method in parallel with conventional splicing assays to test for a splicing enhancer in exon 17 of the human MLH1 gene. Point mutations associated with hereditary non-polyposis colorectal cancer (HNPCC) have previously been found to correlate with exon-skipping in both lymphocytes and tumors from patients. We show that sequences from this exon can replace an ESE from the mouse IgM gene to support RNA splicing in HeLa nuclear extracts. ESE activity was reduced by HNPCC point mutations in codon 659, indicating that their primary effect is on splicing. Surprisingly, the strongest enhancer function mapped to a different region of the exon upstream of this codon. Together, our results indicate that HNPCC point mutations in codon 659 affect an auxillary element that augments the enhancer function to ensure exon inclusion.

Characterization of Mycobacterium tuberculosis Rv3676 (CRPMt), a cyclic AMP receptor protein-like DNA binding protein

Little is known about cyclic AMP (cAMP) function in Mycobacterium tuberculosis, despite its ability to encode 15 adenylate cyclases and 10 cNMP-binding proteins. M. tuberculosis Rv3676, which we have designated CRP(Mt), is predicted to be a cAMP-dependent transcription factor. In this study, we characterized CRP(Mt)'s interactions with DNA and cAMP, using experimental and computational approaches. We used Gibbs sampling to define a CRP(Mt) DNA motif that resembles the cAMP receptor protein (CRP) binding motif model for Escherichia coli. CRP(Mt) binding sites were identified in a total of 73 promoter regions regulating 114 genes in the M. tuberculosis genome, which are being explored as a regulon. Specific CRP(Mt) binding caused DNA bending, and the substitution of highly conserved nucleotides in the binding site resulted in a complete loss of binding to CRP(Mt). cAMP enhanced CRP(Mt)'s ability to bind DNA and caused allosteric alterations in CRP(Mt) conformation. These results provide the first direct evidence for cAMP binding to a transcription factor in M. tuberculosis, suggesting a role for cAMP signal transduction in M. tuberculosis and implicating CRP(Mt) as a cAMP-responsive global regulator.

Chironomus tentans-repressor splicing factor represses SR protein function locally on pre-mRNA exons and is displaced at correct splice sites

Chironomus tentans-repressor splicing factor (Ct-RSF) represses the activation of splicing by SR proteins in vitro. Ct-RSF colocalizes with the Ser-Arg-rich (SR) protein hrp45 in interchromatin granule clusters and coimmunoprecipitates with hrp45 in nuclear extracts. Ct-RSF and hrp45 can also interact directly in vitro. Ct-RSF and hrp45 are recruited together to transcribing genes and associate with growing pre-mRNAs. Ct-RSF and hrp45 colocalize at a large number of gene loci. Injection of anti-Ct-RSF antibodies into nuclei of living cells blocks association of both Ct-RSF and hrp45 with the growing pre-mRNA, whereas binding of U2 small nuclear ribonucleoprotein particle (snRNP) to the pre-mRNA is unaffected. On the intron-rich Balbiani ring (BR) 3 pre-mRNA, hrp45 as well as U1 and U2 snRNPs bind extensively, whereas relatively little Ct-RSF is present. In contrast, the BR1 and BR2 pre-mRNAs, dominated by exon sequences, bind relatively much Ct-RSF compared with hrp45 and snRNPs. Our data suggest that Ct-RSF represses SR protein function at exons and that the assembly of spliceosomes at authentic splice sites displaces Ct-RSF locally.

Global analysis of positive and negative pre-mRNA splicing regulators in Drosophila

To gain insight into splicing regulation, we developed a microarray to assay all annotated alternative splicing events in Drosophila melanogaster and identified the alternative splice events controlled by four splicing regulators: dASF/SF2, B52/SRp55, hrp48, and PSI. The number of events controlled by each of these factors was found to be highly variable: dASF/SF2 strongly affects >300 splicing events, whereas PSI strongly affects only 43 events. Pairwise analysis also revealed many instances of splice site usage affected by multiple factors and provides the framework to understand the network controlling the alternatively spliced mRNA isoforms that compose the Drosophila transcriptome.

Nova regulates brain-specific splicing to shape the synapse

Alternative RNA splicing greatly increases proteome diversity and may thereby contribute to tissue-specific functions. We carried out genome-wide quantitative analysis of alternative splicing using a custom Affymetrix microarray to assess the role of the neuronal splicing factor Nova in the brain. We used a stringent algorithm to identify 591 exons that were differentially spliced in the brain relative to immune tissues, and 6.6% of these showed major splicing defects in the neocortex of Nova2-/- mice. We tested 49 exons with the largest predicted Nova-dependent splicing changes and validated all 49 by RT-PCR. We analyzed the encoded proteins and found that all those with defined brain functions acted in the synapse (34 of 40, including neurotransmitter receptors, cation channels, adhesion and scaffold proteins) or in axon guidance (8 of 40). Moreover, of the 35 proteins with known interaction partners, 74% (26) interact with each other. Validating a large set of Nova RNA targets has led us to identify a multi-tiered network in which Nova regulates the exon content of RNAs encoding proteins that interact in the synapse.

Understanding alternative splicing: towards a cellular code

In violation of the 'one gene, one polypeptide' rule, alternative splicing allows individual genes to produce multiple protein isoforms - thereby playing a central part in generating complex proteomes. Alternative splicing also has a largely hidden function in quantitative gene control, by targeting RNAs for nonsense-mediated decay. Traditional gene-by-gene investigations of alternative splicing mechanisms are now being complemented by global approaches. These promise to reveal details of the nature and operation of cellular codes that are constituted by combinations of regulatory elements in pre-mRNA substrates and by cellular complements of splicing regulators, which together determine regulated splicing pathways.

ASF/SF2-regulated CaMKIIdelta alternative splicing temporally reprograms excitation-contraction coupling in cardiac muscle

The transition from juvenile to adult life is accompanied by programmed remodeling in many tissues and organs, which is key for organisms to adapt to the demand of the environment. Here we report a novel regulated alternative splicing program that is crucial for postnatnal heart remodeling in the mouse. We identify the essential splicing factor ASF/SF2 as a key component of the program, regulating a restricted set of tissue-specific alternative splicing events during heart remodeling. Cardiomyocytes deficient in ASF/SF2 display an unexpected hypercontraction phenotype due to a defect in postnatal splicing switch of the Ca(2+)/calmodulin-dependent kinase IIdelta (CaMKIIdelta) transcript. This failure results in mistargeting of the kinase to sarcolemmal membranes, causing severe excitation-contraction coupling defects. Our results validate ASF/SF2 as a fundamental splicing regulator in the reprogramming pathway and reveal the central contribution of ASF/SF2-regulated CaMKIIdelta alternative splicing to functional remodeling in developing heart.

The splicing regulatory element, UGCAUG, is phylogenetically and spatially conserved in introns that flank tissue-specific alternative exons

Previous studies have identified UGCAUG as an intron splicing enhancer that is frequently located adjacent to tissue-specific alternative exons in the human genome. Here, we show that UGCAUG is phylogenetically and spatially conserved in introns that flank brain-enriched alternative exons from fish to man. Analysis of sequence from the mouse, rat, dog, chicken and pufferfish genomes revealed a strongly statistically significant association of UGCAUG with the proximal intron region downstream of brain-enriched alternative exons. The number, position and sequence context of intronic UGCAUG elements were highly conserved among mammals and in chicken, but more divergent in fish. Control datasets, including constitutive exons and non-tissue-specific alternative exons, exhibited a much lower incidence of closely linked UGCAUG elements. We propose that the high sequence specificity of the UGCAUG element, and its unique association with tissue-specific alternative exons, mark it as a critical component of splicing switch mechanism(s) designed to activate a limited repertoire of splicing events in cell type-specific patterns. We further speculate that highly conserved UGCAUG-binding protein(s) related to the recently described Fox-1 splicing factor play a critical role in mediating this specificity.

Selection of genomic sequences that bind tightly to Ff gene 5 protein: primer-free genomic SELEX

Single-stranded DNA or RNA libraries used in SELEX experiments usually include primer-annealing sequences for PCR amplification. In genomic SELEX, these fixed sequences may form base pairs with the central genomic fragments and interfere with the binding of target molecules to the genomic sequences. In this study, a method has been developed to circumvent these artificial effects. Primer-annealing sequences are removed from the genomic library before selection with the target protein and are then regenerated to allow amplification of the selected genomic fragments. A key step in the regeneration of primer-annealing sequences is to employ thermal cycles of hybridization-extension, using the sequences from unselected pools as templates. The genomic library was derived from the bacteriophage fd, and the gene 5 protein (g5p) from the phage was used as a target protein. After four rounds of primer-free genomic SELEX, most cloned sequences overlapped at a segment within gene 6 of the viral genome. This sequence segment was pyrimidine-rich and contained no stable secondary structures. Compared with a neighboring genomic fragment, a representative sequence from the family of selected sequences had about 23-fold higher g5p-binding affinity. Results from primer-free genomic SELEX were compared with the results from two other genomic SELEX protocols.