Abstract
In gametophytic self-incompatibility systems, many specificities (different 'lock-and-key' combinations) are maintained by frequency-dependent selection for very long evolutionary times. In Solanaceae, trans-specific evolution (the observation that an allele from one species may be more closely related to an allele from another species than to others from the same species) has been taken as an argument for the very old age of specificities. In this work, by determining, for the first time, the age of extant Prunus species, we show that this reasoning cannot be applied to Prunoideae. Furthermore, since our sample size is large (all S-RNase encoding the female component and SFB encoding the male component GenBank sequences), we were able to estimate the age of the oldest Prunus specificities. By doing so, we show that the lower variability levels at the Prunus S-locus, in comparison with Solanaceae, is due to the younger age of Prunus alleles, and not to a difference in silent mutation rates. We show that the ancestor to extant Prunus species harboured at least 102 specificities, in contrast to the maximum of 33 observed in extant Prunus species. Since the number of specificities that can be maintained in a population depends on the effective population size, this observation suggests a bottleneck in Prunus evolutionary history. Loss of specificities may have occurred during this event. Using only information on amino acid sites that determine specificity differences, and a simulation approach, we show that a model that assumes closely related specificities are not preferentially lost during evolution, fails to predict the observed degree of specificity relatedness.

Abstract
Maximum-likelihood methods based on models of codon substitution have been widely used to infer positively selected amino acid sites that are responsible for adaptive changes. Nevertheless, in order to use such an approach, software applications are required to align protein and DNA sequences, infer a phylogenetic tree and run the maximum-likelihood models. Therefore, a significant effort is made in order to prepare input files for the different software applications and in the analysis of the output of every analysis. In this paper we present the ADOPS (Automatic Detection Of Positively Selected Sites) software. It was developed with the goal of providing an automatic and flexible tool for detecting positively selected sites given a set of unaligned nucleotide sequence data. An example of the usefulness of such a pipeline is given by showing, under different conditions, positively selected amino acid sites in a set of 54 Coffea putative S-RNase sequences. ADOPS software is freely available and can be downloaded from http://sing.ei.uvigo.es/ADOPS.

Abstract
In this paper we present the ADOPS (Automatic Detection Of Positively Selected Sites) software that is ideal for research projects involving the analysis of tens of genes. ADOPS is a novel software pipeline that is being implemented with the goal of providing an automatic and flexible tool for detecting positively selected sites given a set of unaligned nucleotide sequence data.

Abstract
Chromosomal inversions can originate from breakage and repair by non-homologous end-joining. Nevertheless, they can also originate from ectopic recombination between transposable elements located on the same chromosome inserted in opposite orientations. Here, we show that a MITE element (DAIBAM), previously involved in the origin of one Drosophila americana polymorphic inversion, is also involved in the origin of one fixed inversion between D. virilis and D. americana and another D. americana polymorphic inversion. Therefore, DAIBAM is responsible for at least 20% of the chromosomal rearrangements that are observed within and between species of the virilis phylad (D. virilis, D. lummei, D. novamexicana and D. americana), having thus played a significant role in the chromosomal evolution of this group of closely related species.

Abstract
Background: Comparative studies using hundreds of sequences can give a detailed picture of the evolution of a given gene family. Nevertheless, retrieving only the sequences of interest from public databases can be difficult, in particular, when working with highly divergent sequences. The difficulty increases substantially when one wants to include in the study sequences from many (or less well studied) species whose genomes are non-annotated or incompletely annotated. Methodology/Principal Findings: In this work we evaluate the usefulness of different approaches of gene retrieval and classification, using the distal-less (DLX) gene family as a test case. Furthermore, we evaluate whether the use of a large number of gene sequences from a wide range of animal species, the use of multiple alternative alignments, and the use of amino acids aligned with high confidence only, is enough to recover the accepted DLX evolutionary history. Conclusions/Significance: The canonical DLX homeobox gene sequence here derived, together with the characteristic amino acid variants here identified in the DLX homeodomain region, can be used to retrieve and classify DLX genes in a simple and efficient way. A program is made available that allows the easy retrieval of synteny information that can be used to classify gene sequences. Maximum likelihood trees using hundreds of sequences can be used for gene identification. Nevertheless, for the DLX case, the proposed DLX evolutionary is not recovered even when multiple alignment algorithms are used.

Abstract
Background: Although homeobox genes have been the subject of many studies, little is known about the main amino acid changes that occurred early in the evolution of genes belonging to different classes. Results: In this study, we report a method for the fast and efficient retrieval of sequences belonging to the ANTP (HOXL and NKL) and PRD classes. Furthermore, we look for diagnostic amino acid residues that can be used to distinguish HOXL, NKL and PRD genes. Conclusion: The reported protein features will facilitate the robust classification of homeobox genes from newly sequenced bilaterian genomes. Nevertheless, in non-bilaterian genomes our findings must be cautiously applied. In principle, as long as a good manually curated data set is available the approach here described can be applied to non-bilaterian organisms as well. Our results help focus experimental studies onto investigating the biochemical functions of key homeodomain residues in different gene classes.