Alternative splicing and protein function

Integrated top-down and bottom-up proteomics mass spectrometry for the characterization of endogenous ribosomal protein heterogeneity

Ribosomes are abundant, large RNA-protein complexes that are the sites of all protein synthesis in cells. Defects in ribosomal proteins (RPs), including proteoforms arising from genetic variations, alternative splicing of RNA transcripts, post-translational modifications and alterations of protein expression level, have been linked to a diverse range of diseases, including cancer and aging. Comprehensive characterization of ribosomal proteoforms is challenging but important for the discovery of potential disease biomarkers or protein targets. In the present work, using E. coli 70S RPs as an example, we first developed a top-down proteomics approach on a Waters Synapt G2 Si mass spectrometry (MS) system, and then applied it to the HeLa 80S ribosome. The results were complemented by a bottom-up approach. In total, 50 out of 55 RPs were identified using the top-down approach. Among these, more than 30 RPs were found to have their N-terminal methionine removed. Additional modifications such as methylation, acetylation, and hydroxylation were also observed, and the modification sites were identified by bottom-up MS. In a HeLa 80S ribosomal sample, we identified 98 ribosomal proteoforms, among which multiple truncated 80S ribosomal proteoforms were observed, the type of information which is often overlooked by bottom-up experiments. Although their relevance to diseases is not yet known, the integration of top-down and bottom-up proteomics approaches paves the way for the discovery of proteoform-specific disease biomarkers or targets.

ApiAP2 Gene-Network Regulates Gametocytogenesis in Plasmodium Parasites

Malaria is a mosquito-borne infectious disease, caused by unicellular Apicomplexan protozoa of the genus Plasmodium. The sexual stage of Plasmodium is one of the most fascinating aspects of the Plasmodium life cycle, yet relatively less explored until now. The production of sexually fit gametocytes through gametocytogenesis is essential to the transmission of the Plasmodium parasite into an anopheline mosquito vector. Understanding how gametocytogenesis is regulated promotes the identification of novel drug targets and also the development of transmission-blocking vaccines that would help reduce the disease burden in endemic areas. Transcriptional regulation in Plasmodium parasites is primarily controlled by a family of twenty-seven Apicomplexan Apetela 2 (ApiAP2) genes which act in a cascade to enable the parasite to progress through its asexual replication as well as gametocytogenesis. Here, we review the latest progress made on members of the ApiAP2 family characterized as key players of the transcriptional machinery of gametocytes. Further, we will highlight the transcriptional regulation network of ApiAP2 genes at each stage of gametocytogenesis.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Alternative splicing dynamics and evolutionary divergence during embryogenesis in wheat species

Among polyploid species with complex genomic architecture, variations in the regulation of alternative splicing (AS) provide opportunities for transcriptional and proteomic plasticity and the potential for generating trait diversities. However, the evolution of AS and its influence on grain development in diploid grass and valuable polyploid wheat crops are poorly understood. To address this knowledge gap, we developed a pipeline for the analysis of alternatively spliced transcript isoforms, which takes the high sequence similarity among polyploid wheat subgenomes into account. Through analysis of synteny and detection of collinearity of homoeologous subgenomes, conserved and specific AS events across five wheat and grass species were identified. A global analysis of the regulation of AS in diploid grass and polyploid wheat grains revealed diversity in AS events not only between the endosperm, pericarp and embryo overdevelopment, but also between subgenomes. Analysis of AS in homoeologous triads of polyploid wheats revealed evolutionary divergence between gene-level and transcript-level regulation of embryogenesis. Evolutionary age analysis indicated that the generation of novel transcript isoforms has occurred in young genes at a more rapid rate than in ancient genes. These findings, together with the development of comprehensive AS resources for wheat and grass species, advance understanding of the evolution of regulatory features of AS during embryogenesis and grain development in wheat.

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

Development and validation of an individual alternative splicing prognostic signature in gastric cancer

Gastric cancer (GC) is a heterogeneous disease with different clinical manifestations and prognoses. Alternative splicing (AS) is a determinant of gene expression and contributes to protein diversity from a rather limited gene transcript in metazoans. AS events are associated with different aspects of cancer biology, including cell proliferation, apoptosis, invasion, etc. Here, we present a comprehensive analysis of the prognostic AS profile in GC. GC-specific AS (GCAS) events were analyzed, and overall survival-associated GCAS (OS-GCAS) events were verified among the genome-wide AS events identified in The Cancer Genome Atlas (TCGA) database. In total, 1,287 GCAS events of 837 genes and 173 OS-GCAS events of 130 genes were identified. The parental genes of OS-GCAS events were significantly enriched in the development of GC. Protein-protein interaction (PPI) and OS-GCAS-associated splicing factor (SF) interaction networks were constructed. Multivariate Cox regression analysis with least absolute shrinkage and selection operator (LASSO) penalty was performed to establish a prognostic risk formula, representing 23 OS-GCAS events. The low-risk group had better OS than the high-risk group and lower immune and stromal scores. Cox proportional hazard regression was applied to generate an AS-clinical integrated prognostic model with a considerable area under the curve (AUC) value in both the training and validation datasets. Our study provides a profile of OS-GCAS events and an AS-clinical nomogram to predict the prognosis of GC.

Isoform function prediction based on bi-random walks on a heterogeneous network

MOTIVATION

Alternative splicing contributes to the functional diversity of protein species and the proteoforms translated from alternatively spliced isoforms of a gene actually execute the biological functions. Computationally predicting the functions of genes has been studied for decades. However, how to distinguish the functional annotations of isoforms, whose annotations are essential for understanding developmental abnormalities and cancers, is rarely explored. The main bottleneck is that functional annotations of isoforms are generally unavailable and functional genomic databases universally store the functional annotations at the gene level.

RESULTS

We propose IsoFun to accomplish Isoform Function prediction based on bi-random walks on a heterogeneous network. IsoFun firstly constructs an isoform functional association network based on the expression profiles of isoforms derived from multiple RNA-seq datasets. Next, IsoFun uses the available Gene Ontology annotations of genes, gene-gene interactions and the relations between genes and isoforms to construct a heterogeneous network. After this, IsoFun performs a tailored bi-random walk on the heterogeneous network to predict the association between GO terms and isoforms, thus accomplishing the prediction of GO annotations of isoforms. Experimental results show that IsoFun significantly outperforms the state-of-the-art algorithms and improves the area under the receiver-operating curve (AUROC) and the area under the precision-recall curve (AUPRC) by 17% and 44% at the gene-level, respectively. We further validated the performance of IsoFun on the genes ADAM15 and BCL2L1. IsoFun accurately differentiates the functions of respective isoforms of these two genes.

AVAILABILITY AND IMPLEMENTATION

The code of IsoFun is available at http://mlda.swu.edu.cn/codes.php? name=IsoFun.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Systematic profiling of alternative splicing in Helicobacter pylori-negative gastric cancer and their clinical significance

Background

Alternative splicing (AS) may cause structural and functional variations in the protein to promote the proliferation of tumor cells. However, there is no comprehensive analysis of the clinical significance of AS in Helicobacter pylori-negative gastric cancer (HP⁻ GC).

Methods

The clinical, gene expression profile data and AS events of 138 HP⁻ GC patients were obtained from the database named the cancer genome atlas. Differently expressed AS (DEAS) events were determined by a comparison of the PSI values between HP⁻ GC samples and adjacent normal samples. Unsupervised clustering analysis, proportional regression and Kaplan–Meier analysis were used to explore the association between clinical data and immune features and to establish two nomograms about the prognosis of HP⁻ GC. Finally, splicing networks were constructed using Cytoscape.

Results

A total of 48141 AS events and 1041 DEAS events were found in HP⁻ GC. Various functions and pathways of DEAS events parent genes were enriched, such as cell-substrate junction, cell leading edge, focal adhension, and AMPK signaling. Seven overall survival (OS)-related and seven disease-free survival (DFS)-related AS events were used to construct the prognostic signatures. Based on the independent prognostic factors, two nomograms were established and showed excellent performance. Then, splicing regulatory networks among the correlations suggested that splicing factors were significantly associated with prognostic DEASs. Finally, the unsupervised clustering analysis revealed that DEAS-based clusters were associated with clinical characteristics, tumor microenvironment, tumor mutation burden, and immune features.

Conclusion

Seven OS-related and seven DFS-related AS events have been found to be correlated with the prognosis of HP⁻ GC and can be used as prognostic factors to establish an effective nomogram.

Alternative Splicing and Protein Diversity: Plants Versus Animals

Plants, unlike animals, exhibit a very high degree of plasticity in their growth and development and employ diverse strategies to cope with the variations during diurnal cycles and stressful conditions. Plants and animals, despite their remarkable morphological and physiological differences, share many basic cellular processes and regulatory mechanisms. Alternative splicing (AS) is one such gene regulatory mechanism that modulates gene expression in multiple ways. It is now well established that AS is prevalent in all multicellular eukaryotes including plants and humans. Emerging evidence indicates that in plants, as in animals, transcription and splicing are coupled. Here, we reviewed recent evidence in support of co-transcriptional splicing in plants and highlighted similarities and differences between plants and humans. An unsettled question in the field of AS is the extent to which splice isoforms contribute to protein diversity. To take a critical look at this question, we presented a comprehensive summary of the current status of research in this area in both plants and humans, discussed limitations with the currently used approaches and suggested improvements to current methods and alternative approaches. We end with a discussion on the potential role of epigenetic modifications and chromatin state in splicing memory in plants primed with stresses.

Alternative splicing and the evolution of phenotypic novelty

Alternative splicing, a mechanism of post-transcriptional RNA processing whereby a single gene can encode multiple distinct transcripts, has been proposed to underlie morphological innovations in multicellular organisms. Genes with developmental functions are enriched for alternative splicing events, suggestive of a contribution of alternative splicing to developmental programmes. The role of alternative splicing as a source of transcript diversification has previously been compared to that of gene duplication, with the relationship between the two extensively explored. Alternative splicing is reduced following gene duplication with the retention of duplicate copies higher for genes which were alternatively spliced prior to duplication. Furthermore, and unlike the case for overall gene number, the proportion of alternatively spliced genes has also increased in line with the evolutionary diversification of cell types, suggesting alternative splicing may contribute to the complexity of developmental programmes. Together these observations suggest a prominent role for alternative splicing as a source of functional innovation. However, it is unknown whether the proliferation of alternative splicing events indeed reflects a functional expansion of the transcriptome or instead results from weaker selection acting on larger species, which tend to have a higher number of cell types and lower population sizes.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

A proteogenomic approach to understand splice isoform functions through sequence and expression-based computational modeling

The products of multi-exon genes are a mixture of alternatively spliced isoforms, from which the translated proteins can have similar, different or even opposing functions. It is therefore essential to differentiate and annotate functions for individual isoforms. Computational approaches provide an efficient complement to expensive and time-consuming experimental studies. The input data of these methods range from DNA sequence, to RNA selection pressure, to expressed sequence tags, to full-length complementary DNA, to exon array, to RNA-seq expression, to proteomic data. Notably, RNA-seq technology generates quantitative profiling of transcript expression at the genome scale, with an unprecedented amount of expression data available for developing isoform function prediction methods. Integrative analysis of these data at different molecular levels enables a proteogenomic approach to systematically interrogate isoform functions. Here, we briefly review the state-of-the-art methods according to their input data sources, discuss their advantages and limitations and point out potential ways to improve prediction accuracies.

Coupling and coordination in gene expression processes with pre-mRNA splicing

RNA processing is a tightly regulated and highly complex pathway which includes transcription, splicing, editing, transportation, translation and degradation. It has been well-documented that splicing of RNA polymerase II medicated nascent transcripts occurs co-transcriptionally and is functionally coupled to other RNA processing. Recently, increasing experimental evidence indicated that pre-mRNA splicing influences RNA degradation and vice versa. In this review, we summarized the recent findings demonstrating the coupling of these two processes. In addition, we highlighted the importance of splicing in the production of intronic miRNA and circular RNAs, and hence the discovery of the novel mechanisms in the regulation of gene expression.

Detecting splicing variants in idiopathic pulmonary fibrosis from non-differentially expressed genes

Idiopathic pulmonary fibrosis (IPF) is an interstitial lung disease of unknown cause that lacks a proven therapy for altering its high mortality rate. Microarrays have been employed to investigate the pathogenesis of IPF, but are presented mostly at the gene-expression level due to technologic limitations. In as much as, alternative RNA splicing isoforms are increasingly identified as potential regulators of human diseases, including IPF, we propose a new approach with the capacity to detect splicing variants using RNA-seq data. We conducted a joint analysis of differential expression and differential splicing on annotated human genes and isoforms, and identified 122 non-differentially expressed genes with a high degree of "switch" between major and minor isoforms. Three cases with variant mechanisms for alternative splicing were validated using qRT-PCR, among the group of genes in which expression was not significantly changed at the gene level. We also identified 35 novel transcripts that were unique to the fibrotic lungs using exon-exon junction evidence, and selected a representative for qRT-PCR validation. The results of our study are likely to provide new insight into the pathogenesis of pulmonary fibrosis and may eventuate in new treatment targets.

Evolution of a potential hormone antagonist following gene splicing during primate evolution

Alternative splicing of genes generates novel mRNAs, leading to the evolution of new functional proteins. Cholecystokinin (CCK) induces the release of pancreatic enzymes and the contraction of the gallbladder to promote the digestion of fat and proteins. CCK activates two G-protein-coupled receptors, CCKA and CCKB. Here, we showed that a CCKsv (splicing variant), originated de novo during Catarrhini evolution by including a portion of intronic sequence of the CCK gene, encodes novel C-terminal peptide sequence followed by a new poly-adenylation signal. CCKsv is expressed in many human tissues and likely a secreted peptide retaining the original signal peptide and the N-terminal proteolytic processing signal, together with novel C-terminal sequences. Although CCKsv cannot activate CCK receptors, it partially inhibits the CRE- or SRF-driven reporter activities stimulated by wide type CCK-8 mediated by both CCK receptors. Co-treatment with CCKsv also partially antagonizes Ewing tumor cell growth stimulated by CCK-8. Our study provides an example of new peptide hormone antagonist evolution in primates.

Identification of a novel isoform of DHRS4 protein with a nuclear localization signal

The DHRS4 gene encodes an NADP(H)-dependent retinol dehydrogenase/reductase (NRDR) and plays an important role in regulating the synthesis of retinoic acid. In the present study, we identified a novel splice RNA variant, designated NRDRA2, of the human DHRS4 gene by RT-PCR, 3' RACE, and 5' RACE. NRDRA2 mRNA lacked exons 4 and 6, and had a shift in the reading frame when compared to DHRS4 mRNA, resulting in loss of the peroxisomal targeting signal of NRDR and gain of a nuclear localization signal in the predicted NRDRA2 protein. Endogenous NRDRA2 protein was identified in the human cervical carcinoma cell line HeLa by immunoprecipitation and mass spectrometric assay. A green fluorescent protein reporter assay showed that NRDRA2 protein mainly localized to the nuclei, confirming the sequence at its C-terminus as a legitimate nuclear localization signal sequence. This study identifies the alternative transcript variant NRDRA2 encoding a subcellular nuclear localized NRDRA2 protein.

Projection of gene-protein networks to the functional space of the proteome and its application to analysis of organism complexity

We consider the problem of biological complexity via a projection of protein-coding genes of complex organisms onto the functional space of the proteome. The latter can be defined as a set of all functions committed by proteins of an organism. Alternative splicing (AS) allows an organism to generate diverse mature RNA transcripts from a single mRNA strand and thus it could be one of the key mechanisms of increasing of functional complexity of the organism's proteome and a driving force of biological evolution. Thus, the projection of transcription units (TU) and alternative splice-variant (SV) forms onto proteome functional space could generate new types of relational networks (e.g. SV-protein function networks, SFN) and lead to discoveries of novel evolutionarily conservative functional modules. Such types of networks might provide new reliable characteristics of organism complexity and a better understanding of the evolutionary integration and plasticity of interconnection of genome-transcriptome-proteome functions.

Stochastic noise in splicing machinery

The number of known alternative human isoforms has been increasing steadily with the amount of available transcription data. To date, over 100 000 isoforms have been detected in EST libraries, and at least 75% of human genes have at least one alternative isoform. In this paper, we propose that most alternative splicing events are the result of noise in the splicing process. We show that the number of isoforms and their abundance can be predicted by a simple stochastic noise model that takes into account two factors: the number of introns in a gene and the expression level of a gene. The results strongly support the hypothesis that most alternative splicing is a consequence of stochastic noise in the splicing machinery, and has no functional significance. The results are also consistent with error rates tuned to ensure that an adequate level of functional product is produced and to reduce the toxic effect of accumulation of misfolding proteins. Based on simulation of sampling of virtual cDNA libraries, we estimate that error rates range from 1 to 10% depending on the number of introns and the expression level of a gene.

Splice-mediated Variants of Proteins (SpliVaP) - data and characterization of changes in signatures among protein isoforms due to alternative splicing

BACKGROUND

It is often the case that mammalian genes are alternatively spliced; the resulting alternate transcripts often encode protein isoforms that differ in amino acid sequences. Changes among the protein isoforms can alter the cellular properties of proteins. The effect can range from a subtle modulation to a complete loss of function.

RESULTS

(i) We examined human splice-mediated protein isoforms (as extracted from a manually curated data set, and from a computationally predicted data set) for differences in the annotation for protein signatures (Pfam domains and PRINTS fingerprints) and we characterized the differences & their effects on protein functionalities. An important question addressed relates to the extent of protein isoforms that may lack any known function in the cell. (ii) We present a database that reports differences in protein signatures among human splice-mediated protein isoform sequences.

CONCLUSION

(i) Characterization: The work points to distinct sets of alternatively spliced genes with varying degrees of annotation for the splice-mediated protein isoforms. Protein molecular functions seen to be often affected are those that relate to: binding, catalytic, transcription regulation, structural molecule, transporter, motor, and antioxidant; and the processes that are often affected are nucleic acid binding, signal transduction, and protein-protein interactions. Signatures are often included/excluded and truncated in length among protein isoforms; truncation is seen as the predominant type of change. Analysis points to the following novel aspects: (a) Analysis using data from the manually curated Vega indicates that one in 8.9 genes can lead to a protein isoform of no "known" function; and one in 18 expressed protein isoforms can be such an "orphan" isoform; the corresponding numbers as seen with computationally predicted ASD data set are: one in 4.9 genes and one in 9.8 isoforms. (b) When swapping of signatures occurs, it is often between those of same functional classifications. (c) Pfam domains can occur in varying lengths, and PRINTS fingerprints can occur with varying number of constituent motifs among isoforms - since such a variation is seen in large number of genes, it could be a general mechanism to modulate protein function. (ii)

DATA

The reported resource (at http://www.bioinformatica.crs4.org/tools/dbs/splivap/) provides the community ability to access data on splice-mediated protein isoforms (with value-added annotation such as association with diseases) through changes in protein signatures.

Alternative translation start sites and hidden coding potential of eukaryotic mRNAs

It is widely suggested that a eukaryotic mRNA typically contains one translation start site and encodes a single functional protein product. However, according to current points of view on translation initiation mechanisms, eukaryotic ribosomes can recognize several alternative translation start sites and the number of experimentally verified examples of alternative translation is growing rapidly. Also, the frequent occurrence of alternative translation events and their functional significance are supported by the results of computational evaluations. The functional role of alternative translation and its contribution to eukaryotic proteome complexity are discussed.

Positive selection in alternatively spliced exons of human genes

Alternative splicing is a well-recognized mechanism of accelerated genome evolution. We have studied single-nucleotide polymorphisms and human-chimpanzee divergence in the exons of 6672 alternatively spliced human genes, with the aim of understanding the forces driving the evolution of alternatively spliced sequences. Here, we show that alternatively spliced exons and exon fragments (alternative exons) from minor isoforms experience lower selective pressure at the amino acid level, accompanied by selection against synonymous sequence variation. The results of the McDonald-Kreitman test suggest that alternatively spliced exons, unlike exons constitutively included in the mRNA, are also subject to positive selection, with up to 27% of amino acids fixed by positive selection.

Sticking together? Falling apart? Exploring the dynamics of the interactome

Advances in techniques for the study of protein-protein interactions have dramatically improved our understanding of the interactome. However, we know little about the dynamics of this complex system. To better understand the dynamics of the interactome, it is important to consider what happens when single proteins are perturbed. Changes in protein abundance and post-translational modification can function as switches in the interactome, affecting protein-complex assembly and function. Changes in protein sequence or a dramatic increase in abundance might cause a promiscuous gain of interactions. These effects are not identical for all proteins and will differ depending on the number and type of interaction partners that a protein has.

Overlapping alternative donor splice sites in the human genome

Over 50% of donor splice sites in the human genome have a potential alternative donor site at a distance of three to six nucleotides. Conservation of these potential sites is determined by the consensus requirements and by its exonic or intronic location. Several hundred pairs of overlapping sites are confirmed to be alternatively spliced as both sites in a pair are supported by a protein, by a full-length mRNA, or by expressed sequence tags (ESTs) from at least two independent clone libraries. Overlapping sites may clash with consensus requirements. Pairs with a site shift of four nucleotides are the most abundant, despite the frameshift in the protein-coding region that they introduce. The site usage in pairs is usually uneven, and the major site is more frequently conserved in other mammalian genomes. Overlapping alternative donor sites and acceptor sites may have different functional roles: alternative splicing of overlapping acceptor sites leads mainly to microvariations in protein sequences; whereas alternative donor sites often lead to frameshifts and thus either yield major differences in the protein sequence and structure, or generate nonsense-mediated decay-inducing mRNA isoforms likely involved in regulated unproductive splicing pathways.

A procedure for identifying homologous alternative splicing events

Background

The study of the functional role of alternative splice isoforms of a gene is a very active area of research in biology. The difficulty of the experimental approach (in particular, in its high-throughput version) leaves ample room for the development of bioinformatics tools that can provide a useful first picture of the problem. Among the possible approaches, one of the simplest is to follow classical protein function annotation protocols and annotate target alternative splice events with the information available from conserved events in other species. However, the application of this protocol requires a procedure capable of recognising such events. Here we present a simple but accurate method developed for this purpose.

Results

We have developed a method for identifying homologous, or equivalent, alternative splicing events, based on the combined use of neural networks and sequence searches. The procedure comprises four steps: (i) BLAST search for homologues of the two isoforms defining the target alternative splicing event; (ii) construction of all possible candidate events; (iii) scoring of the latter with a series of neural networks; and (iv) filtering of the results. When tested in a set of 473 manually annotated pairs of homologous events, our method showed a good performance, with an accuracy of 0.99, a precision of 0.98 and a sensitivity of 0.93. When no candidates were available, the specificity of our method varied between 0.81 and 0.91.

Conclusion

The method described in this article allows the identification of homologous alternative splicing events, with a good success rate, indicating that such method could be used for the development of functional annotation of alternative splice isoforms.

The (in)dependence of alternative splicing and gene duplication

Alternative splicing (AS) and gene duplication (GD) both are processes that diversify the protein repertoire. Recent examples have shown that sequence changes introduced by AS may be comparable to those introduced by GD. In addition, the two processes are inversely correlated at the genomic scale: large gene families are depleted in splice variants and vice versa. All together, these data strongly suggest that both phenomena result in interchangeability between their effects. Here, we tested the extent to which this applies with respect to various protein characteristics. The amounts of AS and GD per gene are anticorrelated even when accounting for different gene functions or degrees of sequence divergence. In contrast, the two processes appear to be independent in their influence on variation in mRNA expression. Further, we conducted a detailed comparison of the effect of sequence changes in both alternative splice variants and gene duplicates on protein structure, in particular the size, location, and types of sequence substitutions and insertions/deletions. We find that, in general, alternative splicing affects protein sequence and structure in a more drastic way than gene duplication and subsequent divergence. Our results reveal an interesting paradox between the anticorrelation of AS and GD at the genomic level, and their impact at the protein level, which shows little or no equivalence in terms of effects on protein sequence, structure, and function. We discuss possible explanations that relate to the order of appearance of AS and GD in a gene family, and to the selection pressure imposed by the environment.

Alarin is a vasoactive peptide

Galanin-like peptide (GALP) is a hypothalamic neuropeptide belonging to the galanin family of peptides. The GALP gene is characterized by extensive differential splicing in a variety of murine tissues. One splice variant excludes exon 3 and results in a frame shift leading to a novel peptide sequence and a stop codon after 49 aa. In this peptide, which we termed alarin, the signal sequence of the GALP precursor peptide and the first 5 aa of the mature GALP are followed by 20 aa without homology to any other murine protein. Alarin mRNA was detected in murine brain, thymus, and skin. In accordance with its vascular localization, the peptide exhibited potent and dose-dependent vasoconstrictor and anti-edema activity in the cutaneous microvasculature, as was also observed with other members of the galanin peptide family. However, in contrast to galanin peptides in general, the physiological effects of alarin do not appear to be mediated via the known galanin receptors. Alarin adds another facet to the surprisingly high-functional redundancy of the galanin family of peptides.

Plant SET domain-containing proteins: structure, function and regulation

Modification of the histone proteins that form the core around which chromosomal DNA is looped has profound epigenetic effects on the accessibility of the associated DNA for transcription, replication and repair. The SET domain is now recognized as generally having methyltransferase activity targeted to specific lysine residues of histone H3 or H4. There is considerable sequence conservation within the SET domain and within its flanking regions. Previous reviews have shown that SET proteins from Arabidopsis and maize fall into five classes according to their sequence and domain architectures. These classes generally reflect specificity for a particular substrate. SET proteins from rice were found to fall into similar groupings, strengthening the merit of the approach taken. Two additional classes, VI and VII, were established that include proteins with truncated/interrupted SET domains. Diverse mechanisms are involved in shaping the function and regulation of SET proteins. These include protein-protein interactions through both intra- and inter-molecular associations that are important in plant developmental processes, such as flowering time control and embryogenesis. Alternative splicing that can result in the generation of two to several different transcript isoforms is now known to be widespread. An exciting and tantalizing question is whether, or how, this alternative splicing affects gene function. For example, it is conceivable that one isoform may debilitate methyltransferase function whereas the other may enhance it, providing an opportunity for differential regulation. The review concludes with the speculation that modulation of SET protein function is mediated by antisense or sense-antisense RNA.

Transcript-based redefinition of grouped oligonucleotide probe sets using AceView: high-resolution annotation for microarrays

BACKGROUND

Extracting biological information from high-density Affymetrix arrays is a multi-step process that begins with the accurate annotation of microarray probes. Shortfalls in the original Affymetrix probe annotation have been described; however, few studies have provided rigorous solutions for routine data analysis.

RESULTS

Using AceView, a comprehensive human transcript database, we have reannotated the probes by matching them to RNA transcripts instead of genes. Based on this transcript-level annotation, a new probe set definition was created in which every probe in a probe set maps to a common set of AceView gene transcripts. In addition, using artificial data sets we identified that a minimal probe set size of 4 is necessary for reliable statistical summarization. We further demonstrate that applying the new probe set definition can detect specific transcript variants contributing to differential expression and it also improves cross-platform concordance.

CONCLUSION

We conclude that our transcript-level reannotation and redefinition of probe sets complement the original Affymetrix design. Redefinitions introduce probe sets whose sizes may not support reliable statistical summarization; therefore, we advocate using our transcript-level mapping redefinition in a secondary analysis step rather than as a replacement. Knowing which specific transcripts are differentially expressed is important to properly design probe/primer pairs for validation purposes. For convenience, we have created custom chip-description-files (CDFs) and annotation files for our new probe set definitions that are compatible with Bioconductor, Affymetrix Expression Console or third party software.

Development of a new set of reference genes for normalization of real-time RT-PCR data of porcine backfat and longissimus dorsi muscle, and evaluation with PPARGC1A

Background

An essential part of using real-time RT-PCR is that expression results have to be normalized before any conclusions can be drawn. This can be done by using one or multiple, validated reference genes, depending on the desired accuracy of the results. In the pig however, very little information is available on the expression stability of reference genes. The aim of this study was therefore to develop a new set of reference genes which can be used for normalization of mRNA expression data of genes expressed in porcine backfat and longissimus dorsi muscle, both representing an economically important part of a pig's carcass. Because of its multiple functions in fat metabolism and muscle fibre type composition, peroxisome proliferative activated receptor γ coactivator 1α (PPARGC1A) is a very interesting candidate gene for meat quality, and was an ideal gene to evaluate our developed set of reference genes for normalization of mRNA expression data of both tissue types.

Results

The mRNA expression stability of 10 reference genes was determined. The expression of RPL13A and SDHA appeared to be highly unstable. After normalization to the geometric mean of the three most stably expressed reference genes (ACTB, TBP and TOP2B), the results not only showed that the mRNA expression of PPARGC1A was significantly higher in each of the longissimus dorsi muscle samples than in backfat (P < 0.05), but also that the expression was significantly higher in the most cranial of the three muscle samples (P < 0.05).

Conclusion

This study provides a new set of reference genes (ACTB, TBP and TOP2B) suitable for normalization of real-time RT-PCR data of backfat and longissimus dorsi muscle in the pig. The obtained PPARGC1A expression results, after application of this set of reference genes, are a first step in unravelling the PPARGC1A expression pattern in the pig and provide a basis for possible selection towards improved meat quality while maintaining a lean carcass.

Comprehensive comparisons of the current human, mouse, and rat RefSeq, Ensembl, EST, and FANTOM3 datasets: Identification of new human genes with specific tissue expression profile

Our understanding of functional genetic elements in the genomes is continuously growing and new entries are entered in various databases on a regular basis. We have here merged the genetic elements in RefSeq, Ensembl, FANTOM3, HINV, and NCBI:s ESTdb using the genome assemblies in order to achieve a comprehensive picture of the current status of the identity and gene number in human, mouse, and rat. The number of human protein coding genes has not increased (25,043) while the increased sequencing of mouse transcripts has provided the considerably higher number of protein coding genes (31,578) in mouse. The results indicate large discrepancies between the datasets, as considerable numbers of unique transcripts can be found in each dataset. Despite the high number of ncRNA (38,129 in mouse) there are also almost 20,000 EST clusters in both mouse and humans with more than one EST that do not overlap any transcript suggesting that several new genetic elements are still to be found. We also demonstrated presence of new genes by identifying new human ones that have specific tissue profiles, using RT-PCR on rat tissues.

PROMPT: a protein mapping and comparison tool

Background

Comparison of large protein datasets has become a standard task in bioinformatics. Typically researchers wish to know whether one group of proteins is significantly enriched in certain annotation attributes or sequence properties compared to another group, and whether this enrichment is statistically significant. In order to conduct such comparisons it is often required to integrate molecular sequence data and experimental information from disparate incompatible sources. While many specialized programs exist for comparisons of this kind in individual problem domains, such as expression data analysis, no generic software solution capable of addressing a wide spectrum of routine tasks in comparative proteomics is currently available.

Results

PROMPT is a comprehensive bioinformatics software environment which enables the user to compare arbitrary protein sequence sets, revealing statistically significant differences in their annotation features. It allows automatic retrieval and integration of data from a multitude of molecular biological databases as well as from a custom XML format. Similarity-based mapping of sequence IDs makes it possible to link experimental information obtained from different sources despite discrepancies in gene identifiers and minor sequence variation. PROMPT provides a full set of statistical procedures to address the following four use cases: i) comparison of the frequencies of categorical annotations between two sets, ii) enrichment of nominal features in one set with respect to another one, iii) comparison of numeric distributions, and iv) correlation of numeric variables. Analysis results can be visualized in the form of plots and spreadsheets and exported in various formats, including Microsoft Excel.

Conclusion

PROMPT is a versatile, platform-independent, easily expandable, stand-alone application designed to be a practical workhorse in analysing and mining protein sequences and associated annotation. The availability of the Java Application Programming Interface and scripting capabilities on one hand, and the intuitive Graphical User Interface with context-sensitive help system on the other, make it equally accessible to professional bioinformaticians and biologically-oriented users. PROMPT is freely available for academic users from .

An expectation-maximization algorithm for probabilistic reconstructions of full-length isoforms from splice graphs

Reconstructing full-length transcript isoforms from sequence fragments (such as ESTs) is a major interest and challenge for bioinformatic analysis of pre-mRNA alternative splicing. This problem has been formulated as finding traversals across the splice graph, which is a directed acyclic graph (DAG) representation of gene structure and alternative splicing. In this manuscript we introduce a probabilistic formulation of the isoform reconstruction problem, and provide an expectation-maximization (EM) algorithm for its maximum likelihood solution. Using a series of simulated data and expressed sequences from real human genes, we demonstrate that our EM algorithm can correctly handle various situations of fragmentation and coupling in the input data. Our work establishes a general probabilistic framework for splice graph-based reconstructions of full-length isoforms.

Induction of revertant fibres in the mdx mouse using antisense oligonucleotides

Background

Duchenne muscular dystrophy is a fatal genetic disorder caused by dystrophin gene mutations that result in premature termination of translation and the absence of functional protein. Despite the primary dystrophin gene lesion, immunostaining studies have shown that at least 50% of DMD patients, mdx mice and a canine model of DMD have rare dystrophin-positive or 'revertant' fibres. Fine epitope mapping has shown that the majority of transcripts responsible for revertant fibres exclude multiple exons, one of which includes the dystrophin mutation.

Methods

The mdx mouse model of muscular dystrophy has a nonsense mutation in exon 23 of the dystrophin gene. We have shown that antisense oligonucleotides (AOs) can induce the removal of this exon, resulting in an in-frame mRNA transcript encoding a shortened but functional dystrophin protein. To emulate one exonic combination associated with revertant fibres, we target multiple exons for removal by the application of a group of AOs combined as a "cocktail".

Results

Exons 19–25 were consistently excluded from the dystrophin gene transcript using a cocktail of AOs. This corresponds to an alternatively processed gene transcript that has been sporadically detected in untreated dystrophic mouse muscle, and is presumed to give rise to a revertant dystrophin isoform. The transcript and the resultant correctly localised smaller protein were confirmed by RT-PCR, immunohistochemistry and western blot analysis.

Conclusion

This work demonstrates the feasibility of AO cocktails to by-pass dystrophin mutation hotspots through multi-exon skipping. Multi-exon skipping could be important in expediting an exon skipping therapy to treat DMD, so that the same AO formulations may be applied to several different mutations within particular domains of the dystrophin gene.