Adaptive radiation in miniature: the minute salamanders of the Mexican highlands (Amphibia: Plethodontidae: Thorius)

The small size and apparent external morphological similarity of the minute salamanders of the genus Thorius have long hindered evolutionary studies of the group. We estimate gene and species trees within the genus using mitochondrial and nuclear DNA from nearly all named and many candidate species and find three main clades. We use this phylogenetic hypothesis to examine patterns of morphological evolution and species coexistence across central and southern Mexico and to test alternative hypotheses of lineage divergence with and without ecomorphological divergence. Sympatric species differ in body size more than expected after accounting for phylogenetic relationship, and morphological traits show no significant phylogenetic signal. Sympatric species tend to differ in a combination of body size, presence or absence of maxillary teeth, and relative limb or tail length, even when they are close relatives. Sister species of Thorius tend to occupy climatically similar environments, which suggests that divergence across climatic gradients does not drive species formation in the genus. Rather than being an example of cryptic species formation, Thorius more closely resembles an adaptive radiation, with ecomorphological divergence that is bounded by organism‐level constraints. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 622–643.     

Discordant mitochondrial and nuclear gene phylogenies in emydid turtles: implications for speciation and conservation

Do phylogenies and branch lengths based on mitochondrial DNA (mtDNA) provide a reasonable approximation to those based on multiple nuclear loci? In the present study, we show widespread discordance between phylogenies based on mtDNA (two genes) and nuclear DNA (nucDNA; six loci) in a phylogenetic analysis of the turtle family Emydidae. We also find an unusual type of discordance involving the unexpected homogeneity of mtDNA sequences across species within genera. Of the 36 clades in the combined nucDNA phylogeny, 24 are contradicted by the mtDNA phylogeny, and six are strongly contested by each data set. Two genera (Graptemys, Pseudemys) show remarkably low mtDNA divergence among species, whereas the combined nuclear data show deep divergences and (for Pseudemys) strongly supported clades. These latter results suggest that the mitochondrial data alone are highly misleading about the rate of speciation in these genera and also about the species status of endangered Graptemys and Pseudemys species. In addition, despite a strongly supported phylogeny from the combined nuclear genes, we find extensive discordance between this tree and individual nuclear gene trees. Overall, the results obtained illustrate the potential dangers of making inferences about phylogeny, speciation, divergence times, and conservation from mtDNA data alone (or even from single nuclear genes), and suggest the benefits of using large numbers of unlinked nuclear loci. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 445–461.     

Increased molecular divergence of two endemic Trachelipus (Isopoda,Oniscidea) species from Greece reveals patterns not congruent with current taxonomy

In the present study, we employed three mitochondrial DNA genetic markers in a phylogenetic analysis aiming at the delineation of the relationships amongst nominal Trachelipus kytherensis populations, as well as between populations of this species and of Trachelipus aegaeus and a new form, occurring syntopically with the latter. Both the phylogenetic analysis and the genetic distances separating populations, show the presence of several distinct and well differentiated clades that undermine the monophyly of T. kytherensis. On the other hand, despite the insular distribution of T. aegaeus populations, their divergence is low and the monophyly of this taxon can be rescued by the inclusion of two more insular populations previously assigned to T. kytherensis. The patterns of genetic divergence among clades are only partially congruent with the geographic distribution of populations. The validity of taxonomic characters used so far in the genus appears to be questionable; therefore, a more comprehensive phylogenetic study at a population level is deemed necessary for understanding the divergence of Trachelipus lineages. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 361–370.     

Molecular patterns of differentiation in canyon treefrogs (Hyla arenicolor): evidence for introgressive hybridization with the Arizona treefrog (H. wrightorum) and correlations with advertisement call differences

Detection of genetic and behavioural diversity within morphologically similar species has led to the discovery of cryptic species complexes. We tested the hypothesis that US populations of the canyon treefrog (Hyla arenicolor) may consist of cryptic species by examining mate‐attraction signals among three divergent clades defined by mtDNA. Using a multi‐locus approach, we re‐analysed phylogenetic relationships among the three clades and a closely related, but morphologically and behaviourally dissimilar species, the Arizona treefrog (H. wrightorum). We found evidence for introgression of H. wrightorum’s mitochondrial genome into H. arenicolor. Additionally, the two‐clade topology based on nuclear data is more congruent with patterns of call variation than the three‐clade topology from the mitochondrial dataset. The magnitude of the call divergence is probably insufficient to promote isolation of the nuclear DNA‐defined clades should they become sympatric, but further divergence in call properties significant in species identification could promote speciation in the future.     

Molecular phylogeny and phylogeography of the Cuban cave-fishes of the genus Lucifuga: evidence for cryptic allopatric diversity

Underground environments are increasingly recognized as reservoirs of faunal diversity. Extreme environmental conditions and limited dispersal ability of underground organisms have been acknowledged as important factors promoting divergence between species and conspecific populations. However, in many instances, there is no correlation between genetic divergence and morphological differentiation. Lucifuga Poey is a stygobiotic fish genus that lives in Cuban and Bahamian caves. In Cuba, it offers a unique opportunity to study the influence of habitat fragmentation on the genetic divergence of stygobiotic species and populations. The genus includes four species and one morphological variant that have contrasting geographical distributions. In this study, we first performed a molecular phylogenetic analysis of the Lucifuga Cuban species using mitochondrial and nuclear markers. The mitochondrial phylogeny revealed three deeply divergent clades that were supported by nuclear and morphological characters. Within two of these main clades, we identified five lineages that are candidate cryptic species and a taxonomical synonymy between Lucifuga subterranea and Lucifuga teresinarum. Secondly, phylogeographic analysis using a fragment of the cytochrome b gene was performed for Lucifuga dentata, the most widely distributed species. We found strong geographical organization of the haplotype clades at different geographic scales that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow. At a larger temporal scale, these processes could also explain the diversification and the distribution of the different species.     

Multilocus phylogeny of a cryptic radiation of Afrotropical long-fingered bats (Chiroptera,Miniopteridae)

The Old World bat family Miniopteridae comprises only the genus Miniopterus, which includes 20 currently recognized species from the Afrotropical realm and 15 species from Eurasia and Australasia. Since 2003, the number of recognized Miniopterus species has grown from 19 to 35, with most newly described species endemic to Madagascar and the Comoros Archipelago. We investigated genetic variation, phylogenetic relationships and clade membership in Miniopterus focusing on Afrotropical taxa. We generated mitochondrial cytochrome-b (cyt-b) and nuclear intron data (five genes) from 352 vouchered individuals collected at 78 georeferenced localities. Including 99 additional mitochondrial sequences from GenBank, we analysed a total of 25 recognized species. Mitochondrial genetic distances among cyt-b-supported clades averaged 9.3%, representing as many as five undescribed species. Multilocus coalescent delimitation strongly supported the genetic isolation of eight of nine tested unnamed clades. A large number of sampled clades in sub-Saharan Africa are distributed wholly or partly in East Africa (nine of 13 clades), suggesting that Miniopterus diversity has been grossly underestimated. Although 25 of 27 cyt-b and 23 of 25 nuclear gene tree lineages from the Afrotropics were strongly supported as monophyletic, a majority of deep nodes were poorly resolved in phylogenetic analyses. Long terminal branches subtending short backbone internodes in the phylogenetic analyses suggest a rapid radiation model of diversification. This hypothesis needs to be tested using more phylogenetically informative data.     

Intraspecific molecular phylogeny and phylogeography of the Meriones meridianus (Rodentia: Cricetidae) complex in northern China reflect the processes of desertification and the Tianshan Mountains uplift

The phylogeographical patterns and demographic history of mitochondrial DNA (cytochrome b, N = 327; D‐loop, N = 252) and nuclear DNA (IRBP gene, N = 235) haplotypes were studied for the Meriones meridianus complex in northern China, a desert‐dwelling gerbil species complex. The phylogenetic analyses, which were performed on the separate and combined (mitochondrial + nuclear) datasets, revealed two divergent clades (Clade A and Clade B) corresponding to distinct geographical regions. Clade A contained the haplotypes found mostly in individuals from the Tianshan Mountains area. Clade B contained haplotypes from populations located in other deserts in northern China. The divergence times indicated that the history of the M. meridianus complex was influenced by the uplift of the Tianshan Mountains and climate‐induced habitat fluctuations. In the Pleistocene, the expansion of forests and grasslands during interglacial period led to the isolation of the M. meridianus complex, which preferred to inhabit deserts. Hence, long geological isolation and the M. meridianus complex adaptation to local ecological conditions led to its genetic divergence. Clade A had long‐lasting demographic stability, most likely because the populations of this clade remained in a stable desert environment for a long time. However, the extension of other deserts and disappearance of palaeolakes during the last glacial period resulted in demographic expansion of Clade B. Furthermore, our genetic data indicated that two subspecies may exist within the M. meridianus complex. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 362–383.     

Endemic evolutionary radiation of Rhagada land snails (Pulmonata: Camaenidae) in a continental archipelago in northern Western Australia

In the Dampier Archipelago of Western Australia's Pilbara Region, several locally endemic, morphologically distinctive species of Rhagada land snails occur, contrasting with the morphologically conservative species with wider distributions on the adjacent mainland. To test alternative origins of this unusual local diversity in a continental archipelago, we examined sequences of the cytochrome oxidase subunit 1 and 16S mitochondrial genes in 22 described species and eight undescribed forms, including all known morphospecies from the Pilbara Region's Dampier Archipelago and adjacent mainland. Phylogenetic analyses consistently resolved four, deep clades within the Pilbara Region, with a mean sequence divergence of 15–18%. All but one of the species from the Dampier Archipelago formed one of the major clades, indicating that the morphological radiation in the archipelago evolved locally, rather than through multiple, relictual mainland lineages. Morphological divergence spanning almost that of the entire genus was within a subclade with sequence divergence < 4%, highlighting the disconnection between morphological diversification and levels of molecular genetic divergence. This in situ morphological radiation in the Dampier Archipelago, which transcends variation seen over much larger distances on the mainland, is unusual for a continental archipelago, and may relate to local heterogeneity of land forms. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 316–327.     

Patterns and processes in the evolutionary history of parrotfishes (Family Labridae)

Phylogenetic reconstruction of the evolutionary relationships among 61 of the 70 species of the parrotfish genera Chlorurus and Scarus (Family Labridae) based on mitochondrial and nuclear gene sequences retrieved 15 well‐supported clades with mid Pliocene/Pleistocene diversification. Twenty‐two reciprocally monophyletic sister‐species pairs were identified: 64% were allopatric, and the remainder were sympatric. Age of divergence was similar for allopatric and sympatric species pairs. Sympatric sister pairs displayed greater divergence in morphology, ecology, and sexually dimorphic colour patterns than did allopatric pairs, suggesting that both genetic drift in allopatric species pairs and ecologically adaptive divergence between members of sympatric pairs have played a role in diversification. Basal species typically have small geographical ranges and are restricted to geographically and ecologically peripheral reef habitats. We found little evidence that a single dominant process has driven diversification, nor did we detect a pattern of discrete, sequential stages of diversification in relation to habitat, ecology, and reproductive biology. The evolution of Chlorurus and Scarus has been complex, involving a number of speciation processes. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 529–557.     

The species identity and biogeography of Blanus (Amphisbaenia: Blanidae) in Lebanon

The genus Blanus Wagler, 1830 represents limbless, burrowing reptiles of the family Blanidae with disjunct circum-Mediterranean distribution. The recently described species Blanus alexandri Sindaco, Kornilios, Sacchi & Lymberakis, 2014 is known from south-eastern Turkey with a presumed occurrence in the Levant and Iraq. We provide here records from Lebanon and confirm the affiliation of Lebanese populations to this species by mitochondrial and nuclear DNA. Blanus alexandri comprises at least seven deeply evolved phylogenetic clades with up to 11% of uncorrected p-distances in their mitochondrial DNA. This suggests a probably older than the Miocene origin of some of these clades. Populations from Lebanon form a different clade with a genetic diversity that is close to populations from southern Turkey.     

Bioacoustic and genetic divergence in a frog with a wide geographical distribution

The study of intraspecific variation of acoustic signals and its relationship with genetic divergence is important for understanding the origin of divergence in communication systems. We studied geographical variation in the acoustic structure of advertisement calls from five populations of the four‐eyed frog, Pleurodema thaul, and its relationship with the genetic divergence among these populations. By analyzing temporal and spectral parameters of the advertisement calls, we report that the signals of northern, central, and southern populations have remarkable differences between them. A phylogeographical analysis from a mitochondrial DNA fragment demonstrated three phylogenetic groups coincident with those found with the bioacoustics analysis. Furthermore, bioacoustic and genetic distances show significant correlations after controlling for geographical distance. These results suggest that behavioural divergence among populations of P. thaul has a phylogenetic basis, supporting three evolutionary units within this species, as well as prompting the exploration of divergence processes in the sound communication system of this species. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 142–155.     

Phylogeography of a pearl oyster (Pinctada maxima) across the Indo‐Australian Archipelago: evidence of strong regional structure and population expansions but no phylogenetic breaks

This study investigates the genetic structure and phylogeography of a broadcast spawning bivalve mollusc, Pinctada maxima, throughout the Indo‐West Pacific and northern Australia. DNA sequence variation of the mitochondrial cytochrome oxidase subunit I (COI) gene was analysed in 367 individuals sampled from nine populations across the Indo‐West Pacific. Hierarchical AMOVA indicated strong genetic structuring amongst populations (Φ_ST = 0.372, P < 0.001); however, sequence divergence between the 47 haplotypes detected was low (maximum 1.8% difference) and no deep phylogenetic divergence was observed. Results suggest the presence of genetic barriers isolating populations of the South China Sea and central Indonesian regions, which, in turn, show patterns of historical separation from northern Australian regions. In P. maxima, historical vicariance during Pleistocene low sea levels is likely to have restricted planktonic larval transport, causing genetic differentiation amongst populations. However, low genetic differentiation is observed where strong ocean currents are present and is most likely due to contemporary larval transport along these pathways. Geographical association with haplotype distributions may indicate signs of early lineage sorting arising from historical population separations, yet an absence of divergent phylogenetic clades related to geography could be the consequence of periodic pulses of high genetic exchange. We compare our results with previous microsatellite DNA analysis of these P. maxima populations, and discuss implications for future conservation management of this species. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 632–646.     

Speciation in the Rana chensinensis species complex and its relationship to the uplift of the Qinghai-Tibetan Plateau

Speciation remains a fundamental issue in biology. Herein, we report an investigation into speciation in the Rana chensinensis species complex using DNA sequence data from one mitochondrial and five nuclear genes. A phylogenetic analysis of the data revealed four major clades in the complex, and each of them was found to likely represent a species, including one cryptic species. Ecological niche models were generated from 19 climatic variables for three of the four major clades, which were represented by widespread sampling, including R. chensinensis, Rana kukunoris and the potential cryptic species. Each clade is associated with a unique ecological unit, and this indicates that ecological divergence probably drove speciation. Ecological divergence is likely related to the late Cenozoic orogenesis of the Qinghai–Tibetan Plateau. In addition, gene flow between species was detected but only in peripheral portions of the ranges of the four major clades, thus likely had little influence on the speciation processes. Discordances between mitochondrial and nuclear genes were also found; the nominal species, R. chensinensis, contains multiple maternal clades, suggesting potential mitochondrial introgression between R. chensinensis and R. kukunoris.     

Phylogeography of two cryptic species of African desert jerboas (Dipodidae: Jaculus)

The lesser Egyptian jerboa Jaculus jaculus is a desert dwelling rodent that inhabits a broad Arabian–Saharan arid zone. Recently, two distant sympatric lineages were described in North‐West Africa, based on morphometric and molecular data, which may correspond to two cryptic species. In the current study, phylogenetic relationships and phylogeographical structure among those lineages and geographical populations from North Africa and the Middle East were investigated. The phylogeographical patterns and genetic diversity of the cytochrome b gene (1110 bp) were addressed on 111 jerboas from 41 localities. We found that the variation in Africa is partitioned into two divergent mitochondrial clades (10.5% divergence relating to 1.65–4.92 Mya) that corresponds to the two cryptic species: J. jaculus and J. deserti. Diversifications within those cryptic species/clades were dated to 0.23–1.13 Mya, suggesting that the Middle Pleistocene climatic change and its environmental consequences affected the evolutionary history of African jerboas. The third distant clade detected, found in the Middle East region, most likely represents a distinct evolutionary unit, independent of the two African lineages. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

Phylogeography of the Australian freshwater turtle Chelodina expansa reveals complex relationships among inland and coastal bioregions

We examined range‐wide mitochondrial phylogeographical structure in the riverine freshwater turtle Chelodina expansa to determine whether this species exhibits deep genetic divergence between coastal and inland hydrological provinces, as seen in co‐distributed freshwater taxa. We sequenced two mitochondrial loci, genealogical relationships were assessed using a network approach, and relationships among biogeographical regions were tested using analyses of molecular variance. Population history was evaluated using neutrality tests, indices of demographic expansion, and mismatch analyses. Twenty‐one haplotypes were recovered across two mitochondrial haplogroups separated by approximately 4% nucleotide divergence. The haplogroups have discrete geographical boundaries but only partially support a hypothesis of deep divergence between coastal and inland bioregions. The first haplogroup comprises populations from the inland Murray‐Darling Basin and from coastal catchments south of the Mary River in south‐east Queensland. The second haplogroup comprises populations from coastal catchments north of the Mary River. Cryptic phylogeographical barriers separating adjacent coastal populations are congruent with those demonstrated for other freshwater taxa and may result from the combined influences of the Conondale Range and alluvial deposits at the mouth of the Mary River. The findings of the present study demonstrate that freshwater taxa commonly display genetic differentiation within a biogeographical region where no boundaries have been recognized, highlighting the need to uncover cryptic microbiogeographical regions to aid conservation of freshwater biota. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 789–805.     

Phylogeographic diversification of antelope squirrels (Ammospermophilus) across North American deserts

We investigated the biogeographic history of antelope squirrels, genus Ammospermophilus, which are widely distributed across the deserts and other arid lands of western North America. We combined range‐wide sampling of all currently recognized species of Ammospermophilus with a multilocus data set to infer phylogenetic relationships. We then estimated divergence times within identified clades of Ammospermophilus using fossil‐calibrated and rate‐calibrated molecular clocks. Lastly, we explored generalized distributional changes of Ammospermophilus since the last glacial maximum using species distribution models, and assessed responses to Quaternary climate change by generating demographic parameter estimates for the three wide‐ranging clades of A. leucurus. From our phylogenetic estimates we inferred strong phylogeographic structure within Ammospermophilus and the presence of three well‐supported major clades. Initial patterns of historical divergence were coincident with dynamic alterations in the landscape of western North America, and the formation of regional deserts during the Late Miocene and Pliocene. Species distribution models and demographic parameter estimates revealed patterns of recent population expansion in response to glacial retreat. When combined with evidence from co‐distributed taxa, the historical biogeography of Ammospermophilus provides additional insight into the mechanisms that impacted diversification of arid‐adapted taxa across the arid lands of western North America. We propose species recognition of populations of the southern Baja California peninsula to best represent our current understanding of evolutionary relationships among genetic units of Ammospermophilus. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 949–967.     

Gene flow and Andean uplift shape the diversification of Gasteracantha cancriformis (Araneae: Araneidae) in Northern South America

The Andean uplift has played a major role in shaping the current Neotropical biodiversity. However, in arthropods other than butterflies, little is known about how this geographic barrier has impacted species historical diversification. Here, we examined the phylogeography of the widespread color polymorphic spider Gasteracantha cancriformis to evaluate the effect of the northern Andean uplift on its divergence and assess whether its diversification occurred in the presence of gene flow. We inferred phylogenetic relationships and divergence times in G. cancriformis using mitochondrial and nuclear data from 105 individuals in northern South America. Genetic diversity, divergence, and population structure were quantified. We also compared multiple demographic scenarios for this species using a model‐based approach (Phrapl ) to determine divergence with or without gene flow. At last, we evaluated the association between genetic variation and color polymorphism. Both nuclear and mitochondrial data supported two well‐differentiated clades, which correspond to populations occurring on opposite sides of the Eastern cordillera of the Colombian Andes. The final uplift of this cordillera was identified as the most likely force that shaped the diversification of G. cancriformis in northern South America, resulting in a cis‐ and trans‐Andean phylogeographic structure for the species. We also found shared genetic variation between the cis‐ and trans‐Andean clades, which is better explained by a scenario of historical divergence in the face of gene flow. This has been likely facilitated by the presence of low‐elevation passes across the Eastern Colombian cordillera. Our work constitutes the first example in which the Andean uplift coupled with gene flow influenced the evolutionary history of an arachnid lineage.     

Contemporary genetic structure of an endemic freshwater turtle reflects Miocene orogenesis of New Guinea

The island of New Guinea lies in one of the most tectonically active regions in the world and has long provided outstanding opportunity for studies of biogeography. Several chelid turtles, of clear Gondwanal origin, occur in New Guinea; all species except one, the endemic Elseya novaeguineae, are restricted to the lowlands south of the Central Ranges. Elseya novaeguineae is found throughout New Guinea. We use mitochondrial and nuclear gene variation among populations of E. novaeguineae throughout its range to test hypotheses of recent extensive dispersal versus more ancient persistence in New Guinea. Its genetic structure bears the signature of Miocene vicariance events. The date of the divergence between a Birds Head (Kepala Burung) clade and clades north and south of the Central Ranges is estimated to be 19.8 Mya [95% highest posterior density (HPD) interval of 13.3–26.8 Mya] and the date between the northern and southern clades is estimated to be slightly more recent at 17.4 Mya (95% HPD interval of 11.0–24.5 Mya). The distribution of this endemic species is best explained by persistent occupation (or early invasion and dispersal) and subsequent isolation initiated by the dramatic landform changes that were part of the Miocene history of the island of New Guinea, rather than as a response to the contemporary landscape of an exceptionally effective disperser. The driving influence on genetic structure appears to have been isolation arising from a combination of: (1) the early uplift of the Central Ranges and establishment of a north‐south drainage divide; (2) development of the Langguru Fold Belt; (3) the opening of Cenderawasih Bay; and (4) the deep waters of the Aru Trough and Cenderawasih Bay that come close to the current coastline to maintain isolation of the Birds Head through periods of sea level minima (?135 m). The dates of divergence of turtle populations north and south of the ranges predate the telescopic uplift of the central ranges associated with oblique subduction of the Australian Plate beneath the Pacific Plate. Their isolation was probably associated with earlier uplift and drainage isolation driven by the accretion of island terranes to the northern boundary of the Australian craton that occurred earlier than the oblique subduction. The opening of Cenderawasih Bay is too recent (6 Mya) to have initiated the isolation of the Birds Head populations from those of the remainder of New Guinea, although its deep waters will have served to sustain the isolation through successive sea level changes. The molecular evidence suggests that the Birds Head docked with New Guinea some time before the Central Ranges emerged as a barrier to turtle dispersal. Overall, deep genetic structure of the species complex reflects events and processes that occurred during Miocene, whereas structure within each clade across the New Guinea landscape relates to Pliocene and Pleistocene times. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 192–208.     

Phylogeography of the Petunia integrifolia complex in southern Brazil

In this study, we evaluated the genetic diversity of the Petunia integrifolia species group using a phylogeographical approach, and attempted to understand better its diversification and taxonomy. Plants from five morphological groups were collected, covering a large part of the geographical distribution of most of the species. Two major clades were found in the phylogenetic tree, and an additional lineage, corresponding to P. inflata, was found in the haplotype network obtained for plastid markers. All three lineages are clearly delimited geographically, but, with the exception of P. inflata, the morphological groups were not genetically distinct. Our results suggest that a population expansion after a size reduction resulted in the establishment of two distinct and allopatric groups c. 0.5 Mya, one group occurring in a geologically ancient area, and the other occurring in areas that were under the influence of a series of marine transgressions/regressions at the end of the Pleistocene. These two clades are evolutionarily significant units with significantly different allele frequencies in their nuclear genome and reciprocal monophyly in maternal, uniparentally inherited markers. All our results suggest that the morphology‐based taxonomy in this group does not reflect its evolutionary history, and revision of its species limits should incorporate the distribution of the genetic diversity. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 174 , 199–213.