Mitochondrial DNA evolution and population history of the Tenerife skink Chalcides viridanus

Recent studies of island lizards have suggested that historical vicariance as a result of volcanism may have played an important role in shaping patterns of within-island genetic diversity. The skink, Chalcides viridanus, shows variation in morphology within the volcanic island of Tenerife. Two mitochondrial DNA (mtDNA) fragments (from the 12S and 16S rRNA regions) were sequenced in individuals from 17 sites to evaluate the relationship between current phylogeography and the geological history of the island. Three main clades were detected. The two most basal clades were restricted to areas representing the ancient precursor islands of Teno and Anaga in the northwest and northeast of Tenerife, respectively. The third clade showed a widespread geographical distribution and provided evidence of a recent rapid expansion after a bottleneck. Within-island cladogenesis appears to have taken place during a recent period of volcanic activity and long after the ancient islands had been united by the eruptions that led to the formation of the Canadas edifice. Evidence of similar biogeographical histories are found in other species in the Canary archipelago, supporting the volcanism scenario as a potentially widespread cause of within-island differentiation in reptiles.     

Palaeo‐islands as refugia and sources of genetic diversity within volcanic archipelagos: the case of the widespread endemic Canarina canariensis (Campanulaceae)

Geographical isolation by oceanic barriers and climatic stability has been postulated as some of the main factors driving diversification within volcanic archipelagos. However, few studies have focused on the effect that catastrophic volcanic events have had on patterns of within‐island differentiation in geological time. This study employed data from the chloroplast (cpDNA haplotypes) and the nuclear (AFLPs) genomes to examine the patterns of genetic variation in Canarina canariensis, an iconic plant species associated with the endemic laurel forest of the Canary Islands. We found a strong geographical population structure, with a first divergence around 0.8 Ma that has Tenerife as its central axis and divides Canarian populations into eastern and western clades. Genetic diversity was greatest in the geologically stable ‘palaeo‐islands’ of Anaga, Teno and Roque del Conde; these areas were also inferred as the ancestral location of migrant alleles towards other disturbed areas within Tenerife or the nearby islands using a Bayesian approach to phylogeographical clustering. Oceanic barriers, in contrast, appear to have played a lesser role in structuring genetic variation, with intra‐island levels of genetic diversity larger than those between‐islands. We argue that volcanic eruptions and landslides after the merging of the palaeo‐islands 3.5 Ma played key roles in generating genetic boundaries within Tenerife, with the palaeo‐islands acting as refugia against extinction, and as cradles and sources of genetic diversity to other areas within the archipelago.     

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.     

Patterns and causes of morphological population differentiation in the Tenerife skink, Chalcides viridanus

Recent studies indicate that differential selection between mesic and xeric habitats is the most plausible explanation of within-island geographic variation in the Gran Ganarian skink, Chalcides sexlineatus. The island of Tenerife shows mesic/xeric heterogeneity similar to that on Gran Canaria, so we tested the prediction of parallel geographic patterns of morphological population differentiation in the Tenerife skink, Chalcides viridanus. Geographic variation was found to be complex; patterns of differentiation in body dimensions and scalation show evidence of both mosaic and latitudinal facets. Using randomization and regression methods, significant relationships were found between a hypothesized xeric/mesic model of population differentiation and geographic variation in female body dimensions, as well as with substantial proportions of the geographic variation in male and female scalation and, to a lesser extent, body dimensions. Matrix comparisons based on the entire morphological distance matrices show that the patterns of variation are not concordant with an alternative hypothesis of historical secondary contact, or with an isolation-by-distance model. It is argued that a substantial component of the geographic variation in body dimensions and scalation of C. viridanus is due to differential selection along an aridity gradient.     

GEOGRAPHIC VARIATION WITHIN AN ISLAND: UNIVARIATE AND MULTIVARIATE CONTOURING OF SCALATION,SIZE, AND SHAPE OF THE LIZARD GALLOTIA GALLOTI

Microgeographic variation of the vegetarian lizard Gallotia galloti within the island of Tenerife is described using univariate analysis, correlation, multiple group principal component analysis, canonical analysis, transects, and contours. The size varies locally in a mosaic pattern while head shape shows a WNW–ESE cline in the south. The scalation (scale and femoral pore counts) has two facets to its geographic variation, both of which are incongruent with the primary variation in the size and shape. The scalation shows categorical variation (stepped cline) between northern and southern populations and also a strong clinal relationship with altitude. The possible causes of this variation are considered and adaptation to current ecological conditions appears to be implicated for at least the altitudinal variation, although the pertinent factors are not obvious. It is apparent that vertebrate populations distributed across small islands do not necessarily offer the advantage of a discrete homogeneous unit for evolutionary studies but can offer the opportunity for studying microgeographic variation.     

Evolution of Synthyris sect. Dissecta (Plantaginaceae) on sky islands in the Northern Rocky Mountains

We examined the phylogeography of Synthyris sect. Dissecta (Plantaginaceae), which is restricted to the Olympic (S. lanuginosa) and Northern Rocky Mountains (S. dissecta and S. canbyi), to infer effects of sky islands and modes of speciation. Sequences of cpDNA trnT-trnL and psbA-trnH intergenic spacers resolved 22 haplotypes among 185 individuals sampled from 16 populations of the three species. Gene flow in the ancestral lineage and random capture of haplotypes in species lineages of sect. Dissecta have resulted in haplotype clades that are not exclusive to species. Nested clade analysis (NCA) indicates that allopatric fragmentation separated Olympic and Northern Rocky Mountain populations, giving rise to the Olympic endemic S. lanuginosa, which is characterized by unique haplotypes consistent with long temporal isolation. Low haplotype and nucleotide diversity in S. canbyi are consistent with newly founded populations experiencing a bottleneck. Furthermore, we infer S. canbyi evolved as a peripheral isolate of S. dissecta. NCA indicated limited migration in S. dissecta with possible isolation by distance. Both isolation on interglacial sky islands and valley glaciers during at least the last glaciation limited gene flow among populations of S. dissecta in different ranges of the Northern Rocky Mountains.     

Morphological evolution in the day gecko Phelsuma sundbergi in the Seychelles: a multivariate study

Geographic variation was investigated in populations of the day gecko Phelsuma sundbergi from 22 islands in the Seychelles, using multivariate ordination procedures. Multiple group principal components analysis was used to negate ontogenetic variation. Seventy-eight characters from three character systems (body proportions, scalation and colour pattern) were analysed from 349 specimens. Three phenetic aggregations of granitic island populations were detected from the northwestern, northeastern and southern island groups. A comparison of the results from the three character systems analysed separately suggests that the separation of the southern form preceded that of the northern forms. Colour pattern characters reflect this closely, while body proportions and scalation characters follow evolution of body size. The populations of Phelsuma sundbergi on the coralline islands are not well differentiated from the races on the granitic islands, and probably represent recent colonizations or introductions by man.     

The dynamics of genetic and morphological variation on volcanic islands

Oceanic archipelagos of volcanic origin have been important in the study of evolution because they provide repeated natural experiments allowing rigorous tests of evolutionary hypotheses. Ongoing volcanism on these islands may, however, affect the evolutionary diversification of species. Analysis of population structure and phylogeographic patterns in island populations can provide insight into evolutionary dynamics on volcanic islands. We analysed genetic and morphological variation in the gecko Tarentola boettgeri on the island of Gran Canaria and compared it with Tarentola delalandii on Tenerife, a neighbouring volcanic island of similar age but distinctly different geological past. Intraspecific divergence of mitochondrial haplotypes indicates long-term persistence of Tarentola on each island, with a phylogeographic signal left by older volcanic events. More recent volcanic eruptions (approximately 0.2 million years ago on Tenerife, approximately 2.2 million years ago on Gran Canaria) have left a signature of population expansion in the population genetic structure, the strength of which depends on the time since the last major volcanic eruption on each island. While these stochastic events have left traces in morphological variation in Tenerife, in Gran Canaria geographical variation was solely associated with environmental variables. This suggests that historically caused patterns in morphology may be overwritten by natural selection within 2 million years.     

Geographic variation in scalation of the lizard Gallotia stehlini within the island of Gran Canaria

The lacertid lizard ( Gallotia stehlini ), an endemic of Gran Canaria, shows no visually obvious geographic variation, yet all seven scalation characters that were examined exhibit significant geographic variation. The number of collar scales, scales along the ventral trunk and femoral pores are correlated with habitat type, while the number of femoral pores is correlated negatively with altitude. Mantel tests were used to compare simultaneously an observed pattern with three hypothesized patterns (habitat type, altitude and proximity). They indicate that, while several individual characters are significantly associated with the putative causal factors of habitat type and altitude, there is no association between an overall scalation distance matrix and habitat type, or altitude when the effect of proximity is removed. Consequently, one should consider the individual characters as well as the multivariate generalized distances. Some of the observed patterns of geographic variation in scalation are very similar to those of the small scincid lizard Chalcides sexlineatus on Gran Canaria and also parallel the altitudinal and latitudinal variation in the scalation of the Tenerife lacertid ( Galotia galloti ). The low level of congruence in patterns of geographic variation in individual characters (i.e. some vary with latitude, some with altitude and one varies with longitude) is consistent with the hypothesis that ecogenetically caused geographic variation may result in lower inter-character congruence than phylogenetically caused geographic variation.     

Recent volcanism and mitochondrial DNA structuring in the lizard Gallotia atlantica from the island of Lanzarote

The phylogeography of the lacertid lizard Gallotia atlantica from the small volcanic island of Lanzarote (Canary Islands) was analysed based on 1075 bp of mitochondrial DNA (mtDNA) sequence (partial cytochrome b and ND2) for 157 individuals from 27 sites (including three sites from neighbouring islets). Levels of sequence divergence were generally low, with the most distant haplotypes separated by only 14 mutational steps. MtDNA divergence appears to coincide with formation of the middle Pleistocene lowland that united formerly separate ancient islands to form the current island of Lanzarote, allowing rejection of a two-island model of phylogeographical structure. There was evidence of large-scale population expansion after island unification, consistent with the colonization of new areas. A nested clade phylogeographical analysis (NCPA) revealed significant phylogeographical structuring. Two-step and higher-level clades each had disjunct distributions, being found to the east and west of a common area with a north-south orientation that extends between coasts in the centre-east of the island (El Jable). Other clades were almost entirely restricted to the El Jable region alone. Bayesian Markov chain Monte Carlo analyses were used to separate ongoing gene flow from historical associations. These supported the NCPA by indicating recent (75,000-150,000 years ago) east-west vicariance across the El Jable region. Lava flows covered El Jable and other parts of the central lowland at this time and likely led to population extinctions and temporary dispersal barriers, although present-day evidence suggests some populations would have survived in small refugia. Expansion of the latter appears to explain the presence of a clade located between the eastern and western components of the disjunct clades. Direct relationships between mtDNA lineages and morphology were not found, although one of two morphological forms on the island has a disjunct distribution that is broadly concordant with east-west components of the phylogeographical pattern. This work demonstrates how recent volcanic activity can cause population fragmentation and thus shape genetic diversity on microgeographical scales.     

Trans-oceanic host dispersal explains high seabird tick diversity on Cape Verde islands

Parasites represent ideal models for unravelling biogeographic patterns and mechanisms of diversification on islands. Both host-mediated dispersal and within-island adaptation can shape parasite island assemblages. In this study, we examined patterns of genetic diversity and structure of Ornithodoros seabird ticks within the Cape Verde Archipelago in relation to their global phylogeography. Contrary to expectations, ticks from multiple, geographically distant clades mixed within the archipelago. Trans-oceanic colonization via host movements probably explains high local tick diversity, contrasting with previous research that suggests little large-scale dispersal in these birds. Although host specificity was not obvious at a global scale, host-associated genetic structure was found within Cape Verde colonies, indicating that post-colonization adaptation to specific hosts probably occurs. These results highlight the role of host metapopulation dynamics in the evolutionary ecology and epidemiology of avian parasites and pathogens.     

The phylogeography of the gecko Cyrtopodion kotschyi (Reptilia: Gekkonidae) in the Aegean archipelago

Cyrtopodion kotschyi is a small gecko, widely distributed on the islands of the Aegean archipelago and the adjacent mainland. We unveiled the phylogeography of the species by using a portion of the cytochrome oxidase I mitochondrial DNA gene from 35 insular and mainland populations. The distinct geographic distribution of the major clades of the phylogenetic tree and its topology suggest a spatial and temporal sequence of phylogenetic separations that coincide with some major paleogeographical separations during the geological history of the Aegean and support a mainly vicariant pattern of differentiation. The separation times and 95% confidence limits among the different clades were estimated according to two different paleogeographical scenarios. The very high interpopulation genetic divergence (up to 20% uncorrected pairwise distances) and the better concordance between paleogeographical and phylogenetic separations for one of the scenarios suggest that species differentiation started during Miocene (about 10 Mya) due to the fragmentation of the united landmass of "Agais" that was Aegean at that time.     

Phylogeography and natural selection in the Tenerife gecko Tarentola delalandii: testing historical and adaptive hypotheses

Combining phylogeographic and matrix correspondence approaches in the analysis of geographical variation provides a fruitful approach to inferring the causes of molecular and morphological evolution within species. Here we present a study on the gecko Tarentola delalandii on the island of Tenerife, Canary Islands, which provides an outstanding model of an exceptionally high degree of phylogeographic differentiation in magnitude and pattern on a small spatial scale. We reconstruct the population history of T. delalandii using phylogeographic information, matrix correspondence tests and estimates of divergence times in conjunction with geological data. It appears that populations differentiated on three precursor islands and secondary contact followed the junction of these islands. The cytochrome b sequence appears to be evolving at least at approximately 1% per million years in this species. Matrix correspondence tests indicate that morphological character systems may reflect ecological selection regimes (colour pattern), history (body dimensions) or both (scalation). The results imply that natural selection can override a historical legacy, but also underline the potential relevance of molecular phylogenetic data for the interpretation of geographical variation in morphology.     

Phylogeography of a morphologically diverse Neotropical montane species, the Common Bush-Tanager (Chlorospingus ophthalmicus)

The Common Bush-Tanager (Chlorospingus ophthalmicus) is distributed in Neotropical cloud-forests from Mexico to Argentina and contains 25 subspecies divided into eight subspecies groups based on biogeography, eye coloration, presence of a postocular spot and chest band. All of Central America is occupied by a single subspecies group; whereas the Andes are believed to be occupied by seven additional subspecies groups. We used five mitochondrial genes to investigate the phylogeography and possible species limits of the ophthalmicus complex. A total of 14 monophyletic lineages were uncovered within the ophthalmicus complex, including three clades currently classified as separate species (C. semifuscus, inornatus and tacarcunae). Divergence estimates for these clades date between 0.8 and 5.2 million years ago (Ma). Contrary to expectations based on morphological diversity, phylogeographic structure was greatest in Mexico and Central America and weakest in the Andes. Morphological and genetic divergences were not significantly correlated and most morphologically defined subspecies groups were not supported. Our evidence suggests the ophthalmicus complex originated in Mexico ca. 6.0 Ma (million years ago) and spread south into the Andes ca. 4.7 Ma before the completion of the Isthmus of Panama. Three genetically divergent lineages of ophthalmicus that formed in the Andes possess a complex checkerboard distribution, with a single lineage represented by disjunct populations from Venezuela and the southern Andes, while intervening populations in Ecuador and Central Peru form two genetically and morphologically divergent lineages.     

Genetic structure, phylogeography and demography of two ground-beetle species endemic to the Tenerife laurel forest (Canary Islands)

The volcanic island of Tenerife (Canary archipelago) was formerly covered at 600-1200 m above sea level on most of its northern side by a cloud forest holding much of the endemic insect fauna. In the most significant surviving patches of this laurel forest at the eastern and western tips of the island occur two forest-specialist, closely related species of Eutrichopus (Coleoptera, Carabidae); here we present data on mitochondrial DNA variation among populations of these species. In total, 116 individuals from 16 localities were sampled and a 638 bp fragment of the cytochrome oxidase subunit II gene was sequenced, obtaining evidence for two distinct evolutionary lineages, in accordance with morphological and biogeographical data. Volcanic events at approximately 0.7 Ma might be responsible for vicariance and the fragmentation of the geographical range of an ancestral species, causing the establishment of two matrilineal lineages. Using nested clade and historical demography analyses we infer past cycles of demographic bottlenecks followed by population expansion, mostly in agreement with the geological time scale of volcanic events. Recent trends, however, refer to fragmentation of the cloud forest due to human intervention.     

Speciation chronology of rockhopper penguins inferred from molecular, geological and palaeoceanographic data

Aim The Southern Ocean is split into several biogeographical provinces between convergence zones that separate watermasses of different temperatures. Recent molecular phylogenies have uncovered a strong phylogeographic structure among rockhopper penguin populations, Eudyptes chrysocome sensu lato, from different biogeographical provinces. These studies suggested a reclassification as three species in two major clades, corresponding, respectively, to warm, subtropical and cold sub‐Antarctic watermasses rather than to geographic proximity. Such a phylogeographic pattern, also observed in plants, invertebrates and fishes of the Southern Ocean, suggests that past changes in the positions of watermasses may have affected the evolutionary history of penguins. We calculated divergence times among various rockhopper penguin clades and calibrated these data with palaeomagmatic and palaeoceanographic events to generate a speciation chronology in rockhopper penguins. Location Southern Ocean. Methods Divergence times between populations were calculated using five distinct mitochondrial DNA loci, and assuming a molecular clock model as implemented in mdiv . The molecular evolution rate of rockhopper penguins was calibrated using the radiochronological age of St Paul Island and Amsterdam Island in the southern Indian Ocean. Separations within other clades were correlated with palaeoceanographic data using this calibrated rate. Results The split between the Atlantic and Indian populations of rockhopper penguins was dated as 0.25 Ma, using the date of emergence of St Paul and Amsterdam islands, and the divergence between sub‐Antarctic and subtropical rockhopper penguins was dated as c. 0.9 Ma (i.e. during the mid‐Pleistocene transition, a major change in the Earth’s climate cycles). Main conclusions The mid‐Pleistocene transition is known to have caused a major southward shift in watermasses in the Southern Ocean, thus changing the environment around the northernmost rockhopper penguin breeding sites. This ecological isolation of northernmost populations may have caused vicariant speciation, splitting the species into two major clades. After the emergence of St Paul and Amsterdam islands in the subtropical Indian Ocean 0.25 Ma, these islands were colonized by penguins from the subtropical Atlantic, 6000 km away, rather than by penguins from the sub‐Antarctic Indian Ocean, 5000 km closer.     

Status and relationships of the extinct giant Canary Island lizard Gallotia goliath (Reptilia: Lacertidae), assessed using ancient mtDNA from its mummified remains

Ancient mitochondrial DNA sequences (378 base pairs of cytochrome b and 368 of 12S rRNA) extracted from a mummified extinct giant lizard, Gallotia goliath , from eastern Tenerife, Canary Islands, were used to assess the species status and relationship of this form within the genus. G. goliath is clearly a member of the G. simonyi group of the western Canary islands (Tenerife, La Gomera, El Hierro and La Palma) and is not closely related to the giant G. stehlini of Gran Canaria. Contrary to recent opinion, it is phylogenetically distinct, within the G. simonyi group, from the extant G. simonyi of El Hierro and also from the recently discovered live G. gomerana on La Gomera and from G. intermedia in north-western Tenerife. It may be the sister taxon of either all the other members of the G. simonyi group or of G. intermedia . The phylogenetic distinctness of G. goliath makes Tenerife unique among oceanic islands in having had one giant and two medium-sized lizard species that were probably substantially herbivorous, the others being G. intermedia and G. galloti . Gallotia shows great community differences on other islands in the Canaries, two having a single small species, one a single giant, and three a giant and a medium-sized form. © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 80 , 659–670.     

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     

Recent origin and adaptive diversification of <Emphasis Type="Italic">Ainsliaea</Emphasis> (Asteraceae) in the Ryukyu Islands: molecular phylogenetic inference using nuclear microsatellite markers

We investigated the phylogeny and phylogeography of four closely related Ainsliaea species (Asteraceae) on the continental Ryukyu Islands of eastern Asia, which consist of two flood-adapted “rheophyte” and two non-rheophyte (inland) species, based on 12 nuclear microsatellite loci. Phylogenetic analyses using 420 individuals from 26 populations showed that rheophytic A. linearis and A. oblonga are genetically distinct. Each species was clustered with the inland species that occur on the same islands, suggesting a different ancestry for the two rheophytes that evolved independently by local adaptation to flooded habitats. The results from the neighbor-joining clustering and principal coordinate analysis (PCoA) indicate that the southern populations of A. macroclinidioides are distinct lineages and ancestral to the northern populations as well as the other diverse species complex in the Ryukyus. These results suggest a pattern of colonization initially from the Asian mainland to the southern islands, followed by the northern islands via land bridges generated during the Quaternary glaciations. After isolation from southern populations, species radiation and regional differentiation within the northern clade occurred possibly via local adaptation and/or geographic isolation of the subdivided island bridge.     

Phylogeography and species biogeography of montane Great Basin stoneflies

Sky islands are ideal systems for determining the effects of climatic oscillations on species distributions and genetic structure. Our study focused on montane stonefly populations in the Great Basin of western North America. We used niche-based distribution modelling, phylogeography and traditional species-based biogeography to test several hypotheses as follows: (i) genetic differentiation among Doroneuria baumanni populations will be independent of hydrologic connectivity (headwater model); (ii) Sky islands were colonized when habitat was more continuous and populations likely experienced multiple expansions and contractions; (iii) Colonization events were coincident with the late Pleistocene and Holocene; and (iv) Shared topography and climate history will result in concordant patterns of genetic differentiation in D. baumanni and occurrences of 32 stonefly species across the region. Overall, Φ(ST) 's and coalescent-based estimates of migration were consistent with the headwater model. Maximum likelihood and Bayesian gene trees identified three major nonoverlapping east-west clades. Distribution modelling indicated more suitable habitat in the Great Basin during the Last Glacial Maximum than at present, but none during the last interglacial period. Demographic analyses showed evidence of population expansion in one of the three major east-west clades. Intra-clade divergence times (60,000-183,000ybp) were well within the late Pleistocene while among-clade divergence times (499.000-719,000ybp) were deeper. Genetic differentiation in D. baumanni and distributions of stonefly species were significantly concordant. These results imply that climatic oscillations have played major roles in shaping the genetic structure and distributions of Great Basin stoneflies, but that divergence among clades occurred much earlier than our late Pleistocence/early Holocene predictions.