Correlated Evolution of Nucleotide Positions within Splice Sites in Mammals

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Conserved and divergent signals in 5' splice site sequences across fungi, metazoa and plants

In eukaryotic organisms the ensemble of 5' splice site sequences reflects the balance between natural nucleotide variability and minimal molecular constraints necessary to ensure splicing fidelity. This compromise shapes the underlying statistical patterns in the composition of donor splice site sequences. The scope of this study was to mine conserved and divergent signals in the composition of 5' splice site sequences. Because 5' donor sequences are a major cue for proper recognition of splice sites, we reasoned that statistical regularities in their composition could reflect the biological functionality and evolutionary history associated with splicing mechanisms. Results: We considered a regularized maximum entropy modeling framework to mine for non-trivial two-site correlations in donor sequence datasets corresponding to 30 different eukaryotes. For each analyzed species, we identified minimal sets of two-site coupling patterns that were able to replicate, at a given regularization level, the observed one-site and two-site frequencies in donor sequences. By performing a systematic and comparative analysis of 5'splice sites we showed that lineage information could be traced from joint di-nucleotide probabilities. We were able to identify characteristic two-site coupling patterns for plants and animals, and propose that they may echo differences in splicing regulation previously reported between these groups.

Deep Splicer: A CNN Model for Splice Site Prediction in Genetic Sequences

Many living organisms have DNA in their cells that is responsible for their biological features. DNA is an organic molecule of two complementary strands of four different nucleotides wound up in a double helix. These nucleotides are adenine (A), thymine (T), guanine (G), and cytosine (C). Genes are DNA sequences containing the information to synthesize proteins. The genes of higher eukaryotic organisms contain coding sequences, known as exons and non-coding sequences, known as introns, which are removed on splice sites after the DNA is transcribed into RNA. Genome annotation is the process of identifying the location of coding regions and determining their function. This process is fundamental for understanding gene structure; however, it is time-consuming and expensive when done by biochemical methods. With technological advances, splice site detection can be done computationally. Although various software tools have been developed to predict splice sites, they need to improve accuracy and reduce false-positive rates. The main goal of this research was to generate Deep Splicer, a deep learning model to identify splice sites in the genomes of humans and other species. This model has good performance metrics and a lower false-positive rate than the currently existing tools. Deep Splicer achieved an accuracy between 93.55% and 99.66% on the genetic sequences of different organisms, while Splice2Deep, another splice site detection tool, had an accuracy between 90.52% and 98.08%. Splice2Deep surpassed Deep Splicer on the accuracy obtained after evaluating C. elegans genomic sequences (97.88% vs. 93.62%) and A. thaliana (95.40% vs. 94.93%); however, Deep Splicer's accuracy was better for H. sapiens (98.94% vs. 97.15%) and D. melanogaster (97.14% vs. 92.30%). The rate of false positives was 0.11% for human genetic sequences and 0.25% for other species' genetic sequences. Another splice prediction tool, Splice Finder, had between 1% and 3% of false positives for human sequences, while other species' sequences had around 4% and 10%.

Alternative splicing patterns of doublesex reveal a missing link between Nix and doublesex in the sex determination cascade of Aedes albopictus

Sexual development in insects is regulated by a complicated hierarchical cascade of sex determination. The primary signals are diverse, whereas the central nexus doublesex (dsx) gene is relatively conserved within the pathway. Aedes (Stegomyia) albopictus is an important vector with an extensive worldwide distribution. We previously reported that Ae. albopictus dsx (Aalbdsx) yields one male- (Aalbdsx^M ) and three female-specific isoforms (Aalbdsx^F1-3 ); however, the spatiotemporal expression profiles and mechanisms regulating sex-specific alternative splicing require further investigation. In this study, we demonstrated that the Aalbdsx^M messenger RNA (mRNA) represents the default pattern when analyzed in human foreskin fibroblasts and HeLa cells. We combined reverse transcription polymerase chain reaction with RNA immunoprecipitation using specific antibodies against tagged Ae. albopictus male-determining factor AalNix and confirmed that AalNix indirectly regulates dsx pre-mRNA and regulates its alternative splicing. During the early embryo stage (0-2 and 4-8 h), maternal dsx^F and default splicing dsx^M were detected in both sexes; the expression of dsx^M then decreased until sufficient AalNix transcripts accumulated in male embryos at 20-24 h. These findings suggest that one or more potential dsx splicing enhancers can shift dsx^M to dsx^F in both sexes; however, the presence of Nix influences the function of this unknown splicing enhancer and ultimately leads to the formation of dsx^M in males. Finally, our results provide important insight into the regulatory mechanism of dsx alternative splicing in the mosquito.

Skeletal and molecular findings in 51 Cleidocranial dysplasia patients from Turkey

Loss or decrease of function in runt-related transcription factor 2 encoded by RUNX2 is known to cause a rare autosomal-dominant skeletal disorder, cleidocranial dysplasia (CCD). Clinical spectrum and genetic findings in 51 CCD patients from 30 unrelated families are herein presented. In a majority of the patients, facial abnormalities, such as delayed fontanel closure (89%), parietal and frontal bossing (80%), metopic groove (77%), midface hypoplasia (94%), and abnormal mobility of shoulders (90%), were recorded following clinical examination. In approximately one-half of the subjects, wormian bone (51%), short stature (43%), bell-shaped thorax (42%), wide pubic symphysis (50%), hypoplastic iliac wing (59%), and chef's hat sign (44%) presented in available radiological examinations. Scoliosis was identified in 28% of the patients. Investigation of RUNX2 revealed small sequence alterations in 90% and gross deletions in 10% of the patients; collectively, 23 variants including 11 novel changes (c.29_30insT, c.203delAinsCG, c.423 + 2delT, c.443_454delTACCAGATGGGAinsG, c.505C > T, c.594_595delCTinsG, c.636_637insC, c.685 + 5G > A, c.1088G > T, c.1281delC, Exon 6-9 deletion) presented high allelic heterogeneity. Novel c.29_30insT is unique in affecting the P1-driven long isoform of RUNX2, which is expected to disrupt the N-terminal region of RUNX2; this was shown in two unrelated phenotypically discordant patients. The clinical findings highlighted mild intra-familial genotype-phenotype correlation in our CCD cohort.

An extended catalogue of tandem alternative splice sites in human tissue transcriptomes

Tandem alternative splice sites (TASS) is a special class of alternative splicing events that are characterized by a close tandem arrangement of splice sites. Most TASS lack functional characterization and are believed to arise from splicing noise. Based on the RNA-seq data from the Genotype Tissue Expression project, we present an extended catalogue of TASS in healthy human tissues and analyze their tissue-specific expression. The expression of TASS is usually dominated by one major splice site (maSS), while the expression of minor splice sites (miSS) is at least an order of magnitude lower. Among 46k miSS with sufficient read support, 9k (20%) are significantly expressed above the expected noise level, and among them 2.5k are expressed tissue-specifically. We found significant correlations between tissue-specific expression of RNA-binding proteins (RBP), tissue-specific expression of miSS, and miSS response to RBP inactivation by shRNA. In combination with RBP profiling by eCLIP, this allowed prediction of novel cases of tissue-specific splicing regulation including a miSS in QKI mRNA that is likely regulated by PTBP1. The analysis of human primary cell transcriptomes suggested that both tissue-specific and cell-type-specific factors contribute to the regulation of miSS expression. More than 20% of tissue-specific miSS affect structured protein regions and may adjust protein-protein interactions or modify the stability of the protein core. The significantly expressed miSS evolve under the same selection pressure as maSS, while other miSS lack signatures of evolutionary selection and conservation. Using mixture models, we estimated that not more than 15% of maSS and not more than 54% of tissue-specific miSS are noisy, while the proportion of noisy splice sites among non-significantly expressed miSS is above 63%.

Pre-mRNA splicing is an important step in the processing of the genomic information during gene expression. During splicing, introns are excised from a gene transcript, and the remaining exons are ligated. Our work concerns one its particular subtype, which involves the so-called tandem alternative splice sites, a group of closely located exon borders that are used alternatively. We analyzed RNA-seq measurements of gene expression provided by the Genotype-Tissue Expression (GTEx) project, the largest to-date collection of such measurements in healthy human tissues, and constructed a detailed catalogue of tandem alternative splice sites. Within this catalogue, we characterized patterns of tissue-specific expression, regulation, impact on protein structure, and evolutionary selection acting on tandem alternative splice sites. In a number of genes, we predicted regulatory mechanisms that could be responsible for choosing one of many tandem alternative splice sites. The results of this study provide an invaluable resource for molecular biologists studying alternative splicing.

Nucleotide spacing distribution analysis for human genome

The distribution of nucleotides spacing in human genome was investigated. An analysis of the frequency of occurrence in the human genome of different sequence lengths flanked by one type of nucleotide was carried out showing that the distribution has no self-similar (fractal) structure. The results nevertheless revealed several characteristic features: (i) the distribution for short-range spacing is quite similar to the purely stochastic sequences; (ii) the distribution for long-range spacing essentially deviates from the random sequence distribution, showing strong long-range correlations; (iii) the differences between (A, T) and (C, G) nucleotides are quite significant; (iv) the spacing distribution displays tiny oscillations.