- Research article
- Open Access
Accelerated evolutionary rates in tropical and oceanic parmelioid lichens (Ascomycota)
BMC Evolutionary Biologyvolume 8, Article number: 257 (2008)
The rate of nucleotide substitutions is not constant across the Tree of Life, and departures from a molecular clock have been commonly reported. Within parmelioid lichens, the largest group of macrolichens, large discrepancies in branch lengths between clades were found in previous studies. Using an extended taxon sampling, we test for presence of significant rate discrepancies within and between these clades and test our a priori hypothesis that such rate discrepancies may be explained by shifts in moisture regime or other environmental conditions.
In this paper, the first statistical evidence for accelerated evolutionary rate in lichenized ascomycetes is presented. Our results give clear evidence for a faster rate of evolution in two Hypotrachyna clades that includes species occurring in tropical and oceanic habitats in comparison with clades consisting of species occurring in semi-arid and temperate habitats. Further we explore potential links between evolutionary rates and shifts in habitat by comparing alternative Ornstein-Uhlenbeck models.
Although there was only weak support for a shift at the base of a second tropical clade, where the observed nucleotide substitution rate is high, overall support for a shift in environmental conditions at cladogenesis is very strong. This suggests that speciation in some lichen clades has proceeded by dispersal into a novel environment, followed by radiation within that environment. We found moderate support for a shift in moisture regime at the base of one tropical clade and a clade occurring in semi-arid regions and a shift in minimum temperature at the base of a boreal-temperate clade.
Differences in nucleotide substitution rates among taxa are a common phenomenon in molecular studies [1–7], and the presence of an exact molecular clock [8, 9] appears to be a rare exception in molecular evolution, if present at all. A number of different causes are invoked to explain differences in evolutionary rates among taxa. Recently, Kay et al.  demonstrated correlations of substitution rates and life histories. In angiosperms, herbaceous plants have substitution rates almost twice as high as woody plants , while no phylogenetic constraint on rates was found in their study among lineages. Previously, higher substitution rates had been found in annual compared to perennial plants [5, 11–14]. These studies in plants agree with previous studies suggesting that shorter generation times are associated with accelerated evolutionary rates [15, 16].
Significant departure from constant rates of a molecular clock have so far only been demonstrated for the kingdom Fungi in Basidiomycota [17, 18] and in four endosymbiotic pyrenomycetes . In these cases accelerated evolutionary rates were significantly associated with mutualism. In one of the largest classes of Ascomycota, Lecanoromycetes [19, 20], significant differences in nucleotide substitution rates have not yet been shown. However, in previous studies [21, 22] we have found remarkable differences in branch lengths between clades in parmelioid lichens (Parmeliaceae, Lecanoromycetes). Taxa of one clade (the Hypotrachyna clade) had consistently longer branches in phylogenetic trees than taxa in other clades, which suggested an accelerated evolutionary rate in that group.
Parmelioid lichens are the largest group within Parmeliaceae , which itself is one of the largest families of lichen-forming fungi and has a worldwide distribution. The Parmelioid group comprises common and well-known species, including taxa such as Parmelia sulcata, Flavoparmelia caperata, and Punctelia borreri, which are frequently used in biomonitoring atmospheric pollution [24, 25]. The group has approximately 1500 taxa [26, 27]. It includes species which are mainly foliose, mostly rhizinate lichens with laminal apothecia and simple, hyaline ascospores . Recent phylogenetic studies have shown that a core group of parmelioid lichens is monophyletic [21–23]. Within this group, seven well-supported clades were found , which show distinct ecological preferences. For example, taxa belonging to the Hypotrachyna clade (c. 300 species) are mostly tropical extending into temperate regions in moist, oceanic habitats [29–31]. In contrast, species of the Xanthoparmelia clade (c. 700 species) [32, 33] have their center of distribution in semi-arid regions.
Lichens are poikilohydric organisms and are able to survive long periods of dry conditions in a dormant stage [34–40]. They lack stomata and any water storage system. Water vapor is lost readily from across the whole surface of lichens. In some lichens, alteration between desiccated and hydrated stages is an almost daily occurrence [41–44]. This has an important ecological consequence, since the duration of physiological activity of lichens is directly dependent on the ability of water in a habitat. Lichens in dry habitats may be physiologically active only for a few hours a day and dormant most of the day [45–49], while lichens growing in moist habitats usually have sufficient water content all day long for physiological activity [50–52].
To further investigate the remarkable differences in branch lengths between species of the Hypotrachyna clade (including tropical species extending into temperate regions in moist, oceanic habitats) and other clades of parmelioid lichens  we obtained sequence data from additional species and/or genes to perform a study focusing on the differences in evolutionary rates within parmelioid lichens and potential causes for rate shifts. Sequences of two nuclear (ITS, nu LSU) and single mitochondrial ribosomal (mtSSU) DNA loci are used to infer phylogenetic relationships. These genes have been successfully used to infer phylogenies in parmelioid lichens [21, 29]. Here we use a phylogenetically based maximum likelihood procedure to identify shifts in rates of nucleotide substitutions. Further we explore the potential links between rates of molecular evolution and ecological parameters by comparing alternative Ornstein-Uhlenbeck models . These models allow us to test the hypothesis that cladogenesis in parmelioid lichens is associated with shifts in habitat that might in turn explain changes in nucleotide substitution rate.
Phylogenetic analyses of the three gene partitions combined included 2022 unambiguously aligned nucleotide, of which 734 were variable and 563 were parsimony informative. The topologies of the Bayesian analysis (Figure 1) and MP results were generally similar and all nodes strongly supported by MP bootstrapping (≥75%) were well supported by Bayesian pp (MP tree not shown). All strongly supported clades are also supported in the relaxed-clock analyses conducted in BEAST (data not shown). Our results were largely consistent with previous phylogenies [21–23, 29, 54], but with the addition of several new taxa and/or an additional gene partition. The species were placed within groups expected by their morphological and chemical characters [21, 28, 30, 31]. However, the placement of some taxa differed slightly from previous analysis. This includes Parmelinella wallichiana, which had an unresolved sister-group relationship to Bulbothrix decurtata and a clade including B. meizospora and B. setschwanensis; Everniastrum, which was resolved as monophyletic in previous analyses; and a few relationships of terminal sister taxa . None of the differences to previous studies involve well supported conflicts. Further, the close relationship of Karoowia saxeti with Xanthoparmelia found in previous studies [21, 32] is confirmed and its placement within that genus is strongly supported for the first time.
Significant Rate Differences
Between major clades of parmelioid lichens considerable branch length differences were evident, but also within clades branch lengths differed (Figure 1). Likelihood ratio tests comparing null models of rate constancy versus alternative models revealed statistically significant departures from rate constancy across the entire tree (Additional file 1). Since some of the well supported clades included only a few OTUs, we restricted our analysis to five clades with sufficient number of taxa (more than six OTUs). Within these five clades selected for the rate difference test, analyses of three of 15 single-gene data sets did not return significant results, meaning that a molecular clock cannot be rejected. This includes the ITS data sets of the Melanohalea- and Xanthoparmelia clades, and the mtSSU data set of the Parmotrema clade.
To determine rate differences between clades and to identify particular clades with deviating rates, two-rate and three-rate models were compared for three selected clades in each analysis using the likelihood ratio test statistics. Significant rate differences were found between the Hypotrachyna clades 1 and 2 and other clades examined (Additional file 1). Rate comparisons of the Melanohalea clade with the Parmotrema and Xanthoparmelia clades, respectively, and comparisons of the latter two clades revealed insignificant results.
The species of the five selected clades live under different ecological conditions (Table 1). On average, the species in the two Hypotrachyna clades occur at higher altitudes. Moreover, these and the species of the Parmotrema clade occupy habitats with higher precipitation and consequently have the highest Emberger Index values. In contrast, species in the Xanthoparmelia and Melanohalea clades occur in habitats with low precipitation and low Emberger's index.
There is moderate support for a shift in moisture regime (represented by precipitation and Emberger's Index) at the base of Hypotrachyna clade 2 and the Xanthoparmelia clade, and strong support for a shift in minimum temperature at the base of the Melanohalea clade. In the local test of alternative Ornstein-Uhlenbeck models (O-U models; the models and our 'local' and 'global' tests are described in the methods section), a shift in precipitation at the base of Clade 2 is supported with Bayes information criterion (BIC) weight = 0.876, and a shift in Emberger Index is supported at BIC weight = 0.842 (Table 2). A shift in Emberger's Index at the base of the Xanthoparmelia clade is similarly supported (BIC weight = 0.895). A shift in minimum temperature is very strongly supported at the base of the Melanohalea clade under both the local and global tests and whether measurement error is accounted for or not (BIC weight > 0.98; Tables 2 &3). These findings argue against a gradual shift in moisture regime in which the covariance among taxa is predicted by time since divergence from their most recent common ancestor. Rather, they support a scenario in which transitions in moisture regime are associated with cladogenetic events. In both the global and local tests of maximum temperature and altitude, the Brownian motion or single-equilibrium O-U model cannot be rejected (Tables 2 &3).
Incorporating estimated measurement error decreases the relative support for the two-equilibrium O-U models in the local test (e.g., BIC weight = 0.567 and 0.276 for precipitation and Emberger's Index respectively for Hypotrachyna clade 2 when measurement error is accounted for; Table 2). In the global test, incorporating measurement error into analysis has little effect on the evidence for a change in moisture regime at the base of clade 2, which is supported at BIC weight = 0.897 to 0.905 either with or without measurement error (Table 3). The model-averaged estimate of α for Emberger's Index is 16.447 when measurement error is incorporated, 19.015 when it is not (Additional files 2, 3). This has the effect of almost erasing the phylogenetic effects in the data. When measurement error is disregarded (treated as zero), there is little support for a Brownian motion model for any characters except minimum temperature (Table 2). In the global analyses of precipitation and Emberger's Index, the 99% confidence interval excludes both the Brownian motion model and the single-moisture-regime O-U model (Table 3, Additional files 2, 3). When measurement error is ignored, there is concordance between the local and global tests.
Our extended sampling revealed phylogenetic estimates congruent with previous analyses [21, 29] including several well-supported clades within parmelioid lichens. However, as in these studies, we failed to resolve with confidence the phylogenetic relationships between these strongly supported main clades. Further, additional sampling of the Hypotrachyna clade supported that it includes two major clades. In contrast to Blanco et al. , but in agreement with Divakar et al. , the sister-group relationship between the two Hypotrachyna clades is not strongly supported, and in the Bayesian analysis using Beast, the two clades do not form a monophyletic group. The lack of confidence in the phylogenetic relationships between the clades is consistent with the scenario of an adaptive radiation at the base of parmelioid clade. However, we cannot exclude the possibility of a lack of power of resolution in the gene partitions used in the analyses.
This study provides clear evidence for a general acceleration in rates of molecular evolution in the two Hypotrachyna clades in comparison with other clades in parmelioid lichens. Although this study presents the first statistical evidence for a significantly accelerated substitution rate in the Hypotrachyna clades, a prior study showed its occurrence without discussing this further . The increased rates of nucleotide substitution seen could result from different factors, such as positive selection on coding regions of the ribosomal DNA, or represent a general acceleration of nucleotide substitutions . Our analyses point to the latter as the most plausible cause.
Positive selection is the process favouring the retention of mutations in a population that are beneficial to the reproductive success of individuals . Positive selection leading to increased fixation of nonsynonymous substitutions has been demonstrated for some parasites [56, 57]. To examine the possibility of this process on coding regions of the ribosomal DNA, we examined the rate increase in a conserved part of the ITS regions (199 bp long) and compared it to the rate increase in the 5.8 S rDNA and nu LSU rDNA. If faster substitution rates were caused by differential selection, there should be no association between faster evolutionary rates and affiliation to the Hypotrachyna clades. However, we found higher rates in the Hypotrachyna clades in all three- gene partitions studied (data not shown). There were no significant differences between the rate increases between the gene partitions, indicating that positive selection on coding regions is not present.
Another possible cause of rate acceleration is an increase in mutation rate . Our results indicate a general acceleration of nucleotide substitutions in the Hypotrachyna clade (see branch lengths, Figure 1). Causes for such an acceleration could be numerous, including shorter generation time , metabolic rate [60, 61], lack of sexuality , or demographic factors [15, 62].
A population bottleneck may result from a significant reduction of population size due to extinction of many genetic lineages. Since mutations have a higher probability of fixation in populations experiencing genetic bottlenecks, this may result in a long branch leading to a clade [6, 63]. In fact, branches leading to Hypotrachyna clade 1 and 2 are long. However, this interpretation fails to account for the general rate acceleration in the Hypotrachyna clades. Populations of species in these two clades are not smaller and species do not have more restricted distributions than in other clades (PKD, unpubl. observations).
The loss of sexual reproduction has been identified as a cause of accelerated evolution in endosymbiotic bacteria [62, 64]. Since the ratio of asexual/sexual species in the Hypotrachyna clades does not differ from other clades (0.29, 0.55 in Hypotrachyna clades and 0.27, 0.5, 1.25 in other clades), this explanation is also not likely.
Two intrinsic biological variables, generation time and metabolic rate, have been interpreted as causes of rate variations in other organisms [60, 61, 65–68]. Since rates of DNA damage are proportional to specific metabolic rate, species with higher metabolic rates should have higher substitution rates given there is a relationship between DNA damage and mutation rate . Unfortunately, there are no ecophysiological measurements available for the lichen clades studied. Hence, we have no data on metabolic rates in these fungi. However, as mentioned above lichenized fungi are poikilohydric organisms and from studies in other groups of lichens, we can assume that species occurring in arid and semi-arid regions are dormant for large parts of their life [45, 49, 69]. For example, Teloschistes capensis (Teloschistaceae) occurs in the same habitats along the Namibian and north-western coast of South Africa as some taxa of the Xanthoparmelia clade studied here. This species has been shown to be physiologically active for only very few hours in the mornings after dew fall . In contrast, species in moist habitats, such as rainforests were shown to be physiologically active throughout the day [51, 52, 70]. Growth measurements in lichens from arid and semi-arid regions revealed that these lichens grow very slowly (0.37 mm/y in Caloplaca aurantia, Teloschistaceae)  and hence reach an age of several hundred years, while some lichens in moist habitats, such as coastal forests of California can increase their length over 35% per year  and have a much shorter generation time. The slow rates of mutations in Xanthoparmelia clade may be related with a strategy adapted to a low metabolic activity due to long dormancy periods in dry conditions. Again, we have no data from the parmelioid lichens studied here and point out that ecophysiological studies need to be extended to these lichens to confirm their general importance.
Given our knowledge of the differences in growth rate and generation time for these poikilohydric lichens in moist vs. arid regions, we hypothesized at the outset of this study that a shift in moisture regime might explain the increased nucleotide substitution rate in Hypotrachyna clades 1 and 2. This might be due either to decreased generation time or increased mutation rate in environments with generally higher precipitation than in the other three clades. The support for O-U models that specify shifts in moisture regime (Emberger's Index and precipitation) at the base of Hypotrachyna clade 2 and minimum temperature at the base of clade 5 (Melanohalea), combined with weak support for the Brownian motion model and single-environment Ornstein-Uhlenbeck model for these variables, suggests that these clades originated by colonization of a novel environment, followed by diversification within those environments. The extreme phylogeny-effacing effects of the high inferred rate of evolution toward equilibrium Emberger index values particularly suggests that shifts in clade-specific moisture regimes govern the evolution of habitat preferences in these lichens. There is support for a change in moisture regime at the base of Hypotrachyna clade 1 (BIC weight = 0.760 without measurement error, 0.803 with error; Table 3), which is comparable to the support for a change at the base of Hypotrachyna clade 2. These two clades have an accelerated substitution rate. However, the support for a change is also approximately the same for the Parmotrema clade (BIC weight = 0.763 without error, 0.806 with error), which does not exhibit elevated substitution rate.
The apparently weak link between moisture regime and substitution rate in this study may be due to a lack of causal relationship. However, it might also be due to an alternative adaptative strategy of the Parmotrema clade. Species in this clade are especially abundant in insular regions  with oceanic climate. Although no ecophysiological data on parmelioid lichens are available, other lichens occurring in similar moisture conditions, were shown to have metabolic (and growth) rates depending on a combination of light intensity and temperature . Most species in the Parmotrema clade may have slow metabolic rates due to ecological conditions related with latitude. In fact, Parmotrema clade species are subtropical or temperate [28, 74], living in regions where, comparing to tropics, light intensity is lower. Wright et al.  showed that evolutionary rates of tropical species are higher than those of subtropical taxa. However, low sample size (inadequate number of taxa to accurately estimate the posterior probability support for environmental shifts associated with shifts in substitution rate) or inadequate characterization of the moisture environment cannot be discarded. Additionally, in our observation (T.L. and A.C.), the moisture available to the lichens in this study is strongly dominated by fog, which is not directly accounted for by either precipitation or Emberger's Index (though it is likely to be correlated). In species occurring in areas with extensive fogs, lichens were shown to be able to use this as the main source of water taken up from water vapour [36, 76]. Future work in the ecology of this group of organisms should focus on characterizing moisture regime more accurately.
Our extended taxon sampling confirmed phylogenetic relationships revealed in previous studies. Here we show significant accelerated evolutionary rates in clades of parmelioid lichens that occur in tropical and oceanic habitats as opposed to those in arid habitats. Comparison of alternative Ornstein-Uhlenbeck models gave moderate support for a shift in moisture regime at the base of one tropical clade and a clade occurring in semi-arid regions and a shift in minimum temperature at the base of a boreal-temperate clade. This finding that cladogenesis may be associated with shifts in environment suggests adaptive radiation as a mechanism of speciation in Parmelioid lichens, a hypothesis that bears testing using finer-grained ecological data.
Data matrices of 128 parmelioid lichens and three outgroup taxa were assembled using sequences of nuclear LSU, ITS, and mitochondrial SSU rDNA sequences. Specimens and sequences used for the molecular analyses are compiled in Additional file 4. The data set includes 11 sequences downloaded from GenBank, 343 from previous publications by us [21, 23, 29, 32, 74, 77], and 39 sequences newly generated for this study. We used three species of the genus Usnea as the out-group since this genus has been shown to be closely related to parmelioid lichens in a previous study .
Small samples prepared from freshly collected and frozen specimens were ground with sterile plastic pestles. Total genomic DNA was extracted using the DNeasy Plant Mini Kit (Qiagen, Hilden) according to the manufacturer's instructions but with slight modifications as described elsewhere . Fungal nu LSU rDNA was amplified using the primers nu-LSU-0155-5' , LR0R (Vilgalys unpublished, http://www.botany.duke.edu/fungi/mycolab), LR5, LR6 and LR7  and mt SSU rDNA using the primers mrSSU1, mrSSU3R , MSU1 and MSU7 .
Amplifications were performed in 50 μl volumes containing 5–25 ng of DNA and a reaction mixture of 5 μl 10× DNA polymerase buffer (Biotools) (containing MgCl2 2 mM, 10 mM Tris-HCl, pH 8.0, 50 mM KCl, 1 mM EDTA, 0.1% Triton X-100), 1 μl of deoxinucleotide triphosphate (dNTPs), containing 10 mM of each base, 2.5 μl of each primer (10 μM) and 1.25 μl of DNA polymerase (1 U/μl). The amplifications were carried out in an automatic thermocycler (Techne Progene) and performed using the same programs as described in . The PCR products were then cleaned using the Bioclean Columns kit (Biotools, Madrid) according to the manufacturer's instructions. The cleaned PCR products were sequenced using the same primers as used for PCR amplification. The ABI Prism™ Dye Terminator Cycle Sequencing Ready reaction kit (Applied Biosystems, Foster City) was used and the following settings were applied: denaturation for 3 min at 94°C, 25 cycles at 96°C for 10 sec, 50°C for 5 sec, and 60°C for 4 min. Sequencing reactions were electrophoresed on a 3730 DNA analyzer (Applied Biosystems, Foster City). Sequence fragments obtained were assembled with SeqMan 4.03 (DNAStar, Madison) and manually adjusted.
The mitochondrial and the ITS data sets contain sequence portions that are highly variable. Standard multiple alignment programs, such as Clustal  become less reliable when sequences show a high degree of divergence. Therefore we employed an alignment procedure that uses a linear Hidden Markov Model (HMM) as implemented in the software SAM (Sequence Alignment and Modelling system)  for separate alignments of the three data sets. Regions that were not aligned with statistical confidence were excluded from the phylogenetic analysis.
The alignments were analysed by maximum parsimony (MP) and a Bayesian approach (B/MCMC) [85, 86]. To test for potential conflict, parsimony bootstrap analyses were performed on each individual data set, and 75% bootstrap consensus trees were examined for conflict .
Maximum parsimony analyses were performed using the program PAUP* . Heuristic searches with 1000 random taxon addition replicates were conducted with TBR branch swapping and MulTrees option in effect, equally weighted characters and gaps treated as missing data. Bootstrapping  was performed based on 2000 replicates with random sequence additions.
The B/MCMC analyses were conducted using the MrBayes 3.1.1 program . The analyses were performed assuming the general time reversible model of nucleotide substitution , including estimation of invariant sites and assuming a discrete gamma distribution with six rate categories (GTR+I+G). The data set was portioned into the three parts (ITS, nu LSU, mt SSU) and each partition was allowed to have its own parameters . No molecular clock was assumed. A run with 3,000,000 generations starting with a random tree and employing 12 simultaneous chains was executed. Every 100th tree was saved into a file. The first 200,000 generations (i.e. the first 2000 trees) were deleted as the "burn in" of the chain. We plotted the log-likelihood scores of sample points against generation time using TRACER 1.0 http://evolve.zoo.ox.ac.uk/software.html?id=tracer to ensure that stationarity was achieved after the first 200,000 generations by checking whether the log-likelihood values of the sample points reached a stable equilibrium value . Of the remaining 56,000 trees (28,000 from each of the parallel runs) a majority rule consensus tree with average branch lengths was calculated using the sumt option of MrBayes. Posterior probabilities were obtained for each clade. Clades with bootstrap support equal or above 75% under MP and posterior probabilities ≥ 0.95 were considered as strongly supported. Phylogenetic trees were visualized using the program Treeview .
Sequence Significance Test for the Presence of Rate Differences
The program Baseml (part of the PAML 3.14 b package;  was used for significance tests for the presence of rate differences, assuming a GTR+G model. The resulting likelihoods were compared using likelihood ratio tests . Two types of analyses were performed and the results of those are listed in Table 1. In a first set of analyses, two models were compared to determine whether observed branch length differences were the result of a significant departure from rate constancy. The null model assumed a molecular clock, while in the alternative model each branch was allowed its own unique rate of molecular evolution. These two models were contrasted for each gene partition separately (ITS, mtSSU, nuLSU) across the whole ingroup of parmelioid lichens, as well as across five selected, well-supported major clades within parmelioid lichens, which are Hypotrachyna 1, Hypotrachyna 2, Melanohalea, Parmotrema and Xanthoparmelia [21, 29] to determine rate differences within these clades. In a second analysis, several tests were performed to determine whether significant departure from rate homogeneity were present between the five selected clades. For these analyses, a pruned tree comprising of the five selected clades of parmelioid lichens was used for comparisons of a null and alternative model. The five clades selected included the two clades (Hypotrachyna clade 1 and 2) with evidently long branches and the Xanthoparmelia clade that has comparably short branches. Comparisons were made between two-rate models in which two selected clades had the same rate but a third clade had a different rate versus a three-rate model in which all three included clades had different rates. These analyses were made on the concatenated data set, and in total 30 analyses were performed.
Phylogenetic Comparative Analyses
Estimates of five environmental variables – Emberger's Index , altitude, precipitation, minimum temperature and maximum temperature – were obtained for one to three specimens of each of 97 species from three online databanks (http://www.ucm.es/info/cif/data/indexcsp.htm, http://climexp.knmi.nl/start.cgi?someone@somewhere, and http://www.worldclimate.com/) and are available in the additional file 4 and 5. The Emberger's index including a combination of simple climate parameters such as the mean of minimal and maximal annual temperature and precipitation was calculated by the formula Q = 100 × P/Tmmax2 - Tmmin2, where P is average annual precipitation, Tmmax is average maximum temperature of warmest month, Tmmin is average minimum temperature of coldest month.
Environmental variables were analyzed in a phylogenetic generalized least squares framework under alternative Ornstein-Uhlenbeck (O-U) models [53, 96]. The O-U process is essentially a Brownian motion model with a pull toward an equilibrium that may be environment-specific or clade specific. The O-U process has been referred to as a "rubber-band" process , because it models random, stochastic processes as well as deterministic evolutionary forces that pull character values toward a central optimum or equilibrium value. Under this model, the parameter α estimates the rate at which the variable being analyzed reaches an equilibrium, and the parameter θ i estimates the equilibrium value of that variable along branches in group i. Analyses were conducted allowing θ to change at the base of each of five strongly supported clades (Figure 1). This models the situation in which environmental conditions shift at the base of a clade, such that the radiation of that clade takes place within a potentially novel set of environmental conditions. Bayes information criterion (BIC) weights were used to assess the support for the hypothesis that each clade radiated under a novel set of conditions relative to each of the five environmental variables. BIC weights are interpreted as Bayesian posterior model probabilities assuming equal prior model probabilities . Strong relative support for a single-equilibrium O-U model or a Brownian motion model (constant variance with respect to phylogeny) indicates that the data do not support a shift in environmental conditions attending cladogenetic events. Emberger's Index and precipitation were log-transformed and altitude was square-root transformed prior to analysis. Transformations had little effect on the conclusions of this paper but improved model fit in trial analyses by 400 to 600 log-likelihood units.
We evaluated the support for alternative hypotheses both locally and globally. In local tests, we compared the fit of three models for each dataset: a two-equilibrium O-U model that allows a change in environment at the base of a single clade; the single-equilibrium O-U model, in which the entire tree comprises a single environmental regime; and the Brownian motion model. In global tests, we compared the fit of the Brownian motion model plus all 5^2 = 32 possible O-U models, allowing changes at all permutations of the five focus clades, and summed BIC weights over these clades to estimate the global support for a change in environment attending the origins of each clade. This approach is analogous to model-averaging and has proven useful as a means of estimating the posterior probability of change occurring along a given phylogenetic branch . Effects of the O-U α parameter on phylogenetic signal were visualized using the Geiger package  in R 2.6.0 (R Development Core Team 2007).
Analyses were conducted on an ultrametric tree estimated using the lognormal relaxed clock model  implemented in an MCMC framework in BEAST v1.4.5 , using the Bayesian consensus as the start tree and a GTR+I+G model of nucleotide substitution, with a total run of 20,000,000 generations. Measurement error was incorporated into analyses by adding an estimate of the squared standard error for each environmental variable to the diagonal of the variance-covariance matrix utilized in generalized least squares computations . Because the small sample sizes (N = 1 to 3 samples per species) increase uncertainty in the estimates of measurement error, the squared standard error was estimated for each environmental variable as the mean variance over all species represented by N = 3 samples, divided by the sample size for each species separately. With such small samples, measurement error in this study represented a substantial component of the total variance in the tree, and analyses incorporating measurement error were consequently compared with analyses in which measurement error was set at zero.
Arbogast BS, Edwards SV, Wakeley J, Beerli P, Slowinski JB: Estimating Divergence Times from Molecular Data on Phylogenetic and Population Genetic Timescales. Annual Review of Ecology and Systematics. 2002, 33: 707-740. 10.1146/annurev.ecolsys.33.010802.150500.
Britten RJ: Rates of DNA-Sequence Evolution Differ between Taxonomic Groups. Science. 1986, 231 (4744): 1393-1398. 10.1126/science.3082006.
Bromham L, Penny D: The Modern Molecular Clock. Nature Reviews Genetics. 2003, 4 (3): 216-224. 10.1038/nrg1020.
Bromham L: Molecular Clocks in Reptiles: Life History Influences Rate of Molecular Evolution. Molecular Biology and Evolution. 2002, 19 (3): 302-309.
Gaut BS, Muse SV, Clark WD, Clegg MT: Relative Rates of Nucleotide Substitution at the Rbcl Locus of Monocotyledonous Plants. Journal of Molecular Evolution. 1992, 35 (4): 292-303. 10.1007/BF00161167.
Woolfit M, Bromham L: Increased Rates of Sequence Evolution in Endosymbiotic Bacteria and Fungi with Small Effective Population Sizes. Molecular Biology and Evolution. 2003, 20 (9): 1545-1555. 10.1093/molbev/msg167.
Langley CH, Fitch WM: An Examination of the Constancy of the Rate of Molecular Evolution. Journal of Molecular Evolution. 1974, 3: 161-167. 10.1007/BF01797451.
Margoliash E: Primary Structure and Evolution of Cytochrome C. Proceedings of the National Academy of Sciences of the United States of America. 1963, 50: 672-679. 10.1073/pnas.50.4.672.
Zuckerkandl E, Pauling L: Evolutionary Divergence and Convergence in Proteins. Evolving Genes and Proteins. Edited by: Bryson V, Vogel H. 1965, New York: Academic Press, 97-166.
Kay KM, Whittall JB, Hodges SA: A Survey of Nuclear Ribosomal Internal Transcribed Spacer Substitution Rates across Angiosperms: An Approximate Molecular Clock with Life History Effects. Bmc Evolutionary Biology. 2006, 6: 9-10.1186/1471-2148-6-36.
Bousquet J, Strauss SH, Doerksen AH, Price RA: Extensive Variation in Evolutionary Rate of Rbcl Gene-Sequences among Seed Plants. Proceedings of the National Academy of Sciences of the United States of America. 1992, 89 (16): 7844-7848. 10.1073/pnas.89.16.7844.
Brunsfeld SJ, Soltis PS, Soltis DE, Gadek PA, Quinn CJ, Strenge DD, Ranker TA: Phylogenetic-Relationships among the Genera of Taxodiaceae and Cupressaceae – Evidence from Rbcl Sequences. Systematic Botany. 1994, 19 (2): 253-262. 10.2307/2419600.
Conti E, Fischbach A, Sytsma KJ: Tribal Relationships in Onagraceae – Implications from Rbcl Sequence Data. Ann Mo Bot Gard. 1993, 80 (3): 672-685. 10.2307/2399853.
Laroche J, Li P, Maggia L, Bousquet J: Molecular Evolution of Angiosperm Mitochondrial Introns and Exons. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94 (11): 5722-5727. 10.1073/pnas.94.11.5722.
Ohta T: The Nearly Neutral Theory of Molecular Evolution. Annual Review of Ecology and Systematics. 1992, 23: 263-286. 10.1146/annurev.es.23.110192.001403.
Gu X, Li W-H: Higher Amino Acid Substitution in Rodents Than in Humans. Molecular Phylogenetics and Evolution. 1992, 1: 211-214. 10.1016/1055-7903(92)90017-B.
Zoller S, Lutzoni F: Slow Algae, Fast Fungi: Exceptionally High Nucleotide Substitution Rate Differences between Lichenized Fungi Omphalina and Their Symbiotic Green Algae Coccomyxa. Molecular Phylogenetics and Evolution. 2003, 29 (3): 629-640. 10.1016/S1055-7903(03)00215-X.
Lutzoni F, Pagel M: Accelerated Evolution as a Consequence of Transitions to Mutualism. Proceedings of the National Academy of Sciences of the United States of America. 1997, 94 (21): 11422-11427. 10.1073/pnas.94.21.11422.
Hibbett DS, Binder M, Bischoff JF, Blackwell M, Cannon PF, Eriksson OE, Huhndorf S, James T, Kirk PM, Lucking R, Lumbsch HT, Lutzoni F, Matheny PB, McLaughlin DJ, Powell MJ, Redhead S, Schoch CL, Spatafora JW, Stalpers JA, Vilgalys R, Aime MC, Aptroot A, Bauer R, Begerow D, Benny GL, Castlebury LA, Crous PW, Dai YC, Gams W, Geiser DM, Griffith GW, Gueidan C, Hawksworth DL, Hestmark G, Hosaka K, Humber RA, Hyde KD, Ironside JE, Koljalg U, Kurtzman CP, Larsson KH, Lichtwardt R, Longcore J, Miadlikowska J, Miller A, Moncalvo JM, Mozley-Standridge S, Oberwinkler F, Parmasto E, Reeb V, Rogers JD, Roux C, Ryvarden L, Sampaio JP, Schussler A, Sugiyama J, Thorn RG, Tibell L, Untereiner WA, Walker C, Wang Z, Weir A, Weiss M, White MM, Winka K, Yao YJ, Zhang N: A Higher-Level Phylogenetic Classification of the Fungi. Mycological Research. 2007, 111: 509-547. 10.1016/j.mycres.2007.03.004.
Lutzoni F, Kauff F, Cox C, McLaughlin D, Celio G, Dentinger B, Padamsee M, Hibbett D, James TY, Baloch E, Grube M, Reeb V, Hofstetter V, Schoch C, Arnold AE, Miadlikowska J, Spatafora J, Johnson D, Hambleton S, Crockett M, Shoemaker R, Sung G-H, Lücking R, Lumbsch T, O'Donnell K, Binder M, Diederich P, Ertz D, Gueidan C, Hansen K, Harris RC, Hosaka K, Lim Y-W, Matheny B, Nishida H, Pfister D, Rogers J, Rossman A, Schmitt I, Sipman H, Stone J, Sugiyama J, Yahr R, Vilgalys R: Assembling the Fungal Tree of Life: Progress, Classification, and Evolution of Subcellular Traits. American Journal of Botany. 2004, 91: 1446-1480. 10.3732/ajb.91.10.1446.
Blanco O, Crespo A, Ree RH, Lumbsch HT: Major Clades of Parmeliold Lichens (Parmeliaceae, Ascomycota) and the Evolution of Their Morphological and Chemical Diversity. Molecular Phylogenetics and Evolution. 2006, 39 (1): 52-69. 10.1016/j.ympev.2005.12.015.
Crespo A, Blanco O, Hawksworth DL: The Potential of Mitochondrial DNA for Establishing Phylogeny and Stabilising Generic Concepts in the Parmelioid Lichens. Taxon. 2001, 50 (3): 807-819. 10.2307/1223708.
Crespo A, Lumbsch HT, Mattsson JE, Blanco O, Divakar PK, Articus K, Wiklund E, Bawingan PA, Wedin M: Testing Morphology-Based Hypotheses of Phylogenetic Relationships in Parmeliaceae (Ascomycota) Using Three Ribosomal Markers and the Nuclear Rpb1 Gene. Molecular Phylogenetics and Evolution. 2007, 44 (2): 812-824. 10.1016/j.ympev.2006.11.029.
Crespo A, Divakar PK, Arguello A, Gasca C, Hawksworth DL: Molecular Studies on Punctelia Species of the Iberian Peninsula, with an Emphasis on Specimens Newly Colonizing Madrid. Lichenologist. 2004, 36: 299-308. 10.1017/S0024282904014434.
Nimis PL, Scheidegger C, Wolseley PA: Monitoring with Lichens Monitoring Lichens. 2002, Dordrecht, The Netherlands: Kluwer Academic Publishers
DePriest PT: Development of Mason E. Hale's List of Epithets in the Parmelioid Genera (Lichen-Forming Ascomycotina): A Bibliographic Review. Bryologist. 1999, 102 (3): 442-461. 10.2307/3244233.
Hale BW, DePriest PT: Mason E. Hale's List of Epithets in the Parmelioid Genera. Bryologist. 1999, 102 (3): 462-544. 10.2307/3244234.
Elix JA: Progress in the Generic Delimitation of Parmelia Sensu-Lato Lichens (Ascomycotina, Parmeliaceae) and a Synoptic Key to the Parmeliaceae. Bryologist. 1993, 96 (3): 359-383. 10.2307/3243867.
Divakar PK, Crespo A, Blanco O, Lumbsch HT: Phylogenetic Significance of Morphological Characters in the Tropical Hypotrachyna Clade of Parmelioid Lichens (Parmeliaceae, Ascomycota). Molecular Phylogenetics and Evolution. 2006, 40 (2): 448-458. 10.1016/j.ympev.2006.03.024.
Hale ME: A Revision of Lichen Genus Hypotrachyna (Parmeliaceae) in Tropical America. Smithsonian Contributions to Botany. 1975, 25: 1-73.
Divakar PK, Upreti DK: Parmelioid Lichens in India (a Revisionary Study). 2005, Dehra Dun: Bishen Singh Mahendra Pal Singh
Blanco O, Crespo A, Elix JA, Hawksworth DL, Lumbsch HT: A Molecular Phylogeny and a New Classification of Parmelioid Lichens Containing Xanthoparmelia-Type Lichenan (Ascomycota: Lecanorales). Taxon. 2004, 53 (4): 959-975. 10.2307/4135563.
Hale ME: A Synopsis of the Lichen Genus Xanthoparmelia (Vainio) Hale (Ascomycotina, Parmeliaceae). Smithsonian Contributions to Botany. 1990, 74: 1-250.
Honegger R: Metabolic Interactions at the Mycobiont-Photobiont Interface in Lichens. The Mycota Ix Fungal Associations. Edited by: Hock Be. 2001, Springer-Verlag, Berlin, Heidelberg, 209-221.
Green TGA, Lange OL, Cowan IR: Ecophysiology of Lichen Photosynthesis: The Role of Water Status and Thallus Diffusion Resistances. Cryptogamic Botany. 1994, 4: 166-178.
Lange OL: Ecophysiology of Photosynthesis – Performance of Poikilohydric Lichens and Homoiohydric Mediterranean Sclerophylls – the 7th Tansley Lecture. Journal of Ecology. 1988, 76 (4): 915-937. 10.2307/2260624.
Beckett RP: Pressure-Volume Analysis of a Range of Poikilohydric Plants Implies the Existence of Negative Turgor in Vegetative Cells. Annals of Botany. 1997, 79 (2): 145-152. 10.1006/anbo.1996.0318.
Green TGA, Lange OL: Photosynthesis in Poikilohydric Plants: A Comparison of Lichens and Bryophytes. Ecophysiology of Photosynthesis. Edited by: Schulze E-D, Caldwell MCe. 1994, Ecological Studies, Springer-Verlag, Berlin, Heidelberg, New York, 319-341.
Kappen L, Valladares F: Opportunistic Growth and Desiccation Tolerance: The Ecological Success of Poikilohydrous Autotrophs. Handbook of Functional Plant Ecology. Edited by: Pugnaire FI, Valladares Fe. 1999, Marcel Dekker, Inc., New York, Basel, 9-80.
Kranner I, Lutzoni F: Evolutionary Consequences of Transition to a Lichen Symbiotic State and Physiological Adaptation to Oxidative Damage Associated with Poikilohydry. Plant Responses to Environmental Stresses: From Phytohormones to Genome Reorganization. Edited by: Lerner HRe. 1999, Marcel Dekker, Inc., New York, Basel, 591-628.
Kershaw KA: Physiological Ecology of Lichens. 1985, Cambridge University Press, Cambridge
Nash TH, (ed.): Lichen Biology. 1996, Cambridge University Press, Cambridge
Sancho LG, Pintado A, Valladares F, Schroeter B, Schlensog M: Photosynthetic Performance of Cosmopolitan Lichens in the Maritime Antarctic. New Species and Novel Aspects in Ecology and Physiology of Lichens in Honour of O L Lange. Edited by: Kappen Le. 1997, Bibliotheca Lichenologica, J. Cramer, Berlin, Stuttgart, 197-210.
del Prado R, Sancho LG: Dew as a Key Factor for the Distribution Pattern of the Lichen Species Teloschistes Lacunosus in the Tabernas Desert (Spain). Flora. 2007, 202 (5): 417-428.
Kappen L, Lange OL, Schulze ED, Evenari M, Buschbom U: Ecophysiological Investigations on Lichens of the Negev Desert .6. Annual Course of the Photosynthetic Production of Ramalina-Maciformis (Del) Bory. Flora. 1979, 168 (1–2): 85-108.
Lange OL: 23 Years of Growth Measurements on the Crustose Lichen Caloplaca-Aurantia in the Central Negev Desert. Israel Journal of Botany. 1990, 39 (4–6): 383-394.
Lange OL, Schulze ED, Kappen L, Buschbom U, Evenari M: Adaptations of Desert Lichens to Drought and Extreme Temperatures. Environmental Physiology of Desert Organisms. 1975, NFH: Dowden, Hutchinson & Ross, Stroudsburg, PA, 20-37.
Lange OL, Meyer A, Zellner H, Ullmann I, Wessels DCJ: Eight Days in the Life of a Desert Lichen: Water Relations and Photosynthesis of Teloschistes Capensis in the Coastal Fog Zone of the Namib Desert. Madoqua. 1990, 17 (1): 17-30.
Lange OL, Meyer A, Ullmann I, Zellner H: Microclimate Conditions, Water-Content and Photosynthesis of Lichens in the Coastal Fog Zone of the Namib Desert – Measurements in the Fall. Flora. 1991, 185 (4): 233-266.
Lange OL, Budel B, Meyer A, Zellner H, Zotz G: Lichen Carbon Gain under Tropical Conditions: Water Relations and Co2 Exchange of Lobariaceae Species of a Lower Montane Rainforest in Panama. Lichenologist. 2004, 36: 329-342. 10.1017/S0024282904014392.
Lange OL, Budel B, Heber U, Meyer A, Zellner H, Green TGA: Temperate Rain-Forest Lichens in New-Zealand – High Thallus Water-Content Can Severely Limit Photosynthetic Co2 Exchange. Oecologia. 1993, 95 (3): 303-313. 10.1007/BF00320981.
Lange OL, Budel B, Meyer A, Zellner H, Zotz G: Lichen Carbon Gain under Tropical Conditions: Water Relations and Co2 Exchange of Three Leptogium Species of a Lower Montane Rainforest in Panama. Flora. 2000, 195 (2): 172-190.
Butler MA, King AA: Phylogenetic Comparative Analysis: A Modeling Approach for Adaptive Evolution. American Naturalist. 2004, 164 (6): 683-695. 10.1086/426002.
Crespo A, Cubero OF: A Molecular Approach to the Circumscription and Evaluation of Some Genera Segregated from Parmelia S. Lat. Lichenologist. 1998, 30: 369-380.
Arbiza L, Dopazo J, Dopazo H: Positive Selection, Relaxation, and Acceleration in the Evolution of the Human and Chimp Genome. Plos Computational Biology. 2006, 2 (4): 288-300. 10.1371/journal.pcbi.0020038.
Haydon DT, Bastos AD, Knowles NJ, Samuel AR: Evidence for Positive Selection in Foot-and-Mouth Disease Virus Capsid Genes from Field Isolates. Genetics. 2001, 157 (1): 7-15.
Urwin R, Holmes EC, Fox AJ, Derrick JP, Maiden MCJ: Phylogenetic Evidence for Frequent Positive Selection and Recombination in the Meningococcal Surface Antigen Porb. Molecular Biology and Evolution. 2002, 19 (10): 1686-1694.
Itoh T, Martin W, Nei M: Acceleration of Genomic Evolution Caused by Enhanced Mutation Rate in Endocellular Symbionts. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (20): 12944-12948. 10.1073/pnas.192449699.
Kohne DE: Evolution of Higher-Organism DNA. Quarterly Reviews of Biophysics. 1970, 3: 327-375.
Martin AP, Palumbi SR: Body Size, Metabolic-Rate, Generation Time, and the Molecular Clock. Proceedings of the National Academy of Sciences of the United States of America. 1993, 90 (9): 4087-4091. 10.1073/pnas.90.9.4087.
Rand DM: Thermal Habit, Metabolic-Rate and the Evolution of Mitochondrial-DNA. Trends in Ecology & Evolution. 1994, 9 (4): 125-131. 10.1016/0169-5347(94)90176-7.
Moran NA: Accelerated Evolution and Muller's Rachet in Endosymbiotic Bacteria. Proceedings of the National Academy of Sciences of the United States of America. 1996, 93 (7): 2873-2878. 10.1073/pnas.93.7.2873.
Rispe C, Moran NA: Accumulation of Deleterious Mutations in Endosymbionts: Muller's Ratchet with Two Levels of Selection. American Naturalist. 2000, 156 (4): 425-441. 10.1086/303396.
Lynch M, Blanchard JL: Deleterious Mutation Accumulation in Organelle Genomes. Genetica. 1998, 103: 29-39. 10.1023/A:1017022522486.
Kisakibaru Y, Matsuda H: Nucleotide Substitution Type Dependence of Generation Time Effect of Molecular Evolution. Japanese Journal of Genetics. 1995, 70 (3): 373-386. 10.1266/jjg.70.373.
Moran NA, Munson MA, Baumann P, Ishikawa H: A Molecular Clock in Endosymbiotic Bacteria Is Calibrated Using the Insect Hosts. Proceedings of the Royal Society of London Series B-Biological Sciences. 1993, 253 (1337): 167-171. 10.1098/rspb.1993.0098.
Moran NA, vonDohlen CD, Baumann P: Faster Evolutionary Rates in Endosymbiotic Bacteria Than in Cospeciating Insect Hosts. Journal of Molecular Evolution. 1995, 41 (6): 727-731. 10.1007/BF00173152.
Moriyama EN: Higher Rates of Nucleotide Substitution in Drosophila Than in Mammals. Japanese Journal of Genetics. 1987, 62 (2): 139-147. 10.1266/jjg.62.139.
Lange OL: Twenty-Three Years of Growth Measurements on the Crustose Lichen Caloplaca Aurantia in the Central Negev Desert. Israel Journal of Botany. 1990, 39: 383-394.
Lange OL, Green TGA, Heber U: Hydration-Dependent Photosynthetic Production of Lichens: What Do Laboratory Studies Tell Us About Field Performance?. Journal of Experimental Botany. 2001, 52 (363): 2033-2042. 10.1093/jexbot/52.363.2033.
Boucher VL, Nash TH: Growth-Patterns in Ramalina-Menziesii in California – Coastal Vs Inland Populations. Bryologist. 1990, 93 (3): 295-302. 10.2307/3243516.
Louwhoff S, Crisp MD: Phylogenetic Analysis of Parmotrema (Parmeliaceae: Lichenized Ascomycotina). Bryologist. 2000, 103 (3): 541-554. 10.1639/0007-2745(2000)103[0541:PAOPPL]2.0.CO;2.
Kappen L: Ecophysiological Relationships in Different Climatic Regions. Crc Handbook of Lichenology. Edited by: Galun M. 1988, Boca Raton: CRC Press Inc., 2: 37-100.
Blanco O, Crespo A, Divakar PK, Elix JA, Lumbsch HT: Molecular Phylogeny of Parmotremoid Lichens (Ascomycota, Parmeliaceae). Mycologia. 2005, 97 (1): 150-159. 10.3852/mycologia.97.1.150.
Wright SD, Gray RD, Gardner RC: Energy and the Rate of Evolution: Inferences from Plant Rdna Substitution Rates in the Western Pacific. Evolution. 2003, 57 (12): 2893-2898.
Lange OL, Redon J: Epiphytic Lichens in the Region of a Chilean Fog Oasis (Fray Jorge) .2. Ecophysiological Characterization of Co2 Exchange and Water Relations. Flora. 1983, 174 (3–4): 245-284.
Blanco O, Crespo A, Divakar PK, Esslinger TL, Hawksworth DL, Lumbsch HT: Melanelixia and Melanohalea, Two New Genera Segregated from Melanelia (Parmeliaceae) Based on Molecular and Morphological Data. Mycological Research. 2004, 108: 873-884. 10.1017/S0953756204000723.
Doring H, Clerc P, Grube M, Wedin M: Mycobiont-Specific Pcr Primers for the Amplification of Nuclear Its and Lsu Rdna from Lichenized Ascomycetes. Lichenologist. 2000, 32: 200-204. 10.1006/lich.1999.0250.
Vilgalys R, Hester M: Rapid Genetic Identification and Mapping of Enzymatically Amplified Ribosomal DNA from Several Cryptococcus Species. Journal of Bacteriology. 1990, 172: 4238-4246.
Zoller S, Scheidegger C, Sperisen C: Pcr Primers for the Amplification of Mitochondrial Small Subunit Ribosomal DNA of Lichen-Forming Ascomycetes. Lichenologist. 1999, 31: 511-516.
Zhou S, Stanosz GR: Primers for Amplification of Mt Ssu Rdna, and a Phylogenetic Study of Botryosphaeria and Associated Anamorphic Fungi. Mycological Research. 2001, 105: 1033-1044. 10.1016/S0953-7562(08)61965-6.
Lumbsch HT, Schmitt L, Lindemuth R, Miller A, Mangold A, Fernandez F, Huhndorf S: Performance of Four Ribosomal DNA Regions to Infer Higher-Level Phylogenetic Relationships of Inoperculate Euascomycetes (Leotiomyceta). Molecular Phylogenetics and Evolution. 2005, 34 (3): 512-524. 10.1016/j.ympev.2004.11.007.
Thompson JD, Higgins DG, Gibson TJ: Clustal-W – Improving the Sensitivity of Progressive Multiple Sequence Alignment through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. Nucleic Acids Research. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.
Karplus K, Barrett C, Hughey R: Hidden Markov Models for Detecting Remote Protein Homologies. Bioinformatics. 1998, 14 (10): 846-856. 10.1093/bioinformatics/14.10.846.
Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP: Evolution – Bayesian Inference of Phylogeny and Its Impact on Evolutionary Biology. Science. 2001, 294 (5550): 2310-2314. 10.1126/science.1065889.
Larget B, Simon DL: Markov Chain Monte Carlo Algorithms for the Bayesian Analysis of Phylogenetic Trees. Molecular Biology and Evolution. 1999, 16 (6): 750-759.
Swofford DL: Paup – a Computer-Program for Phylogenetic Inference Using Maximum Parsimony. Journal of General Physiology. 1993, 102 (6): A9-A9.
Felsenstein J: Confidence-Limits on Phylogenies – an Approach Using the Bootstrap. Evolution. 1985, 39 (4): 783-791. 10.2307/2408678.
Huelsenbeck JP, Ronquist F: Mrbayes: Bayesian Inference of Phylogenetic Trees. Bioinformatics. 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.
Rodriguez F, Oliver JL, Marin A, Medina JR: The General Stochastic-Model of Nucleotide Substitution. Journal of Theoretical Biology. 1990, 142 (4): 485-501. 10.1016/S0022-5193(05)80104-3.
Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian Phylogenetic Analysis of Combined Data. Systematic Biology. 2004, 53 (1): 47-67. 10.1080/10635150490264699.
Page RDM: Treeview: An Application to Display Phylogenetic Trees on Personal Computers. Computer Applied Biosciences. 1996, 12: 357-358.
Yang Z: Paml: A Program Package for Phylogenetic Analysis by Maximum Likelihood. Computer Applied Biosciences. 1997, 13 (5): 555-556.
Felsenstein J: Evolutionary Trees from Gene-Frequencies and Quantitative Characters – Finding Maximum-Likelihood Estimates. Evolution. 1981, 35 (6): 1229-1242. 10.2307/2408134.
Emberger L: La Végétation De La Région Méditérranéenne. Essai D'une Classification Des Groupements Végétaux. Revue de Botanique. 1930, 503: 642-662.
Martins EP, Hansen TF: Phylogenies and the Comparative Method: A General Approach to Incorporating Phylogenetic Information into the Analysis of Interspecific Data. American Naturalist. 1997, 149 (4): 646-667. 10.1086/286013.
Hansen TF, Martins EP: Translating between Microevolutionary Process and Macroevolutionary Patterns: The Correlation Structure of Interspecific Data. Evolution. 1996, 50 (4): 1404-1417. 10.2307/2410878.
Burnham KP, Anderson DR: Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2002, New York: Springer
Hipp AL: Nonuniform Processes of Chromosome Evolution in Sedges (Carex: Cyperaceae). Evolution. 2007, 61 (9): 2175-2194. 10.1111/j.1558-5646.2007.00183.x.
Harmon LJ, Weir J, Brock C, Glor RE, Challenger W: Geiger: Investigating Evolutionary Radiations. Bioinformatics. 2008, 24: 129-131. 10.1093/bioinformatics/btm538.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed Phylogenetics and Dating with Confidence. Plos Biology. 2006, 4 (5): 699-710. 10.1371/journal.pbio.0040088.
Drummond AJ, Rambaut A: Beast: Bayesian Evoluionary Analysis by Sampling Trees. Bmc Evolutionary Biology. 2007, 7: 214-10.1186/1471-2148-7-214.
We thank Rick Ree (Chicago), Ruth Del Prado (Madrid) and Ana Pintado (Madrid) for helpful discussions. Most of the sequencing was performed in the Centro de Genómica y Proteómica del Parque Científico de Madrid, and some sequences were obtained in the Pritzker Laboratory for Molecular Systematics at The Field Museum (Chicago). This work was supported by the Ministry of Education and Science of Spain (CGL 2004-01848/BOS, CGL 2004-20423-E) to AC, and a Juan de la Cierva award to PKD.
HTL carried out the phylogenetic and rate analyses and drafted the manuscript. ALH carried out the phylogenetic comparative analyses and helped draft the manuscript. OB and PKD contributed DNA sequences and helped draft the manuscript. AC initiated and coordinated the work, helped draft the manuscript, and submitted the manuscript.