Summary: Using an appropriate model of amino acid replacement is very important for the study of protein evolution and phylogenetic inference. We have built a tool for the selection of the best-fit model of evolution, among a set of candidate models, for a given protein sequence alignment. Availability: ProtTest is available under the GNU license from http://darwin.uvigo.es Contact:fabascal@uvigo.es

We present TranslatorX, a web server designed to align protein-coding nucleotide sequences based on their corresponding amino acid translations. Many comparisons between biological sequences (nucleic acids and proteins) involve the construction of multiple alignments. Alignments represent a statement regarding the homology between individual nucleotides or amino acids within homologous genes. As protein-coding DNA sequences evolve as triplets of nucleotides (codons) and it is known that sequence similarity degrades more rapidly at the DNA than at the amino acid level, alignments are generally more accurate when based on amino acids than on their corresponding nucleotides. TranslatorX novelties include: (i) use of all documented genetic codes and the possibility of assigning different genetic codes for each sequence; (ii) a battery of different multiple alignment programs; (iii) translation of ambiguous codons when possible; (iv) an innovative criterion to clean nucleotide alignments with GBlocks based on protein information; and (v) a rich output, including Jalview-powered graphical visualization of the alignments, codon-based alignments coloured according to the corresponding amino acids, measures of compositional bias and first, second and third codon position specific alignments. The TranslatorX server is freely available at http://translatorx.co.uk.

A large number of studies in evolutionary biology utilize phylogenetic information obtained from mitochondrial DNA. Researchers place trust in this molecule and expect it generally to be a reliable marker for addressing questions ranging from population genetics to phylogenies among distantly related lineages. Yet, regardless of the phylogenetic method and weighting treatment, individual mitochondrial genes might potentially produce misleading evolutionary inferences and hence might not constitute an adequate representation neither of the entire mitochondrial genome nor of the evolutionary history of the organisms from which they are derived. We investigated the performance of all mitochondrial protein-coding genes to recover two expected phylogenies of tetrapods and mammals. According to these tests, mitochondrial protein-coding genes can be roughly classified into three groups of good (ND4, ND5, ND2, cytb, and COI), medium (COII, COIII, ND1, and ND6), and poor (ATPase 6, ND3, ATPase 8, and ND4L) phylogenetic performers in recovering these expected trees among phylogenetically distant relatives. How general our findings are is unclear. Simple length differences and rate differences between these genes cannot account for their different phylogenetic performance. The phylogenetic performance of these mitochondrial genes might depend on various factors that play a role in determining the probability of discovering the correct phylogeny such as the density of lineage creation events in time, the phylogenetic "depth" of the question, lineage-specific rate heterogeneity, and the completeness of taxa representation.

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication–random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate “WANCY” region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.