Abstract:BackgroundFrameshift translation is an important phenomenon that contributes to the appearance of novel coding DNA sequences (CDS) and functions in gene evolution, by allowing alternative amino acid translations of gene coding regions. Frameshift translations can be identified by aligning two CDS, from a same gene or from homologous genes, while accounting for their codon structure. Two main classes of algorithms have been proposed to solve the problem of aligning CDS, either by amino acid sequence alignment back-translation, or by simultaneously accounting for the nucleotide and amino acid levels. The former does not allow to account for frameshift translations and up to now, the latter exclusively accounts for frameshift translation initiation, not considering the length of the translation disruption caused by a frameshift.ResultsWe introduce a new scoring scheme with an algorithm for the pairwise alignment of CDS accounting for frameshift translation initiation and length, while simultaneously considering nucleotide and amino acid sequences. The main specificity of the scoring scheme is the introduction of a penalty cost accounting for frameshift extension length to compute an adequate similarity score for a CDS alignment. The second specificity of the model is that the search space of the problem solved is the set of all feasible alignments between two CDS. Previous approaches have considered restricted search space or additional constraints on the decomposition of an alignment into length-3 sub-alignments. The algorithm described in this paper has the same asymptotic time complexity as the classical Needleman–Wunsch algorithm.ConclusionsWe compare the method to other CDS alignment methods based on an application to the comparison of pairs of CDS from homologous human, mouse and cow genes of ten mammalian gene families from the Ensembl-Compara database. The results show that our method is particularly robust to parameter changes as compared to existing methods. It also appears to be a good compromise, performing well both in the presence and absence of frameshift translations. An implementation of the method is available at https://github.com/UdeS-CoBIUS/FsePSA.

Pro-Frame: similarity-based gene recognition in eukaryotic DNA sequences with errors

Gene Prediction by the Noise-Assisted MEMD and Wavelet Transform for Identifying the Protein Coding Regions

Aligning coding sequences with frameshift extension penalties

A frame-based representation of genomic sequences for removing errors and rare variant detection in NGS data

Back-translation for discovering distant protein homologies

Alignment of protein-coding sequences with frameshift extension penalties

Using structural and evolutionary information to detect and correct pyrosequencing errors in non-coding RNAs

Interpreting Sequence-Levenshtein distance for determining error type and frequency between two embedded sequences of equal length

Flawed machine-learning confounds coding sequence annotation

GeneScissors: a comprehensive approach to detecting and correcting spurious transcriptome inference owing to RNA-seq reads misalignment.

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

Signal Processing in Sequence Analysis: Advances in Eukaryotic Gene Prediction

Comparative Gene Prediction Based on Gene Structure Conservation.

GeneMarkS-2 : Raising Standards of Accuracy in Gene Recognition

Codon-Based Sequence Alignment for Mutation Analysis by High-Throughput Sequencing

A generalised protein identification method for novel and diverse sequencing technologies

Comparative Exon Prediction Based on Heuristic Coding Region Alignment.

A generalized protein identification method for novel and diverse sequencing technologies

GridWiseTM: The Benefits of a Transformed Energy System

RNA-Seq Mapping Errors When Using Incomplete Reference Transcriptomes of Vertebrates

Fundamental Bounds and Approaches to Sequence Reconstruction from Nanopore Sequencers