A computational survey of candidate exonic splicing enhancer motifs in the model plant <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

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.     

Evolutionarily emerged G tracts between the polypyrimidine tract and 3′ AG are splicing silencers enriched in genes involved in cancer

Background

The 3′ splice site (SS) at the end of pre-mRNA introns has a consensus sequence (Y)_nNYAG for constitutive splicing of mammalian genes. Deviation from this consensus could change or interrupt the usage of the splice site leading to alternative or aberrant splicing, which could affect normal cell function or even the development of diseases. We have shown that the position “N” can be replaced by a CA-rich RNA element called CaRRE1 to regulate the alternative splicing of a group of genes.

Results

Taking it a step further, we searched the human genome for purine-rich elements between the -3 and -10 positions of the 3′ splice sites of annotated introns. This identified several thousand such 3′SS; more than a thousand of them contain at least one copy of G tract. These sites deviate significantly from the consensus of constitutive splice sites and are highly associated with alterative splicing events, particularly alternative 3′ splice and intron retention. We show by mutagenesis analysis and RNA interference that the G tracts are splicing silencers and a group of the associated exons are controlled by the G tract binding proteins hnRNP H/F. Species comparison of a group of the 3′SS among vertebrates suggests that most (~87%) of the G tracts emerged in ancestors of mammals during evolution. Moreover, the host genes are most significantly associated with cancer.

Conclusion

We call these elements together with CaRRE1 regulatory RNA elements between the Py and 3′AG (REPA). The emergence of REPA in this highly constrained region indicates that this location has been remarkably permissive for the emergence of de novo regulatory RNA elements, even purine-rich motifs, in a large group of mammalian genes during evolution. This evolutionary change controls alternative splicing, likely to diversify proteomes for particular cellular functions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1143) contains supplementary material, which is available to authorized users.     

Genomic features defining exonic variants that modulate splicing

Background  

Single point mutations at both synonymous and non-synonymous positions within exons can have severe effects on gene function through disruption of splicing. Predicting these mutations in silico purely from the genomic sequence is difficult due to an incomplete understanding of the multiple factors that may be responsible. In addition, little is known about which computational prediction approaches, such as those involving exonic splicing enhancers and exonic splicing silencers, are most informative.     

Novel KRAS Gene Mutations in Sporadic Colorectal Cancer

Introduction

In this article, we report 7 novel KRAS gene mutations discovered while retrospectively studying the prevalence and pattern of KRAS mutations in cancerous tissue obtained from 56 Saudi sporadic colorectal cancer patients from the Eastern Province.

Methods

Genomic DNA was extracted from formalin-fixed, paraffin-embedded cancerous and noncancerous colorectal tissues. Successful and specific PCR products were then bi-directionally sequenced to detect exon 4 mutations while Mutector II Detection Kits were used for identifying mutations in codons 12, 13 and 61. The functional impact of the novel mutations was assessed using bioinformatics tools and molecular modeling.

Results

KRAS gene mutations were detected in the cancer tissue of 24 cases (42.85%). Of these, 11 had exon 4 mutations (19.64%). They harbored 8 different mutations all of which except two altered the KRAS protein amino acid sequence and all except one were novel as revealed by COSMIC database. The detected novel mutations were found to be somatic. One mutation is predicted to be benign. The remaining mutations are predicted to cause substantial changes in the protein structure. Of these, the Q150X nonsense mutation is the second truncating mutation to be reported in colorectal cancer in the literature.

Conclusions

Our discovery of novel exon 4 KRAS mutations that are, so far, unique to Saudi colorectal cancer patients may be attributed to environmental factors and/or racial/ethnic variations due to genetic differences. Alternatively, it may be related to paucity of clinical studies on mutations other than those in codons 12, 13, 61 and 146. Further KRAS testing on a large number of patients of various ethnicities, particularly beyond the most common hotspot alleles in exons 2 and 3 is needed to assess the prevalence and explore the exact prognostic and predictive significance of the discovered novel mutations as well as their possible role in colorectal carcinogenesis.     

Multisite Phosphorylation of the Guanine Nucleotide Exchange Factor Cdc24 during Yeast Cell Polarization

Background

Cell polarization is essential for processes such as cell migration and asymmetric cell division. A common regulator of cell polarization in most eukaryotic cells is the conserved Rho GTPase, Cdc42. In budding yeast, Cdc42 is activated by a single guanine nucleotide exchange factor, Cdc24. The mechanistic details of Cdc24 activation at the onset of yeast cell polarization are unclear. Previous studies have suggested an important role for phosphorylation of Cdc24, which may regulate activity or function of the protein, representing a key step in the symmetry breaking process.

Methodology/Principal Findings

Here, we directly ask whether multisite phosphorylation of Cdc24 plays a role in its regulation. We identify through mass spectrometry analysis over thirty putative in vivo phosphorylation sites. We first focus on sites matching consensus sequences for cyclin-dependent and p21-activated kinases, two kinase families that have been previously shown to phosphorylate Cdc24. Through site-directed mutagenesis, yeast genetics, and light and fluorescence microscopy, we show that nonphosphorylatable mutations of these consensus sites do not lead to any detectable consequences on growth rate, morphology, kinetics of polarization, or localization of the mutant protein. We do, however, observe a change in the mobility shift of mutant Cdc24 proteins on SDS-PAGE, suggesting that we have indeed perturbed its phosphorylation. Finally, we show that mutation of all identified phosphorylation sites does not cause observable defects in growth rate or morphology.

Conclusions/Significance

We conclude that lack of phosphorylation on Cdc24 has no overt functional consequences in budding yeast. Yeast cell polarization may be more tightly regulated by inactivation of Cdc42 by GTPase activating proteins or by alternative methods of Cdc24 regulation, such as conformational changes or oligomerization.     

An alternative-exon database and its statistical analysis

We compiled a comprehensive database of alternative exons from the literature and analyzed them statistically. Most alternative exons are cassette exons and are expressed in more than two tissues. Of all exons whose expression was reported to be specific for a certain tissue, the majority were expressed in the brain. Whereas the length of constitutive exons follows a normal distribution, the distribution of alternative exons is skewed toward smaller ones. Furthermore, alternative-exon splice sites deviate more from the consensus: their 3' splice sites are characterized by a higher purine content in the polypyrimidine stretch, and their 5' splice sites deviate from the consensus sequence mostly at the +4 and +5 positions. Furthermore, for exons expressed in a single tissue, adenosine is more frequently used at the -3 position of the 3' splice site. In addition to the known AC-rich and purine-rich exonic sequence elements, sequence comparison using a Gibbs algorithm identified several motifs in exons surrounded by weak splice sites and in tissue-specific exons. Together, these data indicate a combinatorial effect of weak splice sites, atypical nucleotide usage at certain positions, and functional enhancers as an important contribution to alternative-exon regulation.     

Exonic splicing regulatory elements skew synonymous codon usage near intron-exon boundaries in mammals

In mammals there is a bias in amino acid usage near splice sites that is explained, in large part, by the high density of exonic splicing enhancers (ESEs) in these regions. Is there a similar bias for the relative use of synonymous codons, and can any such bias be predicted by their abundance in ESEs? Prior reports suggested that such trends may exist. From analysis of human exons, we find that 47 of the 59 codons with at least one synonym show differential usage in the proximity of exon ends, of which 42 remain significant after correction for multiple testing. Within sets of synonymous codons those more preferred near splice sites are generally those that are relatively more abundant within the ESEs. However, the examples given previously appear exceptionally good fits and there exist many exceptions, the usage of lysine's codons being a case in point. Similar results are observed in mouse exons. We conclude that splice regulation impacts on the choice of synonymous codons in mammals, but the magnitude of this effect is less than might at first have been supposed.     

The distribution and mutagenesis of short coding INDELs from 1,128 whole exomes

Background

Identifying insertion/deletion polymorphisms (INDELs) with high confidence has been intrinsically challenging in short-read sequencing data. Here we report our approach for improving INDEL calling accuracy by using a machine learning algorithm to combine call sets generated with three independent methods, and by leveraging the strengths of each individual pipeline. Utilizing this approach, we generated a consensus exome INDEL call set from a large dataset generated by the 1000 Genomes Project (1000G), maximizing both the sensitivity and the specificity of the calls.

Results

This consensus exome INDEL call set features 7,210 INDELs, from 1,128 individuals across 13 populations included in the 1000 Genomes Phase 1 dataset, with a false discovery rate (FDR) of about 7.0%.

Conclusions

In our study we further characterize the patterns and distributions of these exonic INDELs with respect to density, allele length, and site frequency spectrum, as well as the potential mutagenic mechanisms of coding INDELs in humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1333-7) contains supplementary material, which is available to authorized users.     

Highly Conserved Non-Coding Sequences and the 18q Critical Region for Short Stature: A Common Mechanism of Disease?

Background

Isolated growth hormone deficiency (IGHD) and multiple pituitary hormone deficiency (MPHD) are heterogeneous disorders with several different etiologies and they are responsible for most cases of short stature. Mutations in different genes have been identified but still many patients did not present mutations in any of the known genes. Chromosomal rearrangements may also be involved in short stature and, among others, deletions of 18q23 defined a critical region for the disorder. No gene was yet identified.

Methodology/Principal Findings

We now report a balanced translocation X;18 in a patient presenting a breakpoint in 18q23 that was surprisingly mapped about 500 Kb distal from the short stature critical region. It separated from the flanking SALL3 gene a region enriched in highly conserved non-coding elements (HCNE) that appeared to be regulatory sequences, active as enhancers or silencers during embryonic development.

Conclusion

We propose that, during pituitary development, the 18q rearrangement may alter expression of 18q genes or of X chromosome genes that are translocated next to the HCNEs. Alteration of expression of developmentally regulated genes by translocation of HCNEs may represent a common mechanism for disorders associated to isolated chromosomal rearrangements.