Communicated by José Sarukhán, National Autonomous University of Mexico, Mexico D.F., Mexico, September 11, 2003 (received for review May 31, 2003)
Gene duplication is a substrate of evolution. However, the relative importance of positive selection versus relaxation of constraints in the functional divergence of gene copies is still under debate. Plant MADS-box genes encode transcriptional regulators key in various aspects of development and have undergone extensive duplications to form a large family. We recovered 104 MADS sequences from the Arabidopsis genome. Bayesian phylogenetic trees recover type II lineage as a monophyletic group and resolve a branching sequence of monophyletic groups within this lineage. The type I lineage is comprised of several divergent groups. However, contrasting gene structure and patterns of chromosomal distribution between type I and II sequences suggest that they had different evolutionary histories and support the placement of the root of the gene family between these two groups. Site-specific and site-branch analyses of positive Darwinian selection (PDS) suggest that different selection regimes could have affected the evolution of these lineages. We found evidence for PDS along the branch leading to flowering time genes that have a direct impact on plant fitness. Sites with high probabilities of having been under PDS were found in the MADS and K domains, suggesting that these played important roles in the acquisition of novel functions during MADS-box diversification. Detected sites are targets for further experimental analyses. We argue that adaptive changes in MADS-domain protein sequences have been important for their functional divergence, suggesting that changes within coding regions of transcriptional regulators have influenced phenotypic evolution of plants.
Gene duplication provides a substrate for evolution, and understanding the fate of duplicates is fundamental to clarifying mechanisms of genetic redundancy and the link between gene family diversification and phenotypic evolution (1). Several empirical studies have evaluated the roles of duplication in adaptation and diversification (2), but the evolutionary forces at play during functional divergence of duplicates are still under debate (3, 4). Genomic studies are revealing that eukaryotes harbor large families of genes that have arisen during evolution through duplication and have persisted for longer periods of time than expected by classical models (5). Models that incorporate positive selection (6, 7) provide alternative explanations for the persistence of duplicates.
Empirical studies to test models on the fate of duplicates and the evolutionary forces driving their functional divergence will need complete and resolved gene family phylogenies. Several studies suggest that positive Darwinian selection (PDS) might have been important in protein evolution. However, most previous studies have involved few members of a gene family from various species (8, 9). In this article, we annotate, align, and analyze the complete MADS-box gene family of the plant model system Arabidopsis thaliana and provide resolved phylogenies as a basis to infer the role of PDS in protein evolution in this gene family.
The detection of an excess in the ratio of the rate of nonsynonymous (dN) over the synonymous (dS) substitutions (that is dN/dS > 1; dN/dS is also denoted ω) is a nonambiguous indicator of PDS at the coding sequence level. Early studies estimated this ratio as an average over all codon sites within complete or partial sequence stretches and over the entire evolutionary time that separates the sequences compared. This method appears to be conservative because many sites might be under purifying selection because of functional constraint (10). However, in adaptive evolution in developmental regulatory loci, such as the MADS (11) PDS most likely occurs along particular lineages and at specific sites. In such cases, average dN/dS ratios over time and sites might not be significantly greater than 1, even if PDS has occurred.
MADS-box genes are present in plants, animals, and fungi, and previous studies suggested the existence of two main monophyletic lineages (type I and II) among all eukaryotes that probably derived from at least one duplication event before the divergence of plants and animals (14). The trees presented here recover Arabidopsis type II genes as a strongly supported monophyletic lineage, and the type I genes seem to be monophyletic but comprise several divergent sublineages.
MADS-box genes encode transcriptional regulators with diverse functions that could have been key during important events of plant diversification (12, 13). Hence, phylogenetic analyses of MADS-box genes are useful guides for studying their roles in plant evolution. Plant MADS-box gene phylogeny resolution, especially at its basal nodes, has been hindered by incomplete data and by limitations of inference methods (14, 15). Here we show resolved gene phylogenies of the Arabidopsis MADS-box genes.
More than half of the Arabidopsis MADS-box sequences are type I and only share with type II the MADS-box (14). All but one (16) functionally characterized plant MADS-box genes are type II and encode the three floral homeotic functions of the flower development ABC model (17-19). They also encode regulators of flower initiation, flower meristem identity, and various aspects of ovule, fruit, leaf, and root development (11, 20-23). All characterized plant type II MADS-box genes encode proteins that share a stereotypical MIKC structure, with highly conserved MADS and K domains that are putative DNA-binding and protein-protein interaction domains, respectively, and less conserved I and COOH regions.
We show here that sequences of type I and II have contrasting gene structure and chromosome distribution, supporting the idea that these two lineages had different evolutionary histories with a contrasting role of PDS. These contrasting histories also support placing the root of the family tree between the two lineages. Indeed, we found a significant role of PDS at fixing specific residues within the MADS-domain after different duplications in the type I lineage, but in the type II lineage we found evidence of PDS only with branch-site models along specific lineages. We addressed whether PDS played a significant role during the evolution of genes that regulate the transition to flowering, a trait that is clearly linked to plant fitness. Our findings identify target proteins and residues for future functional analyses and suggest that changes in coding sequences of transcriptional regulators, and not only in their regulatory regions, played important roles during phenotypic evolution.
Sequences and Alignment. To detect putative MADS-box genes in the Arabidopsis genome, two tblastn searches were performed on the complete Arabidopsis database (Table 3, which is published as supporting information on the PNAS web site). Sequences were assigned to either type I or II based on exon number and on careful comparisons of their MADS boxes (14, 24). Type I and II were then aligned separately with clustalw (25) launched from bioedit (26) and hand-corrected by using published alignments as guides (refs. 27 and 28; details can be found in Supporting Materials and Methods, which is published as supporting information on the PNAS web site).
Phylogenetic Reconstruction. Bayesian phylogenetic analyses were performed on mrbayes 2.01 (29). All searches were started from a random tree, on four different Markov chains for 2,500,000 generations and saving every 100th tree. At convergence (≈10,000 generations), the first 15,000 trees were discarded and a consensus was built. Posterior Bayesian probabilities were used to evaluate branch support. According to the recommendations of Foster (www.bioinf.org/molsys/data/like.pdf), we used the GTR model with a substitution rate that varies in an intracodon position-specific manner (GTR + SS).
Statistical Tests for Positive Selection. We applied the approach of Yang and coworkers (9, 30) to test for positive selection. First, we ran a test for the existence of sites with dN/dS ratios >1 by using a likelihood ratio test (LRT) to compare a model that does not allow for sites with dN/dS >1 to a model that does. If the LRT was statistically significant, then we identified the sites that were under positive selection. We calculated the posterior probability (PP) that a site was drawn from a given dN/dS class. Sites with PP > 0.5 are reported but we focus on those with a PP > 0.95. The program HYPHY 0.901b (S. L. Kosakovsky-Pond and S. V. Muse, www.hyphy.org) was used. Models tested were M3, M2, and M8 vs. null models M0, M1, and M7, respectively. We used the codon substitution model of Goldman and Yang (31) and 10 classes in the gamma distribution of M7 and M8. To avoid false positives, sites detected by models M3, M2, and M8 were considered as bona fide results only if the same sites were detected with at least two of these models and in both cases the LRT result was significant (32). To further minimize false positives, we performed all analyses on unambiguous and compact alignments (available on request). Thus, our conclusions on the role of PDS are based on very conservative analyses.
Additionally, we performed the branch-site analyses of Yang and Nielsen by using model MB (PAML 3.13, http://abacus.gene.ucl.ac.uk/software/paml.html; refs. 33 and 34) at the basal branch of two clades of flowering-time type II genes [FLC- and SVP-like genes (35-37)]. We hypothesized that functional change would be important precisely at the origin of these clades that evolved a distinct function (details of procedures can be found in refs. 31, 33, and 34). We also performed site-branch analyses for the branch that leads to type I PHERES1 gene.
Gene conversion and concerted evolution may violate the assumptions of site-specific positive selection models. We used geneconv 1.81 (www.math.wustl.edu/~sawyer/geneconv/index.html) and MEGA 2.1 (38) on the alignments tested for PDS.
Annotation, Nomenclature, Gene Structure, and Duplications. The list of 104 MADS-box sequences found in the Arabidopsis genome database (www.ncbi.nlm.nih.gov/blast/Genome/ara.html), given agamous-like (AGL) number, accession number, chromosome location, intron-exon structure, and type of duplication are shown in Table 3.
Type I and II genes have contrasting gene structure and chromosome location. Type I genes have always one or two exons, whereas type II genes have more than five, and typically six to eight, suggesting that these two types of genes have different predisposition to gain or loss of introns or, alternatively, a difference in exon shuffling in the building of both types of genes, perhaps since their origin. All type II genes have a clear MIKC structure, except AGL33 with a very short transcript and AGL30-related sequences that do not have a clear IKC structure but share conserved motifs in their MADS-boxes with the rest of the type II genes. Also, a coiled-coil domain similar to the K domain is inferred for at least AGL104 (data not shown).
We used chromosome map locations to make qualitative inferences on past duplication events during MADS-box gene family evolution. The distribution of type I sequences among the five chromosomes is distinct to that of type II genes. Whereas the former are concentrated in chromosome I and V (χ2 = 26.77; P < 0.001 rejects uniform distribution among chromosomes correcting with chromosome size), type II are uniformly distributed (χ2 = 2.62; P > 0.1) among chromosomes. Also, most type I genes can be traced to intrachromosomal duplications, whereas approximately half of type II genes seem to have originated from interchromosomal duplications. Interestingly, Lynch and Conery (39) have found that recent duplications happened more frequently within than between two chromosomes, suggesting that most type I genes diverged more recently than type II genes (40). Moreover, a survey of the A. thaliana paralogous blocks database (http://wolfe.gen.tcd.ie/athal/dup) indicates that there are more duplicates from the type II group that seem to have persisted than those from type I. Eighteen out of 26 type II sequences that are found in a nonspurious duplicated region had a close paralog in a sister region, whereas among type I genes only 2 out of 22 did. Among type II genes found in nonspurious duplicated blocks, 14 are found in interchromosomal duplications and 12 in intrachromosomal duplications, whereas among type I genes the corresponding numbers are 15 and 7. But when we consider gene pairs of terminal clades of the trees, 13 out of 18 type I gene pairs involve intrachromosomal duplications, whereas in the type II we found this to be true of 6 out of 15 gene pairs (Fig. 1).
A. thaliana MADS-box gene family Bayesian phylogeny. Numbers above or below branches represent Bayesian posterior probabilities of finding a given clade. Branch lengths are proportional to number of nucleotide substitutions. AGL105 was excluded.
MADS-Box Gene Family Phylogeny. To corroborate the monophyletic origin of the two lineages that we had previously proposed (14), we obtained trees that included 103 of the MADS-box sequences found in the A. thaliana genome (Fig. 1). The global phylogeny recovers the two lineages (types I and II) of MADS-box sequences as two monophyletic groups with both alignments used (see supporting information) if the tree root is placed between type I and II genes. However, type I genes are more divergent among them than type II genes (Tables 4 and 5, which are published as supporting information on the PNAS web site). Nonetheless conserved motifs after the MADS suggest that type I sequences are not pseudogenes, an idea supported by the recent characterization of PHERES1 (16).
AGL30 had been incorrectly assigned to type I. However, this and related genes seem to be divergent type II. This is supported by their affiliation to type II-like moss genes that bear K domains as well as by their exon number (27, 41), and by conserved MADS-box motifs with respect to other type II genes. In our global tree (Fig. 1), these genes are resolved in a different position to that in the type II tree (Fig. 2b). Nonetheless, this remaining ambiguity does not affect PDS analyses. We ran PDS analyses for MIKC genes with an alternative type II topology similar to that in Fig. 1 and results recover the same sites with high PP as those obtained with the topology shown in Fig. 2b (data available on request).
Type I and type II A. thaliana MADS-box Bayesian phylogenies. (a) Type I tree polarized with type II sequences. (b) Type II tree polarized with type I sequences. Numbers above or below branches represent posterior probabilities. Branch lengths are proportional to the number of nucleotide substitutions. Boxed letters identify clades in which site-specific tests of positive selection yielded statistically significant LRT results for at least two model comparisons and in which at least one of the models detected sites under PDS with PP > 0.90. Branches underlined with a broken line identify cases in which the branch-site analyses yielded significant PDS results. #, Excluded from the site-specific analyses.
Given the ambiguity of the alignment of sequences after the MADS-boxes of both lineages, we resolved the internal phylogenetic relationships of each lineage separately. For the type I lineage tree, we used type II sequences as outgroups and vice versa for type II lineage. In both cases, tree topologies were very similar if regions downstream of the MADS-boxes of outgroup sequences were assumed to be homologous to those of ingroup sequences or were displaced. In contrast to previously published phylogenies, the trees shown here for both lineages resolve the branching sequence of the monophyletic groups (Fig. 2).
In the type I lineage, several strongly supported monophyletic clades are resolved and these are confirmed in trees with various outgroups (data not shown). It is noteworthy that no type I sequences bear the IKC region typical of lineage II and that several of the previously identified (14) privative amino acids within the putative MADS-domain of this lineage are found in most available sequences. Two main monophyletic groups are distinguished within the type I, suggesting at least one ancestral duplication within this lineage. DNA sequences beyond the MADS-box in the AGL23-like sublineage are conserved within each small monophyletic group resolved within this sublineage but are very divergent among groups (Fig. 2a and Table 4). Within the AGL26-like genes, two groups are resolved [AGL26-like genes themselves and PHERES1-like genes that includes the first type I gene functionally characterized (16)] and within each a very high degree of conservation is found in the putative domains beyond the MADS (Fig. 2a and Table 4). Sisters to these groups, AGL47 and AGL82, are resolved but have divergent putative domains beyond the MADS.
In the type II lineage, which includes all of the MIKC genes functionally characterized up to now, several clades are resolved and all are well supported (Fig. 2b). It is noteworthy that the AGAMOUS clade has AGL12, which had been reported as root-specific (42), as its sister gene. Another important finding is the strong association of AGL79, expressed in roots (data not shown), with AP1, CAL, and FUL, which are well characterized flower development genes (43). Therefore, this tree suggests that not all monophyletic groups resolved include genes with similar expression patterns and functions as previously thought (15, 21), but formal and robust inferences on evolution of MADS-box gene expression and function will have to await more experimental data and the inclusion of genes from additional taxa.
Positive Selection in MADS-Box Gene Evolution. We compared Models M0 and M3 to evaluate whether there had been dN/dS ratio variation among codon positions below each node of type I and II trees from Fig. 2 (Tables 1 and 2). We found rate variation at deep, intermediate, and recent duplications of both type I and II lineages (data available on request).
Secondly, we applied the LRT to compare data fit to models M1 vs. M2 and M7 vs. M8 to address whether PDS promoted divergence of MADS-box genes below nodes and whether the action of selection has been heterogeneous among protein domains codified by these genes, using the trees of Fig. 2 (Tables 1 and 2 and Fig. 3). Below many of the deep nodes of type II tree, model M3 and at least one of M2 or M8 had significant LRT results (Fig. 2b). However, none of the sites with high PP were detected by more than one model comparison with significant LRT results (data available on request). In contrast, a similar analysis for type I lineage reveals several nodes below which specific sites appear to have been under PDS (Tables 1 and 2). For instance, in nodes AH and AL, positions 72-74 appear to have been under PDS with high PP, with position 72 showing five different amino acids for the seven sequences involved (Fig. 2a and Table 4). Positions 73 and 74 are part of the otherwise highly conserved “RQVTF” motif, and in the human serum response factor, position 72 has been shown to be involved in DNA contact (44). Below nodes AH and AR, position 82 was also found to be under PDS in two of the model pairs compared, although model M3 collapsed to only two rate classes. At this position, amino acid diversity is very high. For example, 10 different residues can be found for 20 sequences analyzed below node AR. In contrast, the homologous position of type II genes (position 26 in Table 4) shows only six different amino acids for 45 sequences.
Schematic representation of the distribution of sites under PDS in type I (a) and type II (b) sequences. MADS, I, K, and COOH domains are indicated. In a, the upper row corresponds to branch-site analysis along the branch leading to PHERES1, and the rest of the rows correspond, from top to bottom, to the site analyses performed below the nodes marked with the boxed letters AG, AH, AL, AP, AQ, and AR in Fig. 2a. In b, the upper row corresponds to branch-site analysis along the FLC branch and the lower row corresponds to branch-site analysis along the SVP branch. In a, the upper scale corresponds to amino acid position along the MADS domain and the lower scale corresponds to amino acid position along our alignment. In b, the scale corresponds to amino acid position along our alignment. Sites with PP > 0.70 are included. Shaded regions in b were excluded from PDS to avoid false positives. All sites are listed in Tables 1 and 2.
Strong evidence for positive selection in type I evolution was also found in less variable positions. This is the case for position 123 of node AG and position 58 of nodes AP and AQ. In position 123, there are only two different amino acids, although their distribution suggests that this site mutated twice during the history of descendants of node AG. In position 58, there are only five variable sites out of 14 sequences compared.
In the above analyses, PDS is detected at individual sites only if the average dN rate across lineages is higher than the average dS rate. This is observed mainly when recurrent positive selection occurs. However, PDS may change a few key residues of a protein but only at particular moments during its evolutionary history. In the latter case, detecting a significantly elevated dN would be hard if an average across-lineages estimate is considered. Thus, we applied the branch-site model of Yang and Nielsen (33) to test for PDS affecting individual sites along the branches leading to the flowering-time FLC- and SVP-like genes in type II lineage and along the branch leading to the only functionally characterized type I protein (PHERES1).
Parameter estimates suggest that in the basal branch of the FLC-like lineage, four residues were fixed by PDS with PP > 0.95 (Tables 1 and 2 and Fig. 3). Two of the residues were found within the MADS and correspond to amino acids that in other MADS-domain proteins participate in interactions between subunits (44-46). For example, site 42 is homologous to a site that in the myocyte enhancer factor-2 has been shown to intervene in subunit folding (45). Another site with high PP is 154. This site is within the K-box and has been shown to intervene in AP3/PI dimerization and determines functional specificity in AP1 and AG (47).
PDS seems to have been important also along the branch leading to the SVP-like genes, and we found two sites that appear to have been fixed by PDS (PP > 0.95). Site 16 is within the MADS-domain, and its homologous position in MEF2A plays a role in DNA-protein complex stabilization (46). Position 144 is found within the K-domain. SVP and AGL24 are the only type II proteins that have a lysine at that position.
Branch-site models also detected strong (PP > 0.95; see Tables 1 and 2) support for PDS along the branch leading to PHE1 at sites 92 and 72, which have been reported (44) to be important for α-β folding, and for position 105, which is involved in both dimerization and α-β folding.
Our analyses suggest that gene conversion or concerted evolution has not been prevalent during MADS-box genes evolution and does not bias PDS inferences. The overall modified Nei-Gojobori (48) means of synonymous-nonsynonymous differences were relatively high: 0.18 and 0.09 (transition/transversion ratio of 1.72 and 1.7, respectively) for types I and II, respectively (49). Gene conversion does not seem to have played significant roles in MADS-box evolution because only three conversion events were detected among type II genes, but none of these included genes for which we infer PDS and no conversion events were detected among type I sequences.
We presented an annotated list of 104 MADS-box sequences from the complete A. thaliana genome database. Our phylogenetic analyses provide a resolved evolutionary hypothesis for the A. thaliana MADS-box gene family. This will be a useful reference for establishing orthology relationships, postulating functional hypotheses for uncharacterized MADS-box genes, and evaluating the role of MADS-box genes in plant morphological evolution.
The monophyly of the type II lineage is strongly supported in the present analyses, and type I comprises several sublineages with divergent putative domains after the MADS. However, previous analyses (14), as well as contrasting exon-intron structure and chromosomal distribution between type I and II sequences, still support the placement of the root between the type I and II genes in the tree of the complete gene family (Fig. 1). This tree hence resolves type I and II sequences in two monophyletic lineages. Nonetheless, genes from other plant, animal, and fungal species should be included in future analyses to trace MADS-box gene duplications with respect to taxa divergence and to specifically reevaluate the number of MADS-box gene duplications that occurred before the divergence of plants and animals (14). Such analyses will provide further evidence to reevaluate the monophyly of type I and II lineages.
Gene family structure has to be understood in the context of extensive gene duplications that have occurred in the evolutionary history of A. thaliana (40). Duplications leading to the chromosome stretches identified in the Arabidopsis Genome Initiative occurred 65 million years ago or before (refs. 50 and 51, but see ref. 52). Most retained groups within these stretches belong to type II genes, and duplications among type II seem to have been more ancient than those among type I, as suggested by their differential distribution among chromosomes. Retention due to a balance between genetic drift and mutation (5, 53) would depend on population characteristics (mainly effective population size) and hence would affect sequences of type I and II equally. But contrasting roles of selection between these two lineages could underlie the contrasting retention rates observed between them. The different evolutionary histories could have been determined by the fact that genes from these two lineages were recruited for different functions.
Interestingly, although duplications of type I seem to have occurred more recently than those leading to the type II lineage, type I sequences are more divergent among them in comparison to type II. This finding would suggest that whereas type II genes have been affected by sporadic PDS at the origin of new functions, followed by strong functional constraint, type I genes have been subject to recurrent events of PDS. In turn, this could be an indication that the functional roles of type I genes are overall distinct to those of type II genes. This pattern also allows us to put forward the hypothesis that type I orthologues from other taxa are less conserved than most type II orthologues.
Indeed, PDS analyses presented here suggest that type I and II lineages have been subject to overall contrasting selection regimes. We found that recurrent positive selection could have played a role in fixing specific amino acids after several duplication events during the evolution of type I genes. In contrast, analyses of site models did not provide strong evidence for PDS selection among type II genes. We had to use site-branch models to detect a role for PDS among type II genes. Indeed, we found evidence for PDS along the branches leading to the groups of genes that control flowering time that evolved a new function with respect to most other genes characterized up to now that are involved in cell- or organ-type specification. Indeed, probably by their control of life-history traits, flowering-time genes may have directly impacted plant fitness, and this could also explain the prevalence of positive selection during protein evolution among them.
Sites with high PP of having been fixed by natural selection in both lineages were found mainly in the MADS and K domains. However, our analyses are biased toward these domains because we only focused on conserved stretches that may be unambiguously aligned and excluded most variable domains. Future studies focusing on particular closely related genes for several species will be useful to address the role of PDS within COOH and other divergent domains. Indeed, a recent study showed that regions within the C-terminal domain determine functional specificity in AP3 and PI and may be relevant for floral organ evolution (54). The localization of the sites with high PP identified here suggest a role for PDS in MADS-domain protein diversification through interactions with protein partners and changes in affinity to binding motifs (46, 47, 55). Sites and genes identified here to have been under PDS become interesting targets for functional evaluations.
Evaluations of assumptions and predictions made by models of gene duplication and persistence will require phylogenetically driven analyses of functional and population level data for related genes in different species. The MADS-box gene family might become a good “model family” for such a purpose. For example, our site model analyses did not find that PDS played a role in the divergence of redundant AP1, CAL, and FUL (56). However, population-level data do suggest a role for positive selection in the divergence of these genes (57). More powerful analyses should be performed to rule out false negatives in our analyses due to low gene number (32). Moreover, our conclusions were based on conservative analyses and unambiguously aligned sequences to avoid false positives. Additional tests (data available on request; Fig. 2b) suggest that the role of PDS in MADS-box gene evolution might be more widespread. Other approaches (58) and the inclusion of sequences for additional taxa should be considered when further investigating the role of PDS in MADS evolution.
Our results suggest a role for positive selection during MADS-box evolution in plants. Previous studies have emphasized the role of changes in cis-regulatory regions of transcriptional regulators during plant evolution (59). Fewer recent studies, however, have also demonstrated that the evolution of transcriptional regulators' cDNA sequences played important roles in plant evolution (54). The detection of positive selection in MADS protein sequences that are developmentally important also indicates that changes in cDNA, and not only in the regulatory regions of these genes, have played a role in the evolution of plant body plans.
We shared unpublished results and agreed on AGL numbers with Lucia Colombo and her collaborators. We thank Lorenzo Segovia and Julio Collado for computer time. Discussions with Francisco Vergara-Silva were insightful, and Alejandra Vásquez-Lobo, Rodolfo Salas, and Lev Jardón made comments and helped with the figures. Gary Ditta and Marty Yanofsky provided a preliminary list of MADS-like sequences. The comments of three anonymous reviewers improved this paper. Elizabeth Núñez helped with logistical tasks. This work was supported by Ph.D. fellowships to L.P.M.-C. from Consejo Nacional de Ciencia y Tecnología (CONACYT) and Dirección General de Estudios de Posgrado (DGEP) at the National Autonomous University of Mexico, and by grants from CONACYT, Programa de Apoyo para Proyectos de Investigación e Innovación Tecnológica, the Human Frontiers Science Program, and the University of California-Mexico (to E.R.A.-B.).
↵* To whom correspondence should be addressed. E-mail: ealvarez{at}miranda.ecologia.unam.mx.
Abbreviations: PDS, positive Darwinian selection; LRT, likelihood ratio test; AGL, agamous-like.
