NucAmino: a nucleotide to amino acid alignment optimized for virus gene sequences
- PMID: 28249562
- PMCID: PMC5333393
- DOI: 10.1186/s12859-017-1555-6
NucAmino: a nucleotide to amino acid alignment optimized for virus gene sequences
Abstract
Background: Current nucleotide-to-amino acid alignment software programs were developed primarily for detecting gene exons within eukaryotic genomes and were therefore optimized for speed across long genetic sequences. We developed a nucleotide-to-amino acid alignment program NucAmino optimized for virus sequencing.
Results: NucAmino is an open source program written in the high-level language Go. NucAmino is more likely to align codons flush with a reference sequence's amino acids and can be modified to facilitate the placement of insertions and deletions at specific positions. We compared NucAmino to the nucleotide to amino acid alignment program Local Alignment Program (LAP) using 115,118 human immunodeficiency virus type 1 (HIV-1) protease, reverse transcriptase, and integrase sequences-three genes that are commonly sequenced in clinical laboratories. Discordances between NucAmino and LAP occurred in 512 (16.9%) of the 3,029 sequences containing gaps but in none of 112,910 sequences without gaps. For 242 of the sequences with discordances, NucAmino produced an alignment that was preferable to that found by LAP in that it was more likely to codon align insertions and deletions and to facilitate the placement of an important drug-resistance associated insertion at the position at which most laboratories expect it to occur.
Conclusions: NucAmino is a nucleotide-to-amino acid alignment program with several advantages for clinical laboratories performing virus sequencing compared with older programs designed for gene finding.
Keywords: Drug resistance; HIV-1; Open source; Sequence alignment; Viruses.
Figures
References
-
- Huang X, Zhang J. Methods for comparing a DNA sequence with a protein sequence. Comput Appl Biosci. 1996;12(6):497–506. - PubMed
-
- Winters MA, Coolley KL, Girard YA, Levee DJ, Hamdan H, Shafer RW, Katzenstein DA, Merigan TC. A 6-basepair insert in the reverse transcriptase gene of human immunodeficiency virus type 1 confers resistance to multiple nucleoside inhibitors. J Clin Invest. 1998;102(10):1769–1775. doi: 10.1172/JCI4948. - DOI - PMC - PubMed
-
- Masquelier B, Race E, Tamalet C, Descamps D, Izopet J, Buffet-Janvresse C, Ruffault A, Mohammed AS, Cottalorda J, Schmuck A, et al. Genotypic and phenotypic resistance patterns of human immunodeficiency virus type 1 variants with insertions or deletions in the reverse transcriptase (RT): multicenter study of patients treated with RT inhibitors. Antimicrob Agents Chemother. 2001;45(6):1836–1842. doi: 10.1128/AAC.45.6.1836-1842.2001. - DOI - PMC - PubMed
MeSH terms
Substances
Grants and funding
LinkOut - more resources
Full Text Sources
Other Literature Sources
Research Materials
