The phylogeographical history of the Iberian steppe plant Ferula loscosii (Apiaceae): a test of the abundant-centre hypothesis

Biogeography of Rhaponticoides, an Irano-Turanian element in the Mediterranean flora

Floristic relationships between the Irano-Turanian and Mediterranean regions have been known from old. However, only a few biogeographical analyses based on molecular data have evaluated the history of steppe plants within the Mediterranean basin. Our study aims to contribute to a better understanding of the migratory and diversification processes by reconstructing the biogeography of Rhaponticoides (Cardueae), distributed in the Mediterranean and Irano-Turanian regions. We generated nuclear and plastid sequences that were analyzed by Bayesian inference. We used the resulting phylogeny for dating the diversification of the genus and examining the dispersal pathways. Two clades were recovered, an Irano-Turanian clade and a Mediterranean clade. The origin of the genus was placed in the Anatolian plateau in the Middle Miocene. The genus experienced several diversifications and expansions correlated to the Messinian salinity crisis and the environmental changes in the Pliocene and the Quaternary. Rhaponticoides migrated following two routes reflecting the two souls of the genus: Irano-Turanian taxa colonized the steppes of Eurasia whilst Mediterranean species migrated via eastern and central Mediterranean and North Africa, leaving a trail of species; both pathways ended in the Iberian Peninsula. Our study also confirms that more work is needed to unravel phylogenetic relationships in Rhaponticoides.

Molecular Systematics of Valerianella Mill. (Caprifoliaceae): Challenging the Taxonomic Value of Genetically Controlled Carpological Traits

Valerianella (cornsalad) is a taxonomically complex genus formed by 50–65 annual Holarctic species classified into at least four main sections. Carpological traits (sizes and shapes of achenes and calyx teeth) have been used to characterize its sections and species. However, the potential systematic value of these traits at different taxonomic ranks (from sections to species (and infraspecific taxa)) has not been tested phylogenetically yet. Here, we have assessed the evolutionary systematic value of Valerianella diagnostic carpological traits at different hierarchical ranks and have demonstrated their ability to separate taxa at the sectional level but not at species level for species of several species pairs. A total of 426 individuals (17 species, 4 sections) of Valerianella were analyzed using AFLP and plastid data. Genetic clusters, phylogenetic trees, and haplotype networks support the taxonomic classification of Valerianella at the four studied sectional levels (V. sects. Valerianella, Cornigerae, Coronatae, Platycoelae) but show admixture for ten taxa from five species pairs (V. locusta—V. carinata, V. coronata—V. pumila, V. multidentata—V. discoidea, V. dentata—V. rimosa, V. eriocarpa—V. microcarpa), which are not reciprocally monophyletic. Dating analyses indicate that the Valerianella sections are relatively old (mid-Miocene), while most species diverged in the Pliocene–Pleistocene. A new section Valerianella sect. Stipitae is described to accommodate the highly divergent and taxonomically distinct V. fusiformis type species. Taxonomic treatments that recognize the sectional ranks and that subsume the separate species of each species pair into single species represent a natural classification for Valerianella.

Phylogeny and Historical Biogeography of Veronica Subgenus Pentasepalae (Plantaginaceae): Evidence for Its Origin and Subsequent Dispersal

Simple Summary

The Irano-Turanian phytogeographical region is considered a biodiversity reservoir for adjacent regions. The present phylogeographic study suggests that Veronica subgenus Pentasepalae originated in the Iranian plateau and was dispersed via a North African route to the Mediterranean and the Euro-Siberian regions. These findings highlight the importance of the Iranian plateau as a center of origin for many temperate plant species. Our results also resolve several taxonomic and phylogenetic issues surrounding the Southwest Asian species of this subgenus.

Abstract

Veronica subgenus Pentasepalae is the largest subgenus of Veronica in the Northern Hemisphere with approximately 80 species mainly from Southwest Asia. In order to reconstruct the phylogenetic relationships among the members of V. subgenus Pentasepalae and to test the “out of the Iranian plateau” hypothesis, we applied thorough taxonomic sampling, employing nuclear DNA (ITS) sequence data complimented with morphological studies and chromosome number counts. Several high or moderately supported clades are reconstructed, but the backbone of the phylogenetic tree is generally unresolved, and many Southwest Asian species are scattered along a large polytomy. It is proposed that rapid diversification of the Irano-Turanian species in allopatric glacial refugia and a relatively high rate of extinction during interglacial periods resulted in such phylogenetic topology. The highly variable Asian V. orientalis–V. multifida complex formed a highly polyphyletic assemblage, emphasizing the idea of cryptic speciation within this group. The phylogenetic results allow the re-assignment of two species into this subgenus. In addition, V. bombycina subsp. bolkardaghensis, V. macrostachya subsp. schizostegia and V. fuhsii var. linearis are raised to species rank and the new name V. parsana is proposed for the latter. Molecular dating and ancestral area reconstructions indicate a divergence age of about 9 million years ago and a place of origin on the Iranian Plateau. Migration to the Western Mediterranean region has likely taken place through a North African route during early quaternary glacial times. This study supports the assumption of the Irano-Turanian region as a source of taxa for neighboring regions, particularly in the alpine flora.

Phylogeography and ecological differentiation of strictly Mediterranean taxa: the case of the Iberian endemic Odontites recordonii

PREMISE

Ecological drivers for genetic differentiation in Mediterranean climates are still underexplored. We have used the strictly Mediterranean endemic Odontites recordonii as a model species to address this question. This species is one of the three Iberian representatives of the O. vernus group, which are morphologically similar. Thus, it was additionally necessary to clarify their phylogenetic relationships.

METHODS

We used amplified fragment length polymorphisms to reveal phylogenetic relationships within O. vernus group, and to reconstruct the phylogeographic patterns within O. recordonii. Additionally, ecological niche models were generated to detect refugia along the Quaternary climatic oscillations. And finally, alleles under natural selection were identified, and correlations between allele presences and environmental variables were calculated in order to shed light on the ecological drivers promoting differentiation.

RESULTS

The three species from the O. vernus group were recovered as distinct species. Three genetic groups were found within O. recordonii and a putative refugium was detected for each one. Eighty-one alleles could be under diversifying selection, and 58 alleles showed significant correlations with environmental variables, especially with temperature and precipitation seasonality and summer drought.

CONCLUSIONS

The three Iberian species of the O. vernus group are reciprocal monophyletic taxa. The three genetic groups of O. recordonii could have been restricted to narrow refugia during the Quaternary and displayed present distributions in accordance with bioclimatic conditions. Temperature and precipitation seasonality and the intensity of summer drought are definitory climatic parameters of Mediterranean-type climates, and they could have acted as drivers of genetic differentiation on O. recordonii.

Phylogeography of a gypsum endemic plant across its entire distribution range in the western Mediterranean

PREMISE

Gypsum soils in the Mediterranean Basin house large numbers of edaphic specialists that are adapted to stressful environments. The evolutionary history and standing genetic variation of these taxa have been influenced by the geological and paleoclimatic complexity of this area and the long-standing effect of human activities. However, little is known about the origin of Mediterranean gypsophiles and the factors affecting their genetic diversity and population structure.

METHODS

Using phylogenetic and phylogeographic approaches based on microsatellites and sequence data from nuclear and chloroplast regions, we evaluated the divergence time, genetic diversity, and population structure of 27 different populations of the widespread Iberian gypsophile Lepidium subulatum throughout its entire geographic range.

RESULTS

Lepidium subulatum diverged from its nearest relatives ~3 million years ago, and ITS and psbA/matK trees supported the monophyly of the species. These results suggest that both geological and climatic changes in the region around the Plio-Pleistocene promoted its origin, compared to other evolutionary processes. We found high genetic diversity in both nuclear and chloroplast markers, but a greater population structure in the chloroplast data. These results suggest that while seed dispersal is limited, pollen flow may be favored by the presence of numerous habitat patches that enhance the movement of pollinators.

CONCLUSIONS

Despite being an edaphic endemic, L. subulatum possesses high genetic diversity probably related to its relatively old age and high population sizes across its range. Our study highlights the value of using different markers to fully understand the phylogeographic history of plant species.

A story from the Miocene: Clock-dated phylogeny of Sisymbrium L. (Sisymbrieae, Brassicaceae)

Morphological variability and imprecise generic boundaries have hindered systematic, taxonomical, and nomenclatural studies of Sisymbrium L. (Brassicaceae, Sisymbrieae DC.). The members of this almost exclusively Old-World genus grow mostly on highly porous substrates across open steppe, semidesert, or ruderal habitats in the temperate zone of the Northern Hemisphere and African subtropics. The present study placed the biological history of Sisymbrium L. into time and space and rendered the tribus Sisymbrieae as monotypic. Five nuclear-encoded and three chloroplast-encoded loci of approximately 85% of all currently accepted species were investigated. Several accessions per species covering their whole distribution range allowed for a more representative assessment of intraspecific genetic diversity. In the light of fossil absence, the impact of different secondary calibration methods and taxon sets on time spans was tested, and we showed that such a combinatorial nested dating approach is beneficial. Multigene phylogeny accompanied with a time divergence estimation analysis placed the onset and development of this tribus into the western Irano-Turanian floristic region during the Miocene. Continuous increase in continentality and decrease in temperatures promoted the diversity of the Sisymbrieae, which invaded the open grasslands habitats in Eurasia, Mediterranean, and South Africa throughout the Pliocene and Pleistocene. Our results support the assumption of the Irano-Turanian region as a biodiversity reservoir for adjacent regions.

Geography alone cannot explain Tetranychus truncatus (Acari: Tetranychidae) population abundance and genetic diversity in the context of the center-periphery hypothesis

The center-periphery hypothesis (CPH) states that the genetic diversity, genetic flow, and population abundance of a species are highest at the center of the species' geographic distribution. However, most CPH studies have focused on the geographic distance and have ignored ecological and historical effects. Studies using niche models to define the center and periphery of a distribution and the interactions among geographical, ecological, and historical gradients have rarely been done in the framework of the CPH, especially in biogeographical studies of animal species. Here, we examined the CPH for a widely distributed arthropod, Tetranychus truncatus (Acari: Tetranychidae), in eastern China using three measurements: geographic distance to the center of the distribution (geography), ecological suitability based on current climate data (ecology), and historical climate data from the last glacial maximum (history). We found that the relative abundances of different populations were more strongly related to ecology than to geography and history. Genetic diversity within populations and genetic differentiation among populations based on mitochondrial marker were only significantly related to history. However, the genetic diversity and population differentiation based on microsatellites were significantly related to all three CPH measurements. Overall, population abundance and genetic pattern cannot be explained very well by geography alone. Our results show that ecological gradients explain the variation in population abundance better than geographic gradients and historical factors, and that current and historical factors strongly influence the spatial patterns of genetic variation. This study highlights the importance of examining more than just geography when assessing the CPH.

Long-term survival and successful conservation? Low genetic diversity but no evidence for reduced reproductive success at the north-westernmost range edge of Poa badensis (Poaceae) in Central Europe

Many steppe species reach their (north)westernmost distribution limit in western Central Europe. This also applies to Poa badensis, a rare steppe plant of calcareous rock/sand vegetation. To explore potential differences in reproductive success and genetic composition of peripheral populations, we analysed the absolute (north)westernmost occurrences in Western Germany and populations at the western margin (Eastern Austria) and the centre (Central Hungary) of the Pannonicum, representing a part of the continuous range. Specifically, we discuss the genetic and reproductive constitution of the (north)westernmost exclave and draw conclusions on the species' biogeographical and conservation history in this region. Therefore, we used two independent molecular marker systems (AFLPs, cpDNA sequences) and a set of performance parameters. Overall, lowest regional genetic diversity was found in Western Germany, which is mainly a result of the specific history of two populations. However, this low genetic diversity was not accompanied by reduced reproductive success. The Eastern Austrian populations showed reduced genetic diversity and predominantly reduced performance, interpreted as a consequence of small population sizes. Central Hungarian populations showed the overall highest genetic diversity and comparatively high performance values. We observed high admixture and haplotype sharing between Austrian and Hungarian populations, indicating gene flow among these regions. In contrast, we interpreted the increased population differentiation within, and the clear distinctiveness of the German exclave as a long-term isolation of these (north)westernmost occurrences. Our results, overall, prove the good constitution of these populations and, together with their particular biogeographical history, highlight their conservation value.

No obvious genetic erosion, but evident relict status at the westernmost range edge of the Pontic-Pannonian steppe plant Linum flavum L. (Linaceae) in Central Europe

We investigate patterns of genetic variation along an east–west transect of Central European populations of Linum flavum and interpret the Quaternary history of its peripheral populations, especially those at the westernmost isolated range edge, discussing their migrations and possible relict status. We defined our peripheral transect across three study regions from Central Hungary, eastern Austria to southwestern Germany. Using AFLP fingerprinting and cpDNA sequence variation (rpL16 intron, atpI‐H), we analyzed 267 and 95 individuals, respectively, representing each study region by four populations. Hierarchical AMOVA (AFLPs) indicated significant variation among study regions (12% of total variance) and moderate differentiation between populations (10%). Population differentiation was high at the westernmost range edge (11.5%, Germany), but also in the east (13.4%, Hungary), compared to the Austrian study region (8.6%). Correspondingly, AFLP diversity was highest in the center of the study transect in eastern Austria. CpDNA haplotypes support a pattern of regional structuring with the strongest separation of the westernmost range edge, and some haplotype sharing among Austrian and Hungarian individuals. Equilibrating nucleotide versus haplotype diversity patterns, the highly diverse populations at the Pannonian range edge (Austria) indicate long‐term persistence, while Central Pannonian populations are obviously effected by recent bottlenecks. Intermediate nucleotide, but high haplotype diversity within the westernmost exclave (Swabian Alb), is indicative of a founder bottleneck during its pre‐LGM or early postglacial migration history, followed by sufficient time to accumulate cpDNA variation. The not obviously reduced genetic diversity and distinctiveness of L. flavum at the westernmost range edge suggest a long‐term persistence (relict status) of populations in this region, where the species has survived probably even the Würm glaciation in extra‐Mediterranean refugia. This genetic relict variation represents an important part of the overall genetic diversity found in the western periphery of this steppe plant and highlights the high conservation priority of respective gene pools.

Phylogeography of a tough rock survivor in European dry grasslands

Phylogeographic analyses of plants in Europe have revealed common glacial refugia and migration routes for several trees and herbs with arctic-alpine distributions. The postglacial histories of dry grassland species in central Europe have rarely been analyzed, even though the extremely species-rich habitat is threatened. Sedum album (Crassulaceae) is a common inhabitant of rocky sites in central European dry grasslands. We inferred the phylogeographic history of S. album over its distribution range in Europe. Genetic diversity within and differentiation between 34 S. album populations was examined using AFLP markers. Population isolation was indicated based on the rarity of the fragments and by isolation-by-distance effects. We sequenced the trnL-trnF region in 32 populations and used chloroplast microsatellites to analyze chloroplast haplotype distributions. Two distinct S. album lineages were detected. One lineage was comprised of populations from eastern and central parts of central Europe, and the Apennine Peninsula. A second lineage was comprised of populations from the Iberian Peninsula and western and northern parts of central Europe. Glacial refugia were identified based on the accumulation of ancient chloroplast haplotypes, high diversity of AFLP fragments within populations, and high levels of rare fragments in Liguria, Serbia, the Apennine and Iberian peninsulas. Cryptic refugia were detected in the Czech Republic and Slovakia. Isolation by distance was present all over the distribution range, and it was separately detected in southwestern and central Europe. In western Europe, where a contact zone between the two lineages can be expected, no isolation by distance was detected. Our results suggest migration routes of S. album northeastward from glacial refugia in southern Iberia, northward from the Apennine Peninsula, and northward and westward from the southeastern parts of central Europe. Therefore, central European grasslands were recently colonized by northern cryptic populations and source populations originating in the east and the Apennine Peninsula.

Environmental isolation explains Iberian genetic diversity in the highly homozygous model grass Brachypodium distachyon

BACKGROUND

Brachypodium distachyon (Poaceae), an annual Mediterranean Aluminum (Al)-sensitive grass, is currently being used as a model species to provide new information on cereals and biofuel crops. The plant has a short life cycle and one of the smallest genomes in the grasses being well suited to experimental manipulation. Its genome has been fully sequenced and several genomic resources are being developed to elucidate key traits and gene functions. A reliable germplasm collection that reflects the natural diversity of this species is therefore needed for all these genomic resources. However, despite being a model plant, we still know very little about its genetic diversity. As a first step to overcome this gap, we used nuclear Simple Sequence Repeats (nSSR) to study the patterns of genetic diversity and population structure of B. distachyon in 14 populations sampled across the Iberian Peninsula (Spain), one of its best known areas.

RESULTS

We found very low levels of genetic diversity, allelic number and heterozygosity in B. distachyon, congruent with a highly selfing system. Our results indicate the existence of at least three genetic clusters providing additional evidence for the existence of a significant genetic structure in the Iberian Peninsula and supporting this geographical area as an important genetic reservoir. Several hotspots of genetic diversity were detected and populations growing on basic soils were significantly more diverse than those growing in acidic soils. A partial Mantel test confirmed a statistically significant Isolation-By-Distance (IBD) among all studied populations, as well as a statistically significant Isolation-By-Environment (IBE) revealing the presence of environmental-driven isolation as one explanation for the genetic patterns found in the Iberian Peninsula.

CONCLUSIONS

The finding of higher genetic diversity in eastern Iberian populations occurring in basic soils suggests that these populations can be better adapted than those occurring in western areas of the Iberian Peninsula where the soils are more acidic and accumulate toxic Al ions. This suggests that the western Iberian acidic soils might prevent the establishment of Al-sensitive B. distachyon populations, potentially causing the existence of more genetically depauperated individuals.

Pre-Holocene Origin for the Coronopus navasii Disjunction: Conservation Implications from Its Long Isolation

Integration of unexpected discoveries about charismatic species can disrupt their well-established recovery plans, particularly when this requires coordinate actions among the different governments responsible. The Critically Endangered Coronopus navasii (Brassicaceae) was considered a restricted endemism to a few Mediterranean temporary ponds in a high mountain range of Southeast Spain, until a new group of populations were discovered 500 km North in 2006. Ten years after this finding, its management has not been accommodated due to limited information of the new populations and administrative inertia. In this study, DNA sequences and species distribution models are used to analyse the origin of the C. navasii disjunction as a preliminary step to reassess its recovery plan. Molecular results placed the disjunction during Miocene-Pleistocene (6.30-0.49 Mya, plastid DNA; 1.45-0.03 Mya, ribosomal DNA), which discards a putative human-mediated origin. In fact, the haplotype network and the low gene flow estimated between disjunct areas suggest long-term isolation. Dispersal is the most likely explanation for the disjunction as interpreted from the highly fragmented distribution projected to the past. Particularly, a northward dispersal from Southeast is proposed since C. navasii haplotype network is connected to the sister-group through the southern haplotype. Although the reassessment of C. navasii conservation status is more optimistic under the new extent of occurrence, its long-term survival may be compromised due to the: (1) natural fragmentation and rarity of the species habitat, (2) genetic isolation between the two disjunct areas, and (3) northward shift of suitable areas under future climate change scenarios. Several ex-situ and in-situ conservation measures are proposed for integrating Central East Spanish populations into the on-going recovery plan, which still only contemplates Southeast populations and therefore does not preserve the genetic structure and diversity of the species.

Geographical isolation caused the diversification of the Mediterranean thorny cushion-like Astragalus L. sect. Tragacantha DC. (Fabaceae)

Understanding the origin and evolution of Mediterranean vascular flora within the long-term context of climate change requires a continuous study of historical biogeography supported by molecular phylogenetic approaches. Here we provide new insights into the fascinating but often overlooked diversification of Mediterranean xerophytic plants. Growing in some of the most stressing Mediterranean environments, i.e. coastal and mountainous opened habitats, the circum-Mediterranean Astragalus L. sect. Tragacantha DC. (Fabaceae) gathers several thorny cushion-like taxa. These have been the subjects of recent taxonomical studies, but they have not yet been investigated within a comprehensive molecular framework. Bayesian phylogenetics applied to rDNA ITS sequences reveal that the diversification of A. sect. Tragacantha has roots dating back to the Pliocene, and the same data also indicate an eastern-western split giving rise to the five main lineages that exist today. In addition, AFLP fingerprinting supports an old east-west pattern of vicariance that completely rules out the possibility of a recent eastern origin for western taxa. The observed network of genetic relationships implies that contrary to what is widely claimed in the taxonomic literature, it is range fragmentation, as opposed to a coastal-to-mountain ecological shift, that is likely the main driver of diversification.

The genetic diversity and spatial genetic structure of the Corso-Sardinian endemic Ferula arrigonii Bocchieri (Apiaceae)

Corsica and Sardinia represent major hotspots of plant diversity in the Mediterranean area and are priority regions for conservation due to their high number of endemic plant species. However, information supporting human decision-making on the conservation of these species is still scarce, especially at the genetic level. In this work, the first assessment is reported of the species-wide spatial genetic structure and diversity of Ferula arrigonii Bocchieri, a Corso-Sardinian endemic located in a few coastal sites and on small islands. Nine populations covering the entire natural range of the species were investigated by means of AFLP (amplified fragment length polymorphism) markers. Results indicate that this species is characterised by high levels of genetic polymorphism (92% polymorphic fragments) and of genetic diversity (H(w) = 0.317) and by relatively low differentiation among populations (F(st) = 0.057). PCoA, Bayesian analysis and neighbour-joining clustering were also employed to investigate the genetic structure of this species. Three genetically distinct groups were detected, although with considerable overlap between populations.

Balearic insular isolation and large continental spread framed the phylogeography of the western Mediterranean Cheirolophus intybaceus s.l. (Asteraceae)

Recent Quaternary geological and climate events have shaped the evolutionary histories of plant species in the Mediterranean basin, one of the most important hotspots of biodiversity. Genetic analyses of the western Mediterranean Cheirolophus intybaceus s.l. (Asteraceae) based on AFLP were conducted to establish the relationships between its close species and populations, to reconstruct the phylogeography of the group and to analyse potential unidirectional versus bidirectional dispersals between the Ibero-Provençal belt and the Balearic Islands. AFLP data revealed two main genetic groups, one constituted by the Balearic populations and Garraf (NE Iberia) and the other formed by the remaining mainland populations that were further sub-structured into two geographically separated subgroups (SE + E Iberia and NE Iberia + SW France). Genetic diversity and spatial structure analyses suggested a mid-Pleistocene scenario for the origin of C. intybaceus in southern Iberia, followed by dispersal to the north and a single colonisation event of the Balearic archipelago from the near Dianic NE Iberian area. This hypothesis was supported by paleogeographic data, which showed the existence of terrestrial connections between the continent and the islands during the Middle-Late Pleistocene marine regressions, whereas the more recent single back-colonisation of the mainland from Mallorca might be explained by several hypotheses, such as long-distance dispersal mediated by migratory marine birds or sea currents.

Range expansion of a selfing polyploid plant despite widespread genetic uniformity

BACKGROUND AND AIMS

Ongoing and previous range expansions have a strong influence on population genetic structure of plants. In turn, genetic variation in the new range may affect the population dynamics and the expansion process. The annual Ceratocapnos claviculata (Papaveraceae) has expanded its Atlantic European range in recent decades towards the north and east. Patterns of genetic diversity were investigated across the native range to assess current population structure and phylogeographical patterns. A test was then made as to whether genetic diversity is reduced in the neophytic range and an attempt was made to identify source regions of the expansion.

METHODS

Samples were taken from 55 populations in the native and 34 populations in the neophytic range (Sweden, north-east Germany). Using amplified fragment length polymorphism markers an analysis was made of genetic variation and population structure (Bayesian statistical modelling) and population differentiation was quantified. Pollen/ovule ratio was analysed as a proxy for the breeding system.

KEY RESULTS

Genetic diversity at population level was very low (mean H(e) = 0·004) and two multilocus genotypes dominated large parts of the new range. Population differentiation was strong (F(ST) = 0·812). These results and a low pollen/ovule ratio are consistent with an autogamous breeding system. Genetic variation decreased from the native to the neophytic range. Within the native range, H(e) decreased towards the north-east, whereas population size increased. According to the Bayesian cluster analysis, the putative source regions of the neophytic range are situated in north-west Germany and adjacent regions.

CONCLUSIONS

Ceratocapnos claviculata shows a cline of genetic variation due to postglacial recolonization from putative Pleistocene refugia in south-west Europe. Nevertheless, the species has expanded successfully during the past 40 years to southern Sweden and north-east Germany where it occurs as an opportunistic neophyte. Recent expansion was mainly human-mediated by single long-distance diaspore transport and was facilitated by habitat modification.

Molecular phylogeny and systematics of the highly polymorphic Rumex bucephalophorus complex (Polygonaceae)

Rumex bucephalophorus is a very polymorphic species that has been subjected to various taxonomic studies in which diverse infraspecific taxa have been recognised on the basis of diaspore traits. In this study we used molecular markers (ITS and AFLP) to explore this remarkable diversity, to test previous hypotheses of classification, and attempt to explain biogeographic patterns. Results show that R. bucephalophorus forms a monophyletic group in which diversification began around 4.2 Mya, at the end of Messinian Salinity Crisis. The two molecular markers clearly show a deep divergence separating subsp. bucephalophorus from all other subspecific taxa, among which subsp. canariensis also constitutes a separate and well distinguishable unit. In contrast, subspecies hispanicus and subsp. gallicus constitute a monophyletic group in which three subgroups can be recognised: subsp. hispanicus, subsp. gallicus var. gallicus and subsp. gallicus var. subaegeus. However, these three subgroups are not clearly distinguished genetically or morphologically, so that in formal classification it would be preferable to treat them at the varietal level.