Abstract:Alternative splicing (AS), in higher eukaryotes, is one of the mechanisms of post-transcriptional regulation that generate multiple transcripts from the same gene. One particular mode of AS is the skipping event where an exon may be alternatively excluded or constitutively included in the resulting mature mRNA. Both transcript isoforms from this skipping event site, i.e. in which the exon is either included (inclusion isoform) or excluded (skipping isoform), are typically present in one cell, and maintain a subtle balance that is vital to cellular function and dynamics. However, how the prevailing conditions dictate which isoform is expressed and what biological factors might influence the regulation of this process remain areas requiring further exploration. In this study, we have developed a novel computational method, graph-based exon-skipping scanner (GESS), for de novo detection of skipping event sites from raw RNA-seq reads without prior knowledge of gene annotations, as well as for determining the dominant isoform generated from such sites. We have applied our method to publicly available RNA-seq data in GM12878 and K562 cells from the ENCODE consortium and experimentally validated several skipping site predictions by RT-PCR. Furthermore, we integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. Our computational analysis found that splice sites within the skipping-isoform-dominated group (SIDG) tended to exhibit weaker MaxEntScan-calculated splice site strength around middle, ‘skipping’, exons compared to those in the inclusion-isoform-dominated group (IIDG). We further showed the positional preference pattern of splicing factors, characterized by enrichment in the intronic splice sites immediately bordering middle exons. Finally, our analysis suggested that different epigenetic factors may introduce a variable obstacle in the process of exon–intron boundary establishment leading to skipping events.

Comparative analysis of the nucleotide composition biases in exons and introns of human genes

Relationships among the nucleotide content of human genome sequence, gene structure, and gene expression features (PhD synopsis)

Statistical Properties of Nucleotides in Human Chromosomes 21 and 22

The Contributions of Replication Orientation, Gene Direction, and Signal Sequences to Base-Composition Asymmetries in Bacterial Genomes

Comparing Synonymous Codon Usage of Alternatively Spliced Genes with Non-Alternatively Spliced Genes in Human Genome

New Views On Strand Asymmetry In Insect Mitochondrial Genomes

Lower GC-content in Editing Exons: Implications for Regulation by Molecular Characteristics Maintained by Selection

Codon bias as a factor in regulating expression via translation rate in the human genome

Selection on codon bias in yeast: a transcriptional hypothesis

EXPLORING THE DISTRIBUTION OF SINGLE NUCLEOTIDE POLYMORPHISMS ACROSS HUMAN EXONS AND INTRONS

Human Exonization Through Differential Nucleosome Occupancy

Intragenic spatial patterns of codon usage bias in prokaryotic and eukaryotic genomes

A-to-I Editing Sites Are a Genomically Encoded G: Implications for the Evolutionary Significance and Identification of Novel Editing Sites

Local DNA dynamics shape mutational patterns of mononucleotide repeats in human genomes

From DNA sequence analysis to modeling replication in the human genome

Relationships among Stop Codon Usage Bias, Its Context, Isochores, and Gene Expression Level in Various Eukaryotes

Biased exonization of transposed elements in duplicated genes: A lesson from the TIF-IA gene

Computational Analysis Reveals a Correlation of Exon-Skipping Events with Splicing, Transcription and Epigenetic Factors.

Codon Usage Bias: An Endless Tale

The GC-content at the 5’ends of human protein-coding genes is undergoing mutational decay

Synonymous but not the same: the causes and consequences of codon bias