Dispersal barriers in tropical oceans and speciation in Atlantic and eastern Pacific sea urchins of the genus Echinometra

Echinometra is a pantropical sea urchin made famous through studies of phylogeny, speciation, and genetic structure of the Indo-West Pacific (IWP) species. We sequenced 630 bp of the cytochrome oxidase I (COI) mitochondrial gene to provide comparable information on the eastern Pacific and Atlantic species, using divergence between those separated by closure of the Isthmus of Panama 3.1 million years ago (Ma) to estimate dates for cladogenic events. Most recently (1.27-1. 62 Ma), the Atlantic species E. lucunter and E. viridis diverged from each other, at a time in the Pleistocene that sea levels fell and Caribbean coral speciation and extinction rates were high. An earlier split, assumed to have been coincident with the completion of the Isthmus of Panama, separated the eastern Pacific E. vanbrunti from the Atlantic common ancestor. Transisthmian COI divergence similar to that in the sea urchin genus Eucidaris supports this assumption. The most ancient split in Echinometra occurred between the IWP and the neotropical clades, due to cessation of larval exchange around South Africa or across the Eastern Pacific Barrier. Gene flow within species is generally high; however, there are restrictions to genetic exchange between E. lucunter populations from the Caribbean and those from the rest of the Atlantic. Correlation between cladogenic and vicariant events supports E. Mayr's contention that marine species, despite their high dispersal potential, form by means of geographical separation. That sympatric, nonhybridizing E. lucunter and E. viridis were split so recently suggests, however, that perfection of reproductive barriers between marine species with large populations can occur in less than 1.6 million years (Myr).     

A MULTILOCUS TEST OF SIMULTANEOUS DIVERGENCE ACROSS THE ISTHMUS OF PANAMA USING SNAPPING SHRIMP IN THE GENUS ALPHEUS

The completion of the Panamanian Isthmus is one of the greatest natural experiments in evolution, sending multiple species pairs from a broad range of taxonomic groups on independent evolutionary trajectories. The resulting transisthmian sister species have been used as model systems for examining consequences that accompany cessation of gene flow in formerly panmictic populations. However, variance in pairwise genetic distances of these "geminates" often exceeds expectations, seemingly conflicting with the assumption that separation of populations was contemporaneous with the final closure of the Isthmus. Multilocus datasets and coalescent-based analytical methods can be used to estimate divergence times while accounting for variance in gene divergence that predates isolation, thus removing the need to invoke unequal divergence times. Here we present results from Bayesian analyses of sequence data from seven nuclear and one mitochondrial marker in eight transisthmian species pairs in the snapping shrimp genus Alpheus . Divergence times in two species pairs were shown to occur much earlier than the Isthmus final closure, but much of the variance in pairwise genetic distances from cytochrome oxidase I (COI) was explained when ancestral polymorphisms were accounted for. Results illustrate how coalescent approaches may be more appropriate for dating recent divergences than for estimating ancient speciation events.     

Confrontation of morphological and molecular data: the Praomys group (Rodentia, Murinae) as a case of adaptive convergences and morphological stasis

Phylogenetic relationships in a group of 21 African rodent species designated as the Praomys group (Murinae) were investigated using morphological characters and sequence data from the complete mitochondrial cytochrome b gene and nuclear IRBP gene fragment (840bp). The molecular results confirm the monophyly of the Praomys group, including the species Malacomys verschureni, while the other Malacomys species appear very divergent. The basal relationships within the Praomys group are poorly resolved, suggesting a rapid radiation at about 7-9 million years ago based on genetic divergence rates calibrated from the fossil record. Discrepancies between molecular and morphological results probably reflect of numerous convergences as well as variations in the rates of morphological evolution among lineages. Reconstructions of the ancestral character states suggest a savannah origin for the Praomys group, along with some morphological traits conserved by stasis in savannah taxa. At the same time, forest taxa seem to be characterized by an accelerated morphological evolution, with acquisition of convergent adaptive characters.     

Amphi-panamic geminates of snook (Percoidei: Centropomidae) provide a calibration of the divergence rate in the mitochondrial DNA control region of fishes.

The mitochondrial DNA control region is one of the most frequently utilized sequences for both intra- and interspecific genetic studies of fishes, yet a tenable divergence rate specifically for fish control regions has not been established. We attempted to establish a rate through a comparative study of control region sequences and those of a protein-coding mitochondrial gene region from geminate species of snook (Centropomus) assumed to have been separated by the emergence of the Isthmus of Panama 3.5 million years ago. The divergence rate suggested from the control region alignments between the geminates was markedly higher than rates currently applied in many studies, as was the rate suggested from the ND 5/6 protein coding region alignments. However, the suggested ND 5/6 rate when applied to alignments with the outgroup species was not concordant with the scant centropomid fossil record and therefore the assumed separation time of 3.5 million years seemed implausible. An average control region divergence rate was then estimated based on separation times of snook species derived assuming a divergence rate of 1% per million years for transversion substitutions at third codon positions in the ND 5/6 region. Using these separation times, a tenable average divergence rate for fish control regions of approximately 3.6% per million years +/- 0.46% SE was calculated.     

Evidence for three major clades within the snapping shrimp genus Alpheus inferred from nuclear and mitochondrial gene sequence data

The snapping shrimp genus Alpheus is among the most diverse of caridean shrimps, and analyses of taxa separated by the Isthmus of Panama have been used to estimate rates of molecular evolution. Although seven morphological groups have been informally suggested, no formal phylogenetic analysis of the genus has been previously attempted. Here we infer the phylogenetic relationships within Alpheus using sequence data from two nuclear genes, glucose-6-phosphate isomerase and elongation factor-1alpha, and from the mitochondrial gene cytochrome oxidase I. Three major clades corresponding to previously noted morphological features were identified. Discrepancies between earlier informal morphological groupings and molecular analyses largely consisted of species whose morphologies were not entirely typical of the group to which they had been assigned. The traditional placements of shrimp with highly sessile lifestyles and consequently simplified morphologies were also not supported by molecular analyses. Phylogenies for Alpheus suggest that specialized ecological requirements (e.g., symbiotic associations and estuarine habitats) and modified claw morphologies have evolved independently several times. These new analyses also support the sister species status of transisthmian pairs analyzed previously, although very similar pairs were not always resolved with the more slowly evolving nuclear loci. In addition, six new cryptic species were identified in the course of these studies plus a seventh whose status remains to be determined.     

Albumin evolution and phylogenetic relationships among Greek rodents of the families Arvicolidae and Muridae

The phylogenetic relationships of seven rodent species Microtus atticus, M. thomasi, M. epiroticus (family Arvicolidae) and Mus domesticus, Rattus norvegicus, Apodemus flavicollis and A. mystacinus (family Muridae) have been studied. In order to define these relationships we study the albumin evolution using the micro-complement fixation test (MC'F). No phylogenetic (immunological) distance between M. atticus and M. thomasi was found, a fact which confirms from the biochemical point of view the opinion that the former taxon is a synonym of the latter one. A molecular time scale relating MC'F immunological distances and geological time was established based on the assumption of a rate of 100 amino acid substitutions per–20 million years. The time of divergence between M. epiroticus and M. thomasi was estimated to be 0.5–0.6 million years ago (Pleistocene). Such a recent divergence corroborates the opinion based on morphological and protein electrophoretic criteria according to which Terricola (formerly Pitymys ) must be considered as a subgenus of the genus Microtus and not as a distinct genus Pitymys , as previously had been accepted. Apodemus flavicollis and A. mystacinus were separated about 0.65–0.8 million years ago (Pleistocene). The Rattus norvegicus lineage was separated–12.5 million years ago (end of Miocene), shortly before the Mus and Apodemus divergence. Our data indicate that the common ancestor of Arvicolidae and Muridae lived–25 million years ago (early Miocene). All these results are in agreement with paleontological and some recent DNA-DNA hybridization and electrophoretic data.     

Of mice and men: Fossil-based divergence dates and molecular “clocks”

The magnitude of immunological differences between species has been used to estimate the time of divergence of lineages, especially among primates. Calibration of “clocks” based on immunological differences has relied heavily on data from Pakistan suggest that the lineages leading to Mus and Rattus diverged between 14 and 8 million years ago rather than 30 or more million years ago as suggested by techniques invoking constant rates of molecular evolution. Molecular analyses of pairs of mammal species which appear well documented paleontologically should be undertaken in order to expand our understanding of the tempo and mode of molecular evolution.     

Evolutionary history of the genus Rhamdia (Teleostei: Pimelodidae) in Central America

We constructed phylogenetic hypotheses for Mesoamerican Rhamdia, the only genus of primary freshwater fish represented by sympatric species across Central America. Phylogenetic relationships were inferred from analysis of 1990 base pairs (bp) of mitochondrial DNA (mtDNA), represented by the complete nucleotide sequences of the cytochrome b (cyt b) and the ATP synthase 8 and 6 (ATPase 8/6) genes. We sequenced 120 individuals from 53 drainages to provide a comprehensive geographic picture of Central American Rhamdia systematics and phylogeography. Phylogeographic analysis distinguished multiple Rhamdia mtDNA lineages, and the geographic congruence across evolutionarily independent Rhamdia clades indicated that vicariance has played a strong role in the Mesoamerican diversification of this genus. Phylogenetic analyses of species-level relationships provide strong support for the monophyly of a trans-Andean clade of three evolutionarily equivalent Rhamdia taxa: R. guatemalensis, R. laticauda, and R. cinerascens. Application of fish-based mitochondrial DNA clocks ticking at 1.3-1.5% sequence divergence per million years (Ma), suggests that the split between cis- and trans-Andean Rhamdia extends back about 8 Ma, and the three distinct trans-Andean Rhamdia clades split about 6 Ma ago. Thus the mtDNA divergence observed between cis- and trans-Andean Rhamdia species is too low to support an ancient colonization of Central America in the Late Cretaceous or Paleocene as had been hypothesized in one colonization model for Mesoamerican fishes. Rather the mtDNA data indicate that Rhamdia most likely colonized Central America in the late Miocene or Pliocene, promoting a strong role for the Isthmus of Panamá in the Mesoamerican expansion of this genus. Basal polytomies suggest that both the R. laticauda and R. guatemalensis clades spread rapidly across the Central American landscape, but differences in the average mtDNA genetic distances among clades comprising the two species, indicate that the R. laticauda spread and diversified across Mesoamerica about 1 million years before R. guatemalensis.     

Human and ape molecular clocks and constraints on paleontological hypotheses

Although the relationships of the living hominoid primates (humans and apes) are well known, the relationships of the fossil species, times of divergence of both living and fossil species, and the biogeographic history of hominoids are not well established. Divergence times of living species, estimated from molecular clocks, have the potential to constrain hypotheses of the relationships of fossil species. In this study, new DNA sequences from nine protein-coding nuclear genes in great apes are added to existing datasets to increase the precision of molecular time estimates bearing on the evolutionary history of apes and humans. The divergence of Old World monkeys and hominoids at the Oligocene-Miocene boundary (approximately 23 million years ago) provides the best primate calibration point and yields a time and 95% confidence interval of 5.4 +/- 1.1 million years ago (36 nuclear genes) for the human-chimpanzee divergence. Older splitting events are estimated as 6.4 +/- 1.5 million years ago (gorilla, 31 genes), 11.3 +/- 1.3 million years ago (orangutan, 33 genes), and 14.9 +/- 2.0 million years ago (gibbon, 27 genes). Based on these molecular constraints, we find that several proposed phylogenies of fossil hominoid taxa are unlikely to be correct.     

Deceleration of morphological evolution in a cryptic species complex and its link to paleontological stasis

Morphological stasis or the absence of morphological change is a well-known phenomenon in the paleontological record, yet it is poorly integrated with neontological evidence. Recent evidence suggests that cryptic species complexes may remain morphologically identical due to morphological stasis. Here, we describe a case of long-term stasis in the Stygocapitella cryptic species complex (Parergodrilidae, Orbiniida, Annelida). Using phylogenetic methods and morphological data, we find that rates of morphological evolution in Stygocapitella are significantly slower than in closely related taxa (Nerillidae, Orbiniidae). Assessment of quantitative and qualitative morphology revealed the presence of four morphotypes with only subtle differences, whereas molecular data supports 10 reproductively isolated clades. Notably, estimates for the time of Stygocapitella species divergence range from ∼275 million years to ∼18 million years, including one case of two morphologically similar species that have diverged about 140 million years ago. These findings provide evidence for morphological deceleration and long-term morphological stasis in Stygocapitella, and that speciation is not necessarily accompanied by morphological changes. The deceleration of morphological divergence in Stygocapitella can be potentially linked to niche conservatism and tracking, coupled with the fluctuating dynamics of the interstitial environment, or genetic constraints due to progenetic evolution. Finally, we conclude that failing to integrate speciation without morphological evolution in paleontology may bias estimates of rates of speciation and morphological evolution.     

Reconstruction of phylogenetic relationships within Grapsidae (Crustacea: Brachyura) and comparison of trans-isthmian versus amphi-atlantic gene flow based on mtDNA

Following taxonomic revisions in recent years, the originally large family Grapsidae MacLeay, 1838 has become a relatively small and morphologically homogeneous family in terms of adult and larval morphology. Most available molecular studies including more than one genus of the family have also suggested monophyly of the corresponding taxa. However, no single phylogenetic study has ever included all constituent genera of the Grapsidae. In the current study, a molecular phylogeny based on sequences of the mitochondrial 16S rRNA gene from all eight grapsid genera and 34 species is presented and suggests that up to four genera are not monophyletic. This is mainly due to the polyphyletic nature of the genus Pachygrapsus which can be found in six different lineages of the phylogeny, suggesting that the genus currently does not represent a single evolutionary lineage and is in need of taxonomic revision. Amphi-atlantic and trans-isthmian species pairs or populations in four genera are compared and reveal relatively constant and pronounced divergences across the Panama Isthmus as opposed to moderate divergences across the Atlantic Ocean, thereby suggesting occurrence of gene flow across the Atlantic Ocean during the past three million years.     

Molecular distance and divergence time in carnivores and primates

Numerous studies have used indices of genetic distance between species to reconstruct evolutionary relationships and to estimate divergence time. However, the empirical relationship between molecular-based indices of genetic divergence and divergence time based on the fossil record is poorly known. To date, the results of empirical studies conflict and are difficult to compare because they differ widely in their choice of taxa, genetic techniques, or methods for calibrating rates of molecular evolution. We use a single methodology to analyze the relationship of molecular distance and divergence time in 86 taxa (72 carnivores and 14 primates). These taxa have divergence times of 0.01-55 Myr and provide a graded series of phylogenetic divergences such that the shape of the curve relating genetic distance and divergence time is often well defined. The techniques used to obtain genetic distance estimates include one- and two-dimensional protein electrophoresis, DNA hybridization, and microcomplement fixation. Our results suggest that estimates of molecular distance and divergence time are highly correlated. However, rates of molecular evolution are not constant; rather, in general they decline with increasing divergence time in a linear fashion. The rate of decline may differ according to technique and taxa. Moreover, in some cases the variability in evolutionary rates changes with increasing divergence time such that the accuracy of nodes in a phylogenetic tree varies predictably with time.     

Evolutionary diversification of banded tube-dwelling anemones (Cnidaria; Ceriantharia; Isarachnanthus) in the Atlantic Ocean

The use of molecular data for species delimitation in Anthozoa is still a very delicate issue. This is probably due to the low genetic variation found among the molecular markers (primarily mitochondrial) commonly used for Anthozoa. Ceriantharia is an anthozoan group that has not been tested for genetic divergence at the species level. Recently, all three Atlantic species described for the genus Isarachnanthus of Atlantic Ocean, were deemed synonyms based on morphological simmilarities of only one species: Isarachnanthus maderensis. Here, we aimed to verify whether genetic relationships (using COI, 16S, ITS1 and ITS2 molecular markers) confirmed morphological affinities among members of Isarachnanthus from different regions across the Atlantic Ocean. Results from four DNA markers were completely congruent and revealed that two different species exist in the Atlantic Ocean. The low identification success and substantial overlap between intra and interspecific COI distances render the Anthozoa unsuitable for DNA barcoding, which is not true for Ceriantharia. In addition, genetic divergence within and between Ceriantharia species is more similar to that found in Medusozoa (Hydrozoa and Scyphozoa) than Anthozoa and Porifera that have divergence rates similar to typical metazoans. The two genetic species could also be separated based on micromorphological characteristics of their cnidomes. Using a specimen of Isarachnanthus bandanensis from Pacific Ocean as an outgroup, it was possible to estimate the minimum date of divergence between the clades. The cladogenesis event that formed the species of the Atlantic Ocean is estimated to have occured around 8.5 million years ago (Miocene) and several possible speciation scenarios are discussed.     

Hominid and gelada baboon evolution: Agreement between molecular and fossil time scales

The time of origin of the hominid lineage has long been debated. Macromolecular studies have consistently shown genetic distances between living humans and African apes to be quite small. The molecular clock hypothesis proposes that the time of separation of these lineages is relatively recent (in the range of 4–8 million years ago) and not 15 million years or more ago as usually suggested. Three independent molecular comparisons yield a mean estimate of 4.6 million years for the hominid-African pongid divergence. The relationship of Theropithecusand Papiois a parallel case within Primates of two taxa which are quite similar at the molecular level, but which are usually thought to have separated relatively long ago. The two cases of seeming discordance between different lines of evidence are analogous. Each involves a speciation event which eventually resulted in one substantially derived lineage and one or more relatively unchanged lineages. In each case, claims of the antiquity of the divergence event extend to at least twice the age of the first certain appearance of the more derived lineage in the fossil record. Finally, in each case, the molecular clock model suggests a range of possible divergence times that overlaps with the first appearances of undoubted hominids and Theropithecusin the fossil record. This test involving paleontological evidence supports the molecular clock hypothesis.     

Genetic variation and phylogeography of free-living mouse species (genus Mus) in the Balkans and the Middle East

This work presents a study of the distribution and pattern of variation throughout the ranges of three free-living mouse species of the genus Mus-M. macedonicus, M. spicilegus, and a M. cypriacus - based on sequencing of two segments of the mitochondrial DNA (mtDNA) control region. The study shows a similar level of variability in the three species and suggests their recent population expansion. The highest proportion of variation is found within populations indicating low genetic structuring. Phylogenetic analysis confirms the significant divergence of a mitochondrial lineage of M. macedonicus from Israel, recently described as a new subspecies, M. macedonicus spretoides. Conversely, no genetic hiatus is revealed between European and Asian populations of M. macedonicus macedonicus. Although phylogenetic relationships among M. spicilegus populations could not be unravelled precisely, the results suggest a recent westward expansion of the species. The mtDNA divergence between M. macedonicus and M. spicilegus is 7.3%, suggesting their split between c. 700,000 and 1 million years ago. These dates correspond with a coalescent estimate about 720,000 years ago. On the other hand, M. cypriacus appeared almost twice as divergent from the former species (4.5%) as from the latter (8.8%) suggesting a divergence of c. 430,000-610,000 years ago (coalescent approximately 490,000 years ago) and 830,000-1.2 million years ago (coalescent approximately 780,000 years ago), respectively. Approximate times of population expansion have also been estimated for all taxa and groups of populations. Existence of several glacial refuges and various colonization scenarios are discussed; since all estimated divergence times fall within interglacial periods it seems that climatic oscillations did not play a crucial role in the evolution of the three species.     

Historical biogeography of the livebearing fish genus Poeciliopsis (Poeciliidae: Cyprinodontiformes)

To assess the historical biogeography of freshwater topminnows in the genus Poeciliopsis, we examined sequence variation in two mitochondrial genes, cytochrome b (1140 bp) and NADH subunit 2 (1047 bp). This widespread fish genus is distributed from Arizona to western Colombia, and nearly half of its 21 named species have distributions that border on the geologically active Trans-Mexican Volcanic Belt (TMVB), a region that defines the uplifted plateau (Mesa Central) of Mexico. We used the parametric bootstrap method to test the hypothesis that a single vicariant event associated with the TMVB was responsible for divergence of taxa found to the north and south of this boundary. Because the single-event hypothesis was rejected as highly unlikely, we hypothesize that at least two geological events were responsible for divergence of these species. The first (8-16 million years ago) separated ancestral populations that were distributed across the present TMVB region. A second event (2.8-6.4 million years ago) was associated with northward dispersal and subsequent vicariance of two independent southern lineages across the TMVB. The geological history of this tectonically and volcanically active region is discussed and systematic implications for the genus are outlined.     

Genomic data reveal a protracted window of introgression during the diversification of a neotropical woodcreeper radiation*

The incidence of introgression during the diversification process and the timespan following divergence when introgression is possible are poorly understood in the neotropics where high species richness could provide extensive opportunities for genetic exchange. We used thousands of genome-wide SNPs to infer phylogenetic relationships, calculate ages of splitting, and to estimate the timing of introgression in a widespread avian neotropical genus of woodcreepers. Five distinct introgression events were reconstructed involving taxa classified both as subspecies and species including lineages descending from the basal–most split, dated to 7.3 million years ago. Introgression occurred between just a few hundred thousand to about 2.5 million years following divergence, suggesting substantial portions of the genome are capable of introgressing across taxa boundaries during a protracted time window of a few million years following divergence. Despite this protracted time window, we found that the proportion of the genome introgressing (6–11%) declines with the time of introgression following divergence, suggesting that the genome becomes progressively more immune to introgression as reproductive isolation increases.     

Test for simultaneous divergence using approximate Bayesian computation

Comparative phylogeographic studies often reveal disparate levels of sequence divergence between lineages spanning a common geographic barrier, leading to the conclusion that isolation was nonsynchronous. However, only rarely do researchers account for the expected variance associated with ancestral coalescence and among-taxon variation in demographic history. We introduce a flexible approximate Bayesian computational (ABC) framework that can test for simultaneous divergence (TSD) using a hierarchical model that incorporates idiosyncratic differences in demographic history across taxon pairs. The method is tested across a range of conditions and is shown to be accurate even with single-locus mitochondrial DNA (mtDNA) data. We apply this method to a landmark dataset of putative simultaneous vicariance, eight geminate echinoid taxon pairs thought to have been split by the Isthmus of Panama 3.1 million years ago. The ABC posterior estimates are not consistent with a history of simultaneous vicariance given these data. Subsequent ABC estimates under a constrained model that assumes two divergence times across the eight taxon pairs suggests simultaneous divergence 3.1 million years ago in seven of the taxon pairs and a more recent divergence in the remaining taxon pair. These ABC estimates on the simultaneous divergence of the seven taxon pairs correspond to a DNA substitution rate of approximately 1.59% per lineage per million years at the mtDNA cytochrome oxidase I gene. This ABC framework can easily be modified to analyze single taxon-pair datasets and/or be expanded to include multiple loci, migration, recombination, and other idiosyncratic demographic histories. The flexible aspect of ABC and its built-in evaluation of estimator bias and statistical power has the potential to greatly enhance statistical rigor in phylogeographic studies.     

Timing the eastern Asian-eastern North American floristic disjunction: molecular clock corroborates paleontological estimates

Sequence data of the chloroplast gene rbcL were used to estimate the time of the well-known eastern Asian-eastern North American floristic disjunction. Sequence divergence of rbcL was examined for 22 species of 11 genera (Campsis, Caulophyllum, Cornus, Decumaria, Liriodendron, Menispermum, Mitchella, Pachysandra, Penthorum, Podophyllum, and Phryma) representing a diverse array of flowering plants occurring disjunctly in eastern Asia and eastern North America. Divergence times of putative disjunct species pairs were estimated from synonymous substitutions, using rbcL molecular clocks calibrated for Cornus. Relative rate tests were performed to assess rate constancy of rbcL evolution among lineages. Corrections of estimates of divergence times for each species pair were made based on rate differences of rbcL between Cornus and other species pairs. Results of these analyses indicate that the time of divergence of species pairs examined ranges from 12.56 +/- 4.30 million years to recent (<0.31 million years), with most within the last 10 million years (in the late Miocene and Pliocene). These results suggest that the isolation of most morphologically similar disjunct species in eastern Asia and eastern North America occurred during the global climatic cooling period that took place throughout the late Tertiary and Quaternary. This estimate is closely correlated with paleontological evidence and in agreement with the hypothesis that considers the eastern Asian-eastern North American floristic disjunction to be the result of the range restriction of a once more or less continuously distributed mixed mesophytic forest of the Northern Hemisphere that occurred during the late Tertiary and Quaternary. This implies that in most taxa the disjunction may have resulted from vicariance events. However, long-distance dispersal may explain the disjunct distribution of taxa with low divergence, such as Menispermum.     

Rates of molecular evolution in nuclear genes of east Mediterranean scorpions

Scorpions of the genus Mesobuthus represent a useful terrestrial model system for studying molecular evolution. They are distributed on several Aegean islands and the adjacent mainland, they are believed to have low rates of dispersal, and evolutionary divergence dates of taxa are available based on biogeographic events that separated islands from each other and the mainland. Here, we present data on polymorphism and synonymous (Ks) and non-synonymous (Ka) substitution rates for nine nuclear protein-coding genes of two east Mediterranean scorpion species, Mesobuthus gibbosus and M. cyprius (Buthidae). Levels of polymorphism tend to be lower in populations from islands (mean nucleotide diversity pi = 0.0071 +/- 0.0028) than in mainland populations (mean pi = 0.0201 +/- 0.0085). By using linear regression of genetic divergence versus isolation time, we estimate Ks to be 3.17 +/- 1.54 per (site x 10(9) years), and Ka to be 0.39 +/- 0.94 per (site x 10(9) years). These estimates for both Ks and Ka are considerably lower than for many other invertebrates, such as Drosophila, and may be attributed to scorpions' mammal-like generation times (approximately 2 years) and low metabolic rates. Phylogenetic analysis using maximum likelihood revealed a phylogeny that is congruent with that expected based on biogeographic events and in which divergences at synonymous sites are proportional to the dates that the taxa are believed to have split. Tests of equality of branch lengths for the Cyprus and Crete lineages revealed that Ks-estimates are about the same in both lineages, as expected from the biogeographic events that separated the islands, but Ka was increased in the Cyprus lineage compared to the Cretan lineage.