Morphological differentiation supports the genetic pattern of the geographic structure of Juniperus thurifera (Cupressaceae)

Juniperus thurifera is an important component of woodland communities of dry sites within the West Mediterranean region and is characterised by a strongly disjunctive geographic range. Two subspecies were recognised, subsp. thurifera in Europe and subsp. africana in Africa. The aim of the study was the comparison of phenetic diversity to the pattern of AFLP geographic differentiation of the species described in the literature. The examination of phenetic diversity was based on the biometrical analysis of 17 populations using 12 morphological characters of cone and seed. The differences among populations were analysed using Student’s t test, analysis of discrimination, UPGMA agglomeration and hierarchical analysis of variance. The majority of morphological characters differentiated at a statistically significant level between populations and between J. thurifera subsp. thurifera and subsp. africana. Three groups of populations were detected using multivariate statistical analyses. The first, well separated, is subsp. africana, while the following two concern subsp. thurifera. The morphological differentiation of populations appeared similar to that described on the AFLP. The Gibraltar Straight appeared to be the most important barrier.     

Phylogenetic analysis and differentiation of Veronica subgenus Stenocarpon in the Balkan Peninsula

The Balkan Peninsula is considered the most important refugium for species during the Pleistocene glaciations and today harbours c. 2000 endemic species, but we know surprisingly little about the evolution of taxa in this region. Veronica saturejoides, V. thessalica and V. erinoides are a group of closely related alpine taxa endemic to the Balkan Peninsula. Here, we analyse four DNA regions [the nuclear chalcone synthase intron (CHSi) and ribosomal internal transcribed spacer (ITS) region and the plastid rpoB‐trnC spacer and trnL‐trnL‐trnF region] and amplified fragment length polymorphism (AFLP) fingerprints to provide a phylogenetic hypothesis for the relationships among these taxa. Additionally, we analyse leaf morphological characters used to distinguish the three subspecies of V. saturejoides. The analyses support the distinction of the three subspecies based on previously intuitively suggested characters. Nuclear chalcone synthase intron data indicate that the southern taxa are genetically much more diverse than the more northern V. saturejoides subsp. saturejoides. Phylogenetic relationships inferred from this region and AFLP fingerprints support the monophyly of V. saturejoides. In contrast, plastid DNA regions suggest a closer relationship of V. saturejoides subsp. saturejoides to V. thessalica. The most likely scenario involves introgression into V. saturejoides subsp. saturejoides from V. thessalica. © 2009 The Linnean Society of London, Botanical Journal of the Linnean Society, 2009, 159 , 616–636.     

Genetic and morphological variation in Viola suavis s.l. (Violaceae) in the western Balkan Peninsula: two endemic subspecies revealed

The Balkan Peninsula, with many endemic species, is known as one of the most important speciation and diversification centres in Europe. Here, we present a study of the western Balkan populations of the polymorphic European species, V. suavis s.l., which have been reported under the name V. adriatica, but their taxonomic status and position within the genus have remained uncertain. Viola suavis s.l. and nine close relatives sampled across Europe were subjected to molecular (sequencing of nuclear ribosomal internal transcribed spacers and amplified fragment length polymorphism), karyological and morphometric analyses. Our results revealed the presence of four allopatric, genetically and morphologically differentiated lineages within V. suavis s.l. in Europe, which are suggested here to be recognized at the subspecific rank. Populations from the western Balkans were segregated into two distinct entities: (1) those from north-western Croatia correspond to the previously recognized taxon, V. suavis subsp. adriatica and (2) those from southern Dalmatia (southern Croatia, southern Bosnia and Herzegovina, and south-western Montenegro) are described here as V. suavis subsp. austrodalmatica subsp. nov. The other two lineages of V. suavis s.l., which both harbour blue- and white-flowered morphotypes, occur in central and eastern Europe (V. suavis subsp. suavis) and in north-eastern Spain (plants provisionally treated as V. suavis ‘Spain’). The AFLP and morphological data indicate gene flow between the nominate subspecies and V. suavis subsp. adriatica in a few localities. The distribution of the two western Balkan subspecies is discussed and an identification key to the V. suavis subspecies in Europe is presented.     

Ecotypes and genetic structure of Rhinanthus alectorolophus (Orobanchaceae) in southwestern Germany

European annual species of the genus Rhinanthus often exhibit seasonal ecotypic variation, a phenomenon also known from related genera of hemiparasitic Orobanchaceae. Populations with different flowering times exist, correlated with differences in a number of morphological characters. The present study evaluates the correlation of morphological characters and genetic differentiation of populations of Rhinanthus alectorolophus. Thirty-nine populations of three different subspecies from southwestern Germany were sampled. A total of 798 individuals were used for morphological analyses and 187 of these for AFLP analyses. Principal component analysis showed that morphological variation is mostly continuous. In a discriminant analysis based on morphological characters, only 89.7 % of all individuals were correctly assigned to their previously determined subspecies, indicating that subspecies identification is ambiguous for some populations. Using AFLP data and Bayesian assignment analysis, the sampled individuals could be grouped in three genetic clusters which do not correspond to the three subspecies. Instead, the clustering shows a clear geographic pattern and a Mantel test likewise revealed a significant correlation between genetic and geographic distances. Correlations of genetic distances with differences in morphological characters were weak and mostly insignificant. The results indicate that the subspecies of R. alectorolophus do not form discrete entities and that the character combinations distinguishing them are homoplastic.     

Electrophoretic study of isozyme patterns in some wild populations of Aubrieta columnae Guss. (Cruciferae)

Aubrieta columnae Guss. is currently divided into three subspecies: A. columnae subsp. columnae, A. columnae subsp. italica, both found as isolated and fragmented populations in rocky habitats of Central and Southern Apennines (Italy), and A. columnae subsp. croatica found in the Balkan region. In order to gain information about the degree of genetic variability and to clarify taxonomic relationships among these taxa, we studied the isozyme patterns at 8 marker loci of 376 individuals from 8 populations by means of starch gel electrophoresis. Data analysis by using Wright's F-statistics and UPGMA clustering method was performed. The results show: 1) a general deviation from Hardy–Weinberg equilibrium within each subspecies; 2) a lesser genetic variation in populations occurring in habitats characterized by milder climatic conditions and relatively small seasonal variations; 3) a relatively high degree of differentiation between the three subspecies; 4) the possible common transadriatic origin of A. columnae subsp. italica and A. columnae subsp. croatica; 5) the possible origin of A. columnae subsp. columnae from A. columnae subsp. italica; and 6) that the current taxonomic status of A. columnae may be substantially confirmed, even if the findings are from a limited number of loci explored.     

Extensive gene flow blurs phylogeographic but not phylogenetic signal in Olea europaea L.

Genetic structure and evolutionary patterns of the wild olive tree (Olea europaea L.) were investigated with AFLP fingerprinting data at three geographic levels: (a) phylogenetic relationships of the six currently recognized subspecies in Eurasia and Africa; (b) lineage identification in subsp. europaea of the Mediterranean basin; and (c) phylogeography in the western Mediterranean. Two statistical approaches (Bayesian inference and analysis of molecular variance) were used to analyse the AFLP fingerprints. To determine the congruency and transferability of results across studies previous RAPD and ISSR data were analysed in a similar manner. Comparisons proved that qualitative results were mostly congruent but quantitative values differed, depending on the method of analysis. Neighbour-Joining analysis of AFLP phenotypes supported current classification of subspecies. At a Mediterranean scale no clear cut phylogeographic pattern was recovered, likely due to extensive gene flow between populations of subsp. europaea. Gene flow estimates calculated with conventional F-statistics showed that reproductive barriers separated neither populations nor lineages of O. europaea. Genetic divergence between eastern and western parts of the Mediterranean basin was observed only when geographical and population information were incorporated into the analyses through hierarchical analysis of molecular variance (AMOVA). Within the western Mediterranean, the highest genetic diversity was found in two regions: on both sides of the Strait of Gibraltar and in the Balearic archipelago. Additionally, long-lasting isolation of the northern-most populations of the Iberian Peninsula appeared to be responsible for a significant divergence.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Geographically Related Variation in Epicuticular Wax Traits of Pinus nigra Populations from Southern Carpathians and Central Balkans – Taxonomic Considerations

The chemical composition of epicuticular waxes of nine populations from three Pinus nigra J. F. Arnold subspecies (namely subsp. nigra, subsp. banatica (Borbás ) Novák , and subsp. pallasiana (Lamb .) Holmboe ) from Southern Carpathians and central Balkan Peninsula were analyzed using GC/MS and GC/FID chromatography, and multivariate statistical techniques with respect to biogeography and taxonomy. In the needle waxes, four primary alcohols and 14 n‐alkanes ranging from C₂₁ to C₃₃ were identified, and the most abundant compounds were the four odd‐numbered n‐alkanes C₂₇, C₂₅, C₂₃, and C_29. Multivariate statistical analyses (CDA and CA) have shown existence of three P. nigra groups and suggested clinal differentiation as a mechanism of genetic variation across a geographic area: the first group consisted of the southernmost populations of subsp. pallasiana from Macedonia, the second consisted of the northernmost subsp. banatica populations from Romania, while all populations in Serbia described as three different subspecies (nigra, banatica, and pallasiana) formed the third group together with subsp. nigra population from Bosnia and Herzegovina. According to simple linear regression, geographic latitude and four bioclimatic parameters were moderately correlated with the contents of epicuticular wax compounds that are important in population discrimination, while stepwise multiple regression showed that latitude participated in most of the regression models for predicting the composition of the epicuticular waxes. These results agree with CDA and CA analysis, and confirmed the possibility of recognition of fine geographic differentiation of the analyzed P. nigra populations.     

Composition of Leaf n‐Alkanes in Three Satureja montana L. Subspecies from the Balkan Peninsula: Ecological and Taxonomic Aspects

The composition of the epicuticular leaf n‐alkanes of eight populations of three Satureja montana subspecies (S. montana L. subsp. pisidica (Wettst.) ?ili? , S. montana L. subsp. montana, and S. montana L. subsp. variegata (Host ) P. W. Ball ), from central and western areas of the Balkan Peninsula was characterized by GC‐FID and GC/MS analyses. In the leaf waxes, 15 n‐alkane homologs with chain‐lengths ranging from C₂₁ to C₃₅ were identified. The main n‐alkane in almost all samples was n‐nonacosane (C₂₉), but differences in the contents of three other dominant n‐alkanes allowed separating the coastal from the continental populations. The diversity and variability of the epicuticular‐leaf‐n‐alkane patterns and their relation to different geographic and bioclimatic parameters were analyzed by several statistical methods (principal component, discriminant, and cluster analyses as well as the Mantel test). All tests showed a high correlation between the leaf n‐alkane pattern and the geographical distribution of the investigated populations, confirming the differentiation between S. montana subsp. pisidica and the other two subspecies. The S. montana subsp. variegata and S. montana subsp. montana populations are geographically closer and their differentiation according to the leaf‐n‐alkane patterns was not clear, even though there was some indication of discrimination between them. Moreover, most of the bioclimatic parameters related to temperature were highly correlated with the differentiation of the coastal and the continental populations.     

Hybridizations and genetic relationships among Lindernia species (Scrophulariaceae): L. procumbens and two subspecies of L. dubia

Lindernia procumbens and L. dubia are common annual weeds in flooded rice fields of Japan. Two subspecies of L. dubia, subsp. major and subsp. dubia, are usually recognized in Japan but they are both regarded as synonyms of L. dubia elsewhere. In a cluster analysis based on AFLP, most L. dubia subsp. major formed a separate cluster from L. dubia subsp. dubia although 11% of haplotypes classified using AFLP were not coincident with classification using the shape of leaf bases, which is the commonly used identification trait. Artificial F₁ plants between L. procumbens and L. dubia subsp. major, and between L. procumbens and L. dubia subsp. dubia did not produce seed. Forty percent of capsules produced by F₁ plants from these two subspecies were slimmer and 80% pollen were sterile in slimmer capsules. However, seed number of most F₁ capsules was not different from that of self-fertilized plants, suggesting that there was no complete reproductive isolation between the subspecies. Natural hybridization of these subspecies may have occurred but we are not aware of it because F₁ plants are rare and F₂ plants are indistinguishable from these subspecies.     

The Balkan chamois,an archipelago or a peninsula? Insights from nuclear and mitochondrial DNA

The Balkan chamois (Rupicapra rupicapra balcanica) is widespread on the Balkan Peninsula, along mountain massifs from Croatia in the north to Greece in the south and Bulgaria in the east. Knowledge on the genetic structure of Balkan chamois populations is limited and restricted to local studies. Therefore, the main objective of this study was to use nuclear (16 microsatellites) and mitochondrial (partial 376 base pairs control region) markers to investigate the genetic structure of this chamois subspecies throughout its distribution range and to obtain information on the degree of connectivity of the different (sub)populations. We extracted DNA from bone, dried skin and muscle tissue and successfully genotyped 92 individuals of Balkan chamois and sequenced the partial control region in 44 individuals. The Bayesian analysis suggested 3 genetic clusters and assigned individuals from Serbia and Bulgaria to two separate clusters, while individuals from the other countries belonged to the same cluster. Thirty new haplotypes were obtained from partial mitochondrial DNA sequences, with private haplotypes in all analyzed populations and only two haplotypes shared among populations, indicating the possibility of past translocations. The subspecies genetic composition presented here provides the necessary starting point to assess the conservation status of the Balkan chamois and allows the development of conservation strategies necessary for its sustainable management and conservation.

     

Assessment of Daucus carota L. (Apiaceae) subspecies by chemotaxonomic and DNA content analyses

For a clearer distinction between the four subspecies of Daucus carota native from Portugal (subsp. carota, subsp. maximus, subsp. gummifer and subsp. halophilus), morphological features of the fruits, DNA content analyses by flow cytometry, and chemical characterization of the essential oils were undertaken.We found chemotaxonomic evidences to consider D. carota subsp. maximus as a separate species rather than a subspecies of D. carota. This separation is based on the morphometric analysis of the fruits and in the high levels of asarone present only in the essential oil of the subsp. maximus. The remaining subspecies are difficult to distinguish from each other based on the morphology of the fruits and in DNA content. However, based on the essential oils, it was possible to distinguish the subspecies halophilus from the other two (subsp. gummifer and subsp. carota) because of its high content of elemicin, with the other two having high levels of geranyl acetate.Based on these results, the subspecies maximus is proposed as a different species (Daucus maximus Desf.) and the taxonomic status of other three subspecies (subsp. carota, subsp. gummifer and subsp. halophilus) is maintained. Still, the latter three taxa need to be further studied for a more precise taxonomic characterization.     

Infraspecific classification reflects genetic differentiation in the widespread Petunia axillaris complex: A comparison among morphological,ecological, and genetic patterns of geographic variation

Plants are subjected to natural selection during expansion of their geographic range or when facing changes in environmental conditions, which in turn may affect phenotypic diversity. Studies on geographic phenotypic variation provide insight into the evolutionary processes and have long contributed to better understanding of diversification within and among species. Petunia axillaris is widely distributed in temperate South America, occurring throughout the Pampas region. The current taxonomy of this species recognizes three allopatric subspecies, which occupy nearly adjacent territories, according to floral morphology: P. axillaris subsp. axillaris, P. axillaris subsp. parodii, and P. axillaris subsp. subandina. In this study, we sampled the three P. axillaris subspecies and used both molecular markers (from the plastid and nuclear genomes) and morphological measurements to investigate how genetic diversity and morphological variation correlate to ecological variables and the geographic context. We used different forms of Mantel tests (partial and correlograms) to investigate the geographic distribution patterns in distinct types of similarity/dissimilarity among the populations and their relationships. We also modeled the morphological variation in P. axillaris as a function of the genetic marker frequencies in the populations. We found that the morphological differences leading to the recognition of different subspecies of P. axillaris reflect historical processes of isolation and that adaptation to different ecological conditions faced by each lineage is perhaps not merely a consequence of phenotypic plasticity. These findings suggest that differences between subspecies could represent an incipient speciation stage.     

Differentiation patterns in the E. MediterraneanScutellaria rubicunda group (Lamiaceae)

Strict sterility barriers where found between theScutellaria populations on Sicily and the Greek populations. The material from Samos is separated by a strict sterility barrier from the other Greek populations. In contrast, the populations on mainland Greece and adjacent islands are all ± interfertile. Three species are recognized,S. rubicunda Hornem., endemic to Sicily, andS. brevibracteata subsp.icarica, endemic to Samos and Ikaria. All other populations are referred toS. rupestris with eight subspecies, most of which are endemic to one island or to one mountain. Three new subspecies are described, viz. subsp.rechingeri and subsp.olympica, endemic to mt Vourinos and mt Olympus in North Greece, respectively, and subsp.caroli-henrici, native to the Malea peninsula of Peloponnisos. The phytogeographical connections and genetic differentiation within and between populations are discussed.Dedicated to Prof. DrK. H. Rechinger on the occasion of his 80th birthday.     

Past and present distribution of the genus Capra in Greece

The Alpine ibex Capra ibex and the wild goat Capra aegagrus are the only caprids known to have existed in the Balkan Peninsula. The agrimi is a feral primitive form of the domestic goat Capra hircus but may constitute an island subspecies of C. aegagrus endemic to Crete, Greece. Here, the distribution of C. ibex and C. aegagrus in this southernmost Balkan country was reviewed for the period from Middle Pleistocene until today through a scientific and historical document search. Such knowledge should provide directions to future conservation plans for these species and aid in tracing the origins of agrimi. Fossil locations indicate that both species lived in this country at elevations of similar altitude (≥10–?≤?2,497 m a.s.l.). C. ibex colonised Greece naturally, perhaps, not earlier than 126,000 years ago. It was spread over most of the mainland (36–41° N and 20–24° E) where from it disappeared after 5600 BC due to hunting. C. aegagrus constituted an exotic in Europe. Man introduced it to insular and mainland Greece (between 35–40° N and 19–26° E) from around 8600 and 7000 BC, respectively. It went extinct on the mainland probably before 6000 BC. Only C. ibex should, therefore, be restored to its former Balkan ranges. This review provides further support for the wild origin hypothesis for agrimi.     

Genetic variation in Mediterranean Helichrysum italicum (Asteraceae; Gnaphalieae): do disjunct populations of subsp. microphyllum have a common origin?

The yellow‐flowered everlasting daisy Helichrysum italicum (Asteraceae, Gnaphalieae) is widely distributed in the Mediterranean basin, where it grows in continuous and widespread populations in diverse open habitats. Helichrysum italicum subsp. microphyllum has a disjunct distribution in the Balearic Islands (Majorca and Dragonera), Corsica, Sardinia, Crete and Cyprus. Numerous morphological intermediates between subsp. italicum and subsp. microphyllum are known from Corsica, where the two subspecies co‐occur. The aims of the study were to investigate if subsp. microphyllum has a common origin, constituting an independent gene pool from subsp. italicum, or if the morphological differences between subsp. microphyllum and subsp. italicum have arisen independently in different locations from a common wider gene pool. Our analyses of AFLP, cpDNA sequences and morphological characters show that there is geographic structure to the genetic variation within H. italicum, with eastern and western Mediterranean groups, which do not correspond with the division into subsp. microphyllum and subsp. italicum as currently circumscribed. Local selection on quantitative trait loci provides sufficient explanation for the morphological divergence observed and is consistent with genetic data. Within the western Mediterranean group of the species we found considerable polymorphism in chloroplast DNA sequences among and within some populations. Comparison with chloroplast DNA sequences from other Helichrysum species showed that some chloroplast haplotypes are shared across species.     

Testing the efficiency of nested barriers to dispersal in the Mediterranean high mountain plant Edraianthus graminifolius (Campanulaceae)

Due to strong spatial heterogeneity and limited Pleistocene glaciation, the Balkan Peninsula is a major European biodiversity hot spot. Surprisingly little, however, is known about patterns and processes of intraspecific diversification of its biota in general and of high‐altitude species in particular. A well‐suited system to test hypotheses with respect to various isolating factors acting at different geographic scales and to explore full‐range phylogeographic patterns on the Balkan Peninsula is Edraianthus graminifolius (Campanulaceae), distributed in the western Balkan mountain systems, the southwestern Carpathians and the Apennine Peninsula. To this end, we used a dense population sampling and employed amplified fragment length polymorphism (AFLP) markers and plastid DNA sequences supplemented by ecological niche modelling. The strongest splits were inferred to separate southern and northern Balkan populations from the central ones, from where range extension occurred to the Carpathians and, in more recent times, once or twice to the Apennine Peninsula. The three genetic groups in the western Balkan Peninsula were remarkably congruent among molecular markers, suggesting that the barriers to gene flow acted over long time periods facilitating allopatric differentiation. Each main group of Balkan populations contained genetically and geographically distinct subgroups, which likely are the result of local refugia during warmer periods. Evidently, the topographically highly complex and during the Last Glacial Maximum only locally glaciated Balkan Peninsula is a hot spot of species richness and endemism as well as a sanctuary of intraspecific genetic diversity, even if the underlying causes remain insufficiently understood.     

Genetic diversity,genetic structure and phylogeography of the Iberian endemic Gypsophila struthium (Caryophyllaceae) as revealed by AFLP and plastid DNA sequences: connecting habitat fragmentation and diversification

Iberian gypsum outcrops are highly fragmented and ecologically challenging environments for plant colonization. As gypsophytes occur exclusively in such habitats, they are ideal models for the study of both the effects of habitat fragmentation and selection on population genetic diversity and structure. In this study, we used amplified fragment length polymorphism (AFLP) and plastid DNA sequences to investigate the phylogeographical history of the Iberian plant Gypsophila struthium (Caryophyllaceae), a widespread endemic restricted to Iberian gypsum outcrops. Gypsophila struthium consists of two subspecies that differ in the architecture of their inflorescence and have mostly allopatric ranges. Gypsophila struthium subsp. struthium occurs in central, eastern and south‐eastern Iberia, whereas G. struthium subsp. hispanica occurs in northern and eastern areas. AFLPs revealed low but significant genetic differentiation between the subspecies, probably as a result of a recent diversification during the Pliocene–Pleistocene. In the geographical contact zone between the taxa, the Bayesian analyses revealed populations with mixed ancestries and genetic clusters predominantly of one or the other subspecies, indicating incomplete reproductive barriers between them. Plastid DNA haplotypes revealed strong geographical structure and testified to processes of isolation by distance and continuous range expansion for some haplotype clades. The Bayesian analyses of the population structure of AFLP data and nested clade phylogeographical analysis (NCPA) of plastid haplotypes revealed that the putative ancestral range corresponded to central and eastern populations of G. struthium subsp. struthium, with those lineages contributing through more recent expansion to increased genetic diversity and structure of the south‐eastern and eastern ranges of this subspecies and to the diversification of G. struthium subsp. hispanica in northern and eastern gypsum outcrops. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 654–675.     

Multiple approaches to detect outliers in a genome scan for selection in ocellated lizards (Lacerta lepida) along an environmental gradient

Identification of loci with adaptive importance is a key step to understand the speciation process in natural populations, because those loci are responsible for phenotypic variation that affects fitness in different environments. We conducted an AFLP genome scan in populations of ocellated lizards (Lacerta lepida) to search for candidate loci influenced by selection along an environmental gradient in the Iberian Peninsula. This gradient is strongly influenced by climatic variables, and two subspecies can be recognized at the opposite extremes: L. lepida iberica in the northwest and L. lepida nevadensis in the southeast. Both subspecies show substantial morphological differences that may be involved in their local adaptation to the climatic extremes. To investigate how the use of a particular outlier detection method can influence the results, a frequentist method, DFDIST, and a Bayesian method, BayeScan, were used to search for outliers influenced by selection. Additionally, the spatial analysis method was used to test for associations of AFLP marker band frequencies with 54 climatic variables by logistic regression. Results obtained with each method highlight differences in their sensitivity. DFDIST and BayeScan detected a similar proportion of outliers (3–4%), but only a few loci were simultaneously detected by both methods. Several loci detected as outliers were also associated with temperature, insolation or precipitation according to spatial analysis method. These results are in accordance with reported data in the literature about morphological and life‐history variation of L. lepida subspecies along the environmental gradient.     

Allopolyploid origin of Cardamine silana (Brassicaceae) from Calabria (southern Italy): karyological, morphological and molecular evidence

Endemic Cardamine silana from Calabria (southern Italy) previously reported to be related to C. raphanifolia was found to be hexaploid. Morphological characters and AFLP data were analysed to evaluate the degree of differentiation of C. silana from closely related taxa and to find parental taxa of this polyploid. Cardamine apennina from the C. pratensis group was examined as one putative parent, as indicated in previous studies of nuclear ITS sequences, along with other related taxa based on both cpDNA and ITS sequences. Both multivariate morphometric analyses of quantitative characters and evaluation of qualitative morphological characters showed: (1) closest position of C. silana to two diploids: C. acris from the Balkan Peninsula and C. apennina from Central Italy; (2) good extent of morphological separation of C. silana from related taxa; and (3) within C. acris subspecies, C. acris ssp. vardousiae from Central Greece as closest to C. silana . Neighbour-joining tree and PCoA ordinations of AFLP data, as well as patterns of AFLP bands sharing, corroborated results of multivariate morphometrics. This evidence supports an allopolyploid origin of C. silana , with C. apennina and C. acris as parental taxa. Its origin may be dated to Pleistocene glacial events, because of the presumably wider geographical distributions of its parental taxa during more humid periods at that time.  © 2005 The Linnean Society of London, Botanical Journal of the Linnean Society , 2005, 148 , 101–116.