Phylogeography and monophyly of the swordtail fish species Xiphophorus birchmanni (Cyprinodontiformes, Poeciliidae)

We used sequences of the mitochondria control region to assess the distribution of genetic variation within and among populations of the poeciliid fish species Xiphophorus birchmanni . We collected 122 X. birchmanni samples from 11 sites in three drainage systems comprising the distribution of the species. We found low levels of polymorphism among aligned sequences and low levels of genetic variation within populations but high levels of genetic differentiation among populations. Haplotypes are exclusive to three river drainages (Los Hules, Calabozo and San Pedro). Mantel tests revealed correlations between geographical (both straight-line and river distances) and genetic distance, consistent with an isolation by distance scenario, while nested clade analysis suggested allopatric fragmentation between haplotypes from two of the major drainages, and isolation by distance with restricted gene flow within those drainages. Finally, monophyly of X. birchmanni is strongly supported while the previous hypothesis of the evolutionary origin of this species from X. malinche is not.     

Phylogeny of the Macaronesian endemic Crambe section Dendrocrambe (Brassicaceae) based on internal transcribed spacer sequences of nuclear ribosomal DNA

The 14 species of Crambe L. sect. Dendrocrambe DC. (Brassicaceae) form a monophyletic group endemic to the Canary and Madeira archipelagos. Both parsimony and maximum likelihood analyses of sequence data from the two internal transcribed spacer regions of nuclear ribosomal DNA were used to estimate phylogenetic relationships within this section. These analyses support the monophyly of three major clades. No clade is restricted to a single island, and therefore it appears that inter-island colonization has been the main avenue for speciation in these two archipelagos. The two species endemic to Fuerteventura (C. sventenii) and Madeira (C. fruticosa) comprise a clade, providing the first evidence for a floristic link between the Eastern Canary Islands and the archipelago of Madeira. Both maximum likelihood and weighted parsimony analyses show that this clade is sister to the two other clades, although bootstrap support for this relationship is weak. Parsimony optimizations of ecological zones and island distribution suggest a colonization route from the low-altitude areas of the lowland scrub toward the high-elevation areas of the laurel and pine forests. In addition, Tenerife is likely the ancestral island for species endemic to the five westernmost islands of Gran Canaria, La Gomera, El Hierro, La Palma, and Tenerife.     

The automation and evaluation of nested clade phylogeographic analysis

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.     

Population structure of the banana weevil, an introduced pest in the Canary Islands, studied by RAPD analysis

The banana weevil (BW), Cosmopolites sordidus (Coleoptera: Curculionidae), is one of the most important insect pests of bananas and plantains. The mobility and the origin of BW infestations at the Canary Islands (Tenerife, La Gomera and La Palma) have been analysed using Random Amplified Polymorphic DNA (RAPD) as molecular markers. Populations from Costa Rica, Colombia, Uganda and Madeira were also included for comparison. One hundred and fifteen reproducible bands from eight primers were obtained. The level of polymorphism in the populations from the Canary Islands (40-62%) was in the range of those found in other populations. Nei's genetic distances, pair-wise fixation index (FST) values indicate that the closest populations are Tenerife populations among themselves (Nei's genetic distance=0.054-0.100; FST=0.091-0.157) and Costa Rica and Colombia populations (Nei's genetic distance=0.049; FST=0.113). Our results indicate the existence of BW local biotypes with limited gene flow and affected by genetic drift. These results are compatible with a unique event of colonization at Tenerife; whereas, the outbreaks in La Gomera and La Palma may come from independent introductions. The Madeira population is phylogenetically and geographically closer to the Canary Islands populations, suggesting that it is the most likely source of the insects introduced in the Canary Islands.     

Patterns of spatial genetic structuring in a hydropsychid caddisfly (Cheumatopsyche sp. AV1) from southeastern Australia

We assessed levels of mitochondrial genetic spatial structuring in the hydropsychid caddisfly Cheumatopsyche sp. AV1 in southeastern New South Wales, Australia. No significant spatial structuring was detected within or between catchments using analysis of molecular variance, and nested clade contingency analysis suggested no strong relationship between haplotypes and geographical location, at any clade level. However, tests for association among haplotypes incorporating geographical distance in the nested clade analysis, revealed patterns of historical range expansion and recent restricted gene flow. Most likely, population fragmentation preceded range expansion, although subsequent recontact and gene flow among the previously sundered populations has apparently obscured the geographical signature of the former fragmentation. Taken together, our analyses suggest that a number of populations fragmented during the Pleistocene evolved in isolation for a time and subsequently expanded into secondary contact. Since expansion, there has apparently been substantial (albeit somewhat restricted) dispersal and gene flow of adult female Cheumatopsyche sp. AV1, throughout the study area.     

Radiation of Atlantic goldcrests Regulus regulus spp.: evidence of a new taxon from the Canary Islands

Phylogenetic relationships between goldcrest populations from the Atlantic Islands (Azores and Canary Islands) were investigated by two molecular markers (mitochondrial control region and cytochrome b sequences), and partly by morphology and territorial song. The Azorean goldcrest populations are closely related to European nominate R. r. regulus. Most probably, the Azores were colonized by goldcrests in a single late‐pleistocene invasion, while colonization of the Canary Islands presumably occurred in two steps: An early invasion to Tenerife and La Gomera 1.9–2.3 million years (my) ago and a more recent one to El Hierro and La Palma 1.3–1.8 my ago. Distribution of haplotypes on the Azores suggests a division of R. r. azoricus on São Miguel into an eastern population with close affinities to R. r. sanctaemariae and a western population belonging to the lineage of R. r. inermis on the central and western island group. The Canarian populations are genetically substructured into a northeastern group embracing Tenerife and La Gomera and a second, southwestern group including El Hierro and La Palma. Genetic distances between members of the two Canarian clades range at 3.1–3.4% (TrN distance, control region and cytochrome b). Differentiation between the two groups is also supported by morphology and by territorial song. Substitution rate estimates for the both genes range at approximately the same values of 0.0031 and 0.0044 substitutions per site and lineage per my which roughly corresponds 0.61–0.83% divergence between Regulus lineages per my. Highest local rates occur in island clades of the Azorean and the Canarian population and in R. r. japonensis from the Russian Far East and Japan. However, a general acceleration of a molecular clock in island populations is not evident from the Regulus data set due to extremely low local rate estimates in the Canarian clade of Tenerife and La Gomera. As a taxonomic consequence of the marked differentiation of the two Canarian goldcrest clades the populations from El Hierro and La Palma are described as a taxon new to science and are named Regulus regulus ellenthalerae n. ssp.     

Phylogeography and population genetics of the maize stalk borer Busseola fusca (Lepidoptera, Noctuidae) in sub-Saharan Africa

The population genetics and phylogeography of African phytophagous insects have received little attention. Some, such as the maize stalk borer Busseola fusca, display significant geographic differences in ecological preferences that may be congruent with patterns of molecular variation. To test this, we collected 307 individuals of this species from maize and cultivated sorghum at 52 localities in West, Central and East Africa during the growing season. For all collected individuals, we sequenced a fragment of the mitochondrial cytochrome b. We tested hypotheses concerning the history and demographic structure of this species. Phylogenetic analyses and nested clade phylogeographic analyses (NCPA) separated the populations into three mitochondrial clades, one from West Africa, and two--Kenya I and Kenya II--from East and Central Africa. The similar nucleotide divergence between clades and nucleotide diversity within clades suggest that they became isolated at about the same time in three different refuges in sub-Saharan Africa and have similar demographic histories. The results of mismatch distribution analyses were consistent with the demographic expansion of these clades. Analysis of molecular variance (amova) indicated a high level of geographic differentiation at different hierarchical levels. NCPA suggested that the observed distribution of haplotypes at several hierarchical levels within the three major clades is best accounted for by restricted gene flow with isolation by distance. The domestication of sorghum and the introduction of maize in Africa had no visible effect on the geographic patterns observed in the B. fusca mitochondrial genome.     

Phylogeography of the endangered darkling beetle species of Pimelia endemic to Gran Canaria (Canary Islands)

Phylogenetic and geographical nested clade analysis (NCA) methods were applied to mitochondrial DNA sequences of Pimelia darkling beetles (Coleoptera, Tenebrionidae) endemic to Gran Canaria, an island in the Canary archipelago. The three species P. granulicollis, P. estevezi and P. sparsa occur on the island, the latter with three recognized subspecies. Another species, P. fernandezlopezi (endemic to the island of La Gomera) is a close relative of P. granulicollis based on partial Cytochrome Oxidase I mtDNA sequences obtained in a previous study. Some of these beetles are endangered, so phylogeographical structure within species and populations can help to define conservation priorities. A total of about 700 bp of Cytochrome Oxidase II were examined in 18 populations and up to 75 individuals excluding outgroups. Among them, 22 haplotypes were exclusive to P. granulicollis and P. estevezi and 31 were from P. sparsa. Phylogenetic analysis points to the paraphyly of Gran Canarian Pimelia, as the La Gomera P. fernandezlopezi haplotypes are included in them, and reciprocal monophyly of two species groups: one constituted by P. granulicollis, P. estevezi and P. fernandezlopezi (subgenus Aphanaspis), and the other by P. sparsa'sensu lato'. The two species groups show a remarkably high mtDNA divergence. Within P. sparsa, different analyses all reveal a common result, i.e. conflict between current subspecific taxonomic designations and evolutionary units, while P. estevezi and P. fernandezlopezi are very close to P. granulicollis measured at the mtDNA level. Geographical NCA identifies several cases of nonrandom associations between haplotypes and geography that may be caused by allopatric fragmentation of populations with some cases of restriction of gene flow or range expansion. Analyses of molecular variance and geographical NCA allow definition of evolutionary units for conservation purposes in both species-groups and suggest scenarios in which vicariance caused by geological history of the island may have shaped the pattern of the mitochondrial genetic diversity of these beetles.     

Phylogeography of chamois (Rupicapra spp.) inferred from microsatellites

Evolutionary relationships among populations of chamois (Rupicapra spp.) across their current range from the Caucasus to the Cantabrian Mountains were investigated. The allelic variation in 23 microsatellite loci was assessed in eight geographical populations, recognised as subspecies of the two closely related species R. pyrenaica and R. rupicapra. Analysis of variance in allele frequencies (Fst, statistics) and in repeat numbers (Rst, statistics) showed these data to be highly structured. Two genetic distances between pairs of populations, Ds and (δμ)², were computed and phylogenetic trees were constructed. Similar patterns were produced by the different statistics. All trees indicate a deep divergence between the two recognised species, which is compatible with archaeological data that place their split in the Riss–Würm interglacial period. Genetic distances between pairs of populations are highly correlated with geographical distance. This suggests that the history of the genus during Pleistocene glacial-interglacial periods was dominated by expansions and contractions within limited geographic regions, leading to alternate contact and isolation of contiguous populations. In addition, the alpine barrier has played a substantial role in West–East differentiation.     

Evidence for long-distance dispersal in a sedentary passerine, Gymnorhina tibicen (Artamidae)

Australian magpies (Gymnorhina tibicen) are group-living birds found across much of mainland Australia. Adults commonly remain in a breeding territory until death. Young of the year either remain on the natal (birth) site or are forced by their parents to disperse. Observational studies in south-eastern Australia suggest that most dispersing juveniles settle within 7 km of their natal territory. Therefore, despite potential for considerable gene flow (via flight), social organization predisposes magpies towards local population structuring. In this study, we measured genetic variation at both nuclear (allozyme) and mitochondrial loci and found evidence of substantial gene flow over very large distances (up to 1599 km). Thus, some juvenile magpies may disperse much greater distances than was previously thought. For mtDNA, geographic and genetic distance were strongly correlated, consistent with a pattern of isolation by distance. Therefore, although female gene flow is substantial it is apparently geographically restricted over large distances, in approximately a stepping-stone fashion. We conclude that a strong relationship between gene flow and geographic distance can develop even over large distances if populations have experienced no major historical disturbances to gene flow.     

Phylogeography of a mountain lizard species: an ancient fragmentation process mediated by riverine barriers in the Liolaemus monticola complex (Sauria: Liolaemidae)

Liolaemus monticola is a mountain lizard species, with a widespread distribution from central Chile that displays several highly polymorphic chromosomal races. Our study determined the phylogeographic structuring and relationships among three chromosomal races of L. monticola in Chile. Mitochondrial DNA (mtDNA) sequences of the cytochrome b gene were examined using the following phylogenetic methods: maximum parsimony, maximum likelihood, Bayesian inference and nested clade phylogeographic analyses (NCPAs). These methods revealed two major monophyletic clades (north and south) in the L. monticola species, with non-overlapping geographical locations separated by the Maipo and Yeso rivers (except one hybrid, from a zone of secondary contact). The NCPA showed that a past fragmentation process likely resulted in the separation of the two clades. The southern clade includes all samples of the 'Southern, 2 n  = 34' race; the northern clade is comprised of all remaining derived chromosomal races: the 'Northern, 2 n  = 38–40 and the Multiple Fission, 2 n  = 42–44' races. Our results support the hypothesis of a geographical and genetic split resulting from allopatric processes caused by riparian barriers acting over a long time period. The inferred biogeographical scenario shows that populations have moved from the south to the north using the Andean mountains as the primary corridor for dispersal.     

Phylogeography of skinks (Chalcides) in the Canary Islands inferred from mitochondrial DNA sequences

Mitochondrial DNA (mtDNA) evolution was investigated in skinks of the genus Chalcides found in the Canary Islands ( Ch. sexlineatus, Ch. viridanus and Ch. simonyi ), together with some North African congenerics ( Ch. polylepis and Ch. mionecton ). Several sites were included within islands to cover areas of known within-island geographical variation in morphology. Skinks from the islands of El Hierro and La Gomera appear to be sister taxa. The relationships between this clade and the Tenerife and Gran Canarian skinks were not fully resolved, although the best working hypothesis indicated monophyly with the former, with the latter forming a closely related outgroup. Ch. simonyi from Fuerteventura was more distantly related to the Western Canary Island skinks and did not show close relationships with the North African species Ch. mionecton and Ch . polylepis . Possible colonization sequences for the four most Western Canary Islands were considered. El Hierro appears to have been colonized relatively recently from La Gomera, commensurate with the recent origin of this island, while dispersal between La Gomera and Tenerife and between Gran Canaria and Tenerife or La Gomera appears to have taken place considerably earlier. Substantial within-island haplotype divergence was found in Gran Canaria and Tenerife. This may be a result of recent periods of intense volcanic activity found within these two islands. Lower levels of within-island differentiation are found in La Gomera and El Hierro and may be explained by lower levels of volcanic activity during recent geological history and a more recent colonization, respectively.     

Phylogeography of chamois (Rupicapra spp.) inferred from microsatellites

Evolutionary relationships among populations of chamois (Rupicapra spp.) across their current range from the Caucasus to the Cantabrian Mountains were investigated. The allelic variation in 23 microsatellite loci was assessed in eight geographical populations, recognised as subspecies of the two closely related species R. pyrenaica and R. rupicapra. Analysis of variance in allele frequencies (Fst, statistics) and in repeat numbers (Rst, statistics) showed these data to be highly structured. Two genetic distances between pairs of populations, Ds and (deltamu)(2), were computed and phylogenetic trees were constructed. Similar patterns were produced by the different statistics. All trees indicate a deep divergence between the two recognised species, which is compatible with archaeological data that place their split in the Riss-Würm interglacial period. Genetic distances between pairs of populations are highly correlated with geographical distance. This suggests that the history of the genus during Pleistocene glacial-interglacial periods was dominated by expansions and contractions within limited geographic regions, leading to alternate contact and isolation of contiguous populations. In addition, the alpine barrier has played a substantial role in West-East differentiation.     

Genetic and socio-cultural differentiation in the aborigines of Arnhem Land, Australia

Four tribes of Arnhem Land were surveyed for dermatoglyphics; based on pattern intensity indices, total ridge-counts, and a distance statistic combining the two, it was shown that the tribes can be arranged into western (Tiwi, Gunwinggu) and eastern (“Murngin” and Andilyaugwa) groups. This substantiates observations made on Arnhem Land by linguists and social anthropologists. From a survey of allele frequency traits of blood factors and P.T.C. tasting, distance statistics were computed between the four tribes. These confirmed the relative isolation of the extreme Arnhem Land tribes. Distance statistics were also computed between the four tribes and two Central Australian tribes, the Aranda and Wailbri. The Aranda and “Murngin” were relatively close together agreeing with theories that the Aranda are derived from a not too remote southwards migration from north Arnhem Land, as supported by linguistic data. Correlations between the biological, geographical and linguistic distances were positive, and generally agreed with the expectation that socio-cultural and linguistic barriers are important in regulating gene flow between populations. This study emphasizes the need to consider biological distances in association with ecological and socio-cultural factors.     

Complex population genetic structure in the endemic Canary Island pine revealed using chloroplast microsatellite markers

The Canary archipelago, located on the northwestern Atlantic coast of Africa, is comprised of seven islands aligned from east to west, plus seven minor islets. All the islands were formed by volcanic eruptions and their geological history is well documented providing a historical framework to study colonization events. The Canary Island pine (Pinus canariensis C. Sm.), nowadays restricted to the westernmost Canary Islands (Gran Canaria, Tenerife, La Gomera, La Palma and El Hierro), is considered an old (Lower Cretaceous) relic from an ancient Mediterranean evolutionary centre. Twenty seven chloroplast haplotypes were found in Canary Island pine but only one of them was common to all populations. The distribution of haplotypic variation in P. canariensis suggested the colonization of western Canary Islands from a single continental source located close to the Mediterranean Basin. Present-day populations of Canary Island pine retain levels of genetic diversity equivalent to those found in Mediterranean continental pine species, Pinus pinaster and Pinus halepensis. A hierarchical analysis of variance (AMOVA) showed high differentiation among populations within islands (approximately 19%) but no differentiation among islands. Simple differentiation models such as isolation by distance or stepping-stone colonization from older to younger islands were rejected based on product-moment correlations between pairwise genetic distances and both geographic distances and population-age divergences. However, the distribution of cpSSR diversity within the islands of Tenerife and Gran Canaria pointed towards the importance of the role played by regional Pliocene and Quaternary volcanic activity and long-distance gene flow in shaping the population genetic structure of the Canary Island pine. Therefore, conservation strategies at the population level are strongly recommended for this species.Communicated by D.B. NealeA. Gómez and S.C. González-Martínez as joint authors     

Genetic,morphological, and dietary changes associated with novel habitat colonisation in the Canary Island endemic grasshopper Acrostira bellamyi

1. The large flightless grasshopper Acrostira bellamyi Uvarov, endemic to the island of La Gomera (Canary Islands), inhabits two different environments: the xeric euphorb shrubland, as is typical for congeneric Canarian species, and the humid laurel forest, a novel habitat for the genus. 2. We investigate genetic, morphological, and ecological variation among individuals of A. bellamyi from the two habitats. DNA sequence data were used to evaluate whether grasshoppers from the two environments represent distinct lineages. Morphological and trophic analyses were performed to assess phenotypic differentiation between the two different habitats. 3. Population genetic analyses support the hypothesis that the euphorb shrubland is the ancestral habitat for this species. Female laurel forest specimens are larger than those inhabiting the euphorb shrubland, and some external body parts exhibit significant morphometric differences between the two populations. Diet of shrubland individuals is completely different from that of laurel forest individuals. Although in each habitat they consume the most abundant plants, individuals are able to select food plants, which appear to be explained by their nutrient content. 4. Our results suggest that A. bellamyi has colonised laurel forest from shrubland, and that this habitat shift has resulted in genetic, morphological, and ecological changes, perhaps as an adaptation to this new habitat.     

Population structure and history of a phenotypically variable teiid lizard (Ameiva chrysolaema) from Hispaniola: the influence of a geologically complex island

Ameiva chrysolaema is distributed across the island of Hispaniola in the West Indies. The species is restricted to dry lowlands between major mountain ranges and along the southern and eastern coasts. Phylogenetic and phylogeographic analyses of mtDNA sequence variation from 14 sampling localities identify at least three independent evolutionary lineages, separated from one another by major mountain ranges. Nested clade phylogeographic analysis (NCPA) suggests a complex history of population fragmentation, consistent with geological evidence of seawater incursions into the Azua and Enriquillo basins during the Pliocene/Pleistocene (approximately 1.6 mya). Significantly negative Fu's F(S) values and parameters of mismatch distributions suggest that formerly fragmented populations have recently expanded their ranges. Significantly large average population clade distances (APCD) for two sampling localities in the Azua basin suggest secondary contact at these localities of previously separated populations. The distribution of haplotypes among polymorphic populations of A. chrysolaema suggests that variation in dorsal pattern represents a polymorphism within evolutionary lineages. Ameiva leberi is ecologically indistinguishable from and syntopic with A. chrysolaema. Genetic data suggest that A. leberi is a junior synonym of A. chrysolaema.     

Population structure in the endangered cyprinid fish Pararasbora moltrechti in Taiwan, based on mitochondrial and microsatellite DNAs

The genetic structure of four populations of Pararasbora moltrechti, an endemic species of the Cyprinidae in Taiwan, was investigated based on the genetic variation of mtDNA Cyt b gene and five microsatellite loci. High haplotype diversity (h = 0.92) but low nucleotide diversity (0.004) in mtDNA was detected in this endangered species. In total, 33 haplotypes and four clusters were identified in its mtDNA. Nevertheless, low correspondence was found between geographical division and mtDNA clusters. In contrast, Bayesian cluster analysis of the microsatellite data identified four genetic groups and revealed highly structured populations. Significantly negative Tajima's D statistics and mismatch distribution analyses suggest that P. moltrechti populations may have experienced a demographic expansion. In light of the results of a nested clade analysis of mtDNA haplotypes, we conclude that recent population fluctuations and restricted gene flow played major roles in shaping the spatial genetic structure of P. moltrechti populations.     

Intraspecific Phylogeography of Red Squirrels (Tamiasciurus hudsonicus) in the Central Rocky Mountain Region of North America

We used variation in a portion of the mitochondrial DNA control region to examine phylogeography of Tamiasciurus hudsonicus, a boreal-adapted small mammal in the central Rocky Mountain region. AMOVA revealed that 65.66% of genetic diversity was attributable to variation within populations, 16.93% to variation among populations on different mountain ranges, and 17.41% to variation among populations within mountain ranges. Nested clade analysis revealed two major clades that likely diverged in allopatry during the Pleistocene: a southern clade from southern Colorado and a northern clade comprising northern Colorado, Wyoming, eastern Utah, and eastern Idaho. Historically restricted gene flow as a result of geographic barriers was indicated between populations on opposite sides of the Green River and Wyoming Basin and among populations in eastern Wyoming. In some instances genetic structure indicated isolation by distance.     

Dispersal, habitat differences, and comparative phylogeography of Southeast Asian seahorses (Syngnathidae: Hippocampus)

Four distinct phylogeographical patterns across Southeast Asia were observed for four species of seahorse (genus Hippocampus) with differing ecologies. For all species, genetic differentiation (based on cytochrome b sequence comparisons) was significantly associated with sample site (Phi(ST) = 0.190-0.810, P < 0.0001) and with geographical distance (Mantel's r = 0.37-0.59, P < 0.019). Geographic locations of genetic breaks were inconsistent across species in 7/10 comparisons, although some similarities across species were also observed. The two shallow-water species (Hippocampus barbouri and Hippocampus kuda) have colonized the Sunda Shelf to a lesser degree than the two deeper-water species (Hippocampus spinosissimus and Hippocampus trimaculatus). In all species the presence of geographically restricted haplotypes in the Philippines could indicate past population fragmentation and/or long-distance colonization. A nested clade analysis (NCA) revealed that long-distance colonization and/or fragmentation were likely the dominant forces that structure populations of the two shallow-water species, whereas range expansion and restricted dispersal with isolation by distance were proportionally more important in the history of the two deeper-water species. H. trimaculatus has the most widespread haplotypes [average clade distance (D(c)) of nonsingleton haplotypes = 1169 km], indicating potentially high dispersal capabilities, whereas H. barbouri has the least widespread haplotypes (average D(c) = 67 km) indicating potentially lower dispersal capabilities. Pleistocene separation of marine basins and postglacial flooding of the Sunda Shelf are extrinsic factors likely to have contributed to the phylogeographical structure observed, whereas differences among the species appear to reflect their individual ecologies.