DNAskew: statistical analysis of base compositional asymmetry and prediction of replication boundaries in the genome sequences.
Sueoka and Lobry declared respectively that, in the absence of bias between the two DNA strands for mutation and selection, the base composition within each strand should be A=T and C=G (this state is called Parity Rule type 2, PR2). However, the genome sequences of many bacteria, vertebrates and viruses showed asymmetries in base composition and gene direction. To determine the relationship of base composition skews with replication orientation, gene function, codon usage biases and phylogenetic evolution, in this paper a program called DNAskew was developed for the statistical analysis of strand asymmetry and codon composition bias in the DNA sequence. In addition, the program can also be used to predict the replication boundaries of genome sequences. The method builds on the fact that there are compositional asymmetries between the leading and the lagging strand for replication. DNAskew was written in Perl script language and implemented on the LINUX operating system. It works quickly with annotated or unannotated sequences in GBFF (GenBank flatfile) or fasta format. The source code is freely available for academic use at http://www.epizooty.com/pub/stat/DNAskew.
Ma
Xiang-Ru
XR
Laboratory of Animal Virology, College of Veterinary Medicine, Huazhong Agricultural University, Wuhan 430070, China. hzauvet@public.wh.hb.cn
Xiao
Shao-Bo
SB
Guo
Ai-Zhen
AZ
Lv
Jian-Qing
JQ
Chen
Huan-Chun
HC
D000465:Algorithms; D003062:Codon / Q000235:genetics*; D004261:DNA Replication / Q000235:genetics*; D004269:DNA, Bacterial / Q000235:genetics; D020869:Gene Expression Profiling / Q000379:methods*; D023281:Genomics / Q000379:methods; D008957:Models, Genetic; D015233:Models, Statistical; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*; D012689:Sequence Homology, Nucleic Acid
D003062:Codon; D004269:DNA, Bacterial