Abstract

The purpose of this study is to provide a new focus to contribute, from the perspective of genomic evolution, towards a better understanding of the Valerianaceae evolutionary history. Chromosome numbers were determined by Feulgen staining in 24 populations of 18 species (first count for Valerianella multidentata, –16), and DNA contents were assessed by flow cytometry in 74 populations of 35 species (first assessments in all taxa but Centranthus ruber). A molecular phylogeny based on the trnL-trnF and including 41 new sequences was established, with the first DNA sequence for Centranthus nevadensis, Valeriana rotundifolia, V. saxatilis, Valerianella multidentata, and V. turgida. This work is the first large genome size study devoted to the Valerianaceae, showing a range of DNA amounts from  pg (Valerianella turgida) to  pg (Valeriana officinalis). At the family level, changes in basic chromosome number and genome size coincide with or precede major shifts in the evolutionary history of the group, such as those concerning stamen number and floral symmetry.

1. Introduction

The family Valerianaceae (currently considered within the Caprifoliaceae s.l.; [1]) comprises ca. 400 species of which approximately 200 are included in Valeriana L., the biggest genus of the order Dipsacales. The Valerianaceae are cosmopolitan in their natural distribution, with the exceptions of Australia and the Pacific islands, where they were introduced, and can be nowadays considered as naturalized. Several of their representatives have economic interest as medicinal (e.g., V. officinalis L.), edible (e.g., Valerianella Mill., corn salad or lamb’s lettuce), or ornamental plants (e.g., Centranthus DC.). The family always roused important scientific interest, and its study has broadly benefited from the new molecular techniques, especially the phylogenetic reconstructions. Recent phylogenies based on DNA sequencing [27] have considerably modified the traditional classification of the family [812]. Molecular results recognize only six genera: Centranthus, Fedia Gaertn., Nardostachys DC., Patrinia Juss., Valeriana (including Aretiastrum DC., Astrephia Dufr., Belonanthus Graebn., Phyllactis Pers., Porteria Hook., and Stangea Graebn.), and Valerianella. Plectritis (Lindl.) DC. is nested amongst South American Valeriana [5, 6]. Moreover, Bell [7] suggested a possible further taxonomic treatment of the family considering Fedia as a synonym of Valerianella. New genera may also be described for some of the Valeriana species that do not group with their congenerics in the phylogenetic reconstructions and make the genus Valeriana paraphyletic in its current circumscription. These are Valeriana longiflora Willk., which appears closely related to the genus Centranthus, and also Valeriana celtica L. and V. hardwickii Wall., branched somewhere between Nardostachys and Valerianella plus Fedia clades [4, 5].

Morphology in Valerianaceae is of great interest for its impressive diversity of forms, mainly resulting in adaptations to a wide range of ecological conditions (from the sea board to Andean páramos at 4000 m), and concerning both vegetative and reproductive parts. This diversity has been especially well studied from the inflorescential structure perspective, with a common basic type, the thyrse, and very different forms and levels of complexity [13, 14]. Family Valerianaceae stands out in Dipsacales by presenting four different stamen numbers in a series (mainly four stamens in Patrinia and Nardostachys, three in Plectritis, Valerianella and Valeriana, two in Fedia, and one in Centranthus), while the other families show one or two distinct stamen numbers [15, 16]. These features, along with some cases of strong corolla zygomorphy (in Centranthus, with 4 : 1 petals orientation; in Fedia, with 2 : 3 petals orientation) make the family very attractive for studying the genetic base of floral characters, which has been initiated through an evolutionary-developmental approach by [16, 17]. Evo-devo studies highlight the crucial importance of duplication events in the evolution of genes involved in developmental processes (e.g., [18, 19]). These events often correspond to whole duplication genome ([20], and references therein), making the karyological and cytogenetic data essential for understanding many evolutionary processes, as, for example, the floral morphological changes.

Nevertheless, karyological and cytogenetic data in Valerianaceae are basically limited to chromosome counts. The Valerianaceae exhibit a dysploid series of five basic chromosome numbers [21]: in American Valerianella, in Nardostachys, in Patrinia and Valeriana celtica, in Centranthus, Fedia, Valeriana and Valerianella, and in Valeriana and Valerianella. Polyploidization events are common, some genera being exclusively polyploid, such as Centranthus (tetraploid) or Fedia (tetraploid and hexaploid), while species of Patrinia and Valeriana can exhibit various ploidy levels, from diploid to octoploid. Hence, both polyploidy and dysploidy seem to have played a significant role in the differentiation and evolution of these plants. Cytogenetic data (such as banding, fluorescent in situ hybridization and genome size assessment) would certainly be of a great interest for understanding the evolution of the family. Nuclear DNA amount assessments constitute a fundamental complement to chromosome counts and, in addition, a powerful tool in order to establish the relationships between closely related taxonomical groups (e.g., in Echinops L.; [22], Orobanche L.; [23]), as well as to understand the evolution within related genera from a broader point of view (e.g., Liliaceae; [24], Orchidaceae; [25]). Furthermore, genome size is a useful tool to detect possible hybrid and polyploid origins of taxa (e.g., in Carthamus L.; [26], Nicotiana L.; [27], Artemisia L.; [28]), and intraspecific variation can reveal incipient speciation [29].

The aims of the present work are to: (a) enlarge the knowledge concerning nuclear DNA amounts in the Valerianaceae, to date limited to only one report (“Kentranthus ruber Druce” ; [30]), (b) integrate these results with those of chromosome number and molecular phylogeny, and (c) provide a new focus that could contribute, from the perspective of genome evolution, to verify the hypotheses of previous works on the relationships among the above-mentioned genera.

2. Material and Methods

2.1. Plant Material

Table 1 shows the provenance of the species investigated. The sampling includes representatives of all the six genera considered in the family in regard to the results of molecular phylogenies [46]. Studied plants come from germinated cypselas, collected in the field or obtained from Botanic Gardens. Due to difficulties in germinating seeds of Valeriana representatives, whole individuals were collected in the field and cultivated in the Institut Botànic de Barcelona.

2.2. Karyological and Cytogenetic Analyses
2.2.1. Chromosome Counts and/or Determination of the Ploidy Level

Root tip meristems were obtained either by germinating seeds on wet filter paper in Petri dishes at room temperature, or from plants cultivated in pots in the greenhouse. They were pretreated with 0.002 M 8-hydroxyquinoline for 3 h at 16°C. Material was fixed in absolute ethanol, trichloromethane, and glacial acetic acid (6 : 3 : 1) and stored in the fixative at 4°C during at least two days before processing. Samples were hydrolysed in 1 N HCl for 5 min at 60°C, stained with 1% aqueous aceto-orcein for 1 h minimum, and squashed on slides into a drop of 45% acetic acid-glycerol (9 : 1). Metaphase plates were photographed with a digital camera (Zeiss AxioCam HRm) mounted on a Zeiss Axioplan microscope, and images were analysed with Axio Vision Ac version 4.2.

2.2.2. Nuclear DNA Assessments

Leaf tissue of five individuals for each studied population was chopped in 600  l of LB01 isolation buffer [33] with a razor blade, together with the chosen internal standard. The species Petunia hybrida Vilm. `PxPc6’ (2C = 2.85 pg) and Pisum sativum L. `Express long’ (2C = 8.37 pg) were used as internal standards [34] to cover the range of 2C-values found. Seeds of the standards were provided by the Institut des Sciences du Végétal (CNRS), Gif-sur-Yvette (France). For each individual, two independent samples were extracted and measured the same day. Samples were supplemented with 100  g/ml of ribonuclease A (RNase A, Boehringer) and subsequently stained with 36  l of propidium iodide (1 mg/ml) to a final concentration of 60  g/ml (Sigma-Aldrich Química), kept on ice for 20 min and measured in an Epics XL flow cytometer (Coulter Corporation). The instrument was set up with the standard configuration. Excitation of the sample was performed using a standard 488-nm air-cooled argon-ion laser at 15 mW power. Forward scatter (FSC), side scatter (SSC), and red (620 nm) fluorescence for propidium iodide were acquired. Optical alignment was based on optimized signal from 10-nm fluorescent beads (Immunocheck, Epics Division, Coulter Corporation). Time was used as a control for the stability of the instrument. Red fluorescence was projected onto a 1,024 monoparametric histogram. Gating single cells by their area versus peak fluorescence signal excluded aggregates. Acquisition was automatically stopped at 8,000 nuclei. Measurements were made at the Serveis Cientificotècnics (Universitat de Barcelona).

2.2.3. Effect of Valepotriates on Measurements

Valerianaceae contain, inter alia, valepotriates, a family of chemical compounds of great medicinal interest, thanks to which some representatives are considerably exploited by the pharmaceutical industry. Valepotriates are known DNA intercalators (they do have severe effects on DNA in e.g., PCR amplifications) and might thus possibly influence genome size measurements by altering the hydrodynamic diameter of the DNA. Furthermore, they could constitute endogenous staining inhibitors, which affect results by causing stoichiometric errors [35, and references therein]. Valepotriates are mainly stored in oil vesicles in the roots and rhizomes [36]. Therefore, we performed measurements on some roots and compared them with DNA C-values obtained for leaves of the same individuals to detect a potential effect of these chemical compounds on genome size assessments if any.

2.2.4. Statistical Analyses

ANOVA and LSD test were performed with the Statgraphics Plus 5.1 program (Statistical Graphics Corp., USA).

2.3. Phylogenetic Analyses
2.3.1. DNA Extraction, Amplification and Sequencing

Total genomic DNA was extracted from silica-dried, herbarium voucher or fresh leaves using the Nucleospin Plant extraction kit (Macherey-Nagel, GmbH & Co., Duren, Germany). PCRs were carried out with PTC100 (MJ Research, Inc.) research thermal cyclers in 25  l volume. The whole plastid trnL-trnF region [including the trnL intron, the trnL (UAA) exon, and the intergenic spacer between trnL (UAA) and trnF (GAA)] was amplified and sequenced with the universal primers trnL-c, forward, and trnL-f, reverse, and, in some cases, trnL-d, reverse, and trnL-e, forward [37]. The PCR amplification conditions used were 94°C, 1 min 35 sec; 34x (93°C, 1 min; 58°C, 1 min; 72°C, 1 min); 72°C, 10 min, and storage at 4°C. PCR products were purified with QIAquick PCR Purification Kit (Qiagen, Valencia, California, USA) or with DNA Clean & Concentrator-5 D4003 (Zymo Research, Orange, California, USA). Direct sequencing of the amplified fragments was performed using the Big Dye Terminator Cycle Sequencing v3.1 (PE Biosystems, Foster City, California, USA). Nucleotide sequencing was carried out at the Serveis Cientificotècnics (Universitat de Barcelona) with an ABI PRISM 3700 DNA Analyzer (PE Biosystems, Foster City, California, USA).

2.3.2. Sequence Assembly, Alignment, and Analyses

Nucleotide sequences were assembled and edited using MacClade 4.08 [38]. Bayesian inference (BI) was carried out with MrBayes version 3.1.2 [39]. The most appropriate nucleotide substitution models were chosen with MrModeltest version 2.3 [40]. Four Markov chains were run simultaneously for 1 × 106 generations, and these were sampled every 100 generations. Data from the first 1,000 generations were discarded as the burn-in period, after confirming that likelihood values were stabilized prior to the 1,000th generation. The 50% majority rule consensus trees and posterior probability (PP) of nodes were calculated from the pooled samples. We proceeded to an additional analysis, adding to the matrix the gap data codified with the Barriel method, as previously used for Valerianaceae in [3]. Partitioned dataset analysis was carried out for this dataset including the codified gaps. A gamma-shaped rate variation was stated for the codified data following the manufacturer’s protocol, applying the model selected with MrModeltest for the DNA data.

3. Results

3.1. Chromosome Counts

Chromosome numbers of the populations studied are indicated in Table 1, and metaphase plates are presented in Figure 1. According to our data, we provide the first count for Valerianella multidentata Loscos & Pardo. Ploidy levels ranging from diploid to octoploid were detected although diploid and tetraploid were predominantly reported. Good metaphase plates are quite difficult to obtain for Valerianaceae. In fact, this is certainly the main reason why only few chromosome studies have been carried out in the family, and that especially few chromosome pictures have been published to present.

3.2. Genome Size Assessments

Data on nuclear DNA content are presented in Figure 2 and Table 1. The present data are the first reports on nuclear DNA content in all the genera studied excepting Centranthus and in all the species studied of this genus except for C. ruber (L.) DC., which was to date the only member of the family with a known genome size [30]. Genome size in Valerianaceae varies about 21.3-fold, from 2C = 0.39 pg (Valerianella turgida Betcke, 2x) to 2C = 8.32 pg (Valeriana officinalis, 8x). Although the half peak coefficients of variation (HPCV) are rather high for some species (in all cases but one of the genus Valerianella), their mean value is 4.62% for the target plants and 2.04% for the standards. More in-deep research is suitable in Valerianella to try and obtain measurements with better HPCV.

The difference in genome size between leaves and roots measured in the same population has been addressed for Valeriana officinalis (7) (leaves: 2C of 7.67–8.97 pg and roots: 2C of 8.04–8.14 pg, at 95% confidence interval) and Valerianella dentata (L.) Pollich (2) (leaves: 2C of 0.485–0.495 pg and roots: 2C of 0.464–0.476 pg, at 95% confidence interval). A significant difference has been detected for Valerianella dentata (ANOVA, ), which may indicate a possible effect of valepotriates on genome size assessments in this species, also suggested by the especially high HPCV values found in this annual genus (Table 1).

3.3. Molecular Phylogeny

Both AIC (Akaike Information Criterion) and hLRT (hierarchical Likelihood Ratio Tests) selected GTR+G (General Time Reversible model with gamma distribution) as the best-fit model. The results of the Bayesian analyses are presented in Figure 3. Phylogenetic trees obtained through the analyses of the two different datasets basically led to the same tree topology and supports. Nevertheless, three clades are only significantly supported when codified gap information is taken into account. These are (a) Valeriana clade I (P , Figure 3), (b) Valeriana longiflora Willk. plus Centranthus clade (P , Figure 3), and (c) the grouping of Fedia cornucopiae (L.) Gaertn. AF446986 with Fedia graciliflora Fisch. & C.A.Mey. (P , Figure 3).

4. Discussion

We address here the analysis of genome size variation and its phylogenetic and taxonomical implications in the taxa considered.

4.1. The Basal Grades

They consist of the taxa showing the most ancestral characters within the Valerianaceae, namely, the genera Nardostachys and Patrinia, along with the clade grouping Valeriana celtica and Valeriana saxatilis L. (Valeriana clade III, Figure 3). Shifts in biogeography, flower morphology, and basic chromosome number occurred at this point of the evolutionary history of the family. From this moment on, the distribution area of the Valerianaceae, until that time restricted to Asia (Nardostachys and Patrinia), was enlarged to the other continents. At the same time the number of stamens decreased from four to three. These changes just preceded the change in basic chromosome numbers from to , occurring at the arising of Valerianeae tribe (Figure 3).

The Valeriana clade III ( , Figure 3) constituted by two species from the Alps, V. celtica and V. saxatilis, is branched between Patrinieae and Valerianeae, a position that was previously found for V. celtica [4, 5], and stated for the first time here for V. saxatilis. Although their genome sizes fall within the range of the remaining Valeriana, these two valerian species present as basic chromosome number [41, 42], which had only been detected in Patrinia. Furthermore, V. celtica has a yellow corolla, also like species of Patrinia. Nevertheless, species of clade III differ from Patrinia by a pappus-like Centranthus and Valeriana, and the number of stamens reduced to three as in Valeriana and Valerianella (whereas Nardostachys presents four and occasionally five stamens and Patrinia presents five or less stamens; [43, 44]). Therefore, species of clade III are from morphological, cytogenetic, and phylogenetic points of view in between Patrinieae and Valerianeae. These results raise numerous questions. Should this clade III form a new genus? Could any other Valeriana species be susceptible to join this group? Then, should this group be classified within Valerianeae, or as a new tribe (their inclusion in the Patrinieae being impossible under the monophyly criterion)?

4.2. Fedia Plus Valerianella Clade

These species have the smallest genome sizes known in the family (Table 1) to present. Fedia and Valerianella are exclusively composed of annual herbaceous plants. All the species belonging to these two genera, excepted F. pallescens (Maire) Mathez, are common weeds in waste ground and cultivated land. Valerianella, which comprises around 50 species that are mainly diploid and with three stamens, is widely distributed in Eurasia, Africa and America. The genus Fedia includes only three tetraploid species and one hexaploid subspecies [45], with two stamens, and is restricted to the Mediterranean region. Molecular inferences suggest that Fedia species should have originated from Valerianella (Figure 3; [7]), in other words that autogamous strongly associate to anthropogenic environment (Valerianella) would have given rise to allogamous, polyploid (Fedia).

Genome size in those genera shows that Fedia ( to 0.45 pg) has significantly higher DNA values than Valerianella (from to 0.31 pg) ( ), which corroborates the conclusions of Albach & Greilhuber [46] that allogamy is related to higher DNA values than autogamy. Furthermore, Fedia presents the “allogamy syndrome,” morphological adaptations linked to the breeding system (e.g., zygomorphic corolla, development of a nectar gibbosity, tube elongation, and occurrence of polychroic corolla). The intensity of this syndrome varies between the species and the subspecies of Fedia [47]. The tendency within the genus is toward an increasing genome size and degree of allogamy syndrome, from Fedia pallescens (  pg) to F. graciliflora (  pg) and F. cornucopiae (  pg).

The Valerianella group splits in two major clades (Figure 3), one including V. turgida, V. locusta (L.) Laterr., and the Fedia species. The genus Valerianella is revealed in this study as a relatively homogeneous group respecting genome size within Valerianaceae, with 2C values that range from 0.39 pg (V. turgida) to 0.61 pg (V. vesicaria Moench). Furthermore, some of Valerianella species, such as V. turgida, account for very small genomes, which fall into the range of the smallest genome size records in angiosperms [48], around the same scale than Arabidopsis thaliana (L.) Heynh. (  pg; [49]), for long considered as the lowest angiosperms C-DNA value, nowadays reported in the genus Genlisea (G. margaretae  pg; [50]).

The chromosome number of Valerianella multidentata, , is reported here for the first time. This endemic restricted to a small area of Catalonia, Spain, is very close to V. discoidea Loisel (see [51]; ), which has been in turn related to V. coronata (L.) DC. ([52]; ). All these species belong to section Coronatae. The chromosome number of 2n = 14 is the smaller within the genus. Metaphase plates of V. discoidea (Figure 1(r)) show a chromosome pair much larger than others that may result from a chromosome fusion sparking off the descending diploidy observed in the section.

4.3. Centranthus Plus Valeriana longiflora Clade

Centranthus is a circum-Mediterranean genus of nine species, characterized by flowers with one stamen, related to Valeriana longiflora (P , Figure 3; [5]). All the species show the same chromosome number, , both annuals (C. calcitrapae (L.) Dufr., C. macrosiphon Boiss.) and perennials (the other species), and both widespread (C. calcitrapae, C. ruber) and narrowly distributed species (the remaining). Centranthus is organized in three sections: Section Calcitrapae (C. calcitrapae and C. macrosiphon), sect. Centranthus (represented in our study by C. angustifolius (Mill.) DC., C. lecoqii Jord., C. longiflorus Stev. and C. ruber), and sect. Nervosa (not represented in our study). In this case, section Calcitrapae, annual (  pg), has higher genome size than section Centranthus, perennial (  pg). Annual plants have been usually reported to account for smaller genome sizes than perennials [53], as it is comprehensible that the transcriptional machinery would be more efficient in smaller genomes, in order to quickly complete the life-cycle. Otherwise, exceptions to this rule have been found in different plant groups [22, 54], and this seems to be the trend in the genus Centranthus.

The common occurrence of hybridization events in the overlapping areas of Centranthus species [55], and the fact that some species or infraspecific taxa show intermediate morphological characters, can indicate that several of these taxa could have a hybrid origin. As it has been reported in different plant groups, nuclear DNA content of hybrids corresponds to approximately the mean of both parental genome sizes, or is slightly smaller/larger than expected [26, 27, 5664]. Within the section Calcitrapae, representatives of southern Spain and Morocco have a doubtful taxonomic assignation, somewhere between C. calcitrapae and C. macrosiphon ([55] for southern Spain; e.g., J Mathez, Université Montpellier II, France, “unpubl. res.” for Morocco). One population of Morocco was assessed, C. cf. calcitrapae (5), and the 1Cx value obtained of 0.45 pg turns out to be the average between the means of typical C. calcitrapae (  pg) and C. macrosiphon (  pg), this being consistent with the hypothesis of a hybrid origin for this Moroccan population. Also, C. lecoqii is morphologically intermediate between C. angustifolius and C. ruber [55]. The 1Cx amount of C. lecoqii (0.31 pg), the mean of those of C. angustifolius (0.34 pg) and C. ruber (0.29 pg), is compatible with the hypothesis of a hybrid origin for the former taxon.

4.4. Valeriana montana Clade

This group of alpine plants is constituted by species with hermaphrodite flowers and of restricted area (V. saliunca All. and V. supina Ard. from the Alps, V. pyrenaica L. and V. apula Pourr. from the Pyrenees), two gynodioecious species widely distributed, V. tripteris L. and V. montana L. (from the mountains of Southern Europe). Valeriana rotundifolia Vill. is a gynodioecious taxon of obscure status, present in the Alps and Corsica, and morphologically close to V. tripteris and V. montana [65]. The molecular phylogeny shows a Valeriana tripteris complex constituted by the populations of V. montana, V. rotundifolia, V. supina and V. tripteris, who are forming a robust clade (P , Figure 3), with the exception of V. cf. tripteris (7) that is sister to the remaining ones. This result raises the question of the monophyly of V. tripteris, but also shows a possible interesting biogeographical pattern as the population V. cf. tripteris (7) is the only sequenced that grows at the east of the Alps. Furthermore, this population is not only divergent in terms of DNA sequence, but also for its cytogenetic characters. Although V. cf. tripteris (7) is a tetraploid (Figure 1 ma-b), its genome size  pg ranks low if compared with tetraploid V. montana (2C  pg), and even with diploid V. tripteris (2C  pg).

Intermediate forms between V. tripteris and V. montana, such as V. cf. montana (11-12), V. cf. tripteris (7) or V. rotundifolia, occur throughout the range of distribution, which makes their taxonomic delimitation difficult. This is the reason why these taxa have been sometimes considered to represent a single species with various subspecific entities, and also as two closely related species with their intermediate forms as subspecies (see [65], and references therein). Cryptic hybrid and/or polyploid taxa may largely account for the taxonomic heterogeneity of V. tripteris complex and the genome size variation observed in the group, especially at intraspecific level. However, Briquet & Cavillier [65] pointed out that intermediate forms are found in isolated populations, without V. montana, which makes their hybrid origin improbable and their polyploid origin from V. tripteris maybe more likely. Should this polyploid hypothesis be sustained, the question remains as to why tetraploid V. tripteris are morphologically similar to V. montana. One explanation could be that the speciation of V. montana could have also been induced by a polyploidization event of V. tripteris, which would have arisen previously in evolutionary time. This would explain both the more notable morphological and ecological differentiation between the two species, and in some cases the more important decrease in monoploid genome size with respect to that of the intermediate forms [23].

4.5. Valeriana officinalis Clade (Valeriana Clade II, Figure 3)

Valeriana officinalis and relatives constitute a difficult complex of numerous taxa of specific and infraspecific ranges, with doubtful morphological delimitation. Valeriana officinalis presents an unusual basic chromosome number for the genus of . This basic number is also found in V. wallrothii Kreyer, a species closely related to Valeriana officinalis in the molecular phylogeny [6], but not in other representatives of the group as, for example, V. dioica or V. tuberosa. The populations of V. officinalis measured exhibit different DNA amounts, suggesting four different DNA ploidy levels of 2x, 4x, 6x, and 8x (Table 1). Even so, it would be necessary to determine whether those differences really correspond to different ploidy levels, or if they are partially due to a high intraspecific variability within V. officinalis. This doubt concerns particularly the ploidy levels that are uncommon in the V. officinalis complex, like the hexaploid level (corresponding to 42 chromosomes), known only from V. coreana Briq. [66] and V. transjenisensis Kreyer [67]. The results compiled in the present study allow us to conclude that the complexity within V. officinalis, in which polyploidization events are largely implicated, is much more than previously considered.

5. Concluding Remarks

At the family level, changes in basic chromosome number and genome size coincide with or precede major shifts in the evolutionary history of Valerianaceae. One interesting example is the arising of strong zygomorphic flower in Centranthus and Fedia, which is in both cases consecutive to a polyploidization event. Therefore, cytogenetic studies are essential for understanding the family, and, in this sense, we will follow our effort for providing new data of this type, especially in those genera poorly studied or unknown at different levels (i.e., genome size) such as Plectritis.

Acknowledgments

The authors acknowledge S Siljak-Yakovlev for useful comments for improving the chromosome count protocol, and J Suda and three anonymous reviewers for valuable suggestions on the draft. They thank M Bou Dagher-Kharrat, R Douzet, M Casanovas, C Hidalgo, Odile Hidalgo, JM Martin, N Martin, S Pyke, ML Revel-Hidalgo, K Romashchenko, À Romo and I Soriano for their assistance with collections, N Xena de Enrech and M-B Raymúndez for helpful comments, SC Brown and O Catrice for supplying Petunia hybrida and Pisum sativum used as internal standards, R Álvarez, J Comas, R López and R Martínez for technical support in flow cytometry, M Veny for keeping the collections of living plants, and S Pyke for the amelioration of the English language. The collaboration of the botanical gardens and herbaria listed in Table 1 is also acknowledged. This work was subsidized by the Dirección General de Enseñanza Superior, Spain (Project PB 97/1134), the Ministerio de Ciencia y Tecnología, Spain (projects CGL2007-64839-C02-01 and C02-02/BOS) and the Generalitat de Catalunya (“Ajuts a grups de recerca consolidats” 2005/SGR/00344). O Hidalgo received an MICINN postdoctoral grant from the Ministerio de Ciencia e Innovación, Spain and S Garcia a JAE Doc contract from the CSIC, Spain.