A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

Deciphering 3'ss selection in the yeast genome reveals an RNA thermosensor that mediates alternative splicing

Poor understanding of the spliceosomal mechanisms to select intronic 3' ends (3'ss) is a major obstacle to deciphering eukaryotic genomes. Here, we discern the rules for global 3'ss selection in yeast. We show that, in contrast to the uniformity of yeast splicing, the spliceosome uses all available 3'ss within a distance window from the intronic branch site (BS), and that in 70% of all possible 3'ss this is likely to be mediated by pre-mRNA structures. Our results reveal that one of these RNA folds acts as an RNA thermosensor, modulating alternative splicing in response to heat shock by controlling alternate 3'ss availability. Thus, our data point to a deeper role for the pre-mRNA in the control of its own fate, and to a simple mechanism for some alternative splicing.

Exploring the effect of aminoglycoside guanidinylation on ligands for Tau exon 10 splicing regulatory element RNA

We describe the effect of guanidinylation of the aminoglycoside moiety on acridine-neamine-containing ligands for the stem-loop structure located at the exon 10-5′-intron junction of Tau pre-mRNA, an important regulatory element of tau gene alternative splicing. On the basis of dynamic combinatorial chemistry experiments, ligands that combine guanidinoneamine and two different acridines were synthesized and their RNA-binding properties were compared with those of their amino precursors. Fluorescence titration experiments and UV-monitored melting curves revealed that guanidinylation has a positive effect both on the binding affinity and specificity of the ligands for the stemloop RNA, as well as on the stabilization of all RNA sequences evaluated, particularly some mutated sequences associated with the development of FTDP-17 tauopathy. However, this correlation between binding affinity and stabilization due to guanidinylation was only found in ligands containing a longer spacer between the acridine and guanidinoneamine moieties, since a shorter spacer produced the opposite effect (e.g. lower binding affinity and lower stabilization). Furthermore, spectroscopic studies suggest that ligand binding does not significantly change the overall RNA structure upon binding (circular dichroism) and that the acridine moiety might intercalate near the bulged region of the stem->loop structure (UV-Vis and NMR spectroscopy).

Improving gene annotation using peptide mass spectrometry

Annotation of protein-coding genes is a key goal of genome sequencing projects. In spite of tremendous recent advances in computational gene finding, comprehensive annotation remains a challenge. Peptide mass spectrometry is a powerful tool for researching the dynamic proteome and suggests an attractive approach to discover and validate protein-coding genes. We present algorithms to construct and efficiently search spectra against a genomic database, with no prior knowledge of encoded proteins. By searching a corpus of 18.5 million tandem mass spectra (MS/MS) from human proteomic samples, we validate 39,000 exons and 11,000 introns at the level of translation. We present translation-level evidence for novel or extended exons in 16 genes, confirm translation of 224 hypothetical proteins, and discover or confirm over 40 alternative splicing events. Polymorphisms are efficiently encoded in our database, allowing us to observe variant alleles for 308 coding SNPs. Finally, we demonstrate the use of mass spectrometry to improve automated gene prediction, adding 800 correct exons to our predictions using a simple rescoring strategy. Our results demonstrate that proteomic profiling should play a role in any genome sequencing project.

EGASP: the human ENCODE Genome Annotation Assessment Project

Background: We present the results of EGASP, a community experiment to assess the state-ofthe-art in genome annotation within the ENCODE regions, which span 1% of the human genomesequence. The experiment had two major goals: the assessment of the accuracy of computationalmethods to predict protein coding genes; and the overall assessment of the completeness of thecurrent human genome annotations as represented in the ENCODE regions. For thecomputational prediction assessment, eighteen groups contributed gene predictions. Weevaluated these submissions against each other based on a ‘reference set’ of annotationsgenerated as part of the GENCODE project. These annotations were not available to theprediction groups prior to the submission deadline, so that their predictions were blind and anexternal advisory committee could perform a fair assessment.Results: The best methods had at least one gene transcript correctly predicted for close to 70%of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into accountalternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotidelevel, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programsrelying on mRNA and protein sequences were the most accurate in reproducing the manuallycurated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could beverified.Conclusions: This is the first such experiment in human DNA, and we have followed thestandards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe theresults presented here contribute to the value of ongoing large-scale annotation projects and shouldguide further experimental methods when being scaled up to the entire human genome sequence.

PFKFB3 (6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3)

The human PFKFB3 is composed of 19 exons spanning genomic region about 90,6 Kb (GenBank). Alternative splicing variants have been reported. The main variants corresponding to mRNAs of 4453 bp and 4224 bp for the variant 1 u-PFK2 (NM_004566.3) and variant 2 i-PFK2 (NM_001145443.1), respectively...

Do protein–protein interaction databases identify moonlighting proteins?

One of the most striking results of the human (and mammalian) genomes is the low number of protein-coding genes. To-date, the main molecular mechanism to increase the number of different protein isoforms and functions is alternative splicing. However, a less-known way to increase the number of protein functions is the existence of multifunctional, multitask, or ‘‘moonlighting’’, proteins. By and large, moonlighting proteins are experimentally disclosed by serendipity. Proteomics is becoming one of the very active areas of biomedical research, which permits researchers to identify previously unseen connections among proteins and pathways. In principle, protein–protein interaction (PPI) databases should contain information on moonlighting proteins and could provide suggestions to further analysis in order to prove the multifunctionality. As far as we know, nobody has verified whether PPI databases actually disclose moonlighting proteins. In the present work we check whether well-established moonlighting proteins present in PPI databases connect with their known partners and, therefore, a careful inspection of these databases could help to suggest their different functions. The results of our research suggest that PPI databases could be a valuable tool to suggest multifunctionality.

PFKFB3 (6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3)

The human PFKFB3 is composed of 19 exons spanning genomic region about 90,6 Kb (GenBank). Alternative splicing variants have been reported. The main variants corresponding to mRNAs of 4453 bp and 4224 bp for the variant 1 u-PFK2 (NM_004566.3) and variant 2 i-PFK2 (NM_001145443.1), respectively...

Metal dealing at the origin of the Chordata phylum: The metallothionein system and metal overload response in Amphioxus.

Non-vertebrate chordates, specifically amphioxus, are considered of the utmost interest for gaining insight into the evolutionary trends, i.e. differentiation and specialization, of gene/protein systems. In this work, MTs (metallothioneins), the most important metal binding proteins, are characterized for the first time in the cephalochordate subphylum at both gene and protein level, together with the main features defining the amphioxus response to cadmium and copper overload. Two MT genes (BfMT1 and BfMT2) have been identified in a contiguous region of the genome, as well as several ARE (antioxidant response element) and MRE (metal response element) located upstream the transcribed region. Their corresponding cDNAs exhibit identical sequence in the two lancelet species (B. floridae and B. lanceolatum), BfMT2 cDNA resulting from an alternative splicing event. BfMT1 is a polyvalent metal binding peptide that coordinates any of the studied metal ions (Zn, Cd or Cu) rendering complexes stable enough to last in physiological environments, which is fully concordant with the constitutive expression of its gene, and therefore, with a metal homeostasis housekeeping role. On the contrary, BfMT2 exhibits a clear ability to coordinate Cd(II) ions, while it is absolutely unable to fold into stable Cu (I) complexes, even as mixed species. This identifies it as an essential detoxification agent, which is consequently only induced in emergency situations. The cephalochordate MTs are not directly related to vertebrate MTs, neither by gene structure, protein similarity nor metal-binding behavior of the encoded peptides. The closest relative is the echinoderm MT, which confirm proposed phylogenetic relationships between these two groups. The current findings support the existence in most organisms of two types of MTs as for their metal binding preferences, devoted to different biological functions: multivalent MTs for housekeeping roles, and specialized MTs that evolve either as Cd-thioneins or Cu-thioneins, according to the ecophysiological needs of each kind of organisms.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

An alternative way to model merit good arguments

Besley (1988) uses a scaling approach to model merit good arguments in commodity tax policy. In this paper, I question this approach on the grounds that it produces 'wrong' recommendations--taxation (subsidisation) of merit (demerit) goods--whenever the demand for the (de)merit good is inelastic. I propose an alternative approach that does not suffer from this deficiency, and derive the ensuing first and second best tax rules, as well as the marginal cost expressions to perform tax reform analysis.

Economic impacts of alternative water policy scenarios in the Spanish production system: an input-output analysis

The objective of this paper is to analyse the economic impacts of alternative water policies implemented in the Spanish production system. The methodology uses two versions of the input-output price model: a competitive formulation and a mark-up formulation. The input-output framework evaluates the impact of water policy measures on production prices, consumption prices, intermediate water demand and private welfare. Our results show that a tax on the water used by sectors considerably reduces the intermediate water demand, and increases the production and consumption prices. On the other hand, according to Jevons' paradox, an improvement in technical efficiency, which leads to a reduction in the water requirements of all sectors and an increase in water production, increases the amount of water consumed. The combination of a tax on water and improved technical efficiency takes the pressure off prices and significantly reduces intermediate water demand. JEL Classification: C67 ; D57 ; Q25. Keywords: Production prices; Consumption prices; Water uses; Water policy; Water taxation.

Economic impact of alternative policy measures implemented on the energy activities of the Catalan production system: an input-output analysis

The aim of the paper is to analyse the economic impact of alternative policies implemented on the energy activities of the Catalan production system. Specifically, we analyse the effects of a tax on intermediate energy uses, a reduction in the final production of energy, and a reduction in intermediate energy uses. The methodology involves two versions of the input-output price model: a competitive price formulation and a mark-up price formulation. The input-output price framework will make it possible to evaluate how the alternative measures modify production prices, consumption prices, private welfare, and intermediate energy uses. The empirical application is for the Catalan economy and uses economic data for the year 2001.

Comparació dels sistemes "ribosome skipping" i "splicing acceptor" per a l'expressió de transgenes en adenovirus oncolítics

Report for the scientific sojourn carried out at the University of St. Andrews, United Kingdom, from November 2007 until January 2008. Therapeutic transgene expression is a valuable strategy to counteract the limitations associated with oncolytic adenoviruses. Late phase expression is desirable to avoid early cell death for proper virus production. In this 3 months-collaboration, we have constructed a late expression system based on ribosome skipping downstream fiber protein and compared it with a splicing-based method of late gene expression. Despite expressing high amounts of the transgene when utilizing the ribosome skipping-system, flow cytomety assays indicate a delayed transgene-expression kinetics compared with the splicing-based one. Furthermore, when using the ribosome skipping system not only fiber protein expression is more altered but also viral production. These results suggest splicing-based expression strategy as a more suitable system for expression of transgenes late in the viral life cycle of an oncolytic adenovirus.