Polymorphic microsatellite markers for the Iberian rock lizard species,Iberolacerta cyreni,and cross‐species priming in other Iberolacerta species

Eight polymorphic microsatellite loci are described for the Iberian rock lizard species, Iberolacerta cyreni. Loci were isolated from a partial genomic library that had been enriched for AAAG repeat sequence. Number of alleles per locus ranged from two to 13 in a sample of 24 individuals from the Sierra de Guadarrama (central Spain). Observed heterozygosity ranged from 0.13 to 0.96. At least three loci were amplified and polymorphic in four other Iberian rock lizard species: Iberolacerta monticola, Iberolacerta bonnali, Iberolacerta aranica and Iberolacerta aurelioi. These markers will be used to study mating strategies and reproductive success in I. cyreni.     

Phylogenetic relationships of the European lacertid genera Archaeolacerta and Iberolacerta and their relationships to some other 'Archaeolacertae' (sensu lato) from Near East, derived from mitochondrial DNA sequences

Parts of the mitochondrial genes coding for 12SrRNA and 16SrRNA (together about 960 bp) were sequenced for all Mediterranean species of 'Mountain lizards' of the genera Archaeolacerta ( sensu lato ) and Iberolacerta . All subspecies of the Iberian species Iberolacerta cyreni and I. monticola were included in this study. In addition, samples of Apathya cappadocica and Darevskia rudis were analysed to elucidate the relationships of the European 'Mountain lizards' to their possible relatives in the Near East. Maximum parsimony and neighbour joining analyses lead to the following major conclusions: (i) the monophyly of the genus Iberolacerta is very well supported; (ii) Archaeolacerta bedriagae (the type species of the genus) is most basal with respect to the ingroup taxa. If we accept Iberolacerta as a genus, Archaeolacerta becames paraphyletic. Therefore, we propose to restrict Archaeolacerta to the type species and to treat A. mosorensis and A. oxycephala provisionally as members of the collective genus Lacerta ; (iii) within the genus Iberolacerta three groups were found: a Pyrenean group, an Iberian group and I. horvathi . The relationships among these groups remain unresolved; and (iv) the Peña de Francia lizards, described originally as a subspecies of I. cyreni , are in fact more closely related to I. monticola .     

Seed image analysis and taxonomy of Diplotaxis DC. (Brassicaceae,Brassiceae)

West European Rock lizards, Lacerta (Iberolacerta) have small widely separated ranges in highland areas. Mitochondrial and nuclear DNA sequences corroborate the monophyly of the group and show it is not closely related to any of the other Rock lizards with which it was formerly placed in Archaeolacerta, an assemblage for which there is no evidence of clade status. L. (Iberolacerta) consists of four main units: L. (I.) horvathi of NW Croatia and neighbouring regions; the Pyrenees species, L. (I.) bonnali, L. (I.) aranica and L. (I.) aurelioi; L. (I.) cyreni of the Iberian Sistema Central, with distinctive populations in the Sierras de Bejar, Gredos and Guadarrama; and L. (I.) monticola of the Serra da Estrela of Central Portugal and NW Spain, this unit also contains L. (I.) cyrenimartinezricai of La Peña de Francia, W. Spain and a distinctive population in the Montañas de Sanabria. L. (Iberolacerta) has persisted in some mountain ranges for at least 4.2 ± 1.4 Ma and may have been restricted to mountains by competition from Wall lizards (Podarcis). Its clade status shows it has lost range extensively and has produced few external branches since its initial fragmentation. In contrast, Podarcis diversified about the time L. (Iberolacerta) fragmented, producing a series of widespread lineages that have persisted until the present time. The mainly European subfamily Lacertinae, to which both L. (Iberolacerta) and Podarcis belong, diversified rapidly 13–9 Ma ago, probably largely replacing other lacertid lizards of earlier origin. As another round of replacement started at approximately 9 Ma ago with the spread of Podarcis, this may be a recurrent phenomenon in the evolution of some lizard communities.     

Isolation and characterization of polymorphic microsatellite markers in Iberolacerta monticola, and cross-species amplification in Iberolacerta galani and Zootoca vivipara

Fourteen polymorphic microsatellite loci are described for the Iberian rock lizard, Iberolacerta monticola. Genetic variation in a sample of 20 individuals from Piornedo (northwestern Spain) was quantified both by the number of alleles per locus, which ranged from six to 13, and by the expected frequency of heterozygotes under random mating (heterozygosity), which ranged from 0.761 to 0.902. Single locus and global exclusion probabilities were also computed, and indicate a high power of these markers for paternity assignments and mating system studies of I. monticola. All the analysed loci were also polymorphic in Iberolacerta galani, but only seven in Zootoca vivipara.     

Candidate genes for milk production traits in Portuguese dairy sheep

In this study, the influence of the β-lactoglobulin, prolactin and α_S1-casein variants on milk production traits of three Portuguese sheep breeds, including Serra da Estrela, White Merino and Black Merino, was analyzed. The genetic variants of each marker were identified using PCR-RFLP. β-Lactoglobulin genotype AA was associated with lower milk yield in Serra da Estrela and Merino ewes. This marker also affected milk fat content in Serra da Estrela and protein content in Merino. First results regarding the influence of the sheep prolactin gene on milk production traits are reported. In the Serra da Estrela breed, ewes carrying prolactin genotype AA were significantly associated with lower milk yield, but this influence of prolactin genotypes on milk yield was not detected in Merino. Prolactin genotypes were also significantly associated with milk fat and protein contents in the Serra da Estrela breed. A suggestive effect of the α_S1-casein locus on milk yield was detected in Serra da Estrela, but no associations were found between the variants of this marker with milk fat and protein content.     

The effects of thermal biology and refuge availability on the restricted distribution of an alpine lizard

Aim In an effort to disentangle the ecological processes that confine ectotherms to alpine environments, we studied the thermoregulatory and microhabitat selection behaviours of the rock lizard Iberolacerta cyreni, which is endemic to some mountains of central Spain, and of the wall lizard Podarcis muralis, which is a potential competitor of rock lizards. Location We chose three areas in the Sierra de Guadarrama (central Spain) that differed in their thermal quality [mean deviation of environmental operative temperatures from the lizards’ preferred thermal range (PTR)] and refuge availability: a pine forest (1770 m a.s.l.) in which P. muralis was the only species found, and two mixed shrub and rock sites (1770 and 1900 m a.s.l.) where both species were present. Methods In the field we collected data on refuge availability, sun exposure, body temperature (T_b) and operative temperature (T_e). Thus, we estimated the thermal habitat quality of the areas sampled and the thermoregulation accuracy and effectiveness of both species. Results The pine forest had the lowest thermal quality and refuge availability. The lower‐elevation shrub site offered the best thermal quality, but refuges were much scarcer than at the higher‐elevation site. Both species thermoregulated accurately, because mean deviations of body temperature (T_b) from PTR were considerably smaller than those of T_e. Podarcis muralis had higher T_b values than did I. cyreni, which had similar T_b values at both shrub sites, whereas P. muralis had lower T_b values at higher elevation. Overall, the thermoregulatory effectiveness (extent to which T_b values are closer to the PTR than are T_e values) of both species was similar, but whereas I. cyreni thermoregulated more efficiently at higher elevation, the opposite was true for P. muralis. At the lower‐elevation shrub site, I. cyreni remained closer to refuges than did P. muralis. Main conclusions Our results suggest that the pine forest belt might prevent the expansion of rock lizards towards lower elevations as a result of its low thermal quality and scarcity of refuges, that the thermoregulatory effectiveness of rock lizards in alpine environments depends more on refuge availability than on thermal habitat quality, and that competition with wall lizards is unlikely to explain either the distribution or the thermoregulatory effectiveness of rock lizards.     

Blood Parasite Infection Intensity Covaries with Risk‐Taking Personality in Male Carpetan Rock Lizards (Iberolacerta cyreni)

Identifying evolutionary and developmental mechanisms underlying consistent between‐individual differences in behaviour is the main goal in ‘animal personality studies’. Here, we explored whether activity and risk‐taking varied consistently between individuals and correlated to various – potentially fitness linked – male traits in Carpetan rock lizards (Iberolacerta cyreni). Lizards showed significant consistency within both behaviours, implying the presence of activity and risk‐taking personalities. However, there were no correlation between activity and risk‐taking, neither on the between‐ nor on the within‐individual levels, implying the absence of a behavioural syndrome. We found a strong link between the intensity of blood parasite (Haemogregarinidae) infection and risk‐taking: lizards with higher infection intensity took more risk. While we cannot distinguish cause from causative in the parasite intensity – risk‐taking correlation – our results are in line with the asset protection hypothesis predicting that individuals with lower future reproductive value should focus on the current reproductive event and take higher risk.     

Behavioral buffering of global warming in a cold‐adapted lizard

Alpine lizards living in restricted areas might be particularly sensitive to climate change. We studied thermal biology of Iberolacerta cyreni in high mountains of central Spain. Our results suggest that I. cyreni is a cold‐adapted thermal specialist and an effective thermoregulator. Among ectotherms, thermal specialists are more threatened by global warming than generalists. Alpine lizards have no chance to disperse to new suitable habitats. In addition, physiological plasticity is unlikely to keep pace with the expected rates of environmental warming. Thus, lizards might rely on their behavior in order to deal with ongoing climate warming. Plasticity of thermoregulatory behavior has been proposed to buffer the rise of environmental temperatures. Therefore, we studied the change in body and environmental temperatures, as well as their relationships, for I. cyreni between the 1980s and 2012. Air temperatures have increased more than 3.5°C and substrate temperatures have increased by 6°C in the habitat of I. cyreni over the last 25 years. However, body temperatures of lizards have increased less than 2°C in the same period, and the linear relationship between body and environmental temperatures remains similar. These results show that alpine lizards are buffering the potential impact of the increase in their environmental temperatures, most probably by means of their behavior. Body temperatures of I. cyreni are still cold enough to avoid any drop in fitness. Nonetheless, if warming continues, behavioral buffering might eventually become useless, as it would imply spending too much time in shelter, losing feeding, and mating opportunities. Eventually, if body temperature exceeds the thermal optimum in the near future, fitness would decrease abruptly.     

Chronicle of an extinction foretold: genetic properties of an extremely small population of Iberolacerta monticola

The fragmentation and destruction of natural habitats by human intervention is producing a continuous and inexorable reduction of the size of populations in multitude of species all over the world. Small and isolated populations face higher extinction risks, due to demographic and environmental stochasticity, and also because of several genetic threats, among which inbreeding is considered the most important one. For many of these species, the extinction of a population is an irreversible event, so that determining the immediate importance of these risk factors and understanding their interactions is crucial for conservation plans. Iberolacerta monticola is a small lacertid endemic to the northwestern Iberian Peninsula, distributed mainly across moderate/high altitude mountainous regions. Some populations are found nearly at sea-level, though, in fluvial valleys with relict Atlantic forests, in the severely fragmented western part of its range. One of them has been dramatically reduced over the last 30 years, and presently is on the brink of extinction. Using microsatellite nuclear markers, we obtained different measures of genetic variation at this site, together with demographic and breeding data. Both the level of heterozygosity and the number of alleles per locus indicate that the level of variation in this population is comparatively high, and the average inbreeding coefficient is very low. Individuals appear healthy and long-lived, and are related by a few different lines of descent. These findings are discussed in the context of current theories and experimental evidence of associative overdominance and purging of the genetic load of populations, with special emphasis on the evolutionary potential of recovery of small evolutionary units.     

Phylogenetic relationships and gene flow between sympatric diploid and tetraploid plants ofDactylis glomerata (Gramineae)

Phylogenetic relationships between sympatric, morphologically indistinguishable diploid and tetraploid plants ofDactylis glomerata L. (Gramineae) in Galicia (Spain) were assessed using allozyme markers for 6 distinct systems. The study exploited recent introduction in Galicia and subsequent hybridization of an alien 4xDactylis subspecies possessing distinct allozymes from those of all the native plants. Opportunities for gene exchanges between the ploidies were estimated from in situ observations of flowering, examination of progenies in 2x/4x natural and experimental crosses, and enzyme analyses. Results show a high genetic similarity between the Galician diploids and tetraploids, which possess peculiar alleles in common. Although the ploidy levels usually have distinct flowering periods, interploidal crosses do occasionally occur. Gene flow is likely much more important from the diploid to the tetraploid level. A good genetic intermixing occurs between the Galician and the alien tetraploid entities which have simultaneous flowering. Autopolyploidization of the diploids followed by various rates of hybridization is proposed as one very probable origin of natural tetraploids inDactylis.     

Trophic ecology of Thoropa miliaris (Anura: Cycloramphidae) in two mountain ranges of south-eastern Brazil

Among populations of the same species, the diet composition varies depending on the environment, given the variation in diversity and abundance of available prey. Since Thoropa miliaris is a habitat specialist in reproduction, we tested the hypothesis that this specialization is also verified in its diet. We carried out the work in 10 localities in Serra do Mar and Serra da Mantiqueira, south-eastern Brazil, where we collected frogs and evaluated prey availability. The environment did not provide significant variation in the diet of T. miliaris, and Formicidae was the prey category selectively consumed in greater volume, frequency and abundance in both mountain ranges. Therefore, we can state that T. miliaris has a selective diet on ants.     

Some butterflies do not care much about topography: a single genetic lineage of Erebia euryale (Nymphalidae) along the northern Iberian mountains

The phylogeography of montane species often reveals strong genetic differentiation among mountain ranges. Both classic morphological and genetic studies have indicated distinctiveness of Pyrenean populations of the butterfly Erebia euryale. This hypothesis remained inconclusive until data from the westernmost populations of the distribution area (Cantabrian Mountains) were analysed. In the present study, we set out to describe the population structure of Erebia euryale in western Cantabria, where the species occurs in scattered localities. For this goal, we estimate the genetic diversity and differentiation found in 218 individuals from six Cantabrian (North Spain) localities genotyped by 17 allozyme loci. We also sequence 816 bp of the cytochrome oxidase subunit I mitochondrial gene in 49 individuals from Cantabrian localities and 41 specimens from five other European sites. Mitochondrial data support the recognition of four major genetic groups previously suggested for the European populations based on allozyme polymorphisms. Both mitochondrial and nuclear markers reveal genetic distinctiveness of a single Pyrenean–Cantabrian lineage of E. euryale. The lack of geographical structure and the star‐like topology displayed by the mitochondrial haplotypes indicate a pattern of demographic expansion in northern Iberia, probably related to Upper Pleistocene climatic oscillations. By contrast, within the Pyrenean–Cantabrian lineage, Cantabrian samples are genetically structured in nuclear datasets. In particular, San Isidro is significantly differentiated from the other five populations, which cluster into two groups. We recognize an evolutionary significant unit for Pyrenean–Cantabrian populations of Erebia euryale. Our results also illustrate that the evolutionary history of a species may be shaped by processes undetectable by using mtDNA alone.     

Toxin profiles of natural populations and cultures ofAlexandrium minutum Halim from Galician (Spain) coastal waters

The toxin profiles of three isolates and natural populations of the PSP agentAlexandrium minutum from several Galician rías (NW Spain) was obtained by HPLC. The toxin content of cultures ofA. minutum is dominated by GTX4 (80–90%) and GTX4 (10–15%) with small amounts of GTX3 and GTX2 (less than 3% of each); similar results were obtained for natural populations ofAlexandrium from three different Galician rías, where a mixture ofA. lusitanicum Balech andA. minutum can occur. Important quantitative differences were found between the three isolates, one being highly and two weakly toxic. The results obtained from these isolates and natural populations ofAlexandrium were very similar to those obtained from HPLC analyses of mussels intoxicated during a PSP outbreak in Ría de Ares (Rías Altas) in 1984, confirming thatA. minutum (previously identified asGonyaulax tamarensis Lebour andAlexandrium lusitanicum) was the PSP agent during the toxic outbreak in May 1984. Toxin profiles obtained from natural populations during different PSP outbreaks in different rías and from cultures are fairly consistent and suggest that at least from the toxin point of view,A. lusitanicum andA. minutum are identical, and that the toxin profile ofA. minutum from Galicia can be used as a biochemical marker.     

Plasticity of haemoglobin concentration and thermoregulation in a mountain lizard

The plastic capability of species to cope with the new conditions created by climate change is poorly understood. This is particularly relevant for organisms restricted to high elevations because they are adapted to cold temperatures and low oxygen availability. Therefore, evaluating trait plasticity of mountain specialists is fundamental to understand their vulnerability to environmental change. We transplanted mountain lizards, Iberolacerta cyreni, 800 m downhill to evaluate the plastic response in body condition, thermoregulation traits, haemoglobin level, and haemoparasite load. Initial measurements of body mass, total haemoglobin concentration ([Hb]), hematic parasite intensities, dorsal luminance, and thermoregulatory behaviour were resampled after two and four weeks of acclimation. We also tested whether an anti-parasitic drug reduced haemoparasite intensity. After only two weeks of acclimation to a lower elevation, lizards decreased 42% in [Hb], had 17% less parasite intensities, increased body condition by 25%, and raised by ~3% their mean preferred temperatures and their voluntary thermal maximum. The anti-parasitic treatment had no significant effect on the intensity of hematic parasites, but our results suggest that negative effects of haemoparasites on [Hb] are relaxed at lower elevation. The rapid plastic changes observed in thermal preferences, body condition, [Hb], and parasite intensity of I. cyreni demonstrate a potential plastic response of a mountain specialist. This may be adaptive under the climatic extremes typical of mountain habitats. However, there is uncertainty in whether the observed plasticity can also help overcome long term environmental changes.     

A short-range endemic species from south-eastern Atlantic Rain Forest shows deep signature of historical events: phylogeography of harvestmen Acutisoma longipes (Arachnida: Opiliones)

The Brazilian Atlantic Rain Forest is amongst the most diverse biomes in the world, but the processes that shaped its biodiversity are still poorly understood. We used one mitochondrial and two nuclear markers to evaluate the phylogeographic patterns of the endemic harvestman Acutisoma longipes and its closely related species to investigate the biogeographic history of this biome. The results showed low intrapopulation diversity and strong population structure, suggesting poor dispersion amongst locations. Phylogenetic analyses pointed to three main mitochondrial lineages congruent with the geomorphology of the south-eastern region of Brazil (Serra do Mar, Serra da Mantiqueira, and interior plateau). These older divergences occurred in the middle-Neogene, suggesting that events in this period drove the diversification of the species, but Quaternary events also affected the populations locally. We detected some congruence between A. longipes demographic patterns and the areas of endemism delimited for harvestmen, suggesting that some regions of the distribution could have been more stable in the past (especially in Serra da Mantiqueira). Our findings corroborate that harvestmen are a suitable group for the study of ancient biogeographic events in the Atlantic Rain Forest, even at small-scale ranges. Acutisoma hamatum is here considered as a new junior synonym of A. longipes.     

A multimetric diatom index to assess the ecological status of coastal Galician rivers (NW Spain)

There are many rivers in northwest Spain as a consequence of the mountainous landscape and the granitic geology subjected to Atlantic influences. Water and epilithic diatoms samples were collected at 72 sites in Galicia flowing into the Atlantic Ocean and Cantabrian Sea in summer 2002–2003 and spring 2004. These sites included minimally disturbed sites, defined as reference sites, and impacted sites which were influenced by different human pressures. We used the diatom assemblages to calculate diatom indices using the Omnidia software, but we also developed new metrics based on the similarity of species composition in reference sites. The response of the metrics was tested in relation to physicochemical variables. We developed a diatom multimetric index (MDIAT) as a combination of metric values. The sensitivity of the MDIAT to organic and nutrient stressors supports the use of this index to classify the ecological status of Galician rivers. The MDIAT showed higher correlations with some variables and nutrients than the individual metrics. According to the MDIAT, 69% of the sites sampled in Galician coastal rivers achieve good ecological status. The MDIAT has been developed specifically for Galician granitic rivers (NW Spain), and has been intercalibrated at the European level in the Central Baltic Rivers GIG. Our study validates the application of this multimetric index to evaluate the water quality in coastal Galician rivers.     

Molecular data reveal a complex population genetic structure for Psalidodon scabripinnis (Teleostei: Characidae) in the Atlantic Rainforest,Brazil

Recently renamed, Psalidodon scabripinnis populations of Serra da Mantiqueira, previously known as Astyanax scabripinnis have been deeply studied in the last years. These populations are small and isolated and occur very close to the watershed between Paraíba do Sul River basin and Upper Paraná River basin, in Serra da Mantiqueira region in the Atlantic Rainforest. These conditions arouse the interest in knowing theor genetic conservation status and how they responded to the separation between the two rivers basins. Therefore, we accessed the genetic diversity of five P. scabripinnis populations of this region with microsatellites and mitochondrial data. The results showed a complex structure pattern that doesn’t match the simple basin separation and a reasonably conservation status when compared with other populations of the same family or with similar natural history.     

Investigations into the late Pleistocene and Holocene history of vegetation and climate in Santa Catarina (S Brazil)

Palynological studies have been carried out on three highland peat bogs, and one situated on the Atlantic coastal plain. In the highlands, the late Pleistocene (14,000 - 10,000 uncal B.P.) vegetation was dominated by campos (grassland). Scattered stands of Araucaria forests were preserved in deep valleys. In the region of the sites at Morro da Igreja and Serra do Rio Rastro, the dominance of campos vegetation continued until about 1000 B.P. while at the Serra da Boa Vista site there was an expansion of Atlantic pluvial forest elements followed by Araucaria forests at the beginning of the Holocene. A general expansion of A. angustifolia, clearly related to a change towards an increasingly moist climate, can be dated to the present millenium. On the coastal plains, the late Pleistocene vegetation was dominated by Myrtaceae which were replaced by tropical taxa in the Holocene. The lowland profile (Poço Grande) also covers part of the upper Holocene, where the rich flora of the Atlantic pluvial forests can be characterized by taxa including Alchornea, Urticales and Rapanea. Close to the coring site, there was a repeated alternation between two different dune communities (4840 - 4590 B.P.), followed by a lake stage with aquatic plant succession (4590 - 4265 B.P.), plant communities dominated by Rapanea (4265 - 4230 B.P.) and the spread of Alchornea (4230 - 3525 B.P.). Late Pleistocene climate conditions (14,000 - 10,000 B.P.) can be described as cold and relatively dry, possibly including an equivalent of the Younger Dryas period. In the Holocene, there were changes from a warm and drier climate (10,000 -3000 B.P.) to a cool and more moist regime (ca. 3000 -ca. 1000 B.P.) and finally to a cool and very moist period (from around 1000 B.P.).     

Increased temperature disrupts chemical communication in some species but not others: The importance of local adaptation and distribution

Environmental conditions experienced by a species during its evolutionary history may shape the signals it uses for communication. Consequently, rapid environmental changes may lead to less effective signals, which interfere with communication between individuals, altering life history traits such as predator detection and mate searching. Increased temperature can reduce the efficacy of scent marks released by male lizards, but the extent to which this negative effect is related to specific biological traits and evolutionary histories across species and populations have not been explored. We experimentally tested how increased temperature affects the efficacy of chemical signals of high‐ and low‐altitude populations of three lizard species that differ in their ecological requirements and altitudinal distributions. We tested the behavioral chemosensory responses of males from each species and population to male scent marks that had been incubated at one of two temperatures (cold 16°C or hot 20°C). In high‐altitude populations of a mountain species (Iberolacerta monticola), the efficacy of chemical signals (i.e., latency time and number of tongue flicks) was lower after scent marks had been exposed to a hot temperature. The temperature that scent marks were incubated at did not affect the efficacy of chemical signals in a ubiquitous species (Podarcis muralis) or another mountain species (I. bonalli). Our results suggest that specific ecological traits arising through local adaptation to restricted distributions may be important in determining species vulnerability to climatic change.     

Systematics and evolution of the Dactylorhiza romana/sambucina polyploid complex (Orchidaceae)

The European–Mediterranean–Oriental Dactylorhiza romana/sambucina polyploid complex was studied with regard to genetic and morphological variation patterns. Allozyme and morphometric data were collected from 24 and 19 populations, respectively, initially identified as D. flavescens, D. insularis, D. markusii, D. romana, D. sambucina, and an indeterminate taxon. Genetic distances were calculated and illustrated by an unweighted pair‐group method using arithmetic averages (UPGMA) dendrogram, and principal components analyses (PCAs) were used to summarize morphological variation patterns. Another PCA was performed on combined allozyme and morphometric data. On the basis of the dendrogram and the PCA plots, main groups of populations were delimited, and the probability that each morphological character would distinguish correctly between these groups was estimated. After combining morphometric interpretations with studies of herbarium material and information from the literature, the following taxa were confidently accepted: D. romana ssp. romana, D. romana ssp. guimaraesii (comb. et stat. nov.) , D. romana ssp. georgica, D. sambucina, D. cantabrica (sp. nov.) , and D. insularis. Levels of genetic diversity suggest that D. romana s.s. is the least derived member of the complex. The evolutionary divergence of the diploid species, D. romana and D. sambucina, was probably the outcome of vicariant speciation, whereas D. romana ssp. georgica and D. romana ssp. guimaraesii appear to have evolved from D. romana s.s. through incomplete vicariant and peripheral isolate speciation events, respectively. In some populations of the diploid taxa, a significant deficiency in heterozygotes was found at one to three loci. It is proposed that this pattern may indicate a Wahlund effect, hypothesizing that local populations are subdivided into demes determined by the commonly sympatric occurrence of two distinct colour morphs combined with partial morph constancy of individual pollinators (bumblebees). Several pathways are possible for the origin of the allotriploid D. insularis and the apparently allotetraploid D. cantabrica. A taxonomic revision is provided. © 2006 The Linnean Society of London, Botanical Journal of the Linnean Society, 2006, 152 , 405–434.