Spatial and temporal determinants of genetic structure in Gentianella bohemica

The biennial plant Gentianella bohemica is a subendemic of the Bohemian Massif, where it occurs in seminatural grasslands. It has become rare in recent decades as a result of profound changes in land use. Using amplified fragment length polymorphisms (AFLP) fingerprint data, we investigated the genetic structure within and among populations of G. bohemica in Bavaria, the Czech Republic, and the Austrian border region. The aim of our study was (1) to analyze the genetic structure among populations and to discuss these findings in the context of present and historical patterns of connectivity and isolation of populations, (2) to analyze genetic structure among consecutive generations (cohorts of two consecutive years), and (3) to investigate relationships between intrapopulational diversity and effective population size (N(e)) as well as plant traits. (1) The German populations were strongly isolated from each other (pairwise F(ST)= 0.29-0.60) and from all other populations (F(ST)= 0.24-0.49). We found a pattern of near panmixis among the latter (F(ST)= 0.15-0.35) with geographical distance explaining only 8% of the genetic variance. These results were congruent with a principal coordinate analysis (PCoA) and analysis using STRUCTURE to identify genetically coherent groups. These findings are in line with the strong physical barrier and historical constraints, resulting in separation of the German populations from the others. (2) We found pronounced genetic differences between consecutive cohorts of the German populations (pairwise F(ST)= 0.23 and 0.31), which can be explained by local population history (land use, disturbance). (3) Genetic diversity within populations (Shannon index, H(Sh)) was significantly correlated with N(e) (R(S)= 0.733) and reflected a loss of diversity due to several demographic bottlenecks. Overall, we found that the genetic structure in G. bohemica is strongly influenced by historical periods of high connectivity and isolation as well as by marked demographic fluctuations in declining populations.     

Amplified fragment length polymorphism genetic fingerprinting challenges the taxonomic status of the near-endemic species Chara curta Nolte ex Kütz. (Characeae)

Amplified fragment length polymorphism (AFLP) genetic fingerprinting of 14 accessions of Chara curta and Chara aspera Willd., sampled across a range of habitats and morphologies in Britain, suggests that these taxa are part of the variation within a single species complex. Two primer combinations generating 397 fragments (97% of which were polymorphic), analysed by Jaccard's similarity coefficient and principal co-ordinate analysis, did not recover groups which reflect the current taxonomy. By contrast with the genetic study, a Gower general similarity coefficient and principal co-ordinate analysis of 52 morphological characters recovered the currently recognized species groups. A Mantel test showed no significant correlation between the genetic data and the morphological data, supporting the hypothesis that phenotypic variability in Chara L. is either to some extent environmentally induced or represents developmental stages. Implications for the conservation status of C. curta in Britain are discussed.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 467–476.     

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.     

A brief evolutionary excursion comes to an end: the genetic relationship of British species of Gentianella sect. Gentianella (Gentianaceae)

 There is much taxonomic confusion in the genus Gentianella section Gentianella, five species of which are found in Britain. Gentianella anglica is a British endemic restricted in its distribution to the chalk downland of southern England. It is considered to be threatened because of habitat loss, and due to hybridisation with the closely related, widely distributed G. amarella. The Red Data Book species, G. uliginosa, shares morphological similarities with the other two species but has a very restricted distribution. Individual plants from populations across the range of these three species were analysed using AFLPs and universal cpSSRs to determine the degree of genetic variability within and between them, and to look for evidence of hybridisation. Samples of the two other Gentianella section Gentianella species found in Britain, G. campestris and G. germanica, were also included in the study. The analysis of chloroplast SSRs was not informative for G. anglica, G. amarella and G. uliginosa, while G. campestris and G. germanica differed from each other and the other three species at a single locus. Principal co-ordinate analysis of the AFLP data revealed only three distinct groups: one group contained G. campestris samples, another contained all samples of G. germanica, and a third contained all individuals from the other three species. In mixed populations of G. anglica and G. amarella, individuals of the two species were genetically more similar to each other than they were to individuals of the same species from other populations. It is proposed that G. anglica may be an early flowering form of G. amarella that has been created and maintained as a consequence of former grassland management practices. Received May 20, 2002; accepted October 28, 2002 Published online: March 31, 2003     

Pollen morphology of the Swertiinae (Gentianaceae): phylogenetic implications

Recent molecular surveys of the Swertiinae (Gentianaceae–Gentianeae) revealed unexpected phylogenetic relationships, including polyphyly of the genera Gentianella , Jaeschkea , Lomatogonium and Swertia . To find new non-molecular characters supporting the phylogeny, we examined the exine variation of 73 species of all major lineages of subtribe Swertiinae using environmental scanning electron miscroscopy supplementing older, mainly light microscopical, studies. In contrast to previous studies, we were able to pick out taxa from phylogenetic key positions with particular focus on Swertia . Many distantly related taxa such as parts of Frasera , Gentianopsis , Halenia , Gentianella , Megacodon and several lineages of Swertia share a striate–reticulate or reticulate exine pattern. This is interpreted as the plesiomorphic character state of Swertiinae. There is also considerable variation of derived patterns; for example, different types of microechinate or almost smooth pollen was repeatedly observed in distantly related groups. Another extreme was the ring-shaped reticulation found in a North American species of Gentianopsis . Unfortunately, major relationships as revealed by molecular analyses were rarely supported because of the abundance of the plesiomorphic type and homoplasy even on low taxonomic levels; for example, within Lomatogonium . Exine variation was particularly useful in characterizing independent lineages of Swertia . For example, according to pollen characters and in agreement with other data, the Asian Swertia cuneata is a sister group of a strongly diversified African lineage and Swertia yunnanensis , which is rather aberrant in flower morphology, seems close to parts of Lomatogonium .  © 2008 The Linnean Society of London, Botanical Journal of the Linnean Society , 2008, 157 , 323–341.     

Gagea tisoniana, a new species of Gagea Salisb. sect. Gagea (Liliaceae) from central Italy

Gagea tisoniana Peruzzi et al. sp. nov. is described and its taxonomic relationship is discussed. The new species, which appears to be endemic to central Italy, is close to G. pratensis (Pers.) Dumort., G. pusilla (F. W. Schmidt) Sweet, and other mainly poorly known [i.e. G. succedanea Griseb. et Schenk, G. transversalis Stev., and G. paczoskii (Zapal.) Grossh.] or as yet undescribed taxa from eastern Europe and the eastern Mediterranean area. The recognition of the new taxon G. tisoniana is claimed on the basis of anatomical, karyological, ecological, and phytogeographical data. Maps of the Italian distribution of G. pratensis (2 n  = 60), G. pusilla (2 n  = 24), and G. tisoniana (2 n  = 24) are presented, together with an analytical identification key. Finally, G. pratensis is recorded here for the first time in Tuscany (Monte Cetona).  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 337–347.     

RAPD variation in relation to population size and plant fitness in the rare Gentianella germanica (Gentianaceae)

We investigated the distribution of genetic variation and the relationship between population size and genetic variation in the rare plant Gentianella germanica using RAPD (random amplified polymorphic DNA) profiles. Plants for the analysis were grown from seeds sampled from 72 parent plants in 11 G. germanica populations of different size (40-5000 fruiting individuals). In large populations, seeds were sampled from parents in two spatially distinct subpopulations comparable in area to the total area covered by small populations. Analysis of molecular variance revealed significant genetic variation among populations (P <0.001), while genetic variation among subpopulations was marginally significant (P <0.06). Average molecular variance within subpopulations in large populations did not differ significantly from whole-population values. There was a positive correlation between genetic variation and population size (P <0.01). Genetic variation was also positively correlated with the number of seeds per plant in the field (P <0.02) and the number of flowers per planted seed in a common garden experiment (P <0.051). We conclude that gene flow among natural populations is very limited and that reduced plant fitness in small populations of G. germanica most likely has genetic causes. Management should aim to increase the size of small populations to minimize further loss of genetic variation. Because a large proportion of genetic variation is among populations, even small populations are worth preserving.     

Mating structure and inbreeding and outbreeding depression in the rare plant Gentianella germanica (Gentianaceae)

Isolation and small size of populations as a result of habitat destruction and fragmentation may negatively affect plant fitness through pollinator limitation and increased levels of inbreeding. To increase genetic variation in small populations of rare plants artificial gene flow has been suggested as a management tool. We investigated whether pollinator limitation and inbreeding depression could reduce fitness in Gentianella germanica, an endangered biennial of increasingly fragmented calcareous grasslands in Central Europe. We experimentally excluded pollinators and generated progenies by hand-pollinating flowers with pollen from different distances. G. germanica was highly selfing. Pollinator exclusion strongly reduced seed set, indicating that pollinator limitation could potentially reduce plant fitness. Germination rate as well as number of leaves and rosette size of progeny from 10-m crosses was higher than that of progeny from open pollinations, self-, 1-m, and interpopulation crosses. After 6 mo of growth differences in the number of surviving plants persisted, whereas differences in plant size did not. The results suggest that inbreeding depression may reduce plant performance in G. germanica. Outbreeding depression in the performance of progeny from interpopulation crosses indicates that caution is necessary in using artificial interpopulation gene flow as a management tool.     

A new Pyrenean hybrid Cirsium (Asteraceae) as revealed by morphological and molecular analyses

A new nothosubspecies Cirsium  ×  vivantii is described after a molecular and morphological characterization of eight individuals collected in two separate valleys of the French central Pyrenees. Both hypervariable Amplified Fragment Length Polymorphism (AFLP) markers and nuclear rDNA (ITS) and plastid ( trn L-F, TRNT-L) DNA sequences were analysed. The profiles of these hybrid samples were compared to those of 43 individuals belonging to their presumed parental taxa C. carniolicum ssp. rufescens and C. palustre . A total of 133 AFLP bands were scored from three primer-pair combinations. All 130 AFLP bands that amplified in the hybrid samples were present in either C. carniolicum ssp. rufescens , C. palustre or both taxa, supporting the participation of these plant genomes in the resulting hybrids. Several Additive Polymorphic Sites (APS) detected in the ITS sequences of the hybrid samples also confirmed their derived origins from ribotypes of the two parental taxa. The lack of exclusive AFLP markers and the nonconcerted evolution of rDNA polymorphisms towards either of the parental ribotypes indicated their recent origin. Plastid trn T-L sequences were used to identify C. palustre as the maternal parent of six of these hybrid individuals; either parent could have acted as the plastid genome donor for the other two individuals. The morphological study revealed that all hybrid individuals were morphologically intermediate between their parents showing largely lobed and less spiny basal leaves as in C. carniolicum ssp. rufescens and decurrent leaf bases and pinkish corollas as in C. palustre .  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 154 , 421–434.     

Onosma bulgarica sp. nov. (Boraginaceae–Lithospermeae) found on serpentine in Bulgaria

A new species, Onosma bulgarica (Boraginaceae–Lithospermeae), found in the eastern Rhodope Mountains in Bulgaria is described. It is a typical serpentinitophyte with local distribution and is thus a further addition to the remarkable serpentine flora. The new species belongs to the asterotrichos Onosma species and shows similarities with other endemics distributed on the Balkan Peninsula. Onosma bulgarica is clearly morphologically delimited by its suffruticose dense caespitose habit, very narrow basal and cauline leaves, bracts of lower flowers shorter than calyx and pedicel, corolla pale yellow and glabrous and short anthers. The differences between the new species and related taxa are discussed.     

Spatial arrangement and genetic structure in Gentianella aspera in a regional, local, and temporal context

Gentianella aspera is a biennial plant of various nutrient-poor grasslands that has become rare in the landscapes outside the Alps of eastern Austria. Using AFLP fingerprinting we investigated: (1) effects of spatial structure on genetic structure in a large vineyard population that is confined to the embankments separating the grapevines; (2) temporal variation in genetic diversity and structure in this population; (3) relationships with other populations in a regional context. On the regional scale, moderate isolation by distance among populations was revealed by a Mantel test. Bayesian analysis of population structure indicated three spatially distinct gene pools and an additional one within the vineyard population. Within this population, spatial autocorrelation analysis revealed a positive correlation between genetic and spatial distance up to 50 m. Patterns found by PCoA were not in line with a priori defined subpopulations and indicated substantial gene flow across embankments. AMOVA revealed low differentiation among both the subpopulations that were found on the linear embankments and among two local groups of these subpopulations. We found, however, striking differences in the among-group variation between the 2 years, i.e., between two local groups within the generations and between those groups among generations. This was due to the highly variable larger group of the younger generation, in which an additional gene pool was identified by Bayesian analysis of population structure. Based on these results we discuss scenarios of local and regional dynamics within and among G. aspera populations.     

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.     

AFLP analysis shows high incongruence between genetic differentiation and morphology‐based taxonomy in a widely distributed tortoise

Morphology has traditionally been used to diagnose the taxa of various taxonomic ranks. However, there is growing evidence that morphology is not always able to reveal cryptic taxa, and that pronounced morphological variation could reflect phenotypic plasticity rather than evolutionary divergence. Spur‐thighed tortoises (the Testudo graeca complex), distributed in the western Palaearctic region, are characterized by high morphological variability and complicated taxonomy, which are under debate. Previous molecular studies using mainly mitochondrial DNA (mtDNA) sequences have revealed incongruence between genetic differentiation and morphology‐based taxonomy, suggesting that morphological variability is the result of phenotypic plasticity and stabilizing selection, which masks the true genealogies. In the present study, we used a range‐wide sampling and nuclear Amplified fragment length polymorphism (AFLP) markers to investigate genetic differentiation within the T. graeca complex. We found that spur‐thighed tortoises are differentiated into four geographically well‐defined AFLP groups: Balkans–Middle Eastern, western Mediterranean, Caucasian and central‐eastern Iranian. Compared with the distribution of mtDNA lineages, the groups are largely concordant, although the AFLP markers are less sensitive and distinguish fewer groups than do mtDNA sequences. The AFLP groups show an allopatric or parapatric distribution. The AFLP differentiation conflicts with the previously proposed morphology‐based taxonomy of the complex, suggesting that local adaptation to different environmental conditions may have led to the great extent of morphological variation within the same lineages. We propose a re‐evaluation of the taxa that were confirmed genetically using a thorough morphological analysis corrected for phenotypic plasticity. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ●● , ●●–●●.     

Phylogeny of the pollinating yucca moths, with revision of Mexican species (Tegeticula and Parategeticula; Lepidoptera, Prodoxidae)

The yucca moths ( Tegeticula and Parategeticula ; Lepidoptera, Prodoxidae) are well known for their obligate relationship as exclusive pollinators of yuccas. Revisionary work in recent years has revealed far higher species diversity than historically recognized, increasing the number of described species from four to 20. Based on field surveys in Mexico and examination of collections, we describe five additional species: T. californica Pellmyr sp. nov. , T. tehuacana Pellmyr & Balcázar-Lara sp. nov. , T. tambasi Pellmyr & Balcázar-Lara sp. nov., T. baja Pellmyr & Balcázar-Lara sp. nov. and P. ecdysiastica Pellmyr & Balcázar-Lara sp. nov . Tegeticula treculeanella Pellmyr is identified as a junior synonym of T. mexicana Bastida. A diagnostic key to the adults of all species of the T. yuccasella complex is provided. A phylogeny based on a 2104-bp segment of mitochondrial DNA (mtDNA) in the cytochrome oxidase I and II region supported monophyly of the two pollinator genera, and strongly supported monophyly of the 17 recognized species of the T. yuccasella complex. Most relationships are well supported, but some relationships within a recent and rapidly diversified group of 11 taxa are less robust, and in one case conflicts with a whole-genome data set (amplified fragment length polymorphism, AFLP). The current mtDNA-based analyses, together with previously published AFLP data, provide a robust phylogenetic foundation for future studies of life-history evolution and host interactions in one of the classical models of coevolution and obligate mutualism.  © 2008 The Linnean Society of London, Zoological Journal of the Linnean Society , 2008, 152 , 297–314.     

Taxonomic problems in Carthamus (Asteraceae): RAPD markers and sectional classification

Recent molecular studies have clarified the problem of the generic limits of Carthamus . However, taxonomic problems remain in the sectional classification: section Carthamus was confirmed as a natural group, but molecular phylogenies based on DNA sequences did not resolve the other sections, resulting in a confusing conglomerate. A survey based on RAPD markers has been carried out in this conflicting group of sections, with the aim of evaluating the usefulness of these markers at the section and species levels. On the basis of the results, correlated with morphological and karyological characters, some recommendations are made in relation to the taxonomy of Carthamus . The number of sections is reduced from five to two: section Thamnacanthus is currently classified as a different genus, Phonus ; section Carthamus s.s. is retained in its present delimitation, and section Atractylis is redefined, including current sections Atractylis , Lepidopappus and Odontagnathius . The commonly accepted treatment of some taxa as subspecies is not well founded, and C. alexandrinus , C. creticus , C. tenuis and C. turkestanicus should be considered as independent species. Our results confirm that methodological problems posed by RAPD markers can be avoided by careful laboratory procedures and appropriate data analyses; they also suggest that this kind of marker is useful at low taxonomic levels and is, furthermore, complementary to DNA sequence analysis.  © 2005 The Linnean Society of London, Botanical Journal of the Linnean Society , 2005, 147 , 375–383.     

Little genetic differentiation across Europe between early-flowering and late-flowering populations of the rapidly declining orchid Neotinea ustulata

Amplified fragment length polymorphism (AFLP) was used to characterize genetic diversity of the endangered Burnt Orchid, Neotinea (formerly Orchis ) ustulata . Fingerprinting of Estonian and British populations revealed surprisingly little genetic differentiation between populations but larger amounts of diversity within populations, especially in Britain. The resulting mean F _st value of 0.51 is unusually high for an orchid species. Much of the variation follows a west–east cline across Europe, whereas the much-discussed early- and late-flowering taxa of N. ustulata are considered insufficiently distinct to be viewed as separate subspecies. The later flowering N. ustulata var. aestivalis probably evolved independently on two or three occasions, each time diverging from the earlier flowering nominate race. The identity of the genes underpinning phenology in the species, and the potential selective advantages of phenological divergence, merit further study. Overall genetic diversity within populations is sufficiently high to render impoverishment an unlikely cause of their recent, precipitous decline.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 87 , 13–25.     

Chromosome counts in the genus Festuca section Eskia (Poaceae) in the Iberian Peninsula

Chromosome numbers of taxa of Festuca L. section Eskia Willk. in the Iberian Peninsula are given. The levels of ploidy for five taxa are confirmed. Idiograms and karyotypic formulae of the five taxa are presented for the first time. Two levels of ploidy occur in this section: diploid and tetraploid. One taxon, Festuca elegans ssp. merinoi is tetraploid and two other taxa have diploid and tetraploid populations. The remaining two taxa are solely diploid.  © 2004 The Linnean Society of London, Botanical Journal of the Linnean Society , 2004, 146 , 331–337.     

Taxonomic relationships of the western Asian taxa of Gentiana sect. Pneumonanthe

The taxonomic position of western Asian members of Gentiana sect. Pneumonanthe has long been a matter of conflict. In this paper, the six western Asian species currently recognized as belonging to sect. Pneumonanthe ( G. boissieri , G. calycina , G. freyniana , G. gelida , G. paradoxa , and G. septemfida ) are compared and their relationships established using a morphological phylogenetic analysis. Seed testa and flower and leaf morphological characters were studied and 11 characters were selected for a cladistic analysis. Euro-Siberian and Far Eastern taxa of sect. Pneumonanthe ( G. pneumonanthe and G. scabra ) were used as outgroups. Our results suggest the presence of two morphologically distinct clades within the western Asian gentians: a Septemfida and a Gelida clade. G. calycina and G. freyniana show close affinities to G. boissieri and G. gelida , and are distinct from G. septemfida s.l. Biogeographical aspects of the two groups are discussed.  © 2006 The Linnean Society of London, Botanical Journal of the Linnean Society , 2006, 152 , 197–208.     

Glacial history of Saxifraga paniculata (Saxifragaceae): molecular biogeography of a disjunct arctic-alpine species from Europe and North America

The molecular biogeography of the disjunctly distributed and morphologically highly variable species Saxifraga paniculata Mill. was analysed using amplified fragment length polymorphism (AFLP) and chloroplast microsatellites. The study comprised 77 samples from mountain regions in Europe and North America throughout the complete range of distribution. AFLP data revealed clear genetic differentiation between samples from the Arctic, the Caucasus, and the eastern European mountains. Samples from the Alps were divided into two groups. One group clustered with the samples from central Europe and the Pyrenees, whereas another group with individuals from southern Norway. AFLP diversity was lowest in the Arctic and highest in the Alps. Chloroplast microsatellite analysis revealed eight haplotypes but no unequivocal phylogeographical pattern. However, haplotype diversity was highest in the Alps and central Europe whereas, in the Arctic, only few widespread haplotypes could be found. The results indicate in situ survival of S. paniculata in the Caucasus, the eastern European mountains, and the Alps. The Arctic has presumably been colonized postglacially from North American refugia south of the ice shield. Southern Norway and the Pyrenees have most likely been colonized from two phylogeographically different groups in the Alps. The origin of the central European samples remains ambiguous. In situ survival seems to be as possible as several postglacial recolonization events from the Alps. The obtained molecular data clearly support the subdivision of S. paniculata into three subspecies: ssp. cartilaginea from the Caucasus, ssp. laestadii from northern Norway, Iceland, and North America, and ssp. paniculata from the other geographical regions.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 385–398.