Identification Of Exons Research Articles

Identification of Protein-Coding Regions in DNA Sequences Using A Time-Frequency Filtering Approach

Accurate identification of protein-coding regions (exons) in DNA sequences has been a challenging task in bioinformatics. Particularly the coding regions have a 3-base periodicity, which forms the basis of all exon identification methods. Many signal processing tools and techniques have been applied successfully for the identification task but still improvement in this direction is needed. In this paper, we have introduced a new promising model-independent time-frequency filtering technique based on S-transform for accurate identification of the coding regions. The S-transform is a powerful linear time-frequency representation useful for filtering in time-frequency domain. The potential of the proposed technique has been assessed through simulation study and the results obtained have been compared with the existing methods using standard datasets. The comparative study demonstrates that the proposed method outperforms its counterparts in identifying the coding regions.

DBASS3 and DBASS5: databases of aberrant 3'- and 5'-splice sites

DBASS3 and DBASS5 provide comprehensive repositories of new exon boundaries that were induced by pathogenic mutations in human disease genes. Aberrant 5′- and 3′-splice sites were activated either by mutations in the consensus sequences of natural exon–intron junctions (cryptic sites) or elsewhere (‘de novo’ sites). DBASS3 and DBASS5 currently contain approximately 900 records of cryptic and de novo 3′- and 5′-splice sites that were produced by over a thousand different mutations in approximately 360 genes. DBASS3 and DBASS5 data can be searched by disease phenotype, gene, mutation, location of aberrant splice sites in introns and exons and their distance from authentic counterparts, by bibliographic references and by the splice-site strength estimated with several prediction algorithms. The user can also retrieve reference sequences of both aberrant and authentic splice sites with the underlying mutation. These data will facilitate identification of introns or exons frequently involved in aberrant splicing, mutation analysis of human disease genes and study of germline or somatic mutations that impair RNA processing. Finally, this resource will be useful for fine-tuning splice-site prediction algorithms, better definition of auxiliary splicing signals and design of new reporter assays. DBASS3 and DBASS5 are freely available at http://www.dbass.org.uk/.

Convergent origins and rapid evolution of spliced leader trans-splicing in Metazoa: Insights from the Ctenophora and Hydrozoa

Replacement of mRNA 5' UTR sequences by short sequences trans-spliced from specialized, noncoding, spliced leader (SL) RNAs is an enigmatic phenomenon, occurring in a set of distantly related animal groups including urochordates, nematodes, flatworms, and hydra, as well as in Euglenozoa and dinoflagellates. Whether SL trans-splicing has a common evolutionary origin and biological function among different organisms remains unclear. We have undertaken a systematic identification of SL exons in cDNA sequence data sets from non-bilaterian metazoan species and their closest unicellular relatives. SL exons were identified in ctenophores and in hydrozoan cnidarians, but not in other cnidarians, placozoans, or sponges, or in animal unicellular relatives. Mapping of SL absence/presence obtained from this and previous studies onto current phylogenetic trees favors an evolutionary scenario involving multiple origins for SLs during eumetazoan evolution rather than loss from a common ancestor. In both ctenophore and hydrozoan species, multiple SL sequences were identified, showing high sequence diversity. Detailed analysis of a large data set generated for the hydrozoan Clytia hemisphaerica revealed trans-splicing of given mRNAs by multiple alternative SLs. No evidence was found for a common identity of trans-spliced mRNAs between different hydrozoans. One feature found specifically to characterize SL-spliced mRNAs in hydrozoans, however, was a marked adenosine enrichment immediately 3' of the SL acceptor splice site. Our findings of high sequence divergence and apparently indiscriminate use of SLs in hydrozoans, along with recent findings in other taxa, indicate that SL genes have evolved rapidly in parallel in diverse animal groups, with constraint on SL exon sequence evolution being apparently rare.

Identification of novel exons and transcribed regions by chimpanzee transcriptome sequencing.

BackgroundWe profile the chimpanzee transcriptome by using deep sequencing of cDNA from brain and liver, aiming to quantify expression of known genes and to identify novel transcribed regions.ResultsUsing stringent criteria for transcription, we identify 12,843 expressed genes, with a majority being found in both tissues. We further identify 9,826 novel transcribed regions that are not overlapping with annotated exons, mRNAs or ESTs. Over 80% of the novel transcribed regions map within or in the vicinity of known genes, and by combining sequencing data with de novo splice predictions we predict several of the novel transcribed regions to be new exons or 3' UTRs. For approximately 350 novel transcribed regions, the corresponding DNA sequence is absent in the human reference genome. The presence of novel transcribed regions in five genes and in one intergenic region is further validated with RT-PCR. Finally, we describe and experimentally validate a putative novel multi-exon gene that belongs to the ATP-cassette transporter gene family. This gene does not appear to be functional in human since one exon is absent from the human genome. In addition to novel exons and UTRs, novel transcribed regions may also stem from different types of noncoding transcripts. We note that expressed repeats and introns from unspliced mRNAs are especially common in our data.ConclusionsOur results extend the chimpanzee gene catalogue with a large number of novel exons and 3' UTRs and thus support the view that mammalian gene annotations are not yet complete.

Affy exon tissues: exon levels in normal tissues in human, mouse and rat

Summary: Most genes in human, mouse and rat produce more than one transcript isoform. The Affymetrix Exon Array is a tool for studying the many processes that regulate RNA production, with separate probesets measuring RNA levels at known and putative exons. For insights on how exons levels vary between normal tissues, we constructed the Affy Exon Tissues track from tissue data published by Affymetrix. This track reports exon probeset intensities as log ratios relative to median values across the dataset and renders them as colored heat maps, to yield quick visual identification of exons with intensities that vary between normal tissues.Availability: Affy Exon Tissues track is freely available under the UCSC Genome Browser (http://genome.ucsc.edu/) for human (hg18), mouse (mm8 and mm9), and rat (rn4).Contact: cline@soe.ucsc.eduSupplementary information: Supplementary data are available at Bioinformatics online.

When new exons are born

When new exons are born

Discovery and revision of Arabidopsis genes by proteogenomics

Gene annotation underpins genome science. Most often protein coding sequence is inferred from the genome based on transcript evidence and computational predictions. While generally correct, gene models suffer from errors in reading frame, exon border definition, and exon identification. To ascertain the error rate of Arabidopsis thaliana gene models, we isolated proteins from a sample of Arabidopsis tissues and determined the amino acid sequences of 144,079 distinct peptides by tandem mass spectrometry. The peptides corresponded to 1 or more of 3 different translations of the genome: a 6-frame translation, an exon splice-graph, and the currently annotated proteome. The majority of the peptides (126,055) resided in existing gene models (12,769 confirmed proteins), comprising 40% of annotated genes. Surprisingly, 18,024 novel peptides were found that do not correspond to annotated genes. Using the gene finding program AUGUSTUS and 5,426 novel peptides that occurred in clusters, we discovered 778 new protein-coding genes and refined the annotation of an additional 695 gene models. The remaining 13,449 novel peptides provide high quality annotation (>99% correct) for thousands of additional genes. Our observation that 18,024 of 144,079 peptides did not match current gene models suggests that 13% of the Arabidopsis proteome was incomplete due to approximately equal numbers of missing and incorrect gene models.

Improved identification of conserved cassette exons using Bayesian networks

BackgroundAlternative splicing is a major contributor to the diversity of eukaryotic transcriptomes and proteomes. Currently, large scale detection of alternative splicing using expressed sequence tags (ESTs) or microarrays does not capture all alternative splicing events. Moreover, for many species genomic data is being produced at a far greater rate than corresponding transcript data, hence in silico methods of predicting alternative splicing have to be improved.ResultsHere, we show that the use of Bayesian networks (BNs) allows accurate prediction of evolutionary conserved exon skipping events. At a stringent false positive rate of 0.5%, our BN achieves an improved true positive rate of 61%, compared to a previously reported 50% on the same dataset using support vector machines (SVMs). Incorporating several novel discriminative features such as intronic splicing regulatory elements leads to the improvement. Features related to mRNA secondary structure increase the prediction performance, corroborating previous findings that secondary structures are important for exon recognition. Random labelling tests rule out overfitting. Cross-validation on another dataset confirms the increased performance. When using the same dataset and the same set of features, the BN matches the performance of an SVM in earlier literature. Remarkably, we could show that about half of the exons which are labelled constitutive but receive a high probability of being alternative by the BN, are in fact alternative exons according to the latest EST data. Finally, we predict exon skipping without using conservation-based features, and achieve a true positive rate of 29% at a false positive rate of 0.5%.ConclusionBNs can be used to achieve accurate identification of alternative exons and provide clues about possible dependencies between relevant features. The near-identical performance of the BN and SVM when using the same features shows that good classification depends more on features than on the choice of classifier. Conservation based features continue to be the most informative, and hence distinguishing alternative exons from constitutive ones without using conservation based features remains a challenging problem.

Identify Alternative Splicing Events Based on Position-Specific Evolutionary Conservation

The evolution of eukaryotes is accompanied by the increased complexity of alternative splicing which greatly expands genome information. One of the greatest challenges in the post-genome era is a complete revelation of human transcriptome with consideration of alternative splicing. Here, we introduce a comparative genomics approach to systemically identify alternative splicing events based on the differential evolutionary conservation between exons and introns and the high-quality annotation of the ENCODE regions. Specifically, we focus on exons that are included in some transcripts but are completely spliced out for others and we call them conditional exons. First, we characterize distinguishing features among conditional exons, constitutive exons and introns. One of the most important features is the position-specific conservation score. There are dramatic differences in conservation scores between conditional exons and constitutive exons. More importantly, the differences are position-specific. For flanking intronic regions, the differences between conditional exons and constitutive exons are also position-specific. Using the Random Forests algorithm, we can classify conditional exons with high specificities (97% for the identification of conditional exons from intron regions and 95% for the classification of known exons) and fair sensitivities (64% and 32% respectively). We applied the method to the human genome and identified 39,640 introns that actually contain conditional exons and classified 8,813 conditional exons from the current RefSeq exon list. Among those, 31,673 introns containing conditional exons and 5,294 conditional exons classified from known exons cannot be inferred from RefSeq, UCSC or Ensembl annotations. Some of these de novo predictions were experimentally verified.

Revised genomic structure of the human ghrelin gene and identification of novel exons, alternative splice variants and natural antisense transcripts

BackgroundGhrelin is a multifunctional peptide hormone expressed in a range of normal tissues and pathologies. It has been reported that the human ghrelin gene consists of five exons which span 5 kb of genomic DNA on chromosome 3 and includes a 20 bp non-coding first exon (20 bp exon 0). The availability of bioinformatic tools enabling comparative analysis and the finalisation of the human genome prompted us to re-examine the genomic structure of the ghrelin locus.ResultsWe have demonstrated the presence of an additional novel exon (exon -1) and 5' extensions to exon 0 and 1 using comparative in silico analysis and have demonstrated their existence experimentally using RT-PCR and 5' RACE. A revised exon-intron structure demonstrates that the human ghrelin gene spans 7.2 kb and consists of six rather than five exons. Several ghrelin gene-derived splice forms were detected in a range of human tissues and cell lines. We have demonstrated ghrelin gene-derived mRNA transcripts that do not code for ghrelin, but instead may encode the C-terminal region of full-length preproghrelin (C-ghrelin, which contains the coding region for obestatin) and a transcript encoding obestatin-only. Splice variants that differed in their 5' untranslated regions were also found, suggesting a role of these regions in the post-transcriptional regulation of preproghrelin translation. Finally, several natural antisense transcripts, termed ghrelinOS (ghrelin opposite strand) transcripts, were demonstrated via orientation-specific RT-PCR, 5' RACE and in silico analysis of ESTs and cloned amplicons.ConclusionThe sense and antisense alternative transcripts demonstrated in this study may function as non-coding regulatory RNA, or code for novel protein isoforms. This is the first demonstration of putative obestatin and C-ghrelin specific transcripts and these findings suggest that these ghrelin gene-derived peptides may also be produced independently of preproghrelin. This study reveals several novel aspects of the ghrelin gene and suggests that the ghrelin locus is far more complex than previously recognised.

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury.

BackgroundThe sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome.ResultsBioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury.ConclusionUsing the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

Novel aspects of the sperm-associated antigen 11 (SPAG11) gene organization and expression in cattle (Bos taurus).

Beta-defensins are small cationic peptides exhibiting broad spectrum antimicrobial properties. In humans, many beta-defensin genes are located within a cluster on chromosome 8p23. The sperm associated antigen 11 (SPAG11) gene is contained in this cluster and is unusual among the human beta-defensins due to its complex genomic structure and mRNA splicing pattern. Here we report the genomic organization of the Bos taurus SPAG11 gene located on chromosome 27q1.2, within a cluster of beta-defensin genes. The exon structures of the fused bovine SPAG11 gene and of the mosaic transcripts initiated at both A and B promoters were established, including identification of novel exons and transcripts not previously found in primate or rodent. Evolutionary analysis against primate, rodent, canine, and porcine orthologs was performed. In adult bulls SPAG11C, SPAG11E, and SPAG11U mRNAs were detected predominantly in the male reproductive tract, while SPAG11D transcript was detected in reproductive and nonreproductive tissues and SPAG11V and SPAG11W mRNAs were confined to testis. Differential expression of all six transcripts was observed in tissues from fetal and adult bulls, suggesting that similar mRNA splicing mechanisms govern SPAG11 gene expression during pre- and postnatal development. Immunolocalization of SPAG11C and SPAG11D/E was demonstrated in the epithelium of the epididymis and testis, and SPAG11D in association with epididymal spermatozoa. Recombinant full-length SPAG11D protein was strongly antibacterial, while the SPAG11E C-terminal peptide that contains the beta-defensin motif in its structure was somewhat less potent. Taken together, the results suggest that SPAG11 isoforms perform both immune and reproductive functions in cattle.

Identification, tissue expression, and glucocorticoid responsiveness of alternative first exons of the human glucocorticoid receptor

Transcripts for the human glucocorticoid receptor (NR3C1) are known to contain alternative first exons 1A1, 1A2, and 1A3 from the distal promoter or 1D, 1E, 1B, 1F, 1C, or 1H from the proximal promoter. Here, we report two additional alternative first exons identified by Rapid amplification of cDNA ends (RACE)-PCR. The first, exon 1I, starts approximately 700 bp downstream of the splice donor site of the longest form of exon 1A, 1A3, considerably extending the known distal promoter region with a region containing conserved transcription factor-binding sites as well as a potential glucocorticoid response element (GRE) that differs from the consensus GRE in only two positions. The second, exon 1J, is part of the proximal promoter region and resides between exons 1D and 1E. Since this has been determined by quantitative real-time reverse transcriptase (RT)-PCR, exon 1I is used foremost in cells of the T-lymphocyte lineage. In the T-ALL cell line CEM-C7H2, which is sensitive to glucocorticoid-induced apoptosis, transcripts containing alternative first exons from the distal as well as the proximal promoter regions were markedly autoinduced by glucocorticoid treatment, with more pronounced relative induction in the distal promoter. Neither transcript was autoinduced in the related, resistant cell lines CEM-C1, and CEM-C7R5. In contrast, the glucocorticoid-sensitive PreB697 cell line strongly autoinduced transcripts from the proximal promoter, but not transcripts from the distal promoter, to relevant levels. Therefore, the autoinductive feedback loop implicated in glucocorticoid-induced apoptosis cannot universally rely on the distal promoter of the glucocorticoid receptor.

SAGE2Splice: Unmapped SAGE Tags Reveal Novel Splice Junctions

Serial analysis of gene expression (SAGE) not only is a method for profiling the global expression of genes, but also offers the opportunity for the discovery of novel transcripts. SAGE tags are mapped to known transcripts to determine the gene of origin. Tags that map neither to a known transcript nor to the genome were hypothesized to span a splice junction, for which the exon combination or exon(s) are unknown. To test this hypothesis, we have developed an algorithm, SAGE2Splice, to efficiently map SAGE tags to potential splice junctions in a genome. The algorithm consists of three search levels. A scoring scheme was designed based on position weight matrices to assess the quality of candidates. Using optimized parameters for SAGE2Splice analysis and two sets of SAGE data, candidate junctions were discovered for 5%–6% of unmapped tags. Candidates were classified into three categories, reflecting the previous annotations of the putative splice junctions. Analysis of predicted tags extracted from EST sequences demonstrated that candidate junctions having the splice junction located closer to the center of the tags are more reliable. Nine of these 12 candidates were validated by RT-PCR and sequencing, and among these, four revealed previously uncharacterized exons. Thus, SAGE2Splice provides a new functionality for the identification of novel transcripts and exons. SAGE2Splice is available online at http://www.cisreg.ca.

How Many New Genes Are There?

How Many New Genes Are There?

Identification of human exons overexpressed in tumors through the use of genome and expressed sequence data

Alternative splicing is one of the major sources of the large transcriptional diversity found in human cells. Splicing variants have been shown to be associated with features like spreading and progression in several human tumors. Therefore, such variants may be of great importance as both diagnostic and therapeutic tools. Here, by using a set of criteria regarding the expression pattern of splicing variants and statistical analyses, we were able to screen the genome for exons overexpressed in tumors of specific tissues. However, as in other analyses attempting to identify tumor-associated variants, our list of candidates was seriously inflated with cases of genes differentially expressed in tumors. To exclude these cases and increase the probability of finding bona fide regulated splicing variants, we performed a serial analysis of gene expression (SAGE), excluding those genes that were shown to be upregulated in tumors. This allowed us to predict the overexpression of single exons in specific tumors. Our final group of candidates includes 1,386 exons belonging to 638 genes. Experimental validation of a few candidates in normal tissue, tumor cell lines, and patient samples suggests that most of these candidates are indeed tumor-associated exons. Further functional classification of our candidate genes shows that our final list is slightly inflated with cancer-related genes.

Identification of Novel Exons from Rat–Mouse Comparisons

Exon shuffling, a major mechanism of gene evolution, scrambles existing sequences to create new genes. However, is it possible for an exon to be created from scratch? Here we conduct a survey of rat and mouse genomes and identify 2302 putative rodent-specific exons absent from the human genome. Analysis of rodent transcripts supporting these exons indicates that over half appear to be alternatively spliced in genes orthologous between rodents and human. This study demonstrates the importance of sequencing genomes from multiple species to accurately document the evolution of gene structure.

Another face of the Treacher Collins syndrome (TCOF1) gene: identification of additional exons.

Treacher Collins syndrome (TCS) is characterized by an abnormality in craniofacial development during early embryogenesis. TCS is caused by mutations in the gene TCOF1, which encodes the nucleolar phosphoprotein treacle. Genetic and proteomic characterizations of TCS/treacle are based on the previously reported 26 exons of TCOF1. Here, we report the identification of 231-nucleotide (nt) exon 6A (between exons 6 and 7) and 108-nt exon 16A (between exons 16 and 17). Isoforms with exon 6A are up to 3.7-fold more abundant than alternatively spliced variants without exon 6A, but only minor isoforms contain exon 16A. Exon 6A encodes a peptide sequence containing basic and acidic domains similar to 10 other exons of TCOF1. Unlike the other exons, exon 6A encodes a nuclear localization signal (NLS) which does not, however, alter the nucleolar localization of full-length treacle. The discovery of exons 6A and 16A is relevant to mutational analysis of the TCOF1 gene in TCS patients, and to functional analysis of its gene product.

Vector space classification of DNA sequences

Revisiting the problem of intron–exon identification, we use a principal component analysis (PCA) to classify DNA sequences and present first results that validate our approach. Sequences are translated into document vectors that represent their word content; a principal component analysis then defines Gaussian-distributed sequence classes. The classification uses word content and variation of word usage to distinguish sequences. We test our approach with several data sets of genomic DNA and are able to classify introns and exons with an accuracy of up to 96%. We compare the method with the best traditional coding measure, the non-overlapping hexamer frequency count, and find that the PCA method produces better results. We also investigate the degree of cross-validation between different data sets of introns and exons and find evidence that the quality of a data set can be detected.

CHAPTER 3 - Cancer and Epigenesis: A Developmental Perspective

Traditionally, molecular geneticists have studied the genome, that is, DNA. This has included its sequencing, with identification of promoters, enhancers, introns, exons, and mutations. During the last 10 to 15 years, it has become clear that this study of “naked” DNA imposed major limitations on our understanding of gene regulation, and that DNA must be studied in conjunction with its protein backbone (chromatin). This new vision has shown that chromatin is a very dynamic molecule, and that changes in DNA methylation, in histones (eg, their localization, density, and whether they are acetylated, methylated, phosphorylated and/or ubiquitinated), and in the Polycomb-Trithorax equilibrium are all crucial for the control of gene expression. Furthermore, it was found that like DNA, chromatin-related proteins and DNA methylation are generally passed on unchanged from one cell to its daughters. Thus, gene expression is partly controlled through chromatin modifications that are transmitted from one cell to all its descendants. Such a control is referred to as “epigenetic,” as the DNA sequence is not altered. Hence, these epigenetic “marks” can be erased in the early embryo, rendering these cells totipotent. This chapter presents a basic overview of epigenesis in the control of normal embryonal cell differentiation, including imprinting. Similarities between embryonal and cancer cells are highlighted, and the potential impact of chemotherapy aimed at epigenetic mechanisms is reviewed.