A model of k-mer surprisal to quantify local sequence information content surrounding splice regions
AffiliationCRUK Manchester Institute Cancer Biomarker Centre, The University of Manchester, Mancheste
MetadataShow full item record
AbstractMolecular sequences carry information. Analysis of sequence conservation between homologous loci is a proven approach with which to explore the information content of molecular sequences. This is often done using multiple sequence alignments to support comparisons between homologous loci. These methods therefore rely on sufficient underlying sequence similarity with which to construct a representative alignment. Here we describe a method using a formal metric of information, surprisal, to analyse biological sub-sequences without alignment constraints. We applied our model to the genomes of five different species to reveal similar patterns across a panel of eukaryotes. As the surprisal of a sub-sequence is inversely proportional to its occurrence within the genome, the optimal size of the sub-sequences was selected for each species under consideration. With the model optimized, we found a strong correlation between surprisal and CG dinucleotide usage. The utility of our model was tested by examining the sequences of genes known to undergo splicing. We demonstrate that our model can identify biological features of interest such as known donor and acceptor sites. Analysis across all annotated coding exon junctions in Homo sapiens reveals the information content of coding exons to be greater than the surrounding intron regions, a consequence of increased suppression of the CG dinucleotide in intronic space. Sequences within coding regions proximal to exon junctions exhibited novel patterns within DNA and coding mRNA that are not a function of the encoded amino acid sequence. Our findings are consistent with the presence of secondary information encoding features such as DNA and RNA binding sites, multiplexed through the coding sequence and independent of the information required to define the corresponding amino-acid sequence. We conclude that surprisal provides a complementary methodology with which to locate regions of interest in the genome, particularly in situations that lack an appropriate multiple sequence alignment.
CitationHumphrey S, Kerr A, Rattray M, Dive C, Miller CJ. A model of k-mer surprisal to quantify local sequence information content surrounding splice regions. PeerJ. 2020;8:e10063.
- [Analysis, identification and correction of some errors of model refseqs appeared in NCBI Human Gene Database by in silico cloning and experimental verification of novel human genes].
- Authors: Zhang DL, Ji L, Li YD
- Issue date: 2004 May
- Prediction of mutant mRNA splice isoforms by information theory-based exon definition.
- Authors: Mucaki EJ, Shirley BC, Rogan PK
- Issue date: 2013 Apr
- Compensatory relationship between splice sites and exonic splicing signals depending on the length of vertebrate introns.
- Authors: Dewey CN, Rogozin IB, Koonin EV
- Issue date: 2006 Dec 8
- Finding exonic islands in a sea of non-coding sequence: splicing related constraints on protein composition and evolution are common in intron-rich genomes.
- Authors: Warnecke T, Parmley JL, Hurst LD
- Issue date: 2008
- SplicedFamAlign: CDS-to-gene spliced alignment and identification of transcript orthology groups.
- Authors: Jammali S, Aguilar JD, Kuitche E, Ouangraoua A
- Issue date: 2019 Mar 29