- Research article
- Open Access
Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea(Parmeliaceae, Ascomycota) and Pleistocene population expansions
© Leavitt et al.; licensee BioMed Central Ltd. 2012
- Received: 24 February 2012
- Accepted: 30 August 2012
- Published: 11 September 2012
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
- Divergence times
- 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  and an estimated diversity of greater than 28,000 species . 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 . 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 .
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 .
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 . Otte et al.  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 . 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  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  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 . 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 , 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 . 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.
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
N(# of unique haplotypes)
# of variable sites
# of PI sites
TIM1ef + I + G
TIM1ef + I + G
TPM2 + G
TPM1uf + I
TIM3ef + I
TIM2ef + I + G
TPM2uf + I + G
Age estimates on concatenated and species topologies
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 
Concatenated gene tree analysis
Multilocus species tree analysis
Rate 95% HPD
Rate 95% HPD
2.49 – 4.11
2.67 – 4.17
0.58 – 0.80
0.62 – 0.80
0.42 – 1.05
0.44 – 0.92
1.26 – 2.05
1.36 – 2.08
1.27 – 1.82
1.57 – 1.83
1.08 – 1.75
1.20 – 1.82
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
Stem origin (MYA) concatenated
Stem origin (MYA) coalescent
tmrca extant haplotype (MYA) concatenated
Melanelixia – Melanohalea split
43.5 (95% HPD = 26.9-59.7)
35.2 (95% HPD = 27.0-43.5)
Origin of Melanohalea
30.2 (95% HPD = 23.8- 38.1)
23.7 (95% HPD = 19.4-27.9)
2.7 (95% HPD = 1.0-4.4)
1.0 (95% HPD = 0.5-1.6)
M. elegantula Spain
2.7 (95% HPD = 1.0-4.4)
14.8 (95% HPD = 10.8-19.7)
9.2 (95% HPD = 6.2-12.4)
2.2 (95% HPD = 1.2-3.4)
M. exasperata Spain
5.8 (95% HPD = 3.8-8.6)
1.7 (95% HPD = 0.7-3.0
6.0 (95% HPD = 4.1-8.4)
1.0 (95% HPD = 0.5-1.5)
1.9 (95% HPD = 1.0-3.1)
0.9 (95% HPD = 0.0-1.8)
0.1 (95% HPD = 0.0-0.4)
6.2 (95% HPD = 4.2-8.5)
4.5 (95% HPD = 2.7-6.4)
0.2 (95% HPD = 0.0-0.3)
14.9 (95% HPD = 11.0-20.0)
11.0 (95% HPD = 8.4-13.8)
0.6 (95% HPD = 0.2-1.1)
17.2 (95% HPD = 12.8-22.5)
7.4 (95% HPD = 4.6-10.3)
1.2 (95% HPD = 0.6-2.0)
3.7 (95% HPD = 2.0-5.9)
2.3 (95% HPD = 0.9-3.7)
0.6 (95% HPD = 0.0-1.8)
M. multispora 1
15.9 (95% HPD = 11.7-20.5)
10.6 (95% HPD = 7.9-13.5)
0.9 (95% HPD = 0.5-1.6)
M. multispora 2
7.7 (95% HPD = 5.1-11.1)
5.1 (95% HPD = 2.8-7.3)
0.4 (95% HPD = 0.0-1.1)
M. multispora 3
7.7 (95% HPD = 5.1-11.1)
5.1 (95% HPD = 2.8-7.3)
2.7 (95% HPD = 1.7-3.9)
5.7 (95% HPD = 3.9-7.9)
4.2 (95% HPD = 2.7-5.8)
1.7 (95% HPD = 1.0-2.8)
M. olivaceoides 1
1.9 (95% HPD = 1.0-3.1)
1.7 (95% HPD = 0.7-2.6)
0.3 (95% HPD = 0.0-0.7)
M. olivaceoides 2
1.7 (95% HPD = 0.6-3.4)
1.3 (95% HPD = 0.1-2.4)
0.1 (95% HPD = 0.0-0.3)
2.6 (95% HPD = 1.5-4.0)
0.9 (95% HPD = 0.0-1.8)
1.0 (95% HPD = 0.2-2.3)
6.2 (95% HPD = 4.2-8.5)
4.5 (95% HPD = 2.7-6.4)
2.6 (95% HPD = 1.5-3.9)
8.0 (95% HPD = 5.5-10.6)
6.0 (95% HPD = 4.3-7.8)
1.3 (95% HPD = 0.6-2.3)
1.7 (95% HPD = 0.6-3.4)
1.3 (95% HPD = 0.1-2.4)
1.0 (95% HPD = 0.2-2.2)
M. subolivacea 1
3.9 (95% HPD = 2.7-5.6)
1.1 (95% HPD = 0.7-1.6)
1.6 (95% HPD = 0.9-2.6)
M. subolivacea 2
3.9 (95% HPD = 2.7-5.6)
1.1 (95% HPD = 0.7-1.6)
2.4 (95% HPD = 1.9-4.0)
12.5 (95% HPD = 9.0-16.6)
8.4 (95% HPD = 6.0-11.0)
2.5 (95% HPD = 1.3-4.1)
M. ushuaiensis Argentina
15.7 (95% HPD = 8.5-24.2)
12.9 (95% HPD = 8.0-17.7)
M. ushuaiensis Chile
15.7 (95% HPD = 8.5-24.2)
12.9 (95% HPD = 8.0-17.7)
1.3 (95% HPD = 0.5-2.4)
15.7 (95% HPD = 8.5-24.2)
12.9 (95% HPD = 8.0-17.7)
15.9 (95% HPD = 11.7-20.5)
10.6 (95% HPD = 7.9-12.4)
14.9 (95% HPD = 11.0-20.0)
10.9 (95% HPD = 8.4-13.8)
17.2 (95% HPD = 12.8-22.5)
9.1 (95% HPD = 6.2-12.4)
The divergence of Melanelixia Melanohalea, 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 . 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).
Estimates of genetic diversity for sampled lineages within Melanohalea
Fu and Li'sFs
M. elegantula Spain
M. exasperata Spain
M. multispora 1
M. multispora 2
M. multispora 3*
M. olivaceoides 1
M. olivaceoides 2
M. subolivacea 1
M. subolivacea 2*
M. ushuaiensis Argentina
M. ushuaiensis Chile
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.
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.  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 , 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 , 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 .
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 , 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. , 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 . In general, the distributions are in high concordance with those found in vascular plants and can be widely explained by contemporary ecogeographical factors . 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 .
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  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 , 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.
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 , 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 , 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 . 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. Melanohalea multispora 3 and M. subolivacea 2 are essentially western North American species, commonly found in dry forests . 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; ). 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 .
Nimis  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 ). 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 ). 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 , Alps , and Himalayas . 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 . 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.
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 (http://www.indexfungorum.org/) and seven previously unrecognized species-level lineages , awaiting formal description, were represented in this study. In this study nomenclature of undescribed species-level lineages follows arbitrary names in . 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
All primers used for PCR amplification and sequencing in this study
Forward primer sequence
Annealing temperature (°C)
5’- GTC TCC ACG TAT TCG CAT TCC-3’
We assembled and edited sequences using the program Sequencher v3.1.1 (Gene Codes Corporation, Ann Arbor, MI) and Geneious v5.4 . Sequence identity was confirmed with ‘megaBLAST’ search in GenBank . 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 (http://molevol.cmima.csic.es/castresana/Gblocks_server.html).
In this study we provide estimates of divergence events for species-level lineages presented in . In order to obtain divergence date estimates for Melanohalea, we used two Bayesian approaches implemented in BEAST v1.6.1 . 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. ). 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  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 . 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 . 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 . 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  to assess convergence and effective sampling size (ESS) of sampling parameters. Convergence was also assessed in AWTY  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 .
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 . 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 , to estimate a species tree from lineages defined in . *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 . Models of DNA sequence evolution for each marker were selected with the program jModeltest v0.1 , using the AIC criterion. Implementing an uncorrelated relaxed lognormal clock , 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 , 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  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  and Tajima’s D 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  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  in BEAST v1.6.1  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 . 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  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 . 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.
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).
- 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. 10.1007/s10531-006-9142-6.Google Scholar
- Lucking R, Rivas Plata E, Chaves JL, Umana L, Sipman HJM: How many tropical lichens are there …really?. Biblioth Lichenol. 2009, 100: 399-418.Google Scholar
- Printzen C: Lichen Systematics: The Role of Morphological and Molecular Data to Reconstruct Phylogenetic Relationships. In Progress in Botany 71. 2009, Berlin Heidelberg: Springer, 233-275. 71Google Scholar
- 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. 10.1016/j.ympev.2011.04.025.PubMedGoogle Scholar
- 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. 10.1016/j.ympev.2011.03.020.PubMedGoogle Scholar
- Lumbsch HT, Leavitt SD: Goodbye morphology? A paradigm shift in the delimitation of species in lichenized fungi. Fungal Divers. 2011, 50: 59-72. 10.1007/s13225-011-0123-z.Google Scholar
- 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. 10.1016/j.ympev.2011.05.012.PubMedGoogle Scholar
- Kroken S, Taylor JW: A Gene Genealogical Approach to Recognize Phylogenetic Species Boundaries in the Lichenized Fungus Letharia. Mycologia. 2001, 93: 38-53. 10.2307/3761604.Google Scholar
- 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. 10.3989/ajbm.2225.Google Scholar
- 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. 10.1007/s13225-010-0027-3.Google Scholar
- 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. 10.1007/s13127-011-0060-4.Google Scholar
- 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-10.1371/journal.pone.0028161.PubMedPubMed CentralGoogle Scholar
- 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. 10.3732/ajb.0900064.PubMedGoogle Scholar
- 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. 10.1111/j.1365-2699.2010.02469.x.Google Scholar
- Lücking R, Huhndorf S, Pfister DH, Plata ER, Lumbsch HT: Fungi evolved right on track. Mycologia. 2009, 101: 810-822. 10.3852/09-016.PubMedGoogle Scholar
- Taylor JW, Berbee ML: Dating divergences in the Fungal Tree of Life: review and new analyses. Mycologia. 2006, 98: 838-849. 10.3852/mycologia.98.6.838.PubMedGoogle Scholar
- Berbee ML, Taylor JW: Dating the molecular clock in fungi – how close are we?. Fungal Biol Rev. 2010, 24: 1-16. 10.1016/j.fbr.2010.03.001.Google Scholar
- 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. 10.1111/j.1558-5646.2010.01097.x.PubMedGoogle Scholar
- 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. 10.1016/j.ympev.2007.10.020.PubMedGoogle Scholar
- 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-10.1186/1741-7007-6-6.PubMedPubMed CentralGoogle Scholar
- 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. 10.1146/annurev.ecolsys.33.010802.150500.Google Scholar
- Graur D, Martin W: Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision. Trends Genet. 2004, 20: 80-86. 10.1016/j.tig.2003.12.003.PubMedGoogle Scholar
- Drummond A, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evoly Biol. 2007, 7: 214-10.1186/1471-2148-7-214.Google Scholar
- 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.PubMedGoogle Scholar
- Thorne JL, Kishino H: Divergence time and evolutionary rate estimation with multilocus data. Syst Biol. 2002, 51: 689-702. 10.1080/10635150290102456.PubMedGoogle Scholar
- Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-302. 10.1093/bioinformatics/19.2.301.PubMedGoogle Scholar
- 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-10.1186/1471-2148-8-257.Google Scholar
- 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.Google Scholar
- 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. 10.1016/j.ympev.2005.12.015.PubMedGoogle Scholar
- 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. 10.1016/j.mycres.2007.04.002.PubMedGoogle Scholar
- 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. 10.1016/j.ympev.2006.11.029.PubMedGoogle Scholar
- 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. 10.1016/j.ympev.2006.03.024.PubMedGoogle Scholar
- 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. 10.3732/ajb.0900140.PubMedGoogle Scholar
- 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.Google Scholar
- Crespo A, Lumbsch HT: Cryptic species in lichen-forming fungi. IMA Fungus. 2010, 1: 167-170. 10.5598/imafungus.2010.01.02.09.PubMedPubMed CentralGoogle Scholar
- 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. 10.1016/j.ympev.2010.04.014.PubMedGoogle Scholar
- 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. 10.1111/j.1365-294X.2010.04993.x.PubMedGoogle Scholar
- 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. 10.5248/108.347.Google Scholar
- 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. 10.1017/S0953756204000723.PubMedGoogle Scholar
- Divakar PK, Upreti DK, Elix JA: New species and new records in the lichen family Parmeliaceae (Ascomycotina) from India. Mycotaxon. 2001, 80: 356-362.Google Scholar
- Esslinger TL: A chemosytematic revision of the brown Parmeliae. J Hattori Bot Lab. 1977, 42: 1-211.Google Scholar
- Esslinger TL: A new species of Melanelia from Nepal. Mycotaxon. 1987, 28: 215-217.Google Scholar
- Otte V, Esslinger TL, Litterski B: Global distribution of the European species of the lichen genus Melanelia Essl. J Biogeogr. 2005, 32: 1221-1241. 10.1111/j.1365-2699.2005.01268.x.Google Scholar
- 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. 10.5248/111.65.Google Scholar
- 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. 10.1080/0028825X.1990.10412338.Google Scholar
- Esslinger TL: A new status for brown Parmeliae. Mycotaxon. 1978, 7:Google Scholar
- Purvis OW, James PW: Parmelia Ach., 1803. The Lichen Flora of Great Britain and Ireland. Edited by: Purvis OW, Coppins BJ, Hawksworth DL, James PW, Moore DM. 1992, London: Natural History Museum Publication, 421-437.Google Scholar
- Nimis PL: The lichens of Italy. 1993, Torino: Museo Regionale di Scienze Naturali, Monographia 12Google Scholar
- Seaward MRD: Lichen atlas of the British Isles. 1995, London: British Lichen SocietyGoogle Scholar
- Thell A: Forandringar i utbredningsbilden hos Parmelia elegantula och Parmelia laciniatula in Skåne och Blekinge. Graphis Scripta. 1989, 2: 156-160.Google Scholar
- 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.Google Scholar
- 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. 10.1046/j.1365-294X.2003.01812.x.Google Scholar
- 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. 10.1073/pnas.0809861105.PubMedPubMed CentralGoogle Scholar
- Knowles LL: Tests of Pleistocene speciation in Montane grasshoppers (genus Melanoplus) from sky islands of western North America. Evolution. 2000, 54: 1337-1348.PubMedGoogle Scholar
- 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. 10.1111/j.1365-294X.2005.02645.x.PubMedGoogle Scholar
- Zachos JC, Dickens GR, Zeebe RE: An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature. 2008, 451: 279-283. 10.1038/nature06588.PubMedGoogle Scholar
- Poelt J: Flechtenflora und Eiszeit in Europa. Phyton, Annales Tei Botanicae. 1963, 10: 206-215.Google Scholar
- Nimis PL, Martellos S: On the ecology of sorediate lichens in Italy. Lichenological contributions in honour of G. B. Feige. Bibliotheca Lichenologica 86. Edited by: Jensen M. 2003, Berlin: J Cramer, Berlin & Stuttgart, 393-406.Google Scholar
- 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 manuscriptGoogle Scholar
- Knowles LL, Carstens BC: Estimating a geographically explicit model of population divergence. Evolution. 2007, 61: 477-493. 10.1111/j.1558-5646.2007.00043.x.PubMedGoogle Scholar
- Liu L, Yu L, Kubatko L, Pearl DK, Edwards SV: Coalescent methods for estimating phylogenetic trees. Mol Phylogenet Evol. 2009, 53: 320-328. 10.1016/j.ympev.2009.05.033.PubMedGoogle Scholar
- Edwards SV: Is a new and general theory of molecular systematics emerging?. Evolution. 2009, 63: 1-19. 10.1111/j.1558-5646.2008.00549.x.PubMedGoogle Scholar
- Knowles LL: Estimating Species Trees: Methods of Phylogenetic Analysis When There Is Incongruence cross Genes. Syst Biol. 2009, 58: 463-467. 10.1093/sysbio/syp061. syp061PubMedGoogle Scholar
- Knowles LL, Carstens BC: Delimiting Species without Monophyletic Gene Trees. Syst Biol. 2007, 56: 887-895. 10.1080/10635150701701091.PubMedGoogle Scholar
- Edwards SV, Liu L, Pearl DK: High-Resolution Species Trees without Concatenation. Proc of the Natl Acad Sci USA. 2007, 104: 5936-5941. 10.1073/pnas.0607004104.Google Scholar
- Edwards S, Beerli P: Perspective: Gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies. Evolution. 2000, 54: 1839-1854.PubMedGoogle Scholar
- Van Devender T: The Sonoran desert tortoise. 2002, Tuscon: University of Arizona Press and the Arizona-Sonoran Desert MuseumGoogle Scholar
- 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. 10.1038/386788a0.Google Scholar
- 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. 10.1126/science.1059412.PubMedGoogle Scholar
- Berry EW: The Miocene flora of Idaho. Proc Natl Acad Sci. 1932, 18: 289-291. 10.1073/pnas.18.4.289.PubMedPubMed CentralGoogle Scholar
- Axelrod D: Age and origin of Sonoran desert vegetation. Occasional Papers, California Acedemy of Sciences. 1979, 132: 1-74.Google Scholar
- 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. 10.1130/G21184.1.Google Scholar
- Janis CM: Tertiary mammal evolution in the context of changing climates, vegetation, and tectonic events. Annu Rev Ecol Syst. 1993, 24: 467-500. 10.1146/annurev.es.24.110193.002343.Google Scholar
- Jacobs BF, Kingston JD, Jacobs LL: The origins of grass-dominated ecosystems. Ann Mo Bot Gard. 1999, 86: 590-643. 10.2307/2666186.Google Scholar
- Pagani M, Freeman KH, Arthur MA: Late Miocene Atmospheric CO2 Concentrations and the Expansion of C4 Grasses. Science. 1999, 285: 876-879. 10.1126/science.285.5429.876.PubMedGoogle Scholar
- 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. 10.1639/0007-2745(2005)108[1:PPAAGE]2.0.CO;2.Google Scholar
- Endler JA: Natural selection in the wild. 1986, Princeton: Princeton University PressGoogle Scholar
- Schluter D: The ecology of adaptive radiations. 2000, Oxford: Oxford University PressGoogle Scholar
- 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. 10.1016/j.yqres.2009.09.008.Google Scholar
- 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.Google Scholar
- 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. 10.1111/j.1558-5646.2011.01287.x.PubMedGoogle Scholar
- 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. 10.1046/j.1365-2699.2003.00864.x.Google Scholar
- Thomson JW: The Distribution of Arctic Lichens and Thoughts Concerning Their Origin. Lichenologist. 1995, 27: 411-416.Google Scholar
- 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. 10.1029/JD094iD15p18409.Google Scholar
- 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. 10.1016/0012-821X(86)90024-5.Google Scholar
- Meyers SR, Hinnov LA: Northern Hemisphere glaciation and the evolution of Plio-Pleistocene climate noise. Paleoceanography. 2010, 25: PA3207-Google Scholar
- 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. 10.1073/pnas.95.4.1961.PubMedPubMed CentralGoogle Scholar
- 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. 10.1016/j.ympev.2006.07.002.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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. 10.1111/j.1558-5646.2009.00803.x.PubMedGoogle Scholar
- 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. 10.2307/1313224.Google Scholar
- 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.Google Scholar
- George K: Saalian supercycle, Mindel/Riss interglacial and Milankovitch's dating. Quat Sci Rev. 2005, 24: 1573-1583. 10.1016/j.quascirev.2004.08.023.Google Scholar
- Ho SYW, Shapiro B: Skyline-plot methods for estimating demographic history from nucleotide sequences. Mol Ecol Res. 2011, 11: 423-434. 10.1111/j.1755-0998.2011.02988.x.Google Scholar
- 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. 10.1111/j.1365-2699.2009.02239.x.Google Scholar
- 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. 10.1016/j.ympev.2008.04.035.PubMedGoogle Scholar
- 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-10.1186/1471-2148-11-293.PubMedPubMed CentralGoogle Scholar
- Seifert KA: Progress towards DNA barcoding of fungi. Mol Ecol Res. 2009, 9: 83-89.Google Scholar
- Simon UK, Weiß M: Intragenomic variation of fungal ribosomal genes is higher than previously thought. Mol Biol Evol. 2008, 25: 2251-2254. 10.1093/molbev/msn188.PubMedGoogle Scholar
- 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. 10.1007/s00294-005-0581-5.PubMedGoogle Scholar
- 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. 10.1111/j.1469-8137.2011.03677.x.PubMedGoogle Scholar
- 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. 10.1016/j.mycres.2007.05.006.PubMedGoogle Scholar
- 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. 10.1016/j.plaphy.2007.01.004.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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. 10.1111/j.1558-5646.2009.00685.x.PubMedGoogle Scholar
- 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. 10.1111/j.1365-294X.2006.02880.x.PubMedGoogle Scholar
- Pino-Bodas R, Burgaz AR, Martin MP, Lumbsch HT: Species delimitations in the Cladonia cariosa group (Cladoniaceae, Ascomycota). Lichenologist. 2012, 44: 121-135. 10.1017/S002428291100065X.Google Scholar
- 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. 10.1007/s13127-011-0066-y:1–21. in pressGoogle Scholar
- Gregory-Wodzicki KM: Uplift history of the Central and Northern Andes: A review. Geol Soc Am Bull. 2000, 112: 1091-1105. 10.1130/0016-7606(2000)112<1091:UHOTCA>2.0.CO;2.Google Scholar
- 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. 10.1007/BF00374424.Google Scholar
- 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. 10.1130/0091-7613(1990)018<0354:ERUITH>2.3.CO;2.Google Scholar
- 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. 10.3767/003158509X470602.PubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- White TJ, Bruns TD, Lee S, Taylor J: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. 1990, San Diego: Academic Press, 315-322.Google Scholar
- 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. 10.1016/S0168-1656(99)00154-6.Google Scholar
- 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.Google Scholar
- 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. 10.1046/j.1365-294x.1999.00820.x.PubMedGoogle Scholar
- 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. 10.3732/ajb.1000230.PubMedGoogle Scholar
- 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. 10.3732/ajb.89.4.688.PubMedGoogle Scholar
- Liu YJ, Whelen S, Hall BD: Phylogenetic relationships among ascomycetes: evidence from an RNA polymerse II subunit. Mol Biol Evol. 1999, 16: 1799-1808. 10.1093/oxfordjournals.molbev.a026092.PubMedGoogle Scholar
- Drummond AJ, Ashton B, Cheung M, Heled J, Kearse M, Moir R, Stones-Havas S, Thierer T, Wilson A: Geneious ver. 5.3. 2010, http://www.geneious.com,Google Scholar
- 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:Google Scholar
- Castresana J: Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis. Mol Biol Evol. 2000, 17: 540-552. 10.1093/oxfordjournals.molbev.a026334.PubMedGoogle Scholar
- Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007, 56: 564-577. 10.1080/10635150701472164.PubMedGoogle Scholar
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.PubMedGoogle Scholar
- Stamatakis A, Hoover P, Rougemont J: A Rapid Bootstrap Algorithm for the RAxML Web Servers. Syst Biol. 2008, 57: 758-771. 10.1080/10635150802429642.PubMedGoogle Scholar
- Rambaut A, Charleston M: Tree Edit. Phylogenetic Tree Editor v1.0 alpha 8. 2001, [http://tree.bio.ed.ac.uk/software/treeedit/] webciteGoogle Scholar
- Rambaut A, Drummond A: Tracer [computer program]. 2003Google Scholar
- 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. 10.1093/bioinformatics/btm388.PubMedGoogle Scholar
- Rambault A, Drummond AJ: TreeAnnotator. Version 1.6.1. 2009, Available at: http://www.beastbioedacuk/TreeAnnotatorGoogle Scholar
- Heled J, Drummond AJ: Bayesian Inference of Species Trees from Multilocus Data. Mol Biol Evol. 2010, 27: 570-580. 10.1093/molbev/msp274.PubMedPubMed CentralGoogle Scholar
- Posada D: jModelTest: Phylogenetic Model Averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMedGoogle Scholar
- Drummond A, Ho S, Phillips M, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.PubMedPubMed CentralGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.PubMedGoogle Scholar
- Fu Y-X: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selections. Genet. 1997, 174: 915-925.Google Scholar
- Tajima F: Statistical method for testing the neurtral mutation hypothesis by DNA polymorphism. Genet. 1989, 123: 585-595.Google Scholar
- 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.Google Scholar
- 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. 10.1093/molbev/msi103.PubMedGoogle Scholar
- Donnelly P, Tavare S: Coalescents and genealogical structure under neutrality. Annu Rev Genet. 1995, 29: 401-421. 10.1146/annurev.ge.29.120195.002153.PubMedGoogle Scholar
- Robinson H: Considerations on the evolution of lichens. Phytologia. 1975, 32: 407-413.Google Scholar
- Buschbom J, Barker D: Evolutionary History of Vegetative Reproduction in Porpidia s.l. (Lichen-Forming Ascomycota). Syst Biol. 2006, 55: 471-484. 10.1080/10635150600697465.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.