Open Access Open Badges Research article

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

Steven D Leavitt1*, Theodore L Esslinger2, Pradeep K Divakar3 and H Thorsten Lumbsch1

Author Affiliations

1 The Field Museum, Department of Botany, 1400 S. Lake Shore Drive, Chicago, IL, 60605, USA

2 Department of Biological Sciences #2715, North Dakota State University, PO Box 6050, Stevens Hall, Fargo, ND, 58108-6050, USA

3 Departamento de Biología Vegetal II, Facultad de Farmacia, Universidad Complutense de Madrid, Madrid, 28040, Spain

For all author emails, please log on.

BMC Evolutionary Biology 2012, 12:176  doi:10.1186/1471-2148-12-176

The electronic version of this article is the complete one and can be found online at:

Received:24 February 2012
Accepted:30 August 2012
Published:11 September 2012

© 2012 Leavitt et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



Factors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescent-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies.


We found that divergence estimates, from both concatenated gene tree and coalescent-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages.


Our results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.

Bayesian skyline plot; BEAST; Divergence times; Coalescence; Concatenation; Lichens; Melanohalea; Pleistocene; Reproductive strategy; Species tree


Lichenized fungi form mutualistic relationships with photoautotrophic organisms (photobionts), mainly green algae (Trebouxiophyceae and Ulvophyceae) and/or cyanobacteria. The lichen symbiosis has been highly successful within fungi, especially Ascomycota, with more than 18,000 currently accepted species [1] and an estimated diversity of greater than 28,000 species [2]. Due to the availability of genetic data and analytical improvements, DNA-based approaches play an increasing role in the recognition of diversity in lichenized fungi that would otherwise be impossible to recognize using classical phenotypic characters due to morphological convergence or parallelism [3-11].

In spite of the recent advancements in recognizing diversity in lichen-forming fungi, assessing diversification within a temporal context remains largely unexplored in nearly all groups of these important fungal symbionts, (exceptions include [12-14]). This is largely due to the poor fossil record for lichenized fungi, and also fungi in general, and uncertainties in the interpretation of the few known fossil records [15-17]. However, the timing of speciation events plays a valuable role, complementary to discovering and describing diversity, by addressing biogeographical, climatic, ecological, and other factors associated with diversification and extinction within a temporal context, (e.g. [12,18-20]).

In spite of difficulties in obtaining accurate estimates of divergence times [21,22], recent analytical advances in using relaxed molecular clocks, the inclusion of multiple fossil calibrations, and informative priors on substitution rates have improved accuracy in molecular dating [23-26]. However, a recent study also suggested that methods that fail to incorporate the process of gene lineage coalescence may provide inaccurate estimates of divergence dates [18]. While gene trees, including phylogenies estimated from concatenated sequence data, can overestimate divergence times because they do not correct for genetic divergence that predates speciation, species tree methods incorporating the process of gene lineage coalescence likely provide a more biologically realistic framework for dating divergence events [18].

Within lichen-forming ascomycetes, Parmeliaceae (Lecanorales) constitutes one of the largest and best-studied families [27-32]. Although some lichen-forming fungal lineages have geographically restricted distributions, (e.g. [10,14]), there is mounting evidence that transoceanic dispersal has commonly occurred within some lichenized fungi and has played a major role in diversification within Parmeliaceae [12,33,34].

Unrecognized diversity is common in many phenotype-based species circumscriptions in the family, confounding the current interpretation of ecological and biogeographical patterns [7,9,10,35,36]. Additionally, interpreting biogeographical patterns and factors driving diversification is further complicated by a common pattern of long-distance dispersal in many taxa in Parmeliaceae, (e.g. [12,33,34,37]). Therefore, accurate estimates of divergence times are especially critical for identifying major factors driving diversification. A recent study investigating the origin and diversification in the largest clade of macrolichens, the parmelioid lichens (Parmeliaceae, Ascomycota), provides valuable taxon-specific estimates of substitution rates for three genetic markers based on three calibration points, including two dated fossils from Parmeliaceae [12].

Within Parmeliaceae, the genus Melanohalea O. Blanco et al. includes 22 species based on traditional taxonomy, most of which occur primarily on bark and wood throughout the Holarctic [38-44]. Only four known taxa occur in the southern Hemisphere [45-47]. Many Melanohalea species display broad geographic and ecological distributions, although a limited number of taxa appear to have more restricted ranges [43]. Otte et al. [43] suggested that distribution patterns in Melanohalea are largely determined by contemporary ecogeographical factors, and most species have largely filled their potential distributions in the northern Hemisphere. Furthermore, the distributions of some Melanohalea taxa, including M. elegantula (Zahlbr.) O. Blanco et al. and M. exasperatula (Nyl.) O. Blanco et al., have been linked to eutrophication, air pollution, and other anthropogenic factors [48-52].

While Pleistocene events have been shown to have been important factors driving diversification and affecting distributions in many groups, (e.g. [18,53-55]), recent estimates suggest that the Melanohalea radiation vastly predates Pleistocene glacial cycles [12]. Within many parmelioid genera, major radiations occurred from the Late Oligocene to the early Pliocene, before the climate became cooler, drier, and more seasonal at the end of the Pliocene [12,56]. However, the overall importance of Pleistocene glacial cycles in diversification within Melanohalea is unknown.

Reproductive strategies vary among Melanohalea taxa. Sexual reproduction is restricted to characteristic fungal fruiting bodies (ascomata) producing meiospores (=ascospores), and is common in at least 13 of the 22 described species. Ascospores are dispersed independent of the photosynthesizing partner (photobiont) and require reacquisition of the appropriate photobiont partner in order to re-establish the lichenized condition. In contrast, other species within Melanohalea commonly propagate asexually by means of vegetative diaspores, either isidia or soredia. These specialized vegetative reproductive propagules contain both fungal and algal symbionts, eliminating the need for independent acquisition of the appropriate photobiont partner. The isidiate taxa M. elegantula and M. exasperatula show a remarkable potential for dispersal and may be spreading in some areas [43,48,49]. Poelt [57] hypothesized that lichenized fungi reproducing asexually (via soredia and isidia) are generally more successful in pioneering formerly glaciated areas than forms that reproduce sexually. In contrast, Nimis [58] argued that the distribution of Holarctic lichens is more likely determined by general ecology than by their reproductive strategy alone. The latter argument supports the assumption that the distribution of Melanohalea within the Holarctic today widely reflects their biological constitution, rather than their geographic origin or dispersal strategy [43]. However, population structure and history is poorly understood in most lichen-forming ascomycetes and the role of their reproductive strategy in response to climate fluctuation remains unclear.

Cryptic lineages within phenotypically circumscribed taxa are common in Parmeliaceae [35], and previously unrecognized species-level lineages have now been recognized within six of the phenotype-based Melanohalea species from a broad sampling of populations in the northern Hemisphere [59]. Diversity within this genus is now well-characterized, at least in the northern Hemisphere, providing an excellent study system with which to test the relative influence of Miocene orogeny and climatic conditions and Pleistocene glacial cycles on common lichen-forming fungi.

In this study our goals are twofold: (1) we aim to estimate divergence times in the lichenized genus Melanohalea using both concatenated gene tree and coalescent-based multilocus species tree approaches; and (2) we evaluate the impact of Pleistocene glacial cycles on the population demography among four common sexually reproducing lichen-forming fungal species and two taxa reproducing largely via vegetative diaspores. With estimates of divergence times, we examine the relative roles of Miocene orogeny and climate change and Pleistocene glacial cycles on diversification in the lichenized fungal genus Melanohalea. We are also interested in population demographic changes in common Melanohalea species after the last glacial maximum (LGM), including anthropogenic factors. Here we present estimates of divergence times within Melanohalea and assess population demographic histories in relation to the LGM.


The complete ITS data matrix consisted of 487 sequences and 511 aligned nucleotide position characters (Additional file 1; TreeBase ID: 12364). The six locus data matrix, representing genetic diversity identified from the ITS gene tree, consisted of 138 samples and 3839 aligned nucleotide position characters (Additional file 2; TreeBase ID: 12364). All sequences generated for this study have been deposited in GenBank under accession Nos. JQ812998 – JQ814066. Table 1 summarizes patterns of variation in sampled loci and the best-fit model of evolution selected using the Akaike Information Criterion (AIC) in jModeltest.

Additional file 1. Collection information for all Melanohalea specimens included in the present study.

Format: XLSX Size: 55KB Download fileOpen Data

Additional file 2. Maximum likelihood ITS gene tree of the 487 sampled Melanohalea specimens. Bootstrap support indicated at nodes, and operational taxonomic units in red text indicate 138 Melanohalea specimens selected to represent sampled genetic diversity in multilocus phylogenetic reconstructions.

Format: PDF Size: 380KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Table 1. Genetic variability in Melanohalea of sampled markers used in this study, including: the number of specimens, N , alignment length (number of base pairs); variable and parsimony-informative (PI) sites for each sampled locus; and locus-specific model of evolution identified using the Akaike information criterion in jModeltest

Age estimates on concatenated and species topologies

High posterior effective sample sizes (ESS > 200) were observed for all parameters in the BEAST analyses. The substitution rates in Melanohalea of the six sampled loci (ITS, nuLSU, mtSSU, MCM7, RPB1, and RPB2), estimated under a relaxed clock, are shown in Table 2. The estimated substitution rate in the nuclear ITS was nearly five times higher than the mtSSU and nuLSU, and over two times higher than the sampled protein-coding loci. The protein-coding markers had relatively similar estimated substitution rates.

Table 2. Estimates of substitution rates in Melanohalea from the BEAST analysis of the concatenated six-loci dataset estimated under a relax molecular clock using fixed substitution rates for the nuLSU, mtSSU, and RPB 1 markers [12]

Our calibrated maximum clade credibility chronograms from analyses of the concatenated data matrix (ITS, nuLSU, mtSSU, MCM7, RPB1, RPB2) and multi-locus species-tree approach are shown in Figures 1 and 2, respectively. Divergence times estimated from the species tree approach were between 0.2 and 8.3 million years more recent across the genus, with smaller differences at more recent timescales, i.e. divergences that occurred during the Pleistocene (Figures 1 and 2; Table 3).

thumbnailFigure 1. Dated BEAST maximum clade credibility tree for Melanohalea estimated from concatenated data. The chronogram was estimated from a partitioned data set consisting of six loci (ITS, nuLSU, mtSSU, MCM7, RPB1, RPB2) under a relaxed molecular clock. The divergence times correspond to the mean posterior estimate of their age in millions of years. The bars indicate the 95% HPD interval for the divergence times estimates. The four major Melanohalea clades identified in this study are indicated by ‘1’, ‘2’, ‘3’, and ‘4’ at the corresponding node. Values on branches indicate posterior probability, and only support indices posterior probability values > 0.50 are presented.

thumbnailFigure 2. Dated multilocus coalescent-based species tree for Melanohalea. The chronogram was estimated from a multilocus data (ITS, nuLSU, mtSSU, MCM7, RPB1, RPB2) within a coalescent-based framework in *BEAST. The divergence times correspond to the mean posterior estimate of their age in millions of years. The bars indicate the 95% HPD interval for the divergence times estimates. The four major Melanohalea clades identified in this study are indicated by ‘1’, ‘2’, ‘3’, and ‘4’ at the corresponding node. Values on branches indicate posterior probability, and only support indices posterior probability values > 0.50 are presented.

Table 3. The dates of origin of lineages of Melanohalea from their initial divergence (ancestral splits; stem origin) and the time to the most recent common ancestor (tmrca) of extant haplotypes

The divergence of MelanelixiaMelanohalea, estimated from the concatenated six gene dataset, occurred ca. 43.5 (95% HPD = 26.9-59.7) million years ago (Ma) and is similar to a previous estimate of 41.55 Ma [12]. However, in the multilocus species tree approach this split was estimated to have occurred more than eight million years more recently, ca. 35.2 (95% HPD = 27.0-43.5) Ma (Table 3). Divergence dates estimated in the concatenated multilocus analysis support the initial radiation of Melanohalea during the Oligocene, ca. 30.5 (95% HPD = 23.8-38.1) Ma, while divergence estimates from the multilocus species tree approach suggests a more recent radiation during the Oligocene-Miocene boundary, ca. 23.7 (95% HPD = 19.4-27.9) Ma (Table 3). Both the concatenated and species tree approaches suggest that most diversification of Melanohalea occurred throughout the Miocene and Pliocene, and divergence estimates suggest that diversification during the Pleistocene glacial cycles was not accompanied by speciation in Melanohalea (Figures 1 and 2). Divergence dates from the multilocus concatenated analysis suggests a relatively synchronous initial radiation of the four major clades in Melanohalea with mean diversification times occurring between 14.9 – 17.2 Ma (Figure 1), with no clear temporal pattern in subsequent diversification events. However, dates estimated from the species tree approach provide a somewhat different perspective of diversification, with mean diversification times occurring between 9.1 -12.9 Ma (Table 3). The greatest disparity found in estimates was in the radiation of clade IV, estimated to begin radiating 17.2 (95% HPD = 12.8-22.5) and 9.1 (95% HPD = 6.2-12.4) in the concatenated and species trees approaches, respectively (Figure 2).

In the dated chronogram estimated from concatenated data, species divergences that occurred during Pleistocene glacial cycles were limited to two distinct clades, a clade containing M. subexasperata and M. olivaceoides 2 and a clade with M. gomukhensis, M. olivaceoides 1, and M. poeltii. Dates estimated in the species tree analysis suggest Pleistocene divergence for two additional clades, a split between the two lineages in the M. subolivacea complex and a split between M. elegantula and M. elegantula Spain (Table 3).

Historical demography

Genetic diversity indices (Hd, S, and π) for lineages previously recognized in [59] are summarized in Table 4. Overall, lineages producing sexual fruiting structures (ascomata) showed greater nucleotide and haplotype diversities, relative to lineages that commonly reproduce via vegetative diaspores (isidia/soredia). However, haplotype diversity was relatively high in the isidiate taxon M. laciniatula (0.833), and, in contrast, haplotype diversity was low in the apotheciate lineage M. multispora 1 (0.143).

Table 4. Estimates of genetic diversity for sampled lineages within Melanohalea

Although both Tajima’s D and Fu’s F statistic values were negative in most cases (Table 4), Tajima’s D statistic was significant (P < 0.05) and negative for M. elegantula (−2.137) and M. subelegantula (−1.741); and Fu’s F statistic was significant and negative for M. elegantula (−3.011) and M. olivacea (−2.6561). The mismatch distributions generally showed relatively unimodal patterns for taxa commonly producing vegetative diaspores (M. elegantula, M. exasperatula, M. infumata, and M. laciniatula), and mismatch distributions largely appeared multimodal for taxa producing ascomata (M. exasperata, M. multispora 3, M. septentrionalis, M. subolivacea 1, M. subolivacea 2, and M. trabeculata) (Additional file 3). No tests were carried out for M. exasperata Spain, M. gomukhensis, M. halei, M. elegantula Spain, M. lobulata, M. multispora 2, M. olivaceoides 1, M. olivaceoides 2, M. poeltii, M. subexasperata, M. ushuaiensis Argentina, and M. ushuaiensis Chile, due either to their lack of polymorphisms or small sample sizes.

Additional file 3. Mismatch distributions observed in Melanohalea lineages with a minimum sample size of eight individuals. Black squares: simulated mismatch distributions, black triangles: observed mismatch distributions.

Format: PDF Size: 88KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

The effective population sizes and demographic trends estimated by the Bayesian skyline plot (BSP) analyses for M. elegantula, M. exasperatula, M. multispora 3, M. olivacea, M. septentrionalis, and M. subolivacea indicated population size increases during the Pleistocene for the four sexually reproducing species and two taxa commonly producing vegetative diaspores (Figure 3A-B). However, BSP plot analyses suggest population expansion predated the end of the LGM in all sampled species. Our data suggest that the overall increase of the two taxa commonly producing vegetative diaspores occurred much more recently than the expansions of the lineages strictly reproducing through sexual propagation (i.e. taxa that do not produce vegetative diaspores), with the exception of M. septentrionalis (Figure 3A-B). Based on the mutation rate of 0.033 s/s/Myr for the ITS marker, population growth started at approximately 200,000 years before present (BP) in the sampled isidiate taxa, M. elegantula and M. exasperatula (Figure 3A). In contrast, population expansion for the sexually reproducing taxa M. multispora, M. olivacea, and M. subolivacea appears to have started approximately 750,000 to 100,000 years BP (Figure 3B). However, the BSP for the apotheciate taxon M. septentrionalis suggests a more recent expansion, relative to the other apotheciate taxa (Figure 3B).

thumbnailFigure 3. Bayesian skyline plots for (A) two Melanohalea species commonly propagating via vegetative diaspores (M. elegantula and M. exasperatula ) and (B) four sexually reproducing lineages (M. multispora 3, M. olivacea , M. septentrionalis , and M. subolivacea 2). The solid black line represents the median value for the log of the population size (log Ne) and the grey lines represent the upper and lower 95% credible intervals. The x-axis measures times in thousands of years. Generation times are not known for Melanohalea species and population sizes (y-axis) only represent relative changes.


Concatenated gene tree and multilocus divergence estimates

In this study, we estimated divergence times for the lichen-forming fungal genus Melanohalea using Parmeliaceae-specific substitution rates for three DNA regions. McCormack et al. [18] suggested that estimating divergence times using a multilocus species tree approach is especially relevant for investigating taxonomic groups with poor fossil representation, and here we provide the first implementation of calibrated coalescent-based species-tree divergence time estimates for lichen-forming fungi.

Similar to a previous study [18], we found that divergence times varied greatly among estimates from concatenated gene tree and multilocus species-tree approaches. A comparison of the divergence estimates from the concatenated gene tree versus multilocus coalescent-based species-tree showed a consistent pattern of more recent divergence times estimated from the species trees approach, regardless of node age (Figure 2, Table 3). Although gene trees will necessarily overestimate divergence times in recent speciation events [60], in this study we found the greatest disparity in estimates at deeper relationships, differing by over 8 million years at some nodes (Table 3). Additionally, we found that the multilocus coalescent-based species tree approach consistently provided narrower 95% confidence intervals than the concatenated gene tree approach (Table 3).

Topologies derived from concatenated datasets are not necessarily equivalent to species trees [61-65]. In addition to inferring incorrect relationships and estimates of support, using concatenated gene trees for estimating divergence times will likely produce overestimates for the divergence events [18,66]. In contrast, methods incorporating a full probabilistic coalescent framework, such as the *BEAST analysis implemented herein, have been proposed to provide more biologically realistic estimates because multiple gene trees provide independent realizations of divergence history, accounting for mutational and coalescent stochasticity [18].

Pre-Pleistocene divergence

We estimated that the Melanelixia- Melanohalea split occurred during the Eocene, with stem node age estimates of 41.6 or 35.2 Ma in the concatenated and multilocus species tree approaches, respectively (Figures 1, 2, Table 3). A previous study suggested that the radiation of major parmelioid genera started at the end of the Eocene [12], and our estimates support a similar pattern of divergence for the Melanelixia and Melanohalea split. Our data suggest that the radiation of Melanohalea began either during the Oligocene, or the transition at the Oligocene – Miocene boundary, with crown node age estimates of 30.5 Ma and 23.7 Ma in the concatenated and multilocus species tree approaches, respectively (Figures 1, 2, Table 3). The date from the concatenated analysis is approximately five million years older than that presented in Amo de Paz et al. [12], but this is most likely due to our inclusion of the early diverging South American endemic M. ushuaiensis sensu lato, which was absent from their analysis.

Our results also suggest that most major speciation events within Melanohalea occurred during the Miocene and Pliocene, and were likely the result of complex patterns associated with new habitats formed during Miocene orogeny events, major climatic changes, and global shifts in vegetation. The Miocene was a time of major climatic and vegetative changes worldwide, especially in the northern Hemisphere, including major tectonic activity and orogeny [18,67-70]. Increasing aridity in the middle Miocene (15–8 Ma; [71,72]) resulted in woodlands giving way to more open habitats and grasslands [69,73-75]. However, given the incredible dispersal capacity of some Melanohalea species, and other fungi within Parmeliaceae (e.g. [12,33,34,43]), inferring the geographic origin of lineages and associated temporal biogeographic factors remains challenging.

Distributions of most Melanohalea species are largely determined by contemporary ecogeographical factors, and most species have largely filled their potential distributions in the northern Hemisphere [43]. In general, the distributions are in high concordance with those found in vascular plants and can be widely explained by contemporary ecogeographical factors [76]. Although the distribution of lichenized fungal species is limited by the availability of the appropriate substrata, other unrecognized dispersal barriers appear to have limited the dispersal of some Melanohalea species. For example, the amphiatlantic distribution of M. exasperata and the strictly European distribution of M. laciniatula remain unexplained. Although we currently do not know the origin of the distribution patterns in these two species, incomplete spread into North America, particularly western North America, cannot be excluded [43].

Interestingly, within the morphologically cryptic M. multispora and M. ushuaiensis species complexes, species-level lineages also appear to have diversified largely during the Miocene (Figures 1, 2). The two defined M. olivaceoides sensu lato (s.l.) lineages [59] were not recovered as sister taxa and their most recent common ancestor (MRCA) was estimated over 20 Ma (Figures 1, 2). Similarly, the morphologically similar taxa M. elegantula and M. infumata are found in separate major clades within Melanohalea, with the MRCA estimated to have been over 20 Ma in the dated coalescent-based species tree (Figure 2). These results suggest that in some cases, i.e. M. multispora s.l and M. ushuaiensis s.l, taxonomically diagnostic characters, including number of spores and isidia morphology may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters, among non-sister taxa, including reproductive strategies (i.e. the production of isidia, soredia, and apotheica), suggest the potential for convergent evolution due to similar selective pressures among distinct lineages inhabiting similar environments [7,77,78].

Climatic changes during the Pleistocene are hypothesized as one of the foremost contributors to biological diversity (e.g. [18,53-55]). However, in this study, we show that species divergence during Pleistocene glacial cycles was restricted to only a limited number of clades. Of the morphologically indistinguishable species-level lineages recognized in [59], only M. subolivacea sensu lato and M. elegantula sensu lato appear to have diversified during the Pleistocene. However, additional supporting evidence, such as morphology, geographic range, ecological preference, chemistry, etc., will be required to establish robust species boundaries among the recently derived lineages.

Our data suggests that two other clades also diversified during the Pleistocene, one containing the sorediate taxa M. gomukhensis, M. olivaceoides 1 and the isidiate taxon M. poeltii; and the second containing the apotheciate taxa M. lobulata and M. subexasperata and the sorediate taxon M. olivaceoides 2. In both clades, these species are restricted to Asian highlands (i.e. Tibetan Plateau and northern India) and northern latitudes of western North America. This unique distribution pattern found in two distinct clades impacted by Pleistocene divergence provides circumstantial evidence of transcontinental dispersal from the Asian highlands to the New World through Beringia. Studies have demonstrated that Beringia has served as a biological corridor for the dispersal of a variety of organisms (e.g. [79-83]), and we hypothesize that glacial refugia in Beringia during the Pleistocene have likely played an important role in transcontinental dispersal and diversification for some lichenized fungi, including Melanohalea.

Demographic histories

During the Pleistocene, glacial advances affected the physical and biological environments of the Northern Hemisphere [84-86]. The more traditional methods for inferring population growth used in this study (i.e., Fu’s F, and Tajima’s D, and mismatch distribution) predict population expansion unambiguously for only M. elegantula, although Fu’s F and Tajima’s D were negative but not significant for other taxa, including lineages assessed using BSP analyses (M. exasperatula, M. multispora 3, M. olivacea, M. septentrionalis, and M. subolivacea 2). In our study mismatch distributions generally showed relatively unimodal patterns for taxa commonly producing vegetative diaspores (isidiate taxa), and mismatch distributions largely appeared multimodal for taxa producing ascomata (apotheciate taxa) (Additional file 3). Although a multimodal mismatch distribution is generally interpreted as evidence for population stability [87], BSP analyses of apotheciate taxa with multimodal mismatch distributions suggested Pleistocene population expansion of the four tested lineages. Multimodal distributions that fit sudden-expansion models can also be an indication of infraspecific structuring [88], a pattern consistent with the phylogenetic sub-structuring present in the many apotheciate taxa in the ITS topology (Additional file 2). Additional research investigating infraspecific population structure in fertile Melanohalea species and associated demographic histories will be essential to assess the reliability of the demographic histories inferred in this study.

BSP analyses suggest that population expansions predate the LGM in apotheciate taxa (M. multispora 3, M. olivacea, M. septentrionalis, and M. subolivacea 2), as well as those propagating via vegetative diaspores (M. elegantula and M. exasperatula). Our results indicate that at least some Melanohalea populations were able to expand during Pleistocene glacial cycles, albeit with apparently different demographic histories. Similar to other organisms [54,89,90], these data suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species. However, the relative impact of interglacial periods on population expansion in Melanohalea remains unclear.

With the exception of M. septentrionalis, apotheciate taxa evaluated here appear to have expanded over 700 Ka, suggesting a similar response to Pleistocene climatic factors among both broadly distributed (M. olivacea) and geographically restricted species (M. multispora 3 and M. subolivacea 2). The differences in the timing of population expansions seen between M. septentrionalis and other apotheciate species may be related to ecology as the treeline migrated north [91]. Melanohalea septentrionalis is essentially a boreal species, while outpost populations of M. olivacea south of the boreal zone are markedly more common than in M. septentrionalis[43]. Melanohalea multispora 3 and M. subolivacea 2 are essentially western North American species, commonly found in dry forests [92]. The contemporary distributions patterns and results from the BSP analyses suggest that climatic changes and associated shifts in vegetation in non-glaciated ecosystems during the Pleistocene had a significant impact on population demographics of fertile Melanohalea species.

Although it has been proposed that the contemporary distribution and abundance of common isidiate Melanohalea species is a result of recent expansion after the LGM and subsequent anthropogenic disturbances [48-51,57,76], BSP analyses of vegetatively propagating species (M. elegantula and M. exasperatula) suggest synchronous population expansion ca. 200 Ka, with 95% confidence interval providing strong evidence that this expansion predates the LGM (Figure 3). BSP analyses indicate that the beginning of the population expansion of these taxa coincides with the Saalian Pleniglacial cycle (ca. 182 Ka) and continued the expansion into the Weichselian glacial cycle (beginning ca. 110 Ka; [93]). Both isidiate taxa have broad ecological amplitude and are commonly found in formerly glaciated areas, as well as lower elevation sites, but are also found in other nutrient-rich habitats [43].

Nimis [58] argued that the distribution of Holarctic lichens is likely determined by general ecology rather than by their reproductive strategy alone. The latter argument supports the assumption that the distribution of Melanohalea within the Holarctic today widely reflects their biological constitution, rather than their geographic origin or dispersal strategy. However, contemporary distributions of isidiate taxa, M. elegantula and M. exasperatula, and apotheciate taxa, M. multispora 3, and M. subolivacea 2, largely overlap in western North America, and these species co-occur to varying degrees in some habitats. In spite of the fact that these species share similar distributions and habitats in some regions, BSP analyses suggest that lineages sharing reproductive strategies also share similar demographic histories, with population expansions in apotheciate lineages largely predating expansions in isidiate lineages. Whether these similarities are due to reproductive strategy-dependent fitness advantages, availability of the appropriate photosynthetic partner or fine-scale ecological niches remains unknown and merits future investigation.

Caution must be taken with the results from the BSP analyses presented here due to the potential limitations of using genetic data from the ribosomal tandem repeat and potential violations of the assumption of panmixia in lineages reproducing largely via vegetative diaspores. The selection of genetic loci for inferring demographic histories is not always straightforward (reviewed in [94]). Mitochondrial markers have been the locus of choice in skyline-plot analyses of animal populations [94-97]. However, within Parmeliaceae, mitochondrial markers generally have very low intraspecific variation [28,31,34,36]. Additionally, introns and paralogous copies of other commonly used mitochondrial markers have been found in fungi (reviewed in [98]). Although previous studies have found evidence of multiple ITS sequence types within individuals of several fungal species, including lichenized fungi [99-101], intragenomic ribosomal variation has not been reported in Parmeliaceae. Multiple studies have confirmed the utility of variable markers within the ribosomal cistron (i.e., ITS and IGS) for species and population-level studies in lichenized ascomycetes [5,8,37,102-108]. Furthermore, the ITS region has a high mutation rate and is the most variable marker commonly used in Parmeliaceae [7,8,28,36,37]. In this study we found that the substitution rate was over two times higher than that of the sampled protein-coding loci (MCM7, RPB1, and RPB2) and nearly five times higher than that of the sampled mitochondrial marker (mtSSU).

We are aware of potential limitations of using ribosomal markers to infer demographic histories, but we consider this a good “first pass” marker with which to investigate population demographics in lichenized fungi. While additional studies will be required to confirm the utility of the ITS marker for inferring demographic histories of lichenized ascomycetes using Bayesian skyline analyses, our study provides a valuable “hypothesis generating” approach for assessing historical factors driving population dynamics in lichen-forming fungi.


In this study we provide a comparison of dated divergence time estimates between concatenated gene trees and a calibrated multilocus species tree for the lichen-forming fungal genus Melanohalea. Our results showed a consistent pattern of more recent divergence times estimated from the coalescent-based species tree approach, relative to the concatenated gene tree analysis. Although earlier divergence times were estimated across the species tree topology, relative to dates estimated from the concatenated gene tree, we found the greatest disparity between methods at deeper nodes. This discrepancy clearly has important consequences for understanding factors driving diversification, extinctions, and biogeographic patterns. Given the results of this study, we propose that additional care must be taken when attempting to place divergence events within a temporal context by including estimates from calibrated coalescent-based species tree methods. Results from both divergence dating methods implemented herein indicate that most diversification within Melanohalea occurred during the Miocene and Pliocene, and was likely the result of habitat changes due to mountain uplift, major climatic changes (aridification) and associated shifts in vegetation patterns. This timeframe was also a period of significant mountain uplift, including major phases in the Andes [109], Alps [110], and Himalayas [111]. However, the temporal concordance of these diversification events with major vegetation changes in other geographical regions suggests that a more global environmental driver likely contributed to diversification in lichen-forming fungi in Parmeliaceae [12]. However, Pleistocene glacial cycles appear to have played a limited role in the diversification of some taxa in the Asian highlands and northern latitudes of North America and morphologically two cryptic lineages in western North America and Europe. In addition, we provide evidence for population expansions of common Melanohalea species during the Pleistocene. Our data provide evidence that the abundance of M. elegantula and M. exasperatula is likely not the result of anthropogenic factors. Rather, the population expansion of these common species coincides will the beginning of the Saalian glacial period. Although the timing of population expansions of distinct reproductive strategies predates the LGM, lineages commonly reproducing with vegetative diaspores expanded more recently (ca. 200 Ka) than most strictly sexually reproducing species (> 750 Ka). These data suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species.


Taxon sampling

In the present study we obtained sequence data from a total of 487 Melanohalea specimens (Additional file 1). Overall, 18 of the 22 species listed for the genus in Index Fungorum ( webcite) and seven previously unrecognized species-level lineages [59], awaiting formal description, were represented in this study. In this study nomenclature of undescribed species-level lineages follows arbitrary names in [59]. In order to infer changes in effective population sizes of common Melanohalea species over time, our sampling focused on multiple populations throughout the known distributions of four sexually reproducing species and two taxa commonly producing vegetative diaspores. For the sexually reproducing taxa, we selected two species with distributions restricted to western North America, M. multispora 3 (26 specimens) and M. subolivacea 2 (77); and two species with broader geographic distributions, M. olivacea (35) and M. septentrionalis (31). We sampled the broadly distributed isidiate taxa M. elegantula (100) and M. exasperatula (105) across their known distributions, representing common taxa reproducing vegetatively. Melanelixia californica was used as the outgroup for phylogenetic analyses [12,28,39].

DNA isolation, polymerase chain reaction (PCR) amplification and sequencing

Total genomic DNA was extracted from a small section of thallus material using the Prepease DNA Isolation Kit (USB, Cleveland, Ohio, USA), following the plant leaf extraction protocol. For all samples we sequenced the internal transcribed spacer region (ITS, ~500 bp). For a subset of these samples (Additional file 4), representing sampled genetic diversity identified from ITS sequence data, we amplified fragments from the nuclear ribosomal large subunit (LSU, ~530 bp), mitochondrial small subunit (mtSSU, ~800 bp), and three low-copy protein coding markers, RPB1 (~775 bp), RPB2 (~716 bp) and MCM7 (~514 bp). We developed new taxon-specific primers for the ITS, RPB2, and RPB1 in order to improve PCR specificity and efficiency in some cases. All primers used to amplify and sequence loci used in this study are provided in Table 5. PCR amplifications were conducted in 25 μL reactions. In some cases where standard PCR failed to amplify target loci, we used Ready-To-Go PCR Beads (GE Healthcare) following manufacturer’s recommendations with markedly improved success. PCR cycling parameters for amplifying the ITS and LSU loci followed [31]; cycling parameters for the RPB1, RPB2, and MCM7 fragments followed the MCM7 protocol described by [112]. PCR products were quantified on 1% agarose gel and stained with ethidium bromide. PCR fragments were cleaned using either the PrepEase PCR Purification Kit (USB, Cleveland, Ohio, USA) or ExoSAP-IT (USB, Cleveland, Ohio, USA), following manufacturers’ instructions. Complementary strands were sequenced from cleaned PCR products using the same primers used for amplifications. Sequencing reactions were performed using BigDye v3.1 (Applied Biosystems, Foster City, CA, USA). Products were run on an ABI 3730 automated sequencer according to recommended protocols (Applied Biosystems) at the Pritzker Laboratory for Molecular Systematics at the Field Museum, Chicago, IL, USA.

Additional file 4. Selected specimens representing sampled genetic diversity and GenBank accession numbers for the six sampled loci: nuclear ribosomal internal transcribed spacer region (ITS) and large subunit (nuLSU), mitochondrial small subunit (mtSSU), and protein-coding makers MCM 7, RPB 1, and RPB 2.

Format: DOCX Size: 35KB Download fileOpen Data

Table 5. All primers used for PCR amplification and sequencing in this study

Sequence alignment

We assembled and edited sequences using the program Sequencher v3.1.1 (Gene Codes Corporation, Ann Arbor, MI) and Geneious v5.4 [121]. Sequence identity was confirmed with ‘megaBLAST’ search in GenBank [122]. Sequences were aligned using the program MAFFT v6. For the LSU, MCM7, RPB1, and RPB2 markers we implemented the G-INS-I alignment algorithm, ‘1PAM / K = 2’ scoring matrix, and offset value = 0.8, with the remaining parameters set to default values; for the ITS we used the same G-INS-I alignment algorithm, ‘1PAM / K = 2’ and scoring matrix, and offset value = 0.1; and for the mtSSU we used the E-INS-I alignment algorithm, ‘20PAM / K = 2’ scoring matrix, and offset value = 0.0. For the mtSSU MAFFT alignment we also used the program Gblocks v0.91b [123,124] to remove regions of alignment uncertainty, using options for a “less stringent” selection on the Gblocks web server ( webcite).

Molecular dating

In this study we provide estimates of divergence events for species-level lineages presented in [59]. In order to obtain divergence date estimates for Melanohalea, we used two Bayesian approaches implemented in BEAST v1.6.1 [23]. We estimated divergence times under a gene-tree framework from the concatenated data matrix (i.e., ITS, nuLSU, mtSSU, MCM7, RPB1, and RPB2) and a multi-locus species tree approach.

For the divergence estimates in the concatenated gene tree approach we used a user-specified chronogram as the starting tree, rather than a randomly generated tree (e.g. [12]). To generate the starting topology, we conducted a ML analysis of the six-locus dataset (ITS, nuLSU, mtSSU, MCM7, RPB1, and RPB2) using locus-specific model partitions in RAxML v7.2.7 [125,126], and all loci were treated as separate partitions. We used the GTRGAMMA model, which includes a parameter (Γ) for rate heterogeneity among sites and chose not to include a parameter for estimating the proportion of invariable sites [125,126]. A search combining 200 separate maximum likelihood searches to find the optimal tree was conducted. The ML topology obtained from the RAxML analysis was converted to an ultrametric tree using nonparametric rate smoothing (NPRS) implemented in TreeEdit v.10a10 [127] for the BEAST analyses. In BEAST, the molecular dataset was analyzed with unlinked substitutions models across the loci, a birth-death model was used as a prior for the node heights, and a relaxed clock model (uncorrelated lognormal) for each partition. In the absence of relevant fossil evidence for Melanohalea, we used molecular evolution rates for RPB1, nuLSU, and mtSSU loci recently reported for Parmeliaceae [12]. While we used the Melanohalea-specific rate reported for the RPB1 locus (1.69 × 10-9 s/s/y), Melanohalea-specific rates were not available for the nuLSU and mtSSU loci [12]. Therefore, we used rates estimated for Parmeliaceae (Protoparmelia excluded) for the nuLSU (0.70 × 10-9 s/s/y) and mtSSU (0.69 × 10-9 s/s/y) to estimate the time to the most recent common ancestor (mrca) for all clades [12]. Substitution rates for the ITS, MCM7, and RPB2 markers were co-estimated along the run under a uniform prior, relative to the previously published rates for the nuLSU, mtSSU, and RPB1 loci. The rate in each branch was drawn from either an exponential or lognormal distribution. Two independent analyses were run for 50 million generations and parameter values were sampled every 1000 generations. The output from each analysis was visualized using Tracer v1.5 [128] to assess convergence and effective sampling size (ESS) of sampling parameters. Convergence was also assessed in AWTY [129] to ensure that standard deviations of split frequencies between runs approached zero and visualize split probabilities. Based on these results the first 12.5 million generations were discarded as burn-in, and remaining samples were summarized as a maximum clade credibility tree with mean divergence times and HPD intervals of age estimates using TreeAnnotator v.1.6.1 [130].

Gene trees, including phylogenies estimated from concatenated sequence data, can overestimate divergence times because they do not correct for genetic divergence that predates speciation [18]. Species tree methods incorporating the process of gene lineage coalescence likely provide a more biologically realistic framework for dating divergence events [18]. For comparison to concatenated-based divergence estimates, we thus also estimated divergence dates using a coalescent-based species-tree approach. We used the hierarchical Bayesian model *BEAST, implemented in program BEAST v1.6.1 [131], to estimate a species tree from lineages defined in [59]. *BEAST estimates the species tree directly from the sequence data, and incorporates the coalescent process, uncertainty associated with gene trees, and nucleotide substitution model parameters [131]. Models of DNA sequence evolution for each marker were selected with the program jModeltest v0.1 [132], using the AIC criterion. Implementing an uncorrelated relaxed lognormal clock [133], we selected a birth-death model for the species tree prior; the population size model was set to piecewise linear and constant root. Molecular evolution rates for the nuLSU, mtSSU, and RPB1 loci were as defined above. Default values were used for remaining priors. Two independent MCMC analyses were run for a total of 100 million generations (sampling every 2000 steps and excluding the first 25 million generations of each run as burn-in). Convergence was assessed by ensuring that standard deviations of split frequencies between runs approached zero, visualizing split probabilities in AWTY [129], and comparing summarized tree topologies from separate runs. We assessed the effective sample sizes (ESS) of parameters of interest to ensure sample sizes were all above 200. Posterior probabilities of nodes were computed from sampled trees after burn-in.

Molecular Diversity and population demographics

We used the program DnaSP v4.50 [134] to calculate estimates of genetic diversity (including number of haplotypes, h; haplotypic diversity, Hd; number of polymorphic sites, S; and nucleotide diversity, π) from ITS sequence data for each species. To detect possible departures from a constant population size that could be interpreted as a result of a past demographic expansion, we calculated the Fu’s F statistic [135] and Tajima’s D[136] for each species represented by ≥ 9 individuals. Significant and negative Tajima’s D and Fu’s F statistic values are indicative of possible population expansion, and positive values of these sample statistics provide evidence of a recently bottlenecked population or diversifying selection. These statistics were calculated in DnaSP, and significance was determined using the coalescent process implemented in DnaSP (1000 replicates). Mismatch frequency histograms were also plotted for each species represented by ≥ 9 individuals in DnaSP to determine whether clades exhibited evidence of spatial range expansion or a stationary population history. A smooth bell-shaped curve suggests either population expansion or spatial range expansion, whereas a multimodal curve represents a long history in situ. A raggedness index (RI, sum of the squared difference between neighboring peaks) and sum of squared deviations (SSD) between observed and expected mismatch were calculated in Arlequin v3.5 [137] using a parametric bootstrap approach (1000 replicates) to test for significance.

Although Tajima’s D, Fu’s F and mismatch distributions can provide insight into whether or not populations have undergone changes in population size, they do not provide information about the shape of the population growth over time. Here changes in species-specific demographic histories over time were assessed using Bayesian skyline analysis [138] in BEAST v1.6.1 [23] using the ribosomal ITS marker. A skyline plot is a model of population size fitting a wide range of demographic histories, and if a molecular clock rate is known for the locus in question, the model can be used to put demographic events into a historical context [138]. The Bayesian skyline model uses standard Markov chain Monte Carlo (MCMC) sampling procedures to estimate a posterior distribution of effective population size through time from gene sequences, given a model of sequence evolution. We selected six lineages characterized in [59] from which we have obtained sequences from specimens collected throughout their known distributions to minimize the violation of several simplifying assumptions [94,139]. Sexually reproducing taxa were represented by M. multispora 3 (26 specimens), M. olivacea (35), M. septentrionalis (31), and M. subolivacea 2 (77); and isidiate taxa included M. elegantula (100 specimens) and M. exasperatula (105). Although taxa that have been thought to reproduce primarily via asexual diaspores may violate the assumption of panmictic populations, sexual reproductive structures (apothecia) are occasionally found on some specimens of M. elegantula and M. exasperatula and may play an important role in population-wide gene flow [140]. Furthermore, recent data suggests that reproductive strategies in lichen-forming fungi traditionally considered to reproduce vegetatively are more complex than previously recognized [7,141].

Because systematic rate heterogeneity is not expected in intraspecific data, we used a strict molecular clock model with a fixed rate of 3.30 × 10−9 substitutions/site/year, obtained from the molecular dating analysis described above. Two independent analyses were run for 50 million Monte-Carlo-Markov-Chain generations, sampling parameter values every 1000 generations, using the Bayesian Skyline tree model allowing five discreet changes in population size with linear growth between population size change events was applied with a UPGMA generated tree as the starting point. We used the GTR + G + I model of nucleotide substitution without partitioning the ITS1, 5.8S, and ITS2 regions. We combined two independent runs and all ESS were > 200. We used Tracer v1.5.1 to analyze combined runs for each species and generate skyline plots.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HTL conceived the research and with TLE obtained funding. HTL, SDL, and TLE contributed with the conceptual development of the work. SDL, PKD, TLE, and HTL collected samples in the field. SDL and PKD performed the laboratory work. SDL performed data analyses and drafted the manuscript. All authors read and approved of the final version of the manuscript.


We are indebted to various colleagues for providing valuable material and field assistance, notably: Mónica Adler, Andre Aptroot, E. Berlin, Jarle Bjerke, Curtis Björk, Brian Coppins, Ana Crespo, Yuan-Da Du, Lidia Ferraro, Lawrence Glacy, Trevor Goward, E. Hansen, Ann Henson, Richard Harris, Jason Hollinger, Chris and Donna Howell, Martin Hutten, Charlotte Kellison, Kerry Knudsen, James Lendemer, Ede Leppik, Bruce McCune, Tom Nash, Sergio Pérez-Ortega, Rosa Pérez-Pérez, Tiina Randlane, Uve Rauist, Qiang Ren, Lauri Saa, Ghassan Sheika, Toby Spribille, Larry St. Clair, Ken Sweat, Arne Thell, Erin Tripp, James Walton, Hai-Ying Wang, Li-song Wang, and Zun-Tian Zhao. Kevin Feldheim, Sergio Garcia, Warren Chatwin, Jason Kluga, Nicolas Koutsoubelis, and Jessica Allen made valuable contributions in the lab; and Matthew Nelson provided valuable discussion. We also thank anonymous reviewers who provided valuable comments that improved this manuscript. Support by the National Science Foundation is gratefully acknowledged (“Hidden diversity in parmelioid lichens”, DEB-0949147), and PKD thanks the Ministry of Science and Technology, Spain for Ramon y Cajal start up grant (RYC02007-01576).


  1. Feuerer T, Hawksworth D: Biodiversity of lichens, including a world-wide analysis of checklist data based on Takhtajan’s floristic regions.

    Biodivers Conserv 2007, 16:85-98. Publisher Full Text OpenURL

  2. Lucking R, Rivas Plata E, Chaves JL, Umana L, Sipman HJM: How many tropical lichens are there …really?

    Biblioth Lichenol 2009, 100:399-418. OpenURL

  3. Printzen C: Lichen Systematics: The Role of Morphological and Molecular Data to Reconstruct Phylogenetic Relationships. In Progress in Botany 71. 71st edition. Berlin Heidelberg: Springer; 2009:233-275. OpenURL

  4. Rivas Plata E, Lumbsch HT: Parallel evolution and phenotypic divergence in lichenized fungi: A case study in the lichen-forming fungal family Graphidaceae (Ascomycota: Lecanoromycetes: Ostropales).

    Mol Phylogenet Evol 2011, 61:45-63. PubMed Abstract | Publisher Full Text OpenURL

  5. Leavitt SD, Fankhauser JD, Leavitt DH, Porter LD, Johnson LA, St. Clair LL: Complex patterns of speciation in cosmopolitan “rock posy” lichens – Discovering and delimiting cryptic fungal species in the lichen-forming Rhizoplaca melanophthalma species-complex (Lecanoraceae, Ascomycota).

    Mol Phylogenet Evol 2011, 59:587-602. PubMed Abstract | Publisher Full Text OpenURL

  6. Lumbsch HT, Leavitt SD: Goodbye morphology? A paradigm shift in the delimitation of species in lichenized fungi.

    Fungal Divers 2011, 50:59-72. Publisher Full Text OpenURL

  7. Leavitt SD, Johnson LA, Goward T, St. Clair LL: Species delimitation in taxonomically difficult lichen-forming fungi: An example from morphologically and chemically diverse Xanthoparmelia (Parmeliaceae) in North America.

    Mol Phylogenet Evol 2011, 60:317-332. PubMed Abstract | Publisher Full Text OpenURL

  8. Kroken S, Taylor JW: A Gene Genealogical Approach to Recognize Phylogenetic Species Boundaries in the Lichenized Fungus Letharia.

    Mycologia 2001, 93:38-53. Publisher Full Text OpenURL

  9. Crespo A, Pérez-Ortega S: Cryptic species and species pairs in lichens: A discussion on the relationship between molecular phylogenies and morphological characters.

    An Jard Bot Madr 2009, 66(S1):71-81. Publisher Full Text OpenURL

  10. Divakar PK, Figueras G, Hladun N, Crespo A: Molecular phylogenetic studies reveal an undescribed species within the North American concept of Melanelixia glabra (Parmeliaceae).

    Fungal Divers 2010, 42:47-55. Publisher Full Text OpenURL

  11. Molina M, Del-Prado R, Divakar P, Sánchez-Mata D, Crespo A: Another example of cryptic diversity in lichen-forming fungi: the new species Parmelia mayi (Ascomycota: Parmeliaceae).

    Org Divers & Evol 2011, 11:331-342. PubMed Abstract | Publisher Full Text OpenURL

  12. Amo De Paz G, Cubas P, Divakar PK, Lumbsch HT, Crespo A: Origin and Diversification of Major Clades in Parmelioid Lichens (Parmeliaceae, Ascomycota) during the Paleogene Inferred by Bayesian Analysis.

    PLoS One 2011, 6:e28161. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Otálora MAG, Martínez I, Aragón G, Molina MC: Phylogeography and divergence date estimates of a lichen species complex with a disjunct distribution pattern.

    Am J Bot 2010, 97:216-223. PubMed Abstract | Publisher Full Text OpenURL

  14. Sérusiaux E, Villarreal AJC, Wheeler T, Goffinet B: Recent origin, active speciation and dispersal for the lichen genus Nephroma (Peltigerales) in Macaronesia.

    J Biogeogr 2011, 38:1138-1151. Publisher Full Text OpenURL

  15. Lücking R, Huhndorf S, Pfister DH, Plata ER, Lumbsch HT: Fungi evolved right on track.

    Mycologia 2009, 101:810-822. PubMed Abstract | Publisher Full Text OpenURL

  16. Taylor JW, Berbee ML: Dating divergences in the Fungal Tree of Life: review and new analyses.

    Mycologia 2006, 98:838-849. PubMed Abstract | Publisher Full Text OpenURL

  17. Berbee ML, Taylor JW: Dating the molecular clock in fungi – how close are we?

    Fungal Biol Rev 2010, 24:1-16. Publisher Full Text OpenURL

  18. McCormack JE, Heled J, Delaney KS, Peterson AT, Knowles LL: Calibrating divergence times on species trees versus gene trees: Implications for speciation history of Aphelocoma jays.

    Evolution 2011, 65:184-202. PubMed Abstract | Publisher Full Text OpenURL

  19. Burbrink FT, Fontanella F, Alexander Pyron R, Guiher TJ, Jimenez C: Phylogeography across a continent: The evolutionary and demographic history of the North American racer (Serpentes: Colubridae: Coluber constrictor).

    Mol Phylogenet Evol 2008, 47:274-288. PubMed Abstract | Publisher Full Text OpenURL

  20. Brown J, Rest J, Garcia-Moreno J, Sorenson M, Mindell D: Strong mitochondrial DNA support for a Cretaceous origin of modern avian lineages.

    BMC Biol 2008, 6:6. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Arbogast BS, Edwards SV, Wakeley J, Beerli P, Slowinski JB: Estimating divergence times from molecular data on phylogenetic and population genetic timescales.

    Annu Rev Ecol Syst 2002, 33:707-740. Publisher Full Text OpenURL

  22. Graur D, Martin W: Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision.

    Trends Genet 2004, 20:80-86. PubMed Abstract | Publisher Full Text OpenURL

  23. Drummond A, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees.

    BMC Evoly Biol 2007, 7:214. BioMed Central Full Text OpenURL

  24. Yang Z, Rannala B: Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft Bounds.

    Mol Biol Evol 2006, 23:212-226. PubMed Abstract | Publisher Full Text OpenURL

  25. Thorne JL, Kishino H: Divergence time and evolutionary rate estimation with multilocus data.

    Syst Biol 2002, 51:689-702. PubMed Abstract | Publisher Full Text OpenURL

  26. Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock.

    Bioinformatics 2003, 19:301-302. PubMed Abstract | Publisher Full Text OpenURL

  27. Lumbsch HT, Hipp A, Divakar P, Blanco O, Crespo A: Accelerated evolutionary rates in tropical and oceanic parmelioid lichens (Ascomycota).

    BMC Evol Biol 2008, 8:e257. BioMed Central Full Text OpenURL

  28. Crespo A, Kauff F, Divakar PK, del Prado R, Perez-Ortega S, Amo DePaz G, Ferencova Z, Blanco O, Roca-Valiente B, Nunez-Zapata J, et al.: Phylogenetic generic classification of parmelioid lichens (Parmeliaceae, Ascomycota) based on molecular, morphological and chemical evidence.

    Taxon 2010, 59:1735-1753. OpenURL

  29. Blanco O, Crespo A, Ree RH, Lumbsch HT: Major clades of parmelioid lichens (Parmeliaceae, Ascomycota) and the evolution of their morphological and chemical diversity.

    Mol Phylogenet Evol 2006, 39:52-69. PubMed Abstract | Publisher Full Text OpenURL

  30. Del-Prado R, Ferencová Z, Armas-Crespo V, Amo De Paz G, Cubas P, Crespo A: The arachiform vacuolar body: an overlooked shared character in the ascospores of a large monophyletic group within Parmeliaceae (Xanthoparmelia clade, Lecanorales).

    Mycol Res 2007, 111:685-692. PubMed Abstract | Publisher Full Text OpenURL

  31. Crespo A, Lumbsch HT, Mattsson J-E, 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.

    Mol Phylogenet Evol 2007, 44:812-824. PubMed Abstract | Publisher Full Text OpenURL

  32. Divakar PK, Crespo A, Blanco O, Lumbsch HT: Phylogenetic significance of morphological characters in the tropical Hypotrachyna clade of parmelioid lichens (Parmeliaceae, Ascomycota).

    Mol Phylogenet Evol 2006, 40:448-458. PubMed Abstract | Publisher Full Text OpenURL

  33. Divakar PK, Lumbsch HT, Ferencova Z, Del Prado R, Crespo A: Remototrachyna, a newly recognized tropical lineage of lichens in the Hypotrachyna clade (Parmeliaceae, Ascomycota), originated in the Indian subcontinent.

    Am J Bot 2010, 97:579-590. PubMed Abstract | Publisher Full Text OpenURL

  34. Geml J, Kauff F, Brochmann C, Taylor DL: Surviving climate changes: high genetic diversity and transoceanic gene flow in two arctic–alpine lichens, Flavocetraria cucullata and F. nivalis (Parmeliaceae, Ascomycota).

    J Biogeogr 2010, 37:1529-1542. OpenURL

  35. Crespo A, Lumbsch HT: Cryptic species in lichen-forming fungi.

    IMA Fungus 2010, 1:167-170. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Del-Prado R, Cubas P, Lumbsch HT, Divakar PK, Blanco O, de Paz GA, Molina MC, Crespo A: Genetic distances within and among species in monophyletic lineages of Parmeliaceae (Ascomycota) as a tool for taxon delimitation.

    Mol Phylogenet Evol 2010, 56:125-133. PubMed Abstract | Publisher Full Text OpenURL

  37. Fernández-Mendoza F, Domaschke S, GarcÍA MA, Jordan P, Martin MP, Printzen C: Population structure of mycobionts and photobionts of the widespread lichen Cetraria aculeata.

    Mol Ecol 2011, 20:1208-1232. PubMed Abstract | Publisher Full Text OpenURL

  38. Zhao Z-T, Meng F-G, Li H-M, Wang H-Y: A new species of Melanohalea (Parmeliaceae) from the Tibetan Plateau.

    Mycotaxon 2009, 108:347-352. Publisher Full Text OpenURL

  39. Blanco O, Crespo A, Divakar PK, Esslinger TL, Hawksworth DL, Lumbsch TH: Melanelixia and Melanohalea, two new genera segregated from Melanelia (Parmeliaceae) based on molecular and morphological data.

    Mycol Res 2004, 108:873-884. PubMed Abstract | Publisher Full Text OpenURL

  40. Divakar PK, Upreti DK, Elix JA: New species and new records in the lichen family Parmeliaceae (Ascomycotina) from India.

    Mycotaxon 2001, 80:356-362. OpenURL

  41. Esslinger TL: A chemosytematic revision of the brown Parmeliae.

    J Hattori Bot Lab 1977, 42:1-211. OpenURL

  42. Esslinger TL: A new species of Melanelia from Nepal.

    Mycotaxon 1987, 28:215-217. OpenURL

  43. Otte V, Esslinger TL, Litterski B: Global distribution of the European species of the lichen genus Melanelia Essl.

    J Biogeogr 2005, 32:1221-1241. Publisher Full Text OpenURL

  44. Sun L-Y, Meng F-G, Li H-M, Wang H-Y, Zhao Z-T: A new lichen, Melanohalea subexasperata (Parmeliaceae), from the Tibetan Plateau.

    Mycotaxon 2010, 111:65-69. Publisher Full Text OpenURL

  45. Galloway DJ, Jørgensen PM: Bartlettiella, a new lichen genus from New Zealand, with notes on a new species of Melanelia and a new chemodeme of Bryoria indonesica in New Zealand.

    New Zeal J Bot 1990, 28:5-12. Publisher Full Text OpenURL

  46. Esslinger TL: A new status for brown Parmeliae.

    Mycotaxon 1978., 7 OpenURL

  47. Purvis OW, James PW: Parmelia Ach., 1803. In The Lichen Flora of Great Britain and Ireland. Edited by Purvis OW, Coppins BJ, Hawksworth DL, James PW, Moore DM. London: Natural History Museum Publication; 1992:421-437. OpenURL

  48. Nimis PL: The lichens of Italy. Torino: Museo Regionale di Scienze Naturali, Monographia 12; 1993. OpenURL

  49. Seaward MRD: Lichen atlas of the British Isles. London: British Lichen Society; 1995. OpenURL

  50. Thell A: Forandringar i utbredningsbilden hos Parmelia elegantula och Parmelia laciniatula in Skåne och Blekinge.

    Graphis Scripta 1989, 2:156-160. OpenURL

  51. Otte V, Rätzel S, Kummer V, de Bruyn U: Bemerkenswerte Flechtenfunde aus Brandenburg VI.

    Verhandlungen des Botanischen Vereins von Berlin und Brandenburg 2001, 134:137-154. OpenURL

  52. Printzen C, Ekman S, Tønsberg T: Phylogeography of Cavernularia hultenii: Evidence of slow genetic drift in a widely disjunct lichen.

    Mo Ecol 2003, 12:1473-1486. Publisher Full Text OpenURL

  53. Moyle RG, Filardi CE, Smith CE, Diamond J: Explosive Pleistocene diversification and hemispheric expansion of a “great speciator”.

    Proc Natl Acad Sci 2009, 106:1863-1868. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Knowles LL: Tests of Pleistocene speciation in Montane grasshoppers (genus Melanoplus) from sky islands of western North America.

    Evolution 2000, 54:1337-1348. PubMed Abstract OpenURL

  55. Jaeger JR, Riddle BR, Bradford DF: Cryptic Neogene vicariance and Quaternary dispersal of the red-spotted toad (Bufo punctatus): insights on the evolution of North American warm desert biotas.

    Mol Ecol 2005, 14:3033-3048. PubMed Abstract | Publisher Full Text OpenURL

  56. Zachos JC, Dickens GR, Zeebe RE: An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics.

    Nature 2008, 451:279-283. PubMed Abstract | Publisher Full Text OpenURL

  57. Poelt J: Flechtenflora und Eiszeit in Europa.

    Phyton, Annales Tei Botanicae 1963, 10:206-215. OpenURL

  58. Nimis PL, Martellos S: On the ecology of sorediate lichens in Italy. In Lichenological contributions in honour of G. B. Feige. Bibliotheca Lichenologica 86. Edited by Jensen M. Berlin: J Cramer, Berlin & Stuttgart; 2003:393-406. OpenURL

  59. Leavitt SD, Esslinger TL, Spribille T, Divakar PK, Lumbsch HT: Multilocus phylogeny of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota): Insights on diversity, distributions, and a comparison of species tree and concatenated topologies.

    Mol Phylogenet Evol

    accepted manuscript


  60. Knowles LL, Carstens BC: Estimating a geographically explicit model of population divergence.

    Evolution 2007, 61:477-493. PubMed Abstract | Publisher Full Text OpenURL

  61. Liu L, Yu L, Kubatko L, Pearl DK, Edwards SV: Coalescent methods for estimating phylogenetic trees.

    Mol Phylogenet Evol 2009, 53:320-328. PubMed Abstract | Publisher Full Text OpenURL

  62. Edwards SV: Is a new and general theory of molecular systematics emerging?

    Evolution 2009, 63:1-19. PubMed Abstract | Publisher Full Text OpenURL

  63. Knowles LL: Estimating Species Trees: Methods of Phylogenetic Analysis When There Is Incongruence cross Genes.

    Syst Biol 2009, 58:463-467.


    PubMed Abstract | Publisher Full Text OpenURL

  64. Knowles LL, Carstens BC: Delimiting Species without Monophyletic Gene Trees.

    Syst Biol 2007, 56:887-895. PubMed Abstract | Publisher Full Text OpenURL

  65. Edwards SV, Liu L, Pearl DK: High-Resolution Species Trees without Concatenation.

    Proc of the Natl Acad Sci USA 2007, 104:5936-5941. Publisher Full Text OpenURL

  66. Edwards S, Beerli P: Perspective: Gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies.

    Evolution 2000, 54:1839-1854. PubMed Abstract OpenURL

  67. Van Devender T: The Sonoran desert tortoise. Tuscon: University of Arizona Press and the Arizona-Sonoran Desert Museum; 2002. OpenURL

  68. Ramstein G, Fluteau F, Besse J, Joussaume S: Effect of orogeny, plate motion and land-sea distribution on Eurasian climate change over the past 30 million years.

    Nature 1997, 386:788-795. Publisher Full Text OpenURL

  69. Zachos J, Pagani M, Sloan L, Thomas E, Billups K: Trends in rythms, and aberrations in global climate 65 Ma to present.

    Science 2001, 292:686-693. PubMed Abstract | Publisher Full Text OpenURL

  70. Berry EW: The Miocene flora of Idaho.

    Proc Natl Acad Sci 1932, 18:289-291. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  71. Axelrod D: Age and origin of Sonoran desert vegetation.

    Occasional Papers, California Acedemy of Sciences 1979, 132:1-74. OpenURL

  72. Dunai TJ, López GAG, Juez-Larré J: Oligocene–Miocene age of aridity in the Atacama Desert revealed by exposure dating of erosion-sensitive landforms.

    Geology 2005, 33:321-324. Publisher Full Text OpenURL

  73. Janis CM: Tertiary mammal evolution in the context of changing climates, vegetation, and tectonic events.

    Annu Rev Ecol Syst 1993, 24:467-500. Publisher Full Text OpenURL

  74. Jacobs BF, Kingston JD, Jacobs LL: The origins of grass-dominated ecosystems.

    Ann Mo Bot Gard 1999, 86:590-643. Publisher Full Text OpenURL

  75. Pagani M, Freeman KH, Arthur MA: Late Miocene Atmospheric CO2 Concentrations and the Expansion of C4 Grasses.

    Science 1999, 285:876-879. PubMed Abstract | Publisher Full Text OpenURL

  76. Calvelo S, Stocker-Wörgötter E, Liberatore S, Elix JA: Protousnea (Parmeliaceae, Ascomycota), a Genus Endemic to Southern South America.

    Bryologist 2005, 108:1-15. Publisher Full Text OpenURL

  77. Endler JA: Natural selection in the wild. Princeton: Princeton University Press; 1986. OpenURL

  78. Schluter D: The ecology of adaptive radiations. Oxford: Oxford University Press; 2000. OpenURL

  79. Anderson PM, Lozhkin AV, Solomatkina TB, Brown TA: Paleoclimatic implications of glacial and postglacial refugia for Pinus pumila in western Beringia.

    Quatern Res 2010, 73:269-276. Publisher Full Text OpenURL

  80. Vila R, Bell CD, Macniven R, Goldman-Huertas B, Ree RH, Marshall CR, Bálint Z, Johnson K, Benyamini D, Pierce NE: Phylogeny and palaeoecology of Polyommatus blue butterflies show Beringia was a climate-regulated gateway to the New World.

    Proc R Soc Edinb Biol 2011, 278:2737-2744. OpenURL

  81. Galbreath KE, Cook JA, Eddingsaas AA, DeChaine EG: Diversity and demography in Beringia: Multilocus tests of paleodistribution models reveal the complex history of arctic ground squirrels.

    Evolution 2011, 65:1879-1896. PubMed Abstract | Publisher Full Text OpenURL

  82. Weber WA: The Middle Asian Element in the Southern Rocky Mountain Flora of the western United States: a critical biogeographical review.

    J Biogeogr 2003, 30:649-685. Publisher Full Text OpenURL

  83. Thomson JW: The Distribution of Arctic Lichens and Thoughts Concerning Their Origin.

    Lichenologist 1995, 27:411-416. OpenURL

  84. Ruddiman WF, Kutzbach JE: Forcing of Late Cenozoic Northern Hemisphere Climate by Plateau Uplift in Southern Asia and the American West.

    J Geophys Res Atmospheres 1989, 94:18409-18427. Publisher Full Text OpenURL

  85. Ruddiman WF, Raymo M, McIntyre A: Matuyama 41,000-year cycles: North Atlantic Ocean and northern hemisphere ice sheets.

    Earth Planet Sci Lett 1986, 80:117-129. Publisher Full Text OpenURL

  86. Meyers SR, Hinnov LA: Northern Hemisphere glaciation and the evolution of Plio-Pleistocene climate noise.

    Paleoceanography 2010, 25:PA3207. OpenURL

  87. Harpending HC, Batzer MA, Gurven M, Jorde LB, Rogers AR, Sherry ST: Genetic traces of ancient demography.

    Proc Nat Acad Sci 1998, 95:1961-1967. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  88. Castoe TA, Spencer CL, Parkinson CL: Phylogeographic structure and historical demography of the western diamondback rattlesnake (Crotalus atrox): A perspective on North American desert biogeography.

    Mol Phylogenet Evol 2007, 42:193-212. PubMed Abstract | Publisher Full Text OpenURL

  89. DeChaine EG, Martin AP: Using coalescent simulations to test the impact of Quaternary climate cycles on divergence in an alpine plant-insect association.

    Evolution 2006, 60:1004-1013. PubMed Abstract OpenURL

  90. Galbreath KE, Hafner DJ, Zamudio KR: When cold is better; Climate-driven elevation shifts yield complex patters of diversification and demography in an alpine specialist (American Pika, Ochotona princeps).

    Evolution 2009, 63:2848-2863. PubMed Abstract | Publisher Full Text OpenURL

  91. Clark JS, Fastie C, Hurtt G: Reid's paradox of rapid plant migration-dispersal theory and interpretation of paleoenvironmental records.

    Bioscience 1998, 48:13-24. Publisher Full Text OpenURL

  92. Ahti T: Parmelia olivacea and the allied non-isidiate and non-sorediate corticolous lichens in the northern hemisphere.

    Acta Botanica Fennica 1966, 70:1-68. OpenURL

  93. George K: Saalian supercycle, Mindel/Riss interglacial and Milankovitch's dating.

    Quat Sci Rev 2005, 24:1573-1583. Publisher Full Text OpenURL

  94. Ho SYW, Shapiro B: Skyline-plot methods for estimating demographic history from nucleotide sequences.

    Mol Ecol Res 2011, 11:423-434. Publisher Full Text OpenURL

  95. Buckley TR, Marske K, Attanayake D: Phylogeography and ecological niche modelling of the New Zealand stick insect Clitarchus hookeri (White) support survival in multiple coastal refugia.

    J Biogeogr 2010, 37:682-695. Publisher Full Text OpenURL

  96. Gompert Z, Fordyce JA, Forister ML, Nice CC: Recent colonization and radiation of North American Lycaeides (Plebejus) inferred from mtDNA.

    Mol Phylogenet Evol 2008, 48:481-490. PubMed Abstract | Publisher Full Text OpenURL

  97. Kumar S, Bellis C, Zlojutro M, Melton P, Blangero J, Curran J: Large scale mitochondrial sequencing in Mexican Americans suggests a reappraisal of Native American origins.

    BMC Evol Biol 2011, 11:293. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  98. Seifert KA: Progress towards DNA barcoding of fungi.

    Mol Ecol Res 2009, 9:83-89. OpenURL

  99. Simon UK, Weiß M: Intragenomic variation of fungal ribosomal genes is higher than previously thought.

    Mol Biol Evol 2008, 25:2251-2254. PubMed Abstract | Publisher Full Text OpenURL

  100. Simon DM, Hummel CL, Sheeley SL, Bhattacharya D: Heterogeneity of intron presence or absence in rDNA genes of the lichen species Physcia aipolia and P. stellaris.

    Curr Genet 2005, 47:389-399. PubMed Abstract | Publisher Full Text OpenURL

  101. Kelly LJ, Hollingsworth PM, Coppins BJ, Ellis CJ, Harrold P, Tosh J, Yahr R: DNA barcoding of lichenized fungi demonstrates high identification success in a floristic context.

    New Phytol 2011, 191:288-300. PubMed Abstract | Publisher Full Text OpenURL

  102. Wirtz N, Printzen C, Lumbsch HT: The delimitation of Antarctic and bipolar species of neuropogonoid Usnea (Ascomycota, Lecanorales): a cohesion approach of species recognition for the Usnea perpusilla complex.

    Mycol Res 2008, 112:472-484. PubMed Abstract | Publisher Full Text OpenURL

  103. Brunauer G, Hager A, Grube M, Türk R, Stocker-Wörgötter E: Alterations in secondary metabolism of aposymbiotically grown mycobionts of Xanthoria elegans and cultured resynthesis stages.

    Plant Physiol Biochem 2007, 45:146-151. PubMed Abstract | Publisher Full Text OpenURL

  104. Buschbom J, Mueller GM: Testing "Species Pair" Hypotheses: Evolutionary Processes in the Lichen-Forming Species Complex Porpidia flavocoerulescens and Porpidia melinodes.

    Mol Biol Evol 2006, 23:574-586. PubMed Abstract | Publisher Full Text OpenURL

  105. O'Brien HE, Miadlikowska J, Lutzoni F: Assessing reproductive isolation in highly diverse communities of the lichen-forming fugnal genus Peltigera.

    Evolution 2009, 63:2076-2086. PubMed Abstract | Publisher Full Text OpenURL

  106. Lindblom L, Ekman S: Genetic variation and population differentiation in the lichen-forming ascomycete Xanthoria parietina on the island Storfosna, central Norway.

    Mol Ecol 2006, 15:1545-1559. PubMed Abstract | Publisher Full Text OpenURL

  107. Pino-Bodas R, Burgaz AR, Martin MP, Lumbsch HT: Species delimitations in the Cladonia cariosa group (Cladoniaceae, Ascomycota).

    Lichenologist 2012, 44:121-135. Publisher Full Text OpenURL

  108. Wirtz N, Printzen C, Lumbsch H: Using haplotype networks, estimation of gene flow and phenotypic characters to understand species delimitation in fungi of a predominantly Antarctic Usnea group (Ascomycota, Parmeliaceae).

    Org Divers Evol

    in press

    Publisher Full Text OpenURL

  109. Gregory-Wodzicki KM: Uplift history of the Central and Northern Andes: A review.

    Geol Soc Am Bull 2000, 112:1091-1105. Publisher Full Text OpenURL

  110. Hurford AJ: Cooling and uplift patterns in the Lepontine Alps South Central Switzerland and an age of vertical movement on the Insubric fault line.

    Contributions to Mineralogy and Petrology 1986, 92:413-427. Publisher Full Text OpenURL

  111. Copeland P, Harrison TM: Episodic rapid uplift in the Himalaya revealed by 40Ar/39Ar analysis of detrital K-feldspar and muscovite, Bengal fan.

    Geology 1990, 18:354-357. Publisher Full Text OpenURL

  112. Schmitt I, Crespo A, Divakar PK, Fankhauser JD, Herman-Sackett E, Kalb K, Nelsen MP, Rivas-Plata E, Shimp AD, Widhelm T, Lumbsch HT: New primers for promising single-copy genes in fungal phylogenies and systematics.

    Persoonia 2009, 23:35-40. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  113. Gardes M, Bruns TD: ITS primers with enhanced specificity for basidiomycetes - application to the identification of mycorrhizae and rusts.

    Mol Ecol Notes 1993, 2:113-118. OpenURL

  114. White TJ, Bruns TD, Lee S, Taylor J: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. San Diego: Academic Press; 1990:315-322. OpenURL

  115. Larena I, Salazar O, González V, Julián MC, Rubio V: Design of a primer for ribosomal DNA internal transcribed spacer with enhanced specificity for ascomycetes.

    J Biotechnology 1999, 75:187-194. Publisher Full Text OpenURL

  116. Mangold A, Martín MP, Kalb K, Lücking R, Lumbsch HT: Molecular data show that Topeliopsis (Ascomycota, Thelotremataceae) is polyphyletic.

    Lichenologist 2008, 40:39-46. OpenURL

  117. Zoller S, Lutzoni F, Scheidegger C: Genetic variation within and among populations of the threatened lichen Lobaria pulmonaria in Switzerland and implications for its conservation.

    Mol Ecol 1999, 8:2049-2059. PubMed Abstract | Publisher Full Text OpenURL

  118. Leavitt SD, Johnson L, St. Clair LL: Species delimitation and evolution in morphologically and chemically diverse communities of the lichen-forming genus Xanthoparmelia (Parmeliaceae, Ascomycota) in western North America.

    Am J Bot 2011, 98:175-188. PubMed Abstract | Publisher Full Text OpenURL

  119. Matheny PB, Liu YJ, Ammirati JF, Hall BD: Using RPB1 sequences to improve phylogenetic inference among mushrooms (Inocybe, Agaricales).

    Am J Bot 2002, 89:688-698. PubMed Abstract | Publisher Full Text OpenURL

  120. Liu YJ, Whelen S, Hall BD: Phylogenetic relationships among ascomycetes: evidence from an RNA polymerse II subunit.

    Mol Biol Evol 1999, 16:1799-1808. PubMed Abstract | Publisher Full Text OpenURL

  121. Drummond AJ, Ashton B, Cheung M, Heled J, Kearse M, Moir R, Stones-Havas S, Thierer T, Wilson A:

    Geneious ver. 5.3. 2010. webcite


  122. Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, DiCuccio M, Edgar R, Federhen S, et al.: Database resources of the National Center for Biotechnology Information.

    Nucleic Acids Res 2006., gkl1031 OpenURL

  123. Castresana J: Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis.

    Mol Biol Evol 2000, 17:540-552. PubMed Abstract | Publisher Full Text OpenURL

  124. Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments.

    Syst Biol 2007, 56:564-577. PubMed Abstract | Publisher Full Text OpenURL

  125. Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models.

    Bioinformatics 2006, 22:2688-2690. PubMed Abstract | Publisher Full Text OpenURL

  126. Stamatakis A, Hoover P, Rougemont J: A Rapid Bootstrap Algorithm for the RAxML Web Servers.

    Syst Biol 2008, 57:758-771. PubMed Abstract | Publisher Full Text OpenURL

  127. Rambaut A, Charleston M:

    Tree Edit. Phylogenetic Tree Editor v1.0 alpha 8. 2001.

    [ webcite] webcite


  128. Rambaut A, Drummond A:

    Tracer [computer program]. 2003. OpenURL

  129. Nylander JAA, Wilgenbusch JC, Warren DL, Swofford DL: AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics.

    Bioinformatics 2008, 24:581-583. PubMed Abstract | Publisher Full Text OpenURL

  130. Rambault A, Drummond AJ:

    TreeAnnotator. Version 1.6.1. 2009.

    Available at: http://beastbioedacuk/TreeAnnotator webcite


  131. Heled J, Drummond AJ: Bayesian Inference of Species Trees from Multilocus Data.

    Mol Biol Evol 2010, 27:570-580. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  132. Posada D: jModelTest: Phylogenetic Model Averaging.

    Mol Biol Evol 2008, 25:1253-1256. PubMed Abstract | Publisher Full Text OpenURL

  133. Drummond A, Ho S, Phillips M, Rambaut A: Relaxed phylogenetics and dating with confidence.

    PLoS Biol 2006, 4:e88. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  134. Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data.

    Bioinformatics 2009, 25:1451-1452. PubMed Abstract | Publisher Full Text OpenURL

  135. Fu Y-X: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selections.

    Genet 1997, 174:915-925. OpenURL

  136. Tajima F: Statistical method for testing the neurtral mutation hypothesis by DNA polymorphism.

    Genet 1989, 123:585-595. OpenURL

  137. Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis.

    Evol Bioinform 2005, 1:47-50. OpenURL

  138. Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian Coalescent Inference of Past Population Dynamics from Molecular Sequences.

    Mol Biol Evol 2005, 22:1185-1192. PubMed Abstract | Publisher Full Text OpenURL

  139. Donnelly P, Tavare S: Coalescents and genealogical structure under neutrality.

    Annu Rev Genet 1995, 29:401-421. PubMed Abstract | Publisher Full Text OpenURL

  140. Robinson H: Considerations on the evolution of lichens.

    Phytologia 1975, 32:407-413. OpenURL

  141. Buschbom J, Barker D: Evolutionary History of Vegetative Reproduction in Porpidia s.l. (Lichen-Forming Ascomycota).

    Syst Biol 2006, 55:471-484. PubMed Abstract | Publisher Full Text OpenURL