Adaptive Radiation in Mediterranean Cistus (Cistaceae)

Abstract

Background

Adaptive radiation in Mediterranean plants is poorly understood. The white-flowered Cistus lineage consists of 12 species primarily distributed in Mediterranean habitats and is herein subject to analysis.

Methodology/Principal Findings

We conducted a “total evidence” analysis combining nuclear (ncpGS, ITS) and plastid (trnL-trnF, trnK-matK, trnS-trnG, rbcL) DNA sequences and using MP and BI to test the hypothesis of radiation as suggested by previous phylogenetic results. One of the five well-supported lineages of the Cistus-Halimium complex, the white-flowered Cistus lineage, comprises the higher number of species (12) and is monophyletic. Molecular dating estimates a Mid Pleistocene (1.04±0.25 Ma) diversification of the white-flowered lineage into two groups (C. clusii and C. salviifolius lineages), which display asymmetric characteristics: number of species (2 vs. 10), leaf morphologies (linear vs. linear to ovate), floral characteristics (small, three-sepalled vs. small to large, three- or five-sepalled flowers) and ecological attributes (low-land vs. low-land to mountain environments). A positive phenotype-environment correlation has been detected by historical reconstructions of morphological traits (leaf shape, leaf labdanum content and leaf pubescence). Ecological evidence indicates that modifications of leaf shape and size, coupled with differences in labdanum secretion and pubescence density, appear to be related to success of new species in different Mediterranean habitats.

Conclusions/Significance

The observation that radiation in the Cistus salviifolius lineage has been accompanied by the emergence of divergent leaf traits (such as shape, pubescence and labdanum secretion) in different environments suggets that radiation in the group has been adaptive. Here we argued that the diverse ecological conditions of Mediterranean habitats played a key role in directing the evolution of alternative leaf strategies in this plant group. Key innovation of morphological characteristics is supported by our dated phylogeny, in which a Mediterranean climate establishment (2.8 Ma) predated the adaptive radiation of the white-flowered Cistus.

About the Authors:

Beatriz Guzmán

* E-mail: [email protected]

Affiliation: Real Jardín Botánico - CSIC, Madrid, Spain

María Dolores Lledó

Affiliation: Royal Botanic Gardens, Kew, Richmond, Surrey, United Kingdom

Pablo Vargas

Affiliation: Real Jardín Botánico - CSIC, Madrid, Spain

Introduction

The concept of adaptive radiation implies a rapid ecological diversification, which should be reflected in a greater morphological and/or physiological divergence among species in brief periods of rapid diversification from a single ancestor [1], [2]. Two mechanisms could generate adaptive radiations: (1) extrinsic causes due to new environmental circumstances [3], [4]; (2) intrinsic characters of organisms (key innovation) that allow a taxon to utilize existing niche space in a novel manner [5]. Remoteness and the rich diversity of habitats of island systems help ensure little competition and different environments to test the potential of plant radiations [6], [7]. In contrast to the wealth of studies documenting adaptive radiations in oceanic islands [see 3], [8], [9] and particular mainland habitats [see 10], [11], we have found in literature no study fully focused on the Mediterranean region.

The Mediterranean climatic type, characterized by a strong seasonality (hot dry summers, cool wet winters), occurs in California, South Africa, central Chile, southern Australia, and typically in the Mediterranean Basin [12], [13]. In all five of these areas the native vegetation is a dense scrub characterized by annuals, drought-tolerant deciduous and semi-deciduous malacophyllous species, and woody evergreen sclerophyllous species [14]. Sclerophyllous species are adapted to low water availability during summer by means of small, leathery and dark leaves covered with thick cuticles and small, thick-walled cells [15]. Small leaves and low specific leaf area have been viewed as adaptations to Mediterranean-type climates in many species of evergreen plants [16]. Indeed, sclerophylly is so successful that unrelated genera and families of woody plants converged into similar leaf traits. Two alternative origins have been proposed for the evolution of Mediterranean, woody plants: resprouters corresponding to older lineages (Tertiary with tropical to subtropical conditions) and seeders (such as Cistus) to younger lineages (Quaternary with Mediterranean conditions) [17]. Few studies have, however, addressed the origin of Mediterranean plant groups by means of phylogenetic approaches related to ecological preferences [but see 18].

Significant shrub components in the European-African Mediterranean ecosystems (e.g., “maquis”, “garrigue”) belong to Cistaceae (Tuberaria, Halimium, Cistus). Cistus is a genus of 21 frutescent and suffrutescent shrub species with a predominantly Mediterranean distribution [19], except for five species endemic to the Canary Islands (Table 1). Previous phylogenetic studies revealed the separation of the Cistus-Halimium lineage and identification of two major natural groups: one of purple-flowered Cistus species (hereafter the purple-flowered lineage) and other containing the white-flowered species of Cistus, plus the pinkish-flowered C. parviflorus (hereafter the white-flowered lineage) [20], [21], [22]. Moreover, the white-flowered lineage is divided in two groups: one containing C. clusii and C. munbyi species (hereafter the C. clusii group) and other containing the rest of the white-flowered Cistus species (9), plus C. parviflorus (hereafter the C. salviifolius group) (Fig. 1). Despite the two lineages (the C. clussi and the C. salviifolius groups) are inhabiting the Mediterranean basin, the C. salviifolius group has undergone higher differentiation and displays greater variation in leaf trichome density, size, shape and tissue thickness than do the C. clusii group. These properties influence the resistance to drought stress and solar irradiance [14]. Indeed, ecological analyses of leaf morphological and physiological characters in dry environments [23], [24] appear to be related to speciation of Mediterranean plants.

[Figure omitted. See PDF.]

Figure 1. Phylogenetic hypothesis based on plastid (trnL-F, trnK-matK, trnS-trnG, rbcL) and nuclear (ITS, ncpGS) sequences.

(A) Strict consensus of 104 equally parsimonious trees of 1317 steps (CI = 0.82, RI = 0.80), showing bootstrap support for clades above branches; (B) Bayesian inference tree (50% majority rule consensus tree) showing posterior probabilities above branches.

https://doi.org/10.1371/journal.pone.0006362.g001

[Figure omitted. See PDF.]

Table 1. List of species used in the phylogenetic analysis.

https://doi.org/10.1371/journal.pone.0006362.t001

In this study, we used a molecular phylogenetic approach of DNA sequence data, sampled from both the nuclear (ITS, ncpGS) and the plastid (trnL-trnF, trnK-matK, trnS-trnG, rbcL) genomes, to test the explicit hypothesis of adaptive radiation. We first explored single ancestry in the Cistus-Halimium complex and differentiation in short periods of time by means of phylogenetic and molecular clock analyses [1], [25]. To test evolution in Mediterranean conditions, we chose a lineage exclusive to the Mediterranean basin (C. salviifolius lineage). Phenotype-environment correlation was further conducted to infer the role of ecological and vegetative characteristics [26], [27], [28] involved in speciation of this group.

Results

Phylogenetic analyses

The characteristics of the six sequence data sets are summarized in Table 2. MP analysis using Fitch parsimony resulted in 104 shortest trees of length 1317 steps (Fig. 1) for the combined sequence matrix. The consistency index (CI) for these trees was 0.82 and the retention index (RI) was 0.80. The BI tree displayed similar topology (except for the Halimium umbellatum position) and support values. Plastid and nuclear datasets yielded a similar phylogenetic pattern, although plastid sequences provided a more resolved tree (results not shown). In addition to strong (99% BS, 94 PP) support for the monophyly of the Cistus-Halimium complex, parsimony and Bayesian consensus trees were consistent at different places: (1) Cistus species were not monophyletic; (2) Cistus species were divided in two lineages, one of purple-flowered species (except C. parviflorus) (100% BS, 100 PP) and other of white-flowered species plus C. parviflorus (97% BS, 100 PP); (3) Cistus crispus was the sister-group of the rest of purple-flowered species (100% BS, 100 PP); and (4) a sister-group relationship existed between the C. clusii group (100% BS, 100 PP) and the rest of the white-flowered species plus C. parviflorus (81% BS, 100 PP). Halimium umbellatum appears to be related to the white-flowered lineage in the Bayesian analysis, but not in the MP analysis (Fig. 1).

[Figure omitted. See PDF.]

Table 2. Characteristics of each of the DNA sequence regions used in the phylogenetic analysis of Cistaceae and the white-flowered Cistus.

https://doi.org/10.1371/journal.pone.0006362.t002

Evaluating patterns of trait evolution

The range of interspecific variation in leaf morphology and ecological requirements is shown in Table 3. Character reconstruction of three morphological and three ecological characters mapped on the Bayesian consensus tree (Fig. 2) using MacClade optimization and Bayesian inference to investigate patterns of evolution. The most relevant results from the historical reconstructions are following described:

1. Leaf shape (Fig. 2A). The character state reconstruction showed linear or linear-lanceolate to elliptic leaves as a plesiomorphic state. Ovate-lanceolate and ovate shapes evolved twice in the C. salviifolius lineage.

2. Labdanum secretion (Fig. 2B). The character was equivocal in most of the C. salviifolius lineage because, in part, of missing data from two species (C. mumbyi, C. pouzolzii). A medium percentage (5–10%) of secretion per unit leaf dry weight was, however, traced as the most likely ancestral state.

3. Upper leaf pubescence (Fig. 2C). The character was revealed as very homoplastic within the C. salviifolius lineage. Despite the reconstruction was equivocal tracing the state at some nodes, independent acquisition (up to three times) of a dense tomentum is interpreted. Shifts between glabrous and subglabrous leaves appeared dynamic.

4. Soil (Fig. 2D). The historical reconstruction traced silicolous soils as the ancestral state for the C. salviifolius lineage. It was noteworthy that the only two species inhabiting basic (C. parviflorus) and ultrabasic (C. albanicus) soils within this lineage are sister species.

5. Insolation conditions (Fig. 2E). Character optimization was equivocal reconstructing the ancestral state in the C. salviifolius lineage. Two sister species groups underwent a dramatic change in insolation conditions (C. parviflorus-C. albanicus; C. populifolius-C. pouzolzii). Although ancestral character states were poorly optimised for insolation conditions, reversal to high solar exposure (helioxerophyllous) was unequivocally acquired for C. parviflorus.

6. Environment (Fig. 2F). A high frequency in habitat change was found in the C. salviifolius lineage. Similar environments were shared in a few groups with (C. pouzolzii, C. populifolius) or without (C. psilosepalus, C. parviflorus) sister relationships. In contrast, three habitats were occupied by four closely-related species (C. laurifolius, C. psilosepalus, C. parviflorus, C. albanicus) suggesting a dynamic habitat change in the course of evolution of the C. salviifolius lineage.

[Figure omitted. See PDF.]

Figure 2. Historical patterns of leaf characters and ecological attributes.

Key characteristics are mapped onto the MacClade optimization tree as inferred by the Bayesian analysis: (A) Leaf shape, (B) Leaf labdanum secretion [49], (C) Pubescence of upper leaf surface, (D) Soil requirements, (E) Insolation conditions [77], (F) Environment. Pie charts at nodes represent the posterior probabilities of Bayesian inference character state evolution (Table 2). Node coding above branches in Fig. 2A.

https://doi.org/10.1371/journal.pone.0006362.g002

[Figure omitted. See PDF.]

Table 3. Morphological and environmental characteristics of the white-flowered Cistus lineage. Data were taken from Grosser [75], Martín & Guinea [76], Dansereau [77]*, Warburg [78], Demoly and Montserrat [79], Greuter [80], Gülz et al. [49]** and own observations.

https://doi.org/10.1371/journal.pone.0006362.t003

The BayesTraits analysis of trait evolution was used to test reconstruction uncertainty. Table S2 reports ratedev settings and mean values (±95% confidence intervals) of the log-likelihood and posterior distributions of the rate of coefficients obtained from the reversible jump (RJ) MCMC analysis. The mean of the Bayesian posterior probabilities of each character state at every node (nine nodes) are provided in Table 4 and Fig. 2. The 95% confidence intervals of the posterior probabilities were all lower than±0.004. The Bayesian results mostly supported the MP (MacClade) optimization. Particular points of disagreement between both analyses were: (1) subglabre leaves at the root of the tree in the MP analysis whereas the Bayesian probability (0.52) was higher for glabre leaves; (2) ancestral states at node 2 were reconstructed as linear and subglabre leaves in the MP analysis, while the Bayesian approach estimated a higher probability for linear-lanceolate (0.68) and glabre leaves (0.68) states to be ancestral; (3) the historical reconstruction using the MP optimization traced Quercus suber/ilex and Pinus woodlands as the ancestral state at node 4 (C. ladanifer-C. salviifolius), while Quercus suber/ilex woodlands showed the highest posterior probability (0.38); (4) subheliophyllous condition was ancestral at node 5 in the MP optimization but the submesophyllous condition displayed the highest posterior probability (0.70); (5) scrub vegetation was the ancestral state at node 9 (C. parviflorus-C. albanicus) using the MacClade optimization, while Abies cephalonica woodlands displayed the highest posterior probability (0.42).

[Figure omitted. See PDF.]

Table 4. Mean of posterior probabilities of Bayesian inference character state evolution of successive iterations (9,000,000) by RJ MCMC (see text) for six characters.

https://doi.org/10.1371/journal.pone.0006362.t004

Bayesian analysis of correlated evolution

Table 5 shows the log-Bayes factor calculations and significance following the scale of Bayes factor test presented by Kass & Raftery [29]. The evolution of leaf traits was not closely associated with ancestral changes in environment and insolation conditions. There was evidence against a correlated evolution of insolation conditions to leaf shape (log-Bayes factor = −1.5) and labdanum secretion (log-Bayes factor = −1.8). Additionally, barely evidence against correlated evolution for leaf pubescence and environment has been found (log-Bayes factor = −0.1). In contrast, barely correlated evolution was suggested between three pairs of variables: leaf shape/environment (log-Bayes factor = 0.7), labdanum secretion/environment (log-Bayes factor = 0.8) and leaf pubescence/insolation conditions (log-Bayes factor = 0.6). As already discussed elsewhere for organism radiations [30], estimates of barely correlated evolution have a strong evolutionary significance considering short tree branches.

[Figure omitted. See PDF.]

Table 5. Calculations for log-Bayes factor tests in favour of a dependent model. In the final column, we followed the Bayes factor test [29] in our interpretation of the log-Bayes factor.

https://doi.org/10.1371/journal.pone.0006362.t005

Haplotype analysis of the white-flowered Cistus lineage

Sequence length of the white-flowered Cistus lineage was 417–461 bp for trnL-trnF, 561–585 for tnrS-trnG, 1309–1357 for trnK-matK and 1378–1379 for rbcL (Table 2). The combined data of plastid sequences for 10 species (13 taxa) of the C. salviifolius lineage distinguished only 12 substitution-based haplotypes (Table S3). Haplotypes were exclusive to a single species or subspecies (Table S3), except for one for both C. ladanifer subspp. ladanifer and sulcatus. TCS constructed a single, star-like network (Fig. 3) displaying no loops. This analysis is congruent with a multiple lineage divergence pattern from ancestral haplotypes, as expected in a radiation.

[Figure omitted. See PDF.]

Figure 3. Statistical parsimony network representing relationships of the 12 plastid (trnL-trnF, trnK-matK, rbcL, trnS-trnG) haplotypes of the white-flowered Cistus lineage.

Lines indicate mutation steps (single nucleotide substitutions) and dots (•) represent missing haplotypes (extinct or not found). A star-like shape of the network is congruent with a process of radiation in this group.

https://doi.org/10.1371/journal.pone.0006362.g003

Estimates of divergence times

Results of the dating analysis are shown in Table 6 and Fig. 4. In general, the data indicated a Pliocene-Pleistocene (2.11±0.87 Ma) divergence between the basal-most Halimium and Cistus-Halimium groups, followed by a Pleistocene differentiation of the major clades of the latter group. An ancestor shared by Halimium umbellatum and Cistus appeared to have diverged after the Pliocene-Pleistocene boundary (1.47±0.35 Ma). Short branch lengths may reflect a rapid divergence process in the white-flowered lineage. An early divergent lineage of C. clusii and C. mumbyi (C. clusii lineage) at 1.04±0.25 Ma was followed by differentiation of 10 species (C. salviifolius lineage) in the Mid Pleistocene (0.88±0.22 Ma). The average per-lineage species diversification rate for the C. salviifolius lineage was 1.46–2.44 species per million years.

[Figure omitted. See PDF.]

Figure 4. Phylogenetic chronogram of the Cistus-Halimium complex based on the Bayesian consensus tree.

Fossil calibration points are indicated on the tree. Shaded area delineates the establishment of the Mediterranean climate 2.8 million years ago [50]. Geological timescales are shown both at the top and the bottom. Photographs illustrate diversity in leaf morphology of the white-flowered Cistus species (only subsp. ladanifer of C. ladanifer, subsp. clusii of C. clusii and subsp. populifolius of C. populifolius are illustrated). Species insolation conditions [77] are plotted on the right side of the tree (, helioxerophyllous and subhelioxerophyllous;, subsciophyllous and submesophyllous).

https://doi.org/10.1371/journal.pone.0006362.g004

[Figure omitted. See PDF.]

Table 6. Penalized Likelihood (bootstrapping of 100 trees) molecular clock estimates of ages for constrained and unconstrained nodes.

https://doi.org/10.1371/journal.pone.0006362.t006

Discussion

An adaptive radiation comprises a group of species that inhabit a variety of environments, differ in morphological and other traits important in utilizing these environments, and are descended from a common ancestor that rapidly speciated over a short period of time [1]. Available phylogenetic and ecological evidence suggests that the C. salviifolius lineage of 10 white-flowered species meets the four criteria to strictly test adaptive radiation: common ancestry, rapid speciation, phenotype-environment relationships and trait utility [1], [4].

Monophyly of the white-flowered Cistus lineage is strongly supported irrespective of phylogenetic methods and DNA sequences used (Fig. 1). Accordingly, the 12 white-flowered species form a well-defined natural group and fulfil the common ancestry condition.

The concept of rapid speciation is not very well defined, even though a considerable number of species is needed [1]. Asymmetry between sister clades in their number of descendant species is one of the operational standards to distinguish speciation bursts from stochastic background rates [31]. Compared with the two taxa included in the C. clusii lineage, the remaining taxa (13) form a sister group (C. salviifolius lineage) and can be consider as a significant burst. In fact, asymmetries between both lineages can also be inferred in a temporal pattern. After the split of the most common recent ancestor of the two lineages (1.04±0.25 Ma), a relatively long period of time was necessary to bring about limited (2) extant species in C. clusii lineage, in contrast to the 10 species generated in the C. salviifolius group (Fig. 4). Alternatively, rapid radiation is also interpreted as high rates of differentiation in comparison to those of flowering plants. The estimated rate of diversification in the C. salviifolius lineage was significantly higher (1.46–2.44 species per million years) compared to the median rate of diversification of angiosperm families (0.12 species per million years; with a maximum of 0.39) [32] and to that found in the Andean Valeraniaceae [33], and similar to the explosive radiation described for Andean Lupinus [34]. Rapid diversification in the C. salviifolius lineage was already predicted by a combination of different sources of evidence prior to performing explicit analysis of radiation: (1) lack of resolution and low support values depicted mainly in the parsimony-based tree because of a low number of polymorphisms [20], which is overcome by increasing the number of DNA substitutions (this paper); (2) short branch lengths and low pairwise sequence divergence (Fig. 4, Table 2); (3) low resolution at the core of the haplotype network [35] (Fig. 3).

In addition to evidence for common ancestry and rapid diversification, the fit of the diverse phenotypes observed in a lineage with their environment is necessary in prediction of adaptive radiations [1]. Our character reconstruction suggests that shifts in leaf features allowing the colonization of different habitats have been related with specific speciation events (Fig. 2). Acquisition of diverse leaf features is associated with recent lineage splits, and thus closely related taxa exhibit different leaf morphologies [26]. Our character state optimization reveals that the most common recent ancestors of four sister species diversified in different environmental conditions (Fig. 2D, 2E, 2F) by means of shifts in leaf shape (Fig. 2A) and leaf pubescence (Fig. 2C). In addition, transition in leaf labdanum secretion is observed in the Bayesian inference (Table 4; Fig. 2B). Trends of correlated evolution between leaf traits and at least one ecological trait (environment, insolation conditions) have been found (Table 5). The barely correlated evolution found in three morphological/ecological traits indicates that shifts in environmental conditions must parallel evolutionary changes in Cistus leaf morphology as a whole and not in individual leaf features. Experimental studies testing correlated evolution of all leaf traits should be further performed to analyse compensatory effects (trade-off). Alternatively to sister species approaches, another strong indication of the adaptive value of a trait is when phylogenetically separate, but ecologically similar, species converge or show parallel patterns of variation along similar ecological gradients [36]. Multiple leaf morphological character-states studied across the white-flowered lineage (shape, labdanum secretion, pubescence) have been independently acquired at least twice (Fig. 2).

Evidence that some morphological and/or physiological traits of species are particularly useful is the fourth necessary condition to support the most strict concept of adaptive radiation: trait utility [1]. The adaptive implications of leaf size and shape differences are well documented [37], [38]. In absence of explicit experimental studies (plant translocation, common garden conditions) for all species involved in this adaptive radiation [39], the body of knowledge for particular traits is analysed. Although our six DNA sequence data set rendered certain phylogenetic uncertainty for some sister species relationships because of moderate support (Fig. 1), the most plausible hypothesis allows assessing low character-state reconstruction uncertainty of leaf morphological utility using MacClade optimization and BayesTraits analysis of trait evolution. Leaf size and shape are implicated in important aspects as thermoregulation [40], [41], efficiency of water use [23], [42], photosynthetic potential [43], branching and rooting strategies [44], among others. Moreover, comparative studies have revealed the existence of well-marked ecological and leaf morphological trends [45]. Small leaf size (specifically narrow leaves) are generally favoured under high exposure and/or low water availability as they help to maintain favourable leaf temperature and improves water use efficiency [42], [44]. Indeed, small-leaved species are concentrated at the high exposure end on south-facing slopes in Mediterranean garrigue and Californian chaparral [46]. Although our character reconstruction hypothesis indicates dynamic shifts of leaf shapes, the ancestral state (narrow leaves) appears to have evolved early into linear-lanceolate to elliptic, and then into ovate (plus ovate-lanceolate) leaves independently four times. In fact, helioxerophyllous species (C. clusii, C. munbyi, C. libanotis) show ancestral linear revolute leaves while submesophyllous species with broad, flat leaves inhabit shadier environments (C. laurifolius, C. populifolius, C. psilosepalus) (Table 3, Fig. 2). Leaf shape is not, however, the only phenotypic trait associated with adaptation to dry conditions. Leaf pubescence is reported to be an adaptation to sunnier and hotter environments by reducing transpiration, increasing the probability of water uptake by leaves, maintaining favourable leaf temperature, and protecting against UV-B radiation responsible for photosynthetic inhibition [47], [48]. Accordingly, a combination of leaf trait strategies meets in unrelated species of Cistus. Sister species within the white-flowered Cistus lineage have different leaf traits related to leaf transpiration. For instance, the ovate leaves of C. parviflorus unsuitable for xeric environments are protected by a dense tomentum of stellate hairs. The same is true at a lower extent in C. salviifolius. In addition, leaves can be highly reflective in the visible spectrum by covering the upper surface with labdanum, and then decreasing transpiration [49]. The high leaf secretion of resins (labdanum) in the linear-lanceolate leaves of C. monspeliensis and C. ladanifer may confer a trade-off compared to the narrower leaves of C. clusii, C. mumbyi and C. libanotis, which display linear leaves and lower labdanum concentration (Fig. 2). Further studies are needed to pinpoint whether combination of multiple leaf strategies are equally fit in dry, Mediterranean habitats suffering from dry hot summers and high solar radiation.

In summary, the evolutionary history of the 10 species (13 taxa) of the C. salviifolius lineage fits into utilization of the niche space in a novel manner far after the Mediterranean climate establishment [50]. A Mediterranean Cistus ancestor with linear, medium labdanum content and glabrous or subglabrous leaves may have spawned new lines of evolution exploiting six pre-existing Mediterranean habitats. Multiple leaf strategies were successfully essayed in the course of speciation to occupy particular environments and become part of the dominant element in the Mediterranean scrub. As far as we know, this is the first documented plant group involved in an adaptive radiation process in the Mediterranean region.

Materials and Methods

Sample strategy and DNA sequencing

A total of 36 individuals representing the 21 species of Cistus, one of Fumana, eight of Halimium, one of Helianthemum and one of Tuberaria was sequenced for four plastid (trnL-trnF spacer, trnS-trnG spacer, trnK-matK spacer, rbcL exon) and two nuclear (ITS, ncpGS) DNA regions (Table 1; Table S1) to perform phylogenetic analyses and estimate divergence times of Cistus and related lineages. In addition, a data set comprising only the white-flowered Cistus species (plus C. parviflorus) was used to infer character evolution, correlated evolution and haplotype analyses.

Standard primers were used for amplification of the ITS region [51 for 17SE], [52 for ITS4], the trnL(UAA)-trnF(GAA) [53], the trnK-matK [trnK-3914F and matK-1470R, 54] and the trnS (GCU)-trnG (UCC) [55] spacers. The rbcL exon was amplified in two overlapping segments using the following primer combination: 1F-724R and 636F-1460R [56]. A portion of the glutamine synthetase (ncpGS) was amplified for the first time in 11 Cistus species with the universal primers Gscp687f and Gscp856r [57]. To ensure a homogeneous amplification reaction we design two 24-nucleotide-long primers specific for amplifying and sequencing Cistus species (CIS-687f: 5′GTAGCTGGAATCAACATCAGTGG3′ , CIS-856r: 5′GCTTGTTCAGTGATTCTCTGTCAG3′).

After 1–3 min pretreatment at 94°C, PCR conditions for amplification were: 24–39 cycles of 1 min at 94°C, 30 s-1 min at 48–50–55°C and 1–4 min at 72°C (for details see 19). PCR primers were used for cycle sequencing of the spacers, the rbcL exon and the ncpGS gene while the ITS 5 and ITS 4 [52] primers were used for cycle sequencing the ITS region. Additionally, due to mononucleotide repeat stretches (poly-T, poly-A) the internal primer trnSGpolyTf (5′TTAGATTCTATTTACATTCT3′) was used to sequence the trnS-trnG spacer in the purple-flowered species. Sequenced data were assembled and edited using the program Seqed (Applied Biosystems, California). The limits of the regions were determined by position of flanking primers. IUPAC symbols were used to represent nucleotide ambiguities.

Molecular analyses

Phylogenetic analyses.

Maximum Parsimony (MP) and Bayesian Inference (BI) analyses were performed on a combined molecular data set of trnL-trnF, trnS-trnG, trnK-matK, rbcL, ITS and ncpGS sequences. Sequences were aligned using Clustal X 1.62b [58], with further adjustments by visual inspection. All parsimony analyses were conducted using Fitch parsimony [59] with equal weighting of all characters and of transitions/transversions. Heuristic searches were replicated 1000 times with random taxon–addition sequences, tree–bisection–reconnection (TBR) branch swapping, the options MulTrees and Steepest Descent in effect and holding 10 trees per replicate. Internal support was assessed using 5,000,000 bootstrap (BS) replicates [fast stepwise-addition, 60].

To determine the simplest model of sequence evolution that best fits the sequence data, the Hierarchical Likelihood Ratio Test (hLRT) and Akaike Information Criterion (AIC) were implemented using MrModeltest 1.1b [61], [62] in each data set. A Bayesian Inference analysis (BI) was conducted in MrBayes 3.0b4 [63] using two identical searches with two million generations each (four MCMC, chain temperature = 0.2; sample frequency = 100). In both runs probabilities converged at the same stable value after generation 100,000 approximately. A 50% majority-rule consensus tree was calculated using the sumt command to yield the final Bayesian estimate of phylogeny. We used posterior probability (PP) as an estimate of robustness.

Molecular dating and diversification rates.

Divergence dates were estimated for nodes of the Bayesian consensus tree. To check the constancy of substitution rates we used the Langley and Fitch (LF) test [64]. We rejected the null hypothesis of constant rate (χ2 = 5204.26; d.f. = 34) and, then, divergence times were estimated using the r8S 1.71 program [65] with a Penalized Likelihood (PL) approach. Penalized Likelihood was implemented with the Truncated Newton (TN) algorithm. Initial results were obtained under the following parameters: cvstart = 0.5; cvinc = 0.5; cvnum = 10 with cross-validation enforced to estimate the rate smoothing parameter (measure of the rate variation and autocorrelation of rates from clade to clade). The rate smoothing with the lowest crossvalidation score was selected and the dating procedure was repeated with the following parameters: collapse; num_time_guesses = 5 and num_restarts = 5. Crossvalidation suggested 10 as the best smooth parameter. Branching order and branch lengths from 100 Bayesian trees sampled every 10,000 generations after stationary were analyzed to obtain means and standard deviations of clade ages [34]. To convert relative divergence times into absolute time units we used two maximum-age fossil constraints. Palynological studies identified Helianthemum pollen in Upper Miocene formations (11 Ma) from France [66] and Tuberaria pollen in Pliocene formations (5.3 Ma) from Germany [67].

Species diversification rates, assuming an equal rate of random speciation Yule model, were calculated using the formula SR = [(loge(N)–loge(N0)]/T [34], [68], [69], where N is the total number of extant species in the clade of interest, N0 is the initial species diversity, usually taken as 1, and T is the inferred age of the clade (million years). Upper and lower standard deviations of age estimates were used in calculations of speciation rates.

Character evolution.

Patterns of evolution of six key traits (leaf shape, leaf labdanum secretion, leaf pubescence, soil requirements, insolation conditions, habitat) were explored in the white-flowered Cistus lineage using the Bayesian consensus tree (calculated using the same parameters as above). Optimizations were performed in MacClade 4.06 [70] assuming Fitch Parsimony, equal weighting of all characters, transitions among all states equally probable and treating characters as unordered. Character states were determined from literature and personal observations. Samples of Cistus crispus and C. heterophyllus were used as outgrop sequences.

In addition, to account for values of phylogenetic mapping uncertainty, probabilities of ancestral states for the six traits were estimated individually using the BayesMultiState program, contained in the BayesTraits 1.0 package [71], under Montecarlo Markov Chain (MCMC) method and allowing transitions between character states in both directions. To reduce the autocorrelation of successive samples, 1000 trees were drawn from the distribution of 1.9×106 trees, which equates to sampling every 1900th generation of the chains used in the phylogenetic analysis. As suggested in BayesMultiState manual, to reduce some of the uncertainty and arbitrariness of choosing prior in MCMC studies, we used the hyperprior approach, in concrete the reversible-jump (RJ) hyperprior with a gamma prior (mean and variance seeded from uniform distributions on the interval 0 to 10). Preliminary analyses were run to adjust the ratedev parameter until the acceptance rates of proposed changes was around 20–40%. Using ratedev settings (Table S2), we ran the RJ MCMC analyses for each trait three times independently for 1.0×107 iterations, sampling every hundredth iteration (to produce 90,000 sampled points) and discarding the first 1,000,000 iterations. All runs gave mostly the same results and we report one of them here. We use the “Addnode” command to find the proportion of the likelihood associated with each of the possible states at each node.

Testing correlated evolution.

We modelled correlated evolution of discrete binary traits (leaf shape/insolation conditions, leaf shape/habit, labdanum secretion/insolation conditions, labdanum secretion/habit, leaf pubescence/insolation conditions, leaf pubescence/habit) on 1000 Bayesian trees using the BayesDiscrete program, contained in the BayesTraits 1.0 package [71] and the same parameters described above. The method compares the statistical likelihood of a model in which two binary traits are allowed to evolve independently on the tree, with a model in which the two traits are allowed to evolve in a correlated fashion. Evidence for correlated evolution arises if the dependent or correlated model shows significantly better fit to the data than the independent model. As the independent and dependent models are estimated by MCMC, their goodness of fit is compared using the log-Bayes Factor test: 2*log[harmonic mean(dependent model)] – log[harmonic mean (independent model)].

We used one sample per species of the white-flowered lineage given the monophyly of all species [72]. As binary traits are required, we coded traits as followed: leaf shape, 0 linear to elliptic, 1 ovate-lanceolate to ovate; labdanum secretion, 0 zero to eight percent, 1 nine to fifteen; upper leaf pubescence, 0 glabre to subglabre, 1 dense tomentum; insolation conditions, 0 helioxerophyllous to subheliophyllous, 1 subsciophyllous to submesophyllous; environment, 0 bushy and scrub vegetation, 1 woodlands.

Haplotype data analysis.

Sequences of plastid DNA (trnL-trnF, trnK-matK trnS-trnG and rbcL) were combined to analyze relationships among the white-flowered Cistus (plus C. parviflorus) plastid haplotypes. We used the software TCS 1.21 to infer plastid haplotype ancestry [73]. The program implements a statistical parsimony approach using the algorithm described in Templeton et al. [74] to construct haplotype networks. The maximum number of differences among haplotypes, as a result of single substitutions, was calculated with 95% confidence limits and treating gaps as missing data.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

GenBank accession numbers.

https://doi.org/10.1371/journal.pone.0006362.s001

(0.10 MB DOC)

Table S2.

Bayesian inference of trait evolution of successive iterations of the chain (9,000,000) in the white-flowered Cistus lineage by reversible jump Markov chain Monte Carlo. Means±confidence intervals (95%) of the log-likelihoods (Lh) and rate coefficients are shown.

https://doi.org/10.1371/journal.pone.0006362.s002

(0.07 MB DOC)

Table S3.

List of haplotypes found in 16 species and subspecies of the white-flowered Cistus lineage. Variable sites of the sequences of four plastid DNA regions (trnL-trnF, rbcL, trnK-matK, trnS-trnG) are shown. Nucleotide position for each data set is numbered from the 5′ to the 3′ DNA ends.

https://doi.org/10.1371/journal.pone.0006362.s003

(0.20 MB DOC)

Acknowledgments

We thank E. Cano and N. Rønsted for lab assistance; P. Escobar, A. Fernández, O. Filippi, J. Leralta, J. Martínez and R. G. Page for plant materials; I. Sanmartín for advice on BayesTraits analyses; J. J. Aldasoro for computer assistance.

Author Contributions

Conceived and designed the experiments: PV. Performed the experiments: BG MDL. Analyzed the data: BG PV. Wrote the paper: BG PV.

Citation: Guzmán B, Lledó MD, Vargas P (2009) Adaptive Radiation in Mediterranean Cistus (Cistaceae). PLoS ONE4(7): e6362. https://doi.org/10.1371/journal.pone.0006362

References

1. Schluter D (2000) The ecology of adaptive radiation. New York: Oxford University Press Inc. D. Schluter2000The ecology of adaptive radiation.New YorkOxford University Press Inc

2. Gavrilets S, Losos JB (2009) Adaptive radiation: contrasting theory with data. Science 323: 732–737.S. GavriletsJB Losos2009Adaptive radiation: contrasting theory with data.Science323732737

3. Baldwin BG, Robichaux RH (1995) Historical biogeography and ecology of the Hawaiian silversword alliance (Asteraceae): new molecular phylogenetic perspectives. In: Wagner WL, Funk VA, editors. Hawaiian biogeography: evolution on a hotspot archipelago. Washington, D.C.: Smithsonian Institution Press. pp. 257–287.BG BaldwinRH Robichaux1995Historical biogeography and ecology of the Hawaiian silversword alliance (Asteraceae): new molecular phylogenetic perspectives.WL WagnerVA FunkHawaiian biogeography: evolution on a hotspot archipelagoWashington, D.C.Smithsonian Institution Press257287 In:. In:.

4. Meimberg H, Abele T, Brauchler C, McKay JK, Perez de Paz PL, et al. (2006) Molecular evidence for adaptive radiation of Micromeria Benth. (Lamiaceae) on the Canary Islands as inferred from chloroplast and nuclear DNA sequences and ISSR fingerprint data. Mol Phylogenet Evol 41: 566–578.H. MeimbergT. AbeleC. BrauchlerJK McKayPL Perez de Paz2006Molecular evidence for adaptive radiation of Micromeria Benth. (Lamiaceae) on the Canary Islands as inferred from chloroplast and nuclear DNA sequences and ISSR fingerprint data.Mol Phylogenet Evol41566578

5. Hodges SA (1997) Rapid radiation due to a key innovation in columbines (Ranunculaceae: Aquilegia). In: Givnish TJ, Sytsma KJ, editors. Molecular evolution and adaptive radiation. New York, USA: Cambridge University Press. pp. 391–405.SA Hodges1997Rapid radiation due to a key innovation in columbines (Ranunculaceae: Aquilegia).TJ GivnishKJ SytsmaMolecular evolution and adaptive radiationNew York, USACambridge University Press391405

6. MacArthur RH, Wilson EO (1967) The theory of island biogeography. Princeton, New Jersey: Princeton University Press. RH MacArthurEO Wilson1967The theory of island biogeography.Princeton, New JerseyPrinceton University Press

7. Baldwin BG, Crawford DJ, Francisco-Ortega J, Kim S-C, Sang T, et al. (1998) Molecular phylogenetics insights on the origin and evolution of oceanic islands plants. In: Soltis DE, Soltis PS, Doyle JJ, editors. Molecular systematics of plants II DNA sequencing. Boston/Dordretch/London: Kluwer Academic Publishers. pp. 410–441.BG BaldwinDJ CrawfordJ. Francisco-OrtegaS-C KimT. Sang1998Molecular phylogenetics insights on the origin and evolution of oceanic islands plants.DE SoltisPS SoltisJJ DoyleMolecular systematics of plants II DNA sequencingBoston/Dordretch/LondonKluwer Academic Publishers410441

8. Kim S-C, Crawford DJ, Francisco-Ortega J, Santos-Guerra A (1996) A common origin for woody Sonchus and five related genera in the Macaronesian islands: evidence for extensive radiation. Proc Natl Acad Sci USA 93: 7743–7748.S-C KimDJ CrawfordJ. Francisco-OrtegaA. Santos-Guerra1996A common origin for woody Sonchus and five related genera in the Macaronesian islands: evidence for extensive radiation.Proc Natl Acad Sci USA9377437748

9. Givnish TJ, Millam KC, Mast AR, Patterson TB, Theim TJ, et al. (2008) Origin, adaptive radiation and diversification of the Hawaiian lobeliads (Asterales: Campanulaceae). Proc R Soc Lond [Biol]. TJ GivnishKC MillamAR MastTB PattersonTJ Theim2008Origin, adaptive radiation and diversification of the Hawaiian lobeliads (Asterales: Campanulaceae).Proc R Soc Lond [Biol]

10. Reinthal PN, Meyer A (1997) Molecular phylogenetic tests of speciation models in lake Malawi cichlid fishes. In: Givnish TJ, Sytsma KJ, editors. Molecular evolution and adaptive radiation. Melbourne, Australia: Cambridge University Press. pp. 375–390.PN ReinthalA. Meyer1997Molecular phylogenetic tests of speciation models in lake Malawi cichlid fishes.TJ GivnishKJ SytsmaMolecular evolution and adaptive radiationMelbourne, AustraliaCambridge University Press375390

11. Verboom GA, Linder HP, Stock WD (2004) Testing the adaptive nature of radiation: Growth form and life history divergence in the African grass genus Ehrharta (Poaceae: Ehrhartoideade). Am J Bot 91: 1364–1370.GA VerboomHP LinderWD Stock2004Testing the adaptive nature of radiation: Growth form and life history divergence in the African grass genus Ehrharta (Poaceae: Ehrhartoideade).Am J Bot9113641370

12. Daget P (1977) Le bioclimat méditerranéen caractères generaux, mode de caracterisation. Vegetatio 31: 1–20.P. Daget1977Le bioclimat méditerranéen caractères generaux, mode de caracterisation.Vegetatio31120

13. Nahal I (1981) The Mediterranean climate from a biological viewpoint. In: de Castri F, Goodall DW, Specht RL, editors. Ecosystems of the World 11, Mediterranean-type shrublands. Amsterdam: Elsevier Scientific Publising Company. pp. 63–86.I. Nahal1981The Mediterranean climate from a biological viewpoint.F. de CastriDW GoodallRL SpechtEcosystems of the World 11, Mediterranean-type shrublandsAmsterdamElsevier Scientific Publising Company6386

14. Mooney HA, Dunn EL (1970) Photosynthetic systems of Mediterranean-climate shrubs and trees of California and Chile. Am Nat 104: 447–453.HA MooneyEL Dunn1970Photosynthetic systems of Mediterranean-climate shrubs and trees of California and Chile.Am Nat104447453

15. Read J, Sanson GD (2003) Characterizing sclerophylly: the mechanical properties of a diverse range of leaf types. New Phytol 160: 81–99.J. ReadGD Sanson2003Characterizing sclerophylly: the mechanical properties of a diverse range of leaf types.New Phytol1608199

16. Ackerly DD (2004) Adaptation, niche conservatism, and convergence: comparative studies of leaf evolution in the California chaparral. Am Nat 163: 654–671.DD Ackerly2004Adaptation, niche conservatism, and convergence: comparative studies of leaf evolution in the California chaparral.Am Nat163654671

17. Herrera CM (1992) Historical effects and sorting processes as explanations for contemporary ecological patterns: character syndromes in Mediterranean woody plants. Am Nat 140: 421–446.CM Herrera1992Historical effects and sorting processes as explanations for contemporary ecological patterns: character syndromes in Mediterranean woody plants.Am Nat140421446

18. Pausas JG, Verdu M (2005) Plant persistence traits in fire-prone ecosystems of the Mediterranean basin: A phylogenetic approach. Oikos 109: 196–202.JG PausasM. Verdu2005Plant persistence traits in fire-prone ecosystems of the Mediterranean basin: A phylogenetic approach.Oikos109196202

19. Arrington JM, Kubitzki K (2003) Cistaceae. In: Kubitzki K, editor. The Families and Genera of Vascular Plants IV Flowering Plants Dicotyledons Malvales, Capparales and Non-betalain Caryophyllales. Berlin: Springer. pp. 62–70.JM ArringtonK. Kubitzki2003Cistaceae.K. KubitzkiThe Families and Genera of Vascular Plants IV Flowering Plants Dicotyledons Malvales, Capparales and Non-betalain CaryophyllalesBerlinSpringer6270

20. Guzmán B, Vargas P (2005) Systematics, character evolution, and biogeography of Cistus L. (Cistaceae) based on ITS, trnL-trnF, and matK sequences. Mol Phylogenet Evol 37: 644–660.B. GuzmánP. Vargas2005Systematics, character evolution, and biogeography of Cistus L. (Cistaceae) based on ITS, trnL-trnF, and matK sequences.Mol Phylogenet Evol37644660

21. Guzmán B, Vargas P (2009) Long-distance colonization of the Western Mediterranean by Cistus ladanifer (Cistaceae) despite the absence of special dispersal mechanisms. J Biogeogr 36: 954–968.B. GuzmánP. Vargas2009Long-distance colonization of the Western Mediterranean by Cistus ladanifer (Cistaceae) despite the absence of special dispersal mechanisms.J Biogeogr36954968

22. Guzmán B, Vargas P (2009) Historical biogeography and character evolution of Cistaceae (Malvales) based on analysis of plastid rbcL and trnL-trnF sequences. Org Divers Evol 9: 83–99.B. GuzmánP. Vargas2009Historical biogeography and character evolution of Cistaceae (Malvales) based on analysis of plastid rbcL and trnL-trnF sequences.Org Divers Evol98399

23. Cunningham SA, Summerhayes BA, Westoby M (1999) Evolutionary divergences in leaf structure and chemistry, comparing rainfall and soil nutrient gradients. Ecol Monogr 69: 569–588.SA CunninghamBA SummerhayesM. Westoby1999Evolutionary divergences in leaf structure and chemistry, comparing rainfall and soil nutrient gradients.Ecol Monogr69569588

24. Givnish TJ (1984) Leaf and canopy adaptations in tropical forests. In: Medina E, Mooney HA, Vázquez-Yanes C, editors. Physiological ecology of plants of the wet tropics. The Hague: Junk. pp. 51–84.TJ Givnish1984Leaf and canopy adaptations in tropical forests.E. MedinaHA MooneyC. Vázquez-YanesPhysiological ecology of plants of the wet tropicsThe HagueJunk5184

25. Skelton PW (1993) Adaptive radiation: definition and diagnostic tests. In: Lees DR, Edwards D, editors. Evolutionary patterns and processes. London: Academic Press. pp. 45–58.PW Skelton1993Adaptive radiation: definition and diagnostic tests.DR LeesD. EdwardsEvolutionary patterns and processesLondonAcademic Press4558

26. Givnish TJ, Sytsma KJ, Smith JF, Hahn WJ (1995) Molecular evolution, adaptive radiation, and geographic speciation in Cyanea (Campanulaceae, Lebelioideae). In: Wagner WL, Funk VA, editors. Hawaiian biogeography. Washington: Smithsonian Institution Press. pp. 288–337.TJ GivnishKJ SytsmaJF SmithWJ Hahn1995Molecular evolution, adaptive radiation, and geographic speciation in Cyanea (Campanulaceae, Lebelioideae).WL WagnerVA FunkHawaiian biogeographyWashingtonSmithsonian Institution Press288337

27. Lowrey TK (1995) Phylogeny, adaptive radiation, and biogeography of Hawaiian Tetramolopium (Asteraceae, Astereae). In: Wagner WL, Funk VA, editors. Hawaiian biogeography. Washington: Smithsonian Institution Press. pp. 195–220.TK Lowrey1995Phylogeny, adaptive radiation, and biogeography of Hawaiian Tetramolopium (Asteraceae, Astereae).WL WagnerVA FunkHawaiian biogeographyWashingtonSmithsonian Institution Press195220

28. Robichaux RH, Carr GD, Liebman M, Pearcy RW (1990) Adaptive radiation of the Hawaiian silversword alliance (Compositae-Madiinae): ecological, morphological, and physiological diversity. Ann Missouri Bot Gard 77: 64–72.RH RobichauxGD CarrM. LiebmanRW Pearcy1990Adaptive radiation of the Hawaiian silversword alliance (Compositae-Madiinae): ecological, morphological, and physiological diversity.Ann Missouri Bot Gard776472

29. Kass R, Raftery A (1995) Bayes factors. J Am Stat Assoc 90: 773–795.R. KassA. Raftery1995Bayes factors.J Am Stat Assoc90773795

30. Whitfield JB, Lockhart PJ (2007) Deciphering ancient rapid radiations. Trends Ecol Evol 22: 258–265.JB WhitfieldPJ Lockhart2007Deciphering ancient rapid radiations.Trends Ecol Evol22258265

31. Rüber L, Van Tassell JL, Zardoya R (2003) Rapid speciation and ecological divergence in the American seven-spined Gobies (Gobiidae, Gobiosomatini) inferred from a molecular phylogeny. Evolution 57: 1584–1598.L. RüberJL Van TassellR. Zardoya2003Rapid speciation and ecological divergence in the American seven-spined Gobies (Gobiidae, Gobiosomatini) inferred from a molecular phylogeny.Evolution5715841598

32. Eriksson O, Bremer B (1992) Pollination systems, dispersal modes, life forms, and diversification rates in angiosmperm families. Evolution 46: 258–266.O. ErikssonB. Bremer1992Pollination systems, dispersal modes, life forms, and diversification rates in angiosmperm families.Evolution46258266

33. Bell CD, Donoghue MJ (2005) Phylogeny and biogeography of Valerianaceae (Dipsacales) with special reference to the South American valerians. Org Divers Evol 5: 147–159.CD BellMJ Donoghue2005Phylogeny and biogeography of Valerianaceae (Dipsacales) with special reference to the South American valerians.Org Divers Evol5147159

34. Hughes C, Eastwood R (2006) Island radiation on a continental scale: exceptional rates of plant diversification after uplift of the Andes. Proc Nat Acad Sci 103: 10334–10339.C. HughesR. Eastwood2006Island radiation on a continental scale: exceptional rates of plant diversification after uplift of the Andes.Proc Nat Acad Sci1031033410339

35. Avise JC (2009) Phylogeography: retrospect and prospect. J Biogeogr 36: 3–15.JC Avise2009Phylogeography: retrospect and prospect.J Biogeogr36315

36. Endler JA (1986) Natural selection in the wild. Princeton, New Jersey: Princeton University Press. JA Endler1986Natural selection in the wild.Princeton, New JerseyPrinceton University Press

37. Givnish TJ (1979) On the adaptive significance of leaf form. In: Solbrig OT, Jain S, Johnson GB, Raven PH, editors. Topics in plant population biology. New York: Columbia University Press. pp. 375–407.TJ Givnish1979On the adaptive significance of leaf form.OT SolbrigS. JainGB JohnsonPH RavenTopics in plant population biologyNew YorkColumbia University Press375407

38. Givnish TJ, Montgomery RA, Goldstein G (2004) Adaptive radiation of photosynthetic physiology in the Hawaiian Lobeliads: light regimes, static light responses, and whole-plant compensation points. Am J Bot 91: 228–246.TJ GivnishRA MontgomeryG. Goldstein2004Adaptive radiation of photosynthetic physiology in the Hawaiian Lobeliads: light regimes, static light responses, and whole-plant compensation points.Am J Bot91228246

39. Nuñez-Olivera E, Martinez-Abaigar J, Escudero JC (1996) Adaptability of leaves of Cistus ladanifer to widely varying environmental conditions. Funct Ecol 10: 636–646.E. Nuñez-OliveraJ. Martinez-AbaigarJC Escudero1996Adaptability of leaves of Cistus ladanifer to widely varying environmental conditions.Funct Ecol10636646

40. Gates DM, Alderfer R, Taylor SE (1968) Leaf temperature of desert plants. Science 159: 994–995.DM GatesR. AlderferSE Taylor1968Leaf temperature of desert plants.Science159994995

41. Szwarcbaum I (1982) Influence of leaf morphology and optical properties on leaf temperature and survival in three Mediterranean shrubs. Plant Sci Lett 26: 47–56.I. Szwarcbaum1982Influence of leaf morphology and optical properties on leaf temperature and survival in three Mediterranean shrubs.Plant Sci Lett264756

42. Parkhurst D, Loucks O (1972) Optimal leaf size in relation to environment. J Ecol 60: 505–537.D. ParkhurstO. Loucks1972Optimal leaf size in relation to environment.J Ecol60505537

43. Cunningham GL, Strain BR (1969) An ecological significance of seasonal leaf variability in a desert shrub. Ecology 50: 400–408.GL CunninghamBR Strain1969An ecological significance of seasonal leaf variability in a desert shrub.Ecology50400408

44. Givnish TJ, Vermeji GJ (1976) Sizes and shapes of liane leaves. Am Nat 110: 743–778.TJ GivnishGJ Vermeji1976Sizes and shapes of liane leaves.Am Nat110743778

45. Givnish TJ (1987) Comparative studies of leaf form: assessing the relative roles of selective pressures and phylogenetic constraints. New Phytol 106: (Suppl.)131–160.TJ Givnish1987Comparative studies of leaf form: assessing the relative roles of selective pressures and phylogenetic constraints.New Phytol106(Suppl.)131160

46. Ackerly DD, Knight CA, Weiss SB, Barton K, Starmer KP (2002) Leaf size, specific leaf area and microhabitat distribution of chaparral woody plants: contrasting patterns in species level and community level analyses. Oecologia 130: 449–457.DD AckerlyCA KnightSB WeissK. BartonKP Starmer2002Leaf size, specific leaf area and microhabitat distribution of chaparral woody plants: contrasting patterns in species level and community level analyses.Oecologia130449457

47. Savé R, Biel C, de Herralde F (2000) Leaf pubescence, water relations and chlrorophyll fluorescence in two species of Lotus creticus L. Biol Plant 43: 239–244.R. SavéC. BielF. de Herralde2000Leaf pubescence, water relations and chlrorophyll fluorescence in two species of Lotus creticus L.Biol Plant43239244

48. Ehleringer JR, Clark C (1988) Evolution and adaptation in Encelia (Asteraceae). In: Gottlieb LD, Subodh KJ, editors. Plant evolutionary biology. New York: Chapman and Hall Ltd. pp. 221–248.JR EhleringerC. Clark1988Evolution and adaptation in Encelia (Asteraceae).LD GottliebKJ SubodhPlant evolutionary biologyNew YorkChapman and Hall Ltd221248

49. Gülz PG, Herrmann T, Hangst K (1996) Leaf trichomes in the genus Cistus. Flora 191: 85–104.PG GülzT. HerrmannK. Hangst1996Leaf trichomes in the genus Cistus.Flora19185104

50. Suc JP, Bertini A, Combourieu-Nebout N, Diniz F, Leroy S, et al. (1995) Structure of West Mediterranean vegetation and climate since 5.3 Ma. Acta Zool Cracov 38: 3–16.JP SucA. BertiniN. Combourieu-NeboutF. DinizS. Leroy1995Structure of West Mediterranean vegetation and climate since 5.3 Ma.Acta Zool Cracov38316

51. White TJ, Bruns T, Lee S, Taylor J (1990) Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis M, Geldfand D, Sninsky J, White T, editors. PCR protocols: a guide to methods and applications. San Diego: Academic Press. pp. 315–322.TJ WhiteT. BrunsS. LeeJ. Taylor1990Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics.M. InnisD. GeldfandJ. SninskyT. WhitePCR protocols: a guide to methods and applicationsSan DiegoAcademic Press315322

52. Sun Y, Skinner DZ, Liang GH, Hulbert SH (1994) Phylogenetic analysis of Sorghum and related taxa using Internal Transcribed Spacer of nuclear ribosomal DNA. Theor Appl Genet 89: 26–32.Y. SunDZ SkinnerGH LiangSH Hulbert1994Phylogenetic analysis of Sorghum and related taxa using Internal Transcribed Spacer of nuclear ribosomal DNA.Theor Appl Genet892632

53. Taberlet P, Gielly L, GP , Bouvet J (1991) Universal primers for amplification of three non-coding regions of chloroplast DNA. Pl Mol Biol 17: 1105–1109.P. TaberletL. GiellyGPJ. Bouvet1991Universal primers for amplification of three non-coding regions of chloroplast DNA.Pl Mol Biol1711051109

54. Johnson LA, Soltis DE (1994) matK DNA and phylogenetic reconstruction in Saxifragaceae s. str. Syst Bot 19: 143–156.LA JohnsonDE Soltis1994matK DNA and phylogenetic reconstruction in Saxifragaceae s. str.Syst Bot19143156

55. Hamilton M (1999) Four primers pairs for the amplification of chloroplast intergenic regions with intraspecific variation. Mol Ecol 8: 521–523.M. Hamilton1999Four primers pairs for the amplification of chloroplast intergenic regions with intraspecific variation.Mol Ecol8521523

56. Savolainen V, Chase MW, Hoot SB, Morton CM, Soltis DE, et al. (2000) Phylogenetics of flowering plants based on combined analysis of plastid atpB and rbcL gene sequences. Syst Biol 49: 306–362.V. SavolainenMW ChaseSB HootCM MortonDE Soltis2000Phylogenetics of flowering plants based on combined analysis of plastid atpB and rbcL gene sequences.Syst Biol49306362

57. Emshwiller E, Doyle JJ (1999) Chloroplast-expressed Glutamine synthetase (ncpGS): Potencial utility for phylogenetic studies with an example from Oxalis (Oxalidaceae). Mol Phylogenet Evol 12: 310–319.E. EmshwillerJJ Doyle1999Chloroplast-expressed Glutamine synthetase (ncpGS): Potencial utility for phylogenetic studies with an example from Oxalis (Oxalidaceae).Mol Phylogenet Evol12310319

58. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882.JD ThompsonTJ GibsonF. PlewniakF. JeanmouginDG Higgins1997The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.Nucleic Acids Res2548764882

59. Swofford D (2002) PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland, MA: Sinauer. D. Swofford2002PAUP*.Phylogenetic analysis using parsimony (*and other methods). Version 4Sunderland, MASinauer

60. Mort ME, Soltis PS, Soltis DE, Mabry M (2000) Comparison of three methods for estimating internal support on phylogenetic trees. Syst Biol 49: 160–171.ME MortPS SoltisDE SoltisM. Mabry2000Comparison of three methods for estimating internal support on phylogenetic trees.Syst Biol49160171

61. Nylander JAA (2002) MrModeltest v1.0b. Uppsala: Department of Systematic Zoology. Uppsala University. JAA Nylander2002MrModeltest v1.0b.UppsalaDepartment of Systematic Zoology. Uppsala University

62. Posada D, Crandall KA (1998) Modeltest: testing the model of DNA substitution. Bioinformatics 14: 817–818.D. PosadaKA Crandall1998Modeltest: testing the model of DNA substitution.Bioinformatics14817818

63. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.F. RonquistJP Huelsenbeck2003MrBayes 3: Bayesian phylogenetic inference under mixed models.Bioinformatics1915721574

64. Magallón SA, Sanderson MJ (2005) Angiosperm divergence times: The effect of genes, codon positions, and time constraints. Evolution 59: 1653–1670.SA MagallónMJ Sanderson2005Angiosperm divergence times: The effect of genes, codon positions, and time constraints.Evolution5916531670

65. Sanderson MJ (2002) Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol 19: 101–109.MJ Sanderson2002Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach.Mol Biol Evol19101109

66. Naud G, Suc JP (1975) Contribution à l'étude paléofloristique des Coirons (Ardèche). Première analyses polliniques dans les alluvions sous-basaltique et interbasaltiques de Mirabel (Miocène supérieur). Bull Soc Géol France 17: 820–827.G. NaudJP Suc1975Contribution à l'étude paléofloristique des Coirons (Ardèche). Première analyses polliniques dans les alluvions sous-basaltique et interbasaltiques de Mirabel (Miocène supérieur).Bull Soc Géol France17820827

67. Menke B (1976) Pliozäne und ältestquartäre Sporen- und pollenflora von Schleswig-Holstein. Geologisches Jahrbuch Reihe A 32: 3–197.B. Menke1976Pliozäne und ältestquartäre Sporen- und pollenflora von Schleswig-Holstein.Geologisches Jahrbuch Reihe A323197

68. Kendall DG (1949) Stochastic processes and population growth. J Roy Stat Soc B 11: 230–264.DG Kendall1949Stochastic processes and population growth.J Roy Stat Soc B11230264

69. Moran PA (2000) Estimation methods for evolutive processes. J Roy Stat Soc B 13: 141–146.PA Moran2000Estimation methods for evolutive processes.J Roy Stat Soc B13141146

70. Maddison WP, Maddison DR (1992) MacClade: Analysis of Phylogeny and Character Evolution, version 3.01. Sunderland, MA: Sinauer Associates. WP MaddisonDR Maddison1992MacClade: Analysis of Phylogeny and Character Evolution, version 3.01.Sunderland, MASinauer Associates

71. Pagel M, Meade A (2007) BayesTraits.Version 1.0M. PagelA. Meade2007BayesTraits.Version 1.0

72. Guzmán B (2008) Variabilidad morfológica y nucleotídica en el género Cistus: análisis macro- y microevolutivos. Madrid: Universidad Complutense. PhD Thesis. B. Guzmán2008Variabilidad morfológica y nucleotídica en el género Cistus: análisis macro- y microevolutivos.MadridUniversidad Complutense. PhD Thesis

73. Clement M, Posada D, Crandall KA (2000) TCS: a computer program to estimate gene genealogies. Mol Ecol 9: 1657–1659.M. ClementD. PosadaKA Crandall2000TCS: a computer program to estimate gene genealogies.Mol Ecol916571659

74. Templeton AR, Crandall KA, Sing CF (1992) A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics 132: 619–633.AR TempletonKA CrandallCF Sing1992A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation.Genetics132619633

75. Grosser W (1903) Cistaceae. In: Engler A, editor. Das Pflanzenreich. Berlin: Breitkopf & Härtel. W. Grosser1903Cistaceae.A. EnglerDas PflanzenreichBerlinBreitkopf & Härtel161

76. Martín Bolaños M, Guinea E (1949) M. Martín BolañosE. Guinea1949Jarales y jaras (cistografia hispánica). Madrid: Ministerio de Agricultura. Direccion General de Montes, Caza y Pesca Fluvial. Instituto Forestal de Investigaciones y Experiencias. Jarales y jaras (cistografia hispánica). Madrid: Ministerio de Agricultura. Direccion General de Montes, Caza y Pesca Fluvial. Instituto Forestal de Investigaciones y Experiencias.

77. Dansereau P (1958) Notes sur les Cistes. III: Les conditions de la distribution géographique. C R Seances Soc Biogeogr 304: 22–25.P. Dansereau1958Notes sur les Cistes. III: Les conditions de la distribution géographique.C R Seances Soc Biogeogr3042225

78. Warburg EF (1968) Cistus L. In: Tutin TG, Heywood VH, Burges NA, Moore DM, Valentine DH, editors. Flora Europaea. Cambridge: Cambridge University Press. pp. 282–285.EF Warburg1968Cistus L.TG TutinVH HeywoodNA BurgesDM MooreDH ValentineFlora EuropaeaCambridgeCambridge University Press282285

79. Demoly JP, Montserrat P (1993) Cistus. In: Castroviejo S, Aedo C, Cirujano S, Laínz M, Montserrat P, et al., editors. Flora iberica. Madrid: Consejo de Investigaciones Científicas. pp. 319–337.JP DemolyP. Montserrat1993Cistus.S. CastroviejoC. AedoS. CirujanoM. LaínzP. MontserratFlora ibericaMadridConsejo de Investigaciones Científicas319337

80. Greuter W (1996) Proposal to conserve the name Cistus albanicus (Cistaceae). Taxon 45: 715–716.W. Greuter1996Proposal to conserve the name Cistus albanicus (Cistaceae).Taxon45715716