Modeling dependencies in pre-mRNA splicing signals

Intron exon boundary junctions in human genome have in-built unique structural and energetic signals

Precise identification of correct exon–intron boundaries is a prerequisite to analyze the location and structure of genes. The existing framework for genomic signals, delineating exon and introns in a genomic segment, seems insufficient, predominantly due to poor sequence consensus as well as limitations of training on available experimental data sets. We present here a novel concept for characterizing exon–intron boundaries in genomic segments on the basis of structural and energetic properties. We analyzed boundary junctions on both sides of all the exons (3 28 368) of protein coding genes from human genome (GENCODE database) using 28 structural and three energy parameters. Study of sequence conservation at these sites shows very poor consensus. It is observed that DNA adopts a unique structural and energy state at the boundary junctions. Also, signals are somewhat different for housekeeping and tissue specific genes. Clustering of 31 parameters into four derived vectors gives some additional insights into the physical mechanisms involved in this biological process. Sites of structural and energy signals correlate well to the positions playing important roles in pre-mRNA splicing.

Genome assembly and annotation at the chromosomal level of first Pleuronectidae: Verasper variegatus provides a basis for phylogenetic study of Pleuronectiformes

High quality genome is of great significance for the mining of biological information resources of species. Up to now, the genomic information of several important economic flatfishes has been well explained. All these fishes are eyes on left side-type, and no high-quality genome of eyes on right side-type species has been reported. In this study, we applied a combined strategy involving stLFR and Hi-C technologies to generate sequencing data for constructing the chromosomal genome of Verasper variegates, which belongs to Pleuronectidae with characteristic of eyes on right side. The size of genome of V. variegatus is 556 Mb. More than 97.2% of BUSCO genes were detected, and N50 lengths of the contigs and scaffolds reached 79.8 Kb and 23.8 Mb, respectively, demonstrating the outstanding completeness and sequence continuity of the genome. A total of 22,199 protein-coding genes were predicted in the assembled genome, and more than 95% of those genes could be functionally annotated. Meanwhile, the genomic collinearity, gene family and phylogenetic analyses of similar species in Pleuronectiformes were also investigated and portrayed for metamorphosis and benthic adaptation. Sex related genes mapping has also been achieved at the chromosome level. This study is the first chromosomal level genome of a Pleuronectidae fish (V. variegatus). The chromosomal genome assembly constructed in this work will not only be valuable for conservation and aquaculture studies of the V. variegatus but will also be of general interest in the phylogenetic and taxonomic studies of Pleuronectiformes.

Recognition of alternatively spliced cassette exons based on a hybrid model

Alternative splicing (AS) is an important mechanism of gene regulation that contributes to protein diversity. It is of great significance to recognize different kinds of AS accurately so as to understand the mechanism of gene regulation. Many in silico methods have been applied to detecting AS with vast features, but the result is far from satisfactory. In this paper, we used the features proven to be useful in recognizing AS in previous literature and proposed a hybrid method combining Gene Expression Programming (GEP) and Random Forests (RF) to classify the constitutive exons and cassette exons which is the most common AS phenomenon. GEP will firstly make prediction to the samples of strong signal, and the other samples of weak signal will be distinguished with a more complex classifier based on RF. The experiment result indicates that this method can highly improve the recognition level in this issue.

Haplotype analysis of the promoter region of phosphodiesterase type 8B (PDE8B) in correlation with inactivating PDE8B mutation and the serum thyroid-stimulating hormone levels

BACKGROUND

Human phosphodiesterase (PDE) type 8B (PDE8B) is located at 5q14.1 and is known as the PDE with the highest affinity to cAMP. We recently described a family with bilateral micronodular adrenocortical disease that was apparently caused by an inactivating PDE8B mutation (H305P). As a result of a genome-wide study, a strong association between six polymorphic variants in the PDE8B promoter and serum levels of the thyroid-stimulating hormone (TSH) has been recently reported. Despite an extended analysis of the regions surrounding 5q14.1, no other potential genetic variants that could be responsible for the associated TSH levels were found.

METHODS

In this study, we genotyped by polymerase chain reaction the described six polymorphic variants in the PDE8B promoter in the family with micronodular adrenocortical disease and inactivating PDE8B mutation and analyzed their correlation with individual TSH values in the family members.

RESULTS

We observed complete segregation between the reported association and individual TSH values in the family we studied. Haplotype analysis showed that the haplotype associated with the high TSH levels is different from the one that segregated with H305P, suggesting that the mutation most probably has arisen on an allele independent of the high TSH-associated allele.

CONCLUSIONS

The proposed mechanism by which PDE8B may influence TSH levels is through control of cAMP signaling. Our analysis revealed separate segregation of an inactivating PDE8B allele from the high-TSH-allele and showed low TSH levels in persons who carry an inactivating PDE8B allele. These data suggest that, indeed, PDE8B may be involved in regulation of TSH levels.

Recognition of atypical 5' splice sites by shifted base-pairing to U1 snRNA

Accurate pre-mRNA splicing is critical for gene expression. The 5' splice site (5' ss) — the highly diverse element at the 5' end of introns — is initially recognized via base-pairing to the 5' end of U1 small nuclear RNA (snRNA). However, many natural 5' ss have a very poor match to the consensus sequence, and are predicted to be very weak. Using genetic suppression experiments in human cells, we demonstrate that some atypical 5' ss are actually efficiently recognized by U1, in an alternative base-pairing register that is shifted by one nucleotide. These atypical 5' ss are phylogenetically widespread, and many of them are conserved. Moreover, shifted base-pairing provides an explanation for the effect of a 5' ss mutation associated with pontocerebellar hypoplasia. The unexpected flexibility in 5' ss/U1 base-pairing challenges an established paradigm, and has broad implications for splice-site prediction algorithms and gene-annotation efforts in genome projects.

Fast splice site detection using information content and feature reduction

BACKGROUND

Accurate identification of splice sites in DNA sequences plays a key role in the prediction of gene structure in eukaryotes. Already many computational methods have been proposed for the detection of splice sites and some of them showed high prediction accuracy. However, most of these methods are limited in terms of their long computation time when applied to whole genome sequence data.

RESULTS

In this paper we propose a hybrid algorithm which combines several effective and informative input features with the state of the art support vector machine (SVM). To obtain the input features we employ information content method based on Shannon's information theory, Shapiro's score scheme, and Markovian probabilities. We also use a feature elimination scheme to reduce the less informative features from the input data.

CONCLUSION

In this study we propose a new feature based splice site detection method that shows improved acceptor and donor splice site detection in DNA sequences when the performance is compared with various state of the art and well known methods.

Biased clustered substitutions in the human genome: the footprints of male-driven biased gene conversion

We examined fixed substitutions in the human lineage since divergence from the common ancestor with the chimpanzee, and determined what fraction are AT to GC (weak-to-strong). Substitutions that are densely clustered on the chromosomes show a remarkable excess of weak-to-strong "biased" substitutions. These unexpected biased clustered substitutions (UBCS) are common near the telomeres of all autosomes but not the sex chromosomes. Regions of extreme bias are enriched for genes. Human and chimp orthologous regions show a striking similarity in the shape and magnitude of their respective UBCS maps, suggesting a relatively stable force leads to clustered bias. The strong and stable signal near telomeres may have participated in the evolution of isochores. One exception to the UBCS pattern found in all autosomes is chromosome 2, which shows a UBCS peak midchromosome, mapping to the fusion site of two ancestral chromosomes. This provides evidence that the fusion occurred as recently as 740,000 years ago and no more than approximately 3 million years ago. No biased clustering was found in SNPs, suggesting that clusters of biased substitutions are selected from mutations. UBCS is strongly correlated with male (and not female) recombination rates, which explains the lack of UBCS signal on chromosome X. These observations support the hypothesis that biased gene conversion (BGC), specifically in the male germline, played a significant role in the evolution of the human genome.

Aberrant 5' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

Despite a growing number of splicing mutations found in hereditary diseases, utilization of aberrant splice sites and their effects on gene expression remain challenging to predict. We compiled sequences of 346 aberrant 5′splice sites (5′ss) that were activated by mutations in 166 human disease genes. Mutations within the 5′ss consensus accounted for 254 cryptic 5′ss and mutations elsewhere activated 92 de novo 5′ss. Point mutations leading to cryptic 5′ss activation were most common in the first intron nucleotide, followed by the fifth nucleotide. Substitutions at position +5 were exclusively G>A transitions, which was largely attributable to high mutability rates of C/G>T/A. However, the frequency of point mutations at position +5 was significantly higher than that observed in the Human Gene Mutation Database, suggesting that alterations of this position are particularly prone to aberrant splicing, possibly due to a requirement for sequential interactions with U1 and U6 snRNAs. Cryptic 5′ss were best predicted by computational algorithms that accommodate nucleotide dependencies and not by weight-matrix models. Discrimination of intronic 5′ss from their authentic counterparts was less effective than for exonic sites, as the former were intrinsically stronger than the latter. Computational prediction of exonic de novo 5′ss was poor, suggesting that their activation critically depends on exonic splicing enhancers or silencers. The authentic counterparts of aberrant 5′ss were significantly weaker than the average human 5′ss. The development of an online database of aberrant 5′ss will be useful for studying basic mechanisms of splice-site selection, identifying splicing mutations and optimizing splice-site prediction algorithms.

Extensive in silico analysis of NF1 splicing defects uncovers determinants for splicing outcome upon 5' splice-site disruption

We describe 94 pathogenic NF1 gene alterations in a cohort of 97 Austrian neurofibromatosis type 1 patients meeting the NIH criteria. All mutations were fully characterized at the genomic and mRNA levels. Over half of the patients carried novel mutations, and only a quarter carried recurrent minor-lesion mutations at 16 mutational warm spots. The remaining patients carried NF1 microdeletions (7%) and rare recurring mutations. Thirty-six of the mutations (38%) altered pre-mRNA splicing, and fall into five groups: exon skipping resulting from mutations at authentic splice sites (type I), cryptic exon inclusion caused by deep intronic mutations (type II), creation of de novo splice sites causing loss of exonic sequences (type III), activation of cryptic splice sites upon authentic splice-site disruption (type IV), and exonic sequence alterations causing exon skipping (type V). Extensive in silico analyses of 37 NF1 exons and surrounding intronic sequences suggested that the availability of a cryptic splice site combined with a strong natural upstream 3' splice site (3'ss)is the main determinant of cryptic splice-site activation upon 5' splice-site disruption. Furthermore, the exonic sequences downstream of exonic cryptic 5' splice sites (5'ss) resemble intronic more than exonic sequences with respect to exonic splicing enhancer and silencer density, helping to distinguish between exonic cryptic and pseudo 5'ss. This study provides valuable predictors for the splicing pathway used upon 5'ss mutation, and underscores the importance of using RNA-based techniques, together with methods to identify microdeletions and intragenic copy-number changes, for effective and reliable NF1 mutation detection.

A computational survey of candidate exonic splicing enhancer motifs in the model plant Arabidopsis thaliana

Background

Algorithmic approaches to splice site prediction have relied mainly on the consensus patterns found at the boundaries between protein coding and non-coding regions. However exonic splicing enhancers have been shown to enhance the utilization of nearby splice sites.

Results

We have developed a new computational technique to identify significantly conserved motifs involved in splice site regulation. First, 84 putative exonic splicing enhancer hexamers are identified in Arabidopsis thaliana. Then a Gibbs sampling program called ELPH was used to locate conserved motifs represented by these hexamers in exonic regions near splice sites in confirmed genes. Oligomers containing 35 of these motifs have been shown experimentally to induce significant inclusion of A. thaliana exons. Second, integration of our regulatory motifs into two different splice site recognition programs significantly improved the ability of the software to correctly predict splice sites in a large database of confirmed genes. We have released GeneSplicerESE, the improved splice site recognition code, as open source software.

Conclusion

Our results show that the use of the ESE motifs consistently improves splice site prediction accuracy.

Non-EST-based prediction of novel alternatively spliced cassette exons with cell signaling function in Caenorhabditis elegans and human

To better understand the complex role that alternative splicing plays in intracellular signaling, it is important to catalog the numerous splice variants involved in signal transduction. Therefore, we developed PASE (Prediction of Alternative Signaling Exons), a computational tool to identify novel alternative cassette exons that code for kinase phosphorylation or signaling protein-binding sites. We first applied PASE to the Caenorhabditis elegans genome. In this organism, our algorithm had an overall specificity of ≥76.4%, including 33 novel cassette exons that we experimentally verified. We then used PASE to analyze the human genome and made 804 predictions, of which 308 were found as alternative exons in the transcript database. We experimentally tested 384 of the remaining unobserved predictions and discovered 26 novel human exons for a total specificity of ≥41.5% in human. By using a test set of known alternatively spliced signaling exons, we determined that the sensitivity of PASE is ∼70%. GO term analysis revealed that our exon predictions were found in the introns of known signal transduction genes more often than expected by chance, indicating PASE enriches for splice variants that function in signaling pathways. Overall, PASE was able to uncover 59 novel alternative cassette exons in C. elegans and humans through a genome-wide ab initio prediction method that enriches for exons involved in signaling.

Aberrant 3' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

The frequency distribution of mutation-induced aberrant 3' splice sites (3'ss) in exons and introns is more complex than for 5' splice sites, largely owing to sequence constraints upstream of intron/exon boundaries. As a result, prediction of their localization remains a challenging task. Here, nucleotide sequences of previously reported 218 aberrant 3'ss activated by disease-causing mutations in 131 human genes were compared with their authentic counterparts using currently available splice site prediction tools. Each tested algorithm distinguished authentic 3'ss from cryptic sites more effectively than from de novo sites. The best discrimination between aberrant and authentic 3'ss was achieved by the maximum entropy model. Almost one half of aberrant 3'ss was activated by AG-creating mutations and approximately 95% of the newly created AGs were selected in vivo. The overall nucleotide structure upstream of aberrant 3'ss was characterized by higher purine content than for authentic sites, particularly in position -3, that may be compensated by more stringent requirements for positive and negative nucleotide signatures centred around position -11. A newly developed online database of aberrant 3'ss will facilitate identification of splicing mutations in a gene or phenotype of interest and future optimization of splice site prediction tools.

Characterization and transcriptional analysis of protein tyrosine phosphatase genes and an ankyrin repeat gene of the parasitoid Glyptapanteles indiensis polydnavirus in the parasitized host

Glyptapanteles indiensis (Braconidae, Hymenoptera) is an endoparasitoid of Lymantria dispar, the gypsy moth. Expression of G. indiensis polydnavirus (GiBV)-encoded genes within the pest host results in inhibition of immune response and development and alteration of physiology, enabling successful development of the parasitoid. Here, GiBV genome segment F (segF), an 18·6 kb segment shown to encode nine protein tyrosine phosphatase (PTP) genes and a single ankyrin repeat gene (ank), is analysed. PTPs have presumed function as regulators of signal transduction, while ankyrin repeat genes are hypothesized to function in inhibition of NF-κB signalling in the parasitized host. In this study, transcription of each gene was mapped by 5′- and 3′-RACE (rapid amplification of cDNA ends) and temporal and tissue-specific expression was examined in the parasitized host. For polydnavirus gene prediction in the parasitized host, no available gene prediction parameters were entirely precise. The mRNAs for each GiBV segF gene initiated between 30 and 112 bp upstream of the translation initiation codon. All were encoded in single open reading frames (ORFs), with the exception of PTP9, which was transcribed as a bicistronic message with the adjacent ank gene. RT-PCR indicated that all GiBV segF PTPs were expressed early in parasitization and, for most, expression was sustained over the course of at least 7 days after parasitization, suggesting importance in both early and sustained virus-induced immunosuppression and alteration of physiology. Tissue-specific patterns of PTP expression of GiBV segF genes were variable, suggesting differing roles in facilitating parasitism.

Determinants of the inherent strength of human 5' splice sites

We previously showed that the authentic 5' splice site (5'ss) of the first exon in the human beta-globin gene is intrinsically stronger than a cryptic 5'ss located 16 nucleotides upstream. Here we examined by mutational analysis the contribution of individual 5'ss nucleotides to discrimination between these two 5'ss. Based on the in vitro splicing efficiencies of a panel of 26 wild-type and mutant substrates in two separate 5'ss competition assays, we established a hierarchy of 5'ss and grouped them into three functional subclasses: strong, intermediate, and weak. Competition between two 5'ss from different subclasses always resulted in selection of the 5'ss that belongs to the stronger subclass. Moreover, each subclass has different characteristic features. Strong and intermediate 5'ss can be distinguished by their predicted free energy of base-pairing to the U1 snRNA 5' terminus (DeltaG). Whereas the extent of splicing via the strong 5'ss correlates well with the DeltaG, this is not the case for competition between intermediate 5'ss. Weak 5'ss were used only when the competing authentic 5'ss was inactivated by mutation. These results indicate that extensive complementarity to U1 snRNA exerts a dominant effect for 5'ss selection, but in the case of competing 5'ss with similarly modest complementarity to U1, the role of other 5'ss features is more prominent. This study reveals the importance of additional submotifs present in certain 5'ss sequences, whose characterization will be critical for understanding 5'ss selection in human genes.

Prediction of the prototype of the human Toll-like receptor gene family from the pufferfish, Fugu rubripes, genome

The insect Toll family of proteins and their mammalian counterparts seemingly shared one common ancestor and evolved independently. Here we demonstrated that the prototype of the mammalian-type (M-type) Toll family is shared by the fish and humans. According to the draft of the pufferfish Fugu genome project, the signature Toll-IL-1 receptor homology domain (TIR domain) has been conserved during evolution. FuguTLR2, 3, 5, 7, 8 and 9 members correspond structurally to respective mammalian TLRs. One Fugu TLR showed equally high amino acid identity to human TLR1, 6 and 10, and we named it FuguTLR1. Fugu rubripes has genes for TLR21 and 22, which are unique to fish. One possible interpretation of these findings is that TLR1, 2, 3, 4, 5, 7, 8, 9, 21 and 22 existed in the ancestral genome common to fish and mammals, and that TLR4 was lost in the fish lineage, while TLR21 and 22 were lost in the mammalian lineage. Strikingly, a solitary ascidian, Halocynthia roretzi, has only a few Toll-like proteins, which, like Caenorhabditis elegans Toll, represent primitive ones before the expansion of the Toll family. Therefore, the expansion of TLR genes should have occurred earlier than fish, but not C. intestinalis, separated evolutionarily from mammals. These results infer that the appearance of the M-type innate system was completed before or concomitant with the appearance of acquired immunity. We interpret the present data to mean that the differences of TLRs identified in this study between fishes and humans may be rather peripheral, partially due to selection pressure exerted by pathogens in distinct environments.

SHORT INTEGUMENTS1/SUSPENSOR1/CARPEL FACTORY, a Dicer homolog, is a maternal effect gene required for embryo development in Arabidopsis

The importance of maternal cells in controlling early embryogenesis is well understood in animal development, yet in plants the precise role of maternal cells in embryogenesis is unclear. We demonstrated previously that maternal activity of the SIN1 (SHORT INTEGUMENTS1) gene of Arabidopsis is essential for embryo pattern formation and viability, and that its postembryonic activity is required for several processes in reproductive development, including flowering time control and ovule morphogenesis. Here, we report the cloning of SIN1, and demonstrate its identity to the CAF (CARPEL FACTORY) gene important for normal flower morphogenesis and to the SUS1 (SUSPENSOR1) gene essential for embryogenesis. SIN1/SUS1/CAF has sequence similarity to the Drosophila melanogaster gene Dicer, which encodes a multidomain ribonuclease specific for double-stranded RNA, first identified by its role in RNA silencing. The Dicer protein is essential for temporal control of development in animals, through the processing of small RNA hairpins that in turn inhibit the translation of target mRNAs. Structural modeling of the wild-type and sin1 mutant proteins indicates that the RNA helicase domain of SIN1/SUS1/CAF is important for function. The mRNA was detected in floral meristems, ovules, and early embryos, consistent with the mutant phenotypes. A 3.3-kb region 5' of the SIN1/SUS1/CAF gene shows asymmetric parent-of-origin activity in the embryo: It confers transcriptional activation of a reporter gene in early embryos only when transmitted through the maternal gamete. These results suggest that maternal SIN1/SUS1/CAF functions early in Arabidopsis development, presumably through posttranscriptional regulation of specific mRNA molecules.

A computational analysis of sequence features involved in recognition of short introns

Splicing of short introns by the nuclear pre-mRNA splicing machinery is thought to proceed via an "intron definition" mechanism, in which the 5' and 3' splice sites (5'ss, 3'ss, respectively) are initially recognized and paired across the intron. Here, we describe a computational analysis of sequence features involved in recognition of short introns by using available transcript data from five eukaryotes with complete or nearly complete genomic sequences. The information content of five different transcript features was measured by using methods from information theory, and Monte Carlo simulations were used to determine the amount of information required for accurate recognition of short introns in each organism. We conclude: (i) that short introns in Drosophila melanogaster and Caenorhabditis elegans contain essentially all of the information for their recognition by the splicing machinery, and computer programs that simulate splicing specificity can predict the exact boundaries of approximately 95% of short introns in both organisms; (ii) that in yeast, the 5'ss, branch signal, and 3'ss can accurately identify intron locations but do not precisely determine the location of 3' cleavage in every intron; and (iii) that the 5'ss, branch signal, and 3'ss are not sufficient to accurately identify short introns in plant and human transcripts, but that specific subsets of candidate intronic enhancer motifs can be identified in both human and Arabidopsis that contribute dramatically to the accuracy of splicing simulators.