VESPA: Very large-scale Evolutionary and Selective Pressure Analyses

Abstract

Background. Large-scale molecular evolutionary analyses of protein coding sequences requires a number of preparatory inter-related steps from finding gene families, to generating alignments and phylogenetic trees and assessing selective pressure variation. Each phase of these analyses can represent significant challenges, particularly when working with entire proteomes (all protein coding sequences in a genome) from a large number of species. Methods. We present VESPA, software capable of automating a selective pressure analysis using codeML in addition to the preparatory analyses and summary statistics. VESPA is written in python and Perl and is designed to run within a UNIX environment. Results. We have benchmarked VESPA and our results show that the method is consistent, performs well on both large scale and smaller scale datasets, and produces results in line with previously published datasets. Discussion. Large-scale gene family identification, sequence alignment, and phylogeny reconstruction are all important aspects of large-scale molecular evolutionary analyses. VESPA provides flexible software for simplifying these processes along with downstream selective pressure variation analyses. The software automatically interprets results from codeML and produces simplified summary files to assist the user in better understanding the results. VESPA may be found at the following website: http://www.mol-evol.org/VESPA.

Comparative genomics provides new insights into the remarkable adaptations of the African wild dog (Lycaon pictus)

“…To further investigate candidate genes detected to be under positive selection, protein structure information from the relevant gene available in the literature was used. Nucleotide sequences of orthologous genes in FASTA format were translated into amino acids using “vespa.py translate” from the VEPSA tool 95 . The AWD protein sequences were then aligned to human annotated versions of the orthologous genes within Geneious v11.1.1 96 to identify amino acid changes that occur within a given protein domain (see Fig.…”

Comparative Genomics of Synechococcus elongatus Explains the Phenotypic Diversity of the Strains

“…Ortholog assignments from Mauve were further refined using Pfam and COG category analysis in eggNOG-mapper ( 77 , 78 ). These ortholog assignments were checked against homology groups created through reciprocal nucleotide BLAST search using Vespa ( 79 ). Homology groups were translated and realigned using MAFFT ( 80 ) in Geneious Prime, and the percentages of identical residues for the nucleotide and amino acid alignments were reported.…”

Improving orthologous signal and model fit in datasets addressing the root of the animal phylogeny

“…A pair of alignments with a MetAl score of <0.15 were judged to be in mutual agreement, whereas a pair with a MetAl score of >0.15 were judged to be mutually discordant. 3 If any pair of alignments were judged to be mutually discordant by MetAl, the columnbased normalized mean distance of amino acid similarity across all sites in all three alignments was calculated using norMD, and the alignment with the highest similarity score was selected for orthogroup tree reconstruction (Thompson et al, 2001;Muller et al, 2010;Webb et al, 2017). Otherwise if all three alignment methods were in mutual agreement, the best-fit alignment was randomly selected.…”