Transoceanic Dispersal and Subsequent Diversification on Separate Continents Shaped Diversity of the Xanthoparmelia pulla Group (Ascomycota)

Transcription

1 Transoceanic Dispersal and Subsequent Diversification on Separate Continents Shaped Diversity of the Xanthoparmelia pulla Group (Ascomycota) Guillermo Amo de Paz 1, Paloma Cubas 1, Ana Crespo 1, John A. Elix 2, H. Thorsten Lumbsch 3 * 1 Departamento de Biología Vegetal II, Facultad de Farmacia, Universidad Complutense de Madrid, Madrid, Spain, 2 Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia, 3 Department of Botany, The Field Museum, Chicago, Illinois, United States of America Abstract In traditional morphology-based concepts many species of lichenized fungi have world-wide distributions. Molecular data have revolutionized the species delimitation in lichens and have demonstrated that we underestimated the diversity of these organisms. The aim of this study is to explore the phylogeography and the evolutionary patterns of the Xanthoparmelia pulla group, a widespread group of one of largest genera of macrolichens. We used a dated phylogeny based on nuits and nulsu rdna sequences and performed an ancestral range reconstruction to understand the processes and explain their current distribution, dating the divergence of the major lineages in the group. An inferred age of radiation of parmelioid lichens and the age of a Parmelia fossil were used as the calibration points for the phylogeny. The results show that many species of the X. pulla group as currently delimited are polyphyletic and five major lineages correlate with their geographical distribution and the biosynthetic pathways of secondary metabolites. is the area where the X. pulla group radiated during the Miocene times, and currently is the region with the highest genetic, morphological and chemical diversity. From this center of radiation the different lineages migrated by long-distance dispersal to others areas, where secondary radiations developed. The ancestral range reconstruction also detected that a secondary lineage migrated from Australia to South America via long-distance dispersal and subsequent continental radiation. Citation: Amo de Paz G, Cubas P, Crespo A, Elix JA, Lumbsch HT (2012) Transoceanic Dispersal and Subsequent Diversification on Separate Continents Shaped Diversity of the Xanthoparmelia pulla Group (Ascomycota). PLoS ONE 7(6): e doi: /journal.pone Editor: Zhengguang Zhang, Nanjing Agricultural University, China Received April 17, 2012; Accepted May 25, 2012; Published June 20, 2012 Copyright: ß 2012 Amo de Paz et al. 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 author and source are credited. Funding: This study was supported by a grant from the Spanish Ministerio de Ciencia e Innovación (CGL ) and a grant by the National Science Foundation (USA) (DEB ). The Universidad Complutense de Madrid is thanked for providing all complementary facilities and funds (SYSTEMOL research group UCM-BSCH GR35/10A). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * tlumbsch@fieldmuseum.org Introduction Methods for delimiting species, the fundamental taxonomic unit, have always fascinated evolutionary biologists [1 3]. Understanding the circumscription of species is important for biological and ecological studies and for conservation issues. However, the main challenge is to recognize species in organisms with relatively simple morphologies. In lichenized fungi traditional species circumscriptions are based on phenotypic characters, such as thallus and ascomatal morphology or chemical characters. However, there is a growing body of evidence from molecular studies that the traditional morphology-based species circumscriptions are insufficient to represent the diversity in lichenized ascomycetes [4 26]. A number of DNA sequence-based phylogenetic studies revealed the presence of distinct lineages within currently delimited species. Subsequent, detailed studies often revealed previously overlooked morphological subtleties or chemical differences among those clades and authors often refer to these as semi-cryptic species [8]. On a par with the phenotypic-based species circumscription, researchers often accepted wide distribution ranges for species occurring on different continents. This was at least partially due to a common belief by lichenologists in the everything is everywhere hypothesis [27,28] applied to fungi, as discussed elsewhere [8,29]. In several cases molecular data assisted in a better understanding of the biogeography of lichen-forming fungi where taxa were shown to represent different species on different continents, e.g. the Leptogium furfuraceum aggr. on different continents [18], Melanelixia glabra s. lat. in Europe and North America [14], Parmelina quercina s. lat. on different continents [4], Physcia aipolia aggr. in Europe and Australia [9], Xanthoparmelia spp. in North America and Australia [15,16]. In the Leptogium furfuraceum aggr. complex sister-group relationships were found between populations from the same hemispheres which were incongruent with previous classifications based on morphological differences [18], and the dated phylogeny indicated that the species had migrated via transoceanic dispersal to different continents. Here we report another case of a group of lichenized fungi where transoceanic dispersal to different continents is correlated with the phylogenetic lineages. The group studied here is the Xanthoparmelia pulla group which belongs to the family Parmeliaceae. This family represents one of the largest families of lichenized fungi [30,31]. The main clade of the family is the parmelioid clade with almost 2000 species [32] currently classified in 27 genera, with Xanthoparmelia being the largest with over 800 accepted species [33]. The species in this genus characteristically occur on siliceous rocks or soil, predominantly in arid to subarid PLoS ONE 1 June 2012 Volume 7 Issue 6 e39683

2 regions, with a center of distribution in the southern hemisphere. The genus is characterized by having cell wall polysaccharides of the Xanthoparmelia-type, small ascospores with an arachiform vacuolar body [34], and the presence of a pored epicortex [33,35]. It has been hypothesized that the genus diversified in a rapid radiation following a shift towards drier habitats at the base of the Xanthoparmelia clade [36] leading to the high current diversity. Previously, the Xanthoparmelia pulla group has been classified within the separate genus Neofuscelia based on the different cortical chemistry (having melanoid pigments and lacking usnic acid or atranorin, characteristic of the majority of Xanthoparmelia species) [37,38]. A subsequent molecular study showed that the genus Neofuscelia was polyphyletic, with its clades scattered within Xanthoparmelia [33]. Consequently, the genus Neofuscelia was reduced to synonymy with Xanthoparmelia, as have other genera previously distinguished by cortical chemistry or growth form [33,39 43]. The Xanthoparmelia pulla group is a monophyletic clade within the complete Xanthoparmelia clade, that includes the former Esslinger s Xanthoparmelia pulla species and other related species [44]. Although the clades were largely incongruent with the current species circumscription, we found a correlation of the main clades of X. pulla group with their geographical distribution and chemical profile (Fig. 1, 2). Clade 1 includes specimens from California, Macaronesia and areas around the Mediterranean basin, all of which contain depsides and depsidones derived from the orcinol pathway or with aliphatic acids; clade 2 includes specimens from Australia (subclade 2.1) and South America (subclade 2.2) with depsides and depsidones derived from the orcinol pathway; clade 3 derives from n specimens containing olivetoric acid; clade 4 specimens with hypostictic acid from California (subclade 4.1) and South America (subclade 4.2); and clade 5 specimens with physodic acid (orcinol depsidones) from. The species delimitation within the Xanthoparmelia pulla group is currently based on a combination of morphological and chemical characters (Table 1). The morphological characters include the color of the lower surface, shape of the lobes, attachment to the substrate, and presence of vegetative propagules while the chemical differences pertain to upper cortical and medullary secondary metabolites. A number of the currently accepted species have a wide distribution spanning several continents. To address the species delimitation in this group and to test the hypothesis of widely distributed species we have generated a data set using two loci (nuits rdna, nulsu rdna) from specimens collected on different continents. The molecular data were used to perform phylogenetic reconstructions in a maximum likelihood (ML) and Bayesian (B/MCMC) framework. We have also estimated the timing of the diversification events leading to the main clades found in our study to discriminate between vicariance and longdistance dispersal as possible explanations for the current distribution patterns. A Bayesian-based approach of ancestral range reconstruction was used to identify potential areas in which the group and major clades within the group originated. Results Phylogenetic analyses One hundred sixty-eight DNA sequences of ITS and nulsu rdna of 88 representative specimens of Xanthoparmelia were assembled. One hundred forty of these sequences were newly generated in this study. The specimens included 25 currently accepted species in the Xanthoparmelia pulla group, four unassigned specimens, and six samples of four species as outgroup. A data matrix of 1283 unambiguously aligned characters, with 454 characters in the ITS and 829 characters in the nulsu rdna was used for phylogenetic analyses. The data set included 1081 constant characters. The general time-reversible model with a gamma distribution and invariant model of rate heterogeneity (GTR+I+G) was employed for analyses of the single-loci and concatenated data sets. Since no strongly supported conflicts between the two single-locus ML phylogenetic trees were detected, a combined data set was analyzed. In the B/MCMC analysis of the combined data set, the likelihood parameters in the sample had the following averaged values for the partitioned data set (6 standard deviation): base frequencies p(a) = 0.25 (61.54E-4), p(c) = 0.24 (61.42E-4), p(g) = 0.28 (61.58E-4), p(t) = 0.23 (61.51E-4); rate matrix r(ac) = 4.42 (61.43E-4), r(ag) = 0.23 (68.35E-4), r(at) = 9.53 (62.19E-4), r(cg) = 4.82 (61.45E-4), r(ct) = 0.54 (69.26E-4), r(gt) = 2.98 (61.11E-4) and the gamma shape parameter a = 0.21 (63.86E-4). The likelihood parameters in the sample had a mean likelihood of LnL = (60.49), while the ML tree had a likelihood of LnL = The phylogenetic estimates of the ML and B/MCMC analyses were congruent, hence only the ML tree (Fig. 1) is shown here. Specimens of the Xanthoparmelia pulla group form a strongly supported monophyletic group with five main, mostly wellsupported, clades (Fig. 1). The clades do not agree with the current species circumscription, with 11 species being polyphyletic, five of them with specimens from different continents entering different clades. For example, all the Northern Hemisphere specimens identified as X. luteonotata, X. pulla, X. delisei or X. glabrans belong to clade 1 while all the Australian specimens of the same species belong to clade 2.1. Similarly, specimens of X. imitatrix from South America belong to clade 2.2 and are not directly related to the n specimen. Estimates of divergence times and ancestral range reconstructions A Bayesian phylogenetic tree was dated to estimate the age of the X. pulla group and its main clades. The results of the divergence time analysis are summarized in Fig. 2, and the whole parmelioid tree is shown as supp. mat. (Fig. S1). The Xanthoparmelia pulla group started to diversify around Ma ( Ma), the age of the crown node of clade 1 was estimated at 5.31 Ma ( Ma), the ancestor of clade 2 around 8.10 Ma ( Ma), and the crown of clade 4 around 3.44 Ma ( Ma). The results of the ancestral range reconstruction analyses are summarized in Figs. 1 and 2. This established that was the most likely origin of the X. pulla group, with a marginal probability of 0.745, indicating localized uncertainty. The four other areas explored (South America, Australia, California and the Mediterranean basin) were rejected with probabilities below For the base of clade 1, the Mediterranean basin was reconstructed as the most likely ancestral range with a marginal probability of 0.555, but California could not be rejected (probability of 0.104). For clade 2, Australia was recovered as the most likely ancestral area with marginal probability of 0.672; while South America, the other area from which specimens of this clade occur, was rejected as potential ancestral area (p,0.05); similarly, California, the Mediterranean basin and were also rejected as ancestral areas. South America was found to be the most likely origin for clade 4 (which also includes specimens from California and South America) with a marginal probability of 0.48, although neither California nor were rejected (probabilities of and 0.083, respectively). Australia and the Mediterranean basin were rejected as ancestral ranges for clade 4. PLoS ONE 2 June 2012 Volume 7 Issue 6 e39683

3 Table 1. Main differences of the species of Xanthoparmelia pulla group studied in this paper [38,44,47,59,60]. Characters and distribution Species Morphology of lobes Isidia Lower surface Chemistry Distribution X. atroviridis 1 2 mm broad, subirregular, contiguous to imbricate X. caliginosa (23.5) mm broad, subirregular, contiguous to imbricate X. delisei 1 4 mm broad, sublinear to irregular, flat to slightly concave or convex, becoming laciniate, often imbricate and entangled X. fissurina 1 3 mm broad, contiguous to imbricate or entangled X. glabrans mm broad, sublinear to linear-elongate, imbricate to loosely entangled. X. imitatrix mm broad, sublinear to linear-elongate, imbricate to laciniate entangled, rarely developing subfruticose branches X. lineella mm broad, linear and dichotomously branched and entangled X. loxodes (0.5-)1 3(25) mm broad, subirregular to sublinear, contiguous to entangled. X. luteonotata (0.5-)1 3 mm broad, sublinear to irregular, discrete to imbricate, rarely developing subfruticose branches X. pokornyi 1 2 mm broad, sublinear to linear, discrete to loosely imbricate or entangled X. perrugata 1 3 (25) mm broad, sublinear to linear-elongate, discrete to imbricate or entangled. X. pseudoglabrans mm broad, subirregular to sublinear, imbricate to entangled X. pulla 1 3 (25) mm broad, sublinear to linear-elongate, discrete to imbricate or entangled Sparse to crowded and areolate. Isidia erumpent, globose to cylindrical, (20.8) mm tall. Sparsely to densely isidiate, isidia more or less spherical and distinctly pustular, erumpent Black, moderately concolorous, up to 0.4 mm long. moderately rhizinate, rhizines more or less concolorous, up to 0.6 mm long. often paler near the apices, moderately to densely rhizinate, the rhizines simple and concolorous with the lower surface, to 1 mm long Pale tan to pale brown, moderately rhizinate, rhizines concolorous, to 1 mm long black, dull, slightly rugulose, moderately rhizinate; rhizines black, simple or fasciculate, to 1 mm long sparsely to moderately simple, to 1.5 mm long Black, sparsely rhizinate, rhizines concolorous, to 1 mm long smooth to somewhat rugulose, moderately concolorous, to 1mmlong Pale tan to pale brown, moderately to densely simple, to 0.5 mm long Pale tan to brown, moderately to sparsely concolorous or darkening, to 1 (21.5) mm long moderately to densely simple, to 1.5 mm long. Black; moderately rhizinate or rhizines rather parse, concolorous with the lower surface moderately to densely simple, to 1.5 mm long Medulla: Hypoconstictic, hypostictic, hyposalazinic acids; Cortex: HNO 3 + violet Medulla: Olivetoric acid. Cortex: HNO 3 + Medulla: glomelliferic, glomellic, perlatolic acids; 6 gyrophoric acid. Cortex: HNO 3 + Medulla: hypostictic, hypoconstictic, hyposalazinic acids, unknown compounds. Cortex: HNO 3 + Medulla: alectoronic acid; 6 a-collatolic and gyrophoric acids. Cortex: HNO 3 + Europe, Asia, Africa, Australia, Macaronesia, South America and South America Australia, Europe, Africa, South America, New Zealand Medulla: physodic acid; 64-Omethylphysodic and alectoronic America, New Zealand Australia, Africa, South acids. Cortex: HNO 3 + Medulla: physodic acid; Cortex: HNO 3 + Medulla: glomelliferic, glomellic Europe, North Africa, and perlatolic acids; 6gyrophoric Asia, North America, acid. Cortex: HNO 3 + Macaronesia Medulla: 6 divaricatic and stenosporic acids; 6 gyrophoric acid. Cortex: HNO 3 + Medulla: stenosporic acid; 6 gyrophoric and divaricatic acids. Cortex: HNO 3 + Medulla: divaricatic acid; 6 stenosporic, oxostenosporic, gyrophoric, lecanoric acids. Cortex: HNO 3 + Medulla: alectoronic acid; 6 a-collatolic acid. Cortex: HNO 3 - Medulla: stenosporic acid; 6 divaricatic, gyrophoric, perlatolic, 4-O-demethylstenosporic acid, oxostenosporic acids. Cortex: HNO 3 + Australia, Europe, Africa, New Zealand Europe, Asia Europe, North Africa, Australia, Asia Europe, Australia, New Zealand, Africa PLoS ONE 3 June 2012 Volume 7 Issue 6 e39683

4 Table 1. Cont. Characters and distribution Species Morphology of lobes Isidia Lower surface Chemistry Distribution X. pulloides 1 2 mm broad, subirregular to sublinear, contiguous to subimbricate X. quintarioides (23) mm broad, strongly convex and shortflabellate, discrete but close to more or less contiguous X. ryssolea 1 3 mm broad, linear elongate, subterete, convex. X. squamans 1 2 mm broad, sublinear, imbricate to contiguous X. subhosseana 1 2 mm broad, subirregular, contiguous to slightly imbricate X. subimitatrix mm broad, sublinear to subirregular, discrete to subimbricate X. subincerta mm broad, flat, sublinear, more or less imbricate X. subprolixa 1 3 (25) mm broad, sublinear to linear-elongate, discrete to imbricate or entangled. Sparsely to densely isidiate. Isidia pustular, erumpent Isidia cylindrical, simple or densely branched, mm tall. Apices syncorticate X. torulosa mm broad, sublinear to sublirregular, imbricate; laciniae at periphery and within thallus, 6 subfruticose, sublinear to elongate, mm broad. X. verisidiosa 1 3 mm broad, irregular to sublinear, flat, short and rounded, contiguous to imbricate X. verrucella mm wide, irregular to sublinear, flat, imbricate to entangled. X. verruculifera 1 2 mm broad, subirregular to sublinear, contiguous to imbricate doi: /journal.pone t001 Sparsely to densely isidiate. Isidia cylindrical, simple or becoming densely branched, mm tall. Apices syncorticate Moderate to densely isidiate. Isidia cylindrical, simple or becoming branched, to 1 mm tall. Apices syncorticate Sparcely to densely isidiate. Isidia pustular, erumpent. Black, moderately concolorous, to 0.5 mm long Tan to pale brown, sparsely to moderately rhizinate, the rhizines short and hapterate Pale yellow-brown to red-brown, canaliculate, sparsely concolorous, to 0.6 mm long. moderately to sparsely concolorous, to 1 mm long moderately rhizinate, rhizines concolorous, to 0.6 mm long Pale tan to brown, moderately rhizinate, rhizines simple, brown or often blackened, to 0.8 mm long Black, moderately simple, black, to 0.3 mm long often paler at apices, moderately to densely simple, to 1.5 mm long Black, moderately to densely rhizinate; rhizines simple or occasionally tufted, slender. Black, moderately to sparsely rhizinate, rhizines simple, black to 0.4 mm long Black, sparsely to moderately rhizinate, rhizines simple, simple, black, to 0.4 mm long. moderately rhizinate, rhizines concolorous, to 0.8 mm long Medulla: constipatic and protoconstipatic acids; 6 gyrophoric acid. Cortex: HNO 3 + Medulla: hypostictic, hypoconstictic, cryptostictic acids; 6 hyposalazinic acid. Cortex: HNO 3 + Medulla: stenosporic acid; 6 gyrophoric, oxostenosporic, divaricatic acids. HNO 3 + Medulla: hypostictic, hypoconstictic, hyposalazinic acids. HNO 3 + Medulla: hypostictic, hyposalazinic, hypoconstictic acids. Cortex: HNO 3 + Medulla: physodic acid and alectoronic acids. Cortex: HNO 3 + Medulla: glomelliferonic acid; 6 loxodellonic and glomellonic acids. Cortex: HNO 3 + Medulla: divaricatic acid; 6 stenosporic, nordivaricatic acids. Cortex: HNO 3 + Medulla: divaricatic acid; 6 nordivaricatic, stenosporic acids. Cortex: HNO 3 + Medulla: alectoronic and a-collatolic acids. Cortex: HNO 3 + Medulla: divaricatic acid; 6 stenosporic acid. Cortex: HNO 3 + Medulla: divaricatic acid; 6 stenosporic and gyrophoric acids. Cortex: HNO 3 + Macaronesia, Asia Europe, Asia, South America, New Zealand, North America, South America, New Zealand, Australia. Australia, Australia, New Zealand Australia Australia, New Zealand, Australia, New Zealand, North Africa, Europe, North America Discussion Understanding the diversity and delimiting species in lichenized fungi has been a long standing challenge and current studies using molecular data have dramatically changed our ability to distinguish species in this group [8,23,45]. The Xanthoparmelia pulla group is a good example for illustrating the difficulties in in distinguishing species by morphology due to the remarkable PLoS ONE 4 June 2012 Volume 7 Issue 6 e39683

5 PLoS ONE 5 June 2012 Volume 7 Issue 6 e39683

6 Figure 1. Phylogenetic relationships of the Xanthoparmelia pulla group based on nuits and nulsu rdna sequences. Topology based on maximum-likelihood (ML) analyses. Posterior probabilities and bootstrap values are indicated on each branch. Branches with posterior probabilities under Bayesian analysis equal or above 0.95 and/or bootstrap values equal or above 70% under MP are in bold. Medullary compounds, and results of the divergence time estimation and ancestral range reconstruction analyses are shown. Black arrow head indicates the polyphyletic species. doi: /journal.pone g001 plasticity of morphological characters in this group. Consequently, secondary metabolites have played an important role in delimiting species in this group [44,46,47]. Following the current classification using a combination of vegetative morphology and secondary chemistry, a number of species have broad geographical distributions spanning several continents. Here we have used molecular data to investigate the current classification within the group and attempt to explain their distribution. We used likelihood-based and Bayesian approaches to investigate the evolutionary origin of the group and timing of speciation events. Hopefully such data will reveal evolutionary patterns so we may develop a framework for their taxonomic classification which better reflects the phylogenetic relationships in the X. pulla group. Our results clearly indicate that the species as currently delimited are polyphyletic (Fig. 1). This is consistent with results from other studies of Xanthoparmelia species believed to occur on different continents which were subsequently found to represent distinct lineages [15,16]. Further, similar patterns have been found in other groups of lichenized fungi [4,14,18,48]. The ancestral range reconstruction points to as the most likely origin of the X. pulla group. Although there is a certain degree of uncertainty in the reconstruction (marginal probability of 0.745), the analysis rejected other areas as potential ancestral areas for the group. Interestingly, has the highest morphological and chemical diversity within the group and the specimens studied here belong to different, unrelated lineages (Fig. 1). n specimens containing olivetoric acid cluster in clade 3 and those with physodic acid in clade 5. The phylogenetic relationships of other specimens from with hypostictic acid, physodic acid or other orcinol depsides and depsidones are still unresolved. The n specimens show remarkable morphological variability, including subcrustose and foliose species. Further, many Xanthoparmelia species occur in arid climates and the diversification of X. pulla group occurred around Ma ( ). At this time the Cape Region underwent a major aridification [49], which may be responsible for the rapid radiation and current richness of the Cape flora. Thus, it is likely that the X. pulla group originated in around the same time. Unlike most Cape region elements in flowering plants, the species of the X. pulla group subsequently extended their distribution by transoceanic dispersal. Within the X. pulla group, the five lineages identified are characterized by the presence of different substance classes and in some cases they diverged secondarily in different geographical areas. The correlation between chemical pathways and the lineages found in molecular studies has also been found in Pertusariaceae among lichenized fungi [50 52]. Clade 1 includes specimens containing orcinol depsides and depsidones that occur in California and the Mediterranean basin. Neither internally supported subclades nor a geographical pattern was found within this clade, and specimens with different phenotypical characters from different geographical areas are Figure 2. Schematic map showing the relationships between phylogeny, medullary compounds, ancestral range and divergence times estimation for the Xanthoparmelia pulla group. doi: /journal.pone g002 PLoS ONE 6 June 2012 Volume 7 Issue 6 e39683

7 Table 2. Specimens used in this study with country of collection, voucher information and GenBank accession numbers. GenBank accession no. Species Country Herbarium acc. no. nuits nulsu Xanthoparmelia adhaerens 1 MAF-Lich HM HM X. adhaerens 2 MAF-Lich HM HM X. atroviridis 1 MAF-Lich JQ JQ X. atroviridis 1 MAF-Lich JQ JQ X. atroviridis 2 MAF-Lich JQ JQ X. atroviridis 2 MAF-Lich JQ JQ X. atroviridis 3 MAF-Lich JQ JQ X. caliginosa 1 MAF-Lich JQ JQ X. caliginosa 2 MAF-Lich JQ X. caliginosa 3 MAF-Lich JQ JQ X. caliginosa 4 MAF-Lich JQ X. caliginosa 5 MAF-Lich JQ JQ X. delisei 1 Turkey MAF-Lich JQ JQ X. delisei 2 Australia MAF-Lich 7432 AY AY X. delisei 3 Spain MAF-Lich 7659 AY AY X. delisei 4 Turkey MAF-Lich JQ JQ X. delisei 5 Turkey MAF-Lich JQ JQ X. fissurina MAF-Lich JQ JQ X. glabrans 1 Australia CANB JQ JQ X. glabrans 2 Australia CANB JQ JQ X. glabrans 3 Australia MAF-Lich 7665 AY AY X. glabrans 4 Australia CANB JQ JQ X. glabrans 5 Spain MAF-Lich 9912 AY AY X. glabrans 6 Turkey MAF-Lich JQ JQ X. glabrans 7 Morocco MAF-Lich JQ JQ X. imitatrix 1 Chile MAF-Lich JQ JQ X. imitatrix 2 Chile MAF-Lich JQ JQ X. imitatrix 3 Chile MAF-Lich JQ JQ X. imitatrix 4 Chile MAF-Lich JQ JQ X. imitatrix 5 Chile MAF-Lich JQ JQ X. imitatrix 6 Chile MAF-Lich JQ JQ X. imitatrix 7 MAF-Lich JQ JQ X. lineella MAF-Lich JQ JQ X. loxodes 1 Spain MAF-Lich 7072 AY AY X. loxodes 2 Spain MAF-Lich 6206 AY AY X. luteonotata 1 Spain MAF-Lich JQ JQ X. luteonotata 2 Australia CANB JQ X. luteonotata 3 Australia CANB JQ JQ X. luteonotata 4 Spain MAF-Lich JQ JQ X. mougeotii 1 Spain MAF-Lich 6802 AY37006 AY X. mougeotii 2 Spain MAF-Lich 9916 AY AY X. perrugata Spain MAF-Lich JQ X. pokornyi 1 Spain MAF-Lich 6052 AY AY X. pokornyi 2 Spain MAF-Lich 9908 AY AY X. pokornyi 3 Turkey MAF-Lich JQ JQ X. pokornyi 4 Spain MAF-Lich JQ X. pokornyi 5 Turkey MAF-Lich JQ JQ X. pokornyi 6 Turkey MAF-Lich JQ JQ X. pokornyi 7 Turkey MAF-Lich JQ JQ PLoS ONE 7 June 2012 Volume 7 Issue 6 e39683

8 Table 2. Cont. GenBank accession no. Species Country Herbarium acc. no. nuits nulsu X. pseudoglabrans MAF-Lich JQ JQ X. pulla 1 Spain MAF-Lich JQ X. pulla 2 Australia CANB JQ JQ X. pulla 3 Australia CBG JQ JQ X. pulla 5 Spain MAF-Lich 6794 AY AJ X. pulloides 1 Spain MAF-Lich JQ JQ X. pulloides 2 Spain MAF-Lich 6784 AY AY X. quintarioides MAF-Lich JQ JQ X. ryssolea 1 Turkey MAF-Lich JQ JQ X. ryssolea 2 Turkey MAF-Lich JQ JQ X. ryssolea 3 Spain MAF-Lich JQ X. sp. 1 MAF-Lich JQ JQ X. sp. 2 MAF-Lich JQ JQ X. sp. 3 MAF-Lich JQ X. sp. 4 MAF-Lich JQ JQ X. squamans 1 Chile MAF-Lich JQ JQ X. squamans 2 Chile MAF-Lich JQ JQ X. squamans 3 Chile MAF-Lich JQ JQ X. subhosseana 1 USA MAF-Lich JQ JQ X. subhosseana 2 Chile MAF-Lich JQ JQ X. subimitatrix Chile MAF-Lich JQ JQ X. subincerta 1 Australia CANB JQ JQ X. subincerta 2 Australia MAF-Lich 7494 AY AY X. subprolixa 1 Australia MAF-Lich 7667 AY AY X. subprolixa 2 Australia CANB JQ JQ X. tegeta Australia MAF-Lich 7523 AY AY X. torulosa 1 Australia CANB JQ JQ X. torulosa 2 Australia CANB JQ JQ X. verisidiosa 1 Australia CANB JQ JQ X. verisidiosa 2 Australia CANB JQ JQ X. verrucella 1 Australia CANB JQ JQ X. verrucella 2 Australia CANB JQ JQ X. verrucella 3 MAF-Lich JQ JQ X. verruculifera 1 USA MAF-Lich JQ JQ X. verruculifera 2 USA MAF-Lich JQ JQ X. verruculifera 3 USA MAF-Lich JQ JQ X. verruculifera 4 USA MAF-Lich JQ JQ X. verruculifera 5 Spain MAF-Lich JQ X. xanthomelaena Australia MAF-Lich HM HM New sequences are in bold. doi: /journal.pone t002 intermingled. In fact, shifting between orcinol depsides and depsidones can occur by one-step transformations [46]. By contrast, in other genera (e.g. Melanelixia, Parmelina, Leptogium) with similar disjunct distributions (North America and the Mediterranean basin), the geographical distribution correlates with the clades found in the molecular study. In the X. pulla group this pattern was not found, possibly due to insufficient sampling or absence of a phylogenetic signal in the markers used. This might be due to the slower evolutionary rates of lichenized fungi from arid and subarid regions compared to oceanic parmelioid lichens [36]. The diversification age of clade 1 was estimated at 5.31 Ma ( Ma), at the end of the Miocene, which was a geological period when numerous groups radiated in arid conditions. The Mediterranean region was suggested as the ancestral area of this clade by the ancestral range analysis, but the result was poorly supported. PLoS ONE 8 June 2012 Volume 7 Issue 6 e39683

9 Clade 2 also contains species with orcinol depsides and depsidones and comprises two disjunct lineages, one occurring in Australia (clade 2.1) and the other in South America (2.2). The estimated age for clade 2 (8.10 Ma; Ma) rules out the possibility of vicariance, since the breakup of Australia, Antarctica and South America occurred between Ma ago [53]. The ancestral range reconstruction points to Australia as the ancestral area of clade 2. This would be consistent with long distance dispersal from Australia to South America, a phenomenon frequently found in many plants groups [54]. Within the Australian clade several strongly supported lineages are not consistent with the current species delimitation of the group, indicating that the phenotypical characters used to distinguish species in the group have limited phylogenetic validity. Similar disparities between phylogenetic relationships and current species delimitations were found within the yellow Xanthoparmelia species from western North American [55,56]. Clade 4 includes specimens containing hypostictic acid from California and South America. Here again all the South American specimens form a monophyletic group. Specimens of X. subhosseana occurring in different continents are not closely related. The most likely ancestral origin of clade 4 is South America (marginal probability of 0.48), although neither California nor could be rejected. Our study indicates that the X. pulla group started to radiate during the Miocene in, where the highest diversity of this group is found. From this region, different lineages with distinct secondary metabolites belonging to different chemical pathways were dispersed to other regions, where they experienced rapid and more recent radiations. In some cases our results showed that the sympatric species of the X. pulla group in an area belong to distantly related groups. For example, the Californian X. pulla flora includes species from clade 1 and clade 4.1, the latter most probably having migrated from South America. Indeed our study indicated that the current taxonomic circumscription of species in the group does not agree with the evolutionary hypotheses inferred by molecular markers. The incongruence of phenotype-based classification and molecular phylogeny is a challenge for the classification of these fungi. Additional studies will be needed to determine whether the lineages found here represent cryptic species or whether new phenotypical characters can be found to distinguish these distinct lineages (as has been found in some other Parmeliaceae [10,57]). Future research should address how such parallel evolution of phenotypical characters in lichenized fungi could be explained in order to provide a better framework to test the adaptive value of these characters [58]. Our results here have important implications for conservation and ecological issues, since species were found to have much more restricted distribution than previously thought. Materials and Methods Taxon sampling Eighty two specimens of the Xanthoparmelia pulla group from California, the Mediterranean basin, Macaronesia, South America, Australia and were used for the phylogenetic study. The specimens were identified following the current species delimitations [38,44,47,59,60]. Chemical constituents were identified using thin layer chromatography (TLC) [61 64], and gradient-elution high performance liquid chromatography (HPLC) [65]. The major medullary compounds were classified into four major groups based on their chemical structure: 1) Orcinol depsides: olivetoric, divaricatic, stenosporic and glomelliferic acids. 2) Orcinol depsidones, physodic, alectoronic and glomelliferonic acids. 3) b-orcinol depsidones: stictic acid (only present in the outgroup) and hypostictic acid. 4) aliphatic acids: constipatic acid. All necessary permits were obtained for the described field studies. Collecting permits in Australia were all obtained by J.A. Elix (ca. 50 permit numbers for each states and over several years) and in Chile by W. Quilhot. For European locations specific permission was not required, since the locations were neither in privatelyowned or protected areas. The field studies did not involve endangered or protected species. Molecular study Total DNA was extracted from frozen lobes of thalli crushed with sterile glass pestles, using the DNeasy Plant Mini Kit (Qiagen) following the manufacturer s instructions and modifications of Crespo et al. [66]. The following primers were used: ITS1-LM [67] and ITS2-KL [68] for nuits rdna, and LR0R and LR5 [69] for nulsu rdna. For each amplification we used a reaction mixture of 25 ml, containing: 2.5 ml of 10x DNA polymerase buffer (including MgCl 2 2mM) (Biotools), 1.25 ml of each primer, 0.75 ml of DNA polymerase (1U/mL), 0.5 ml of dntps containing 10 mm of each base (Biotools), 5 ml of DNA (third elution of DNA extraction) and 13.5 ml dh 2 O. Amplifications were carried out in an automatic thermocycler (Techne Progene 3000) with the following steps: an initial denaturation at 94uC for 5 min; 35 cycles of 94uC for 1 min, 58uC (nuits rdna) or 56uC (nulsu rdna) for 1 min, and 72uC for 1.5 min; a final extension at 72uC for 5 min. PCR products were cleaned with DNA Purification Kit (Flavorgen) and sequenced with the same primers using the ABI Prism Dye Terminator Cycle Sequencing Ready reaction kit (Applied Biosystems) with the following program: initial denaturation at 94uC for 3 min, 25 cycles at 96uC for 10s, 50uC for 5s and 60uC for 4 min. Sequencing reactions were electrophoresed on a 3730 DNA analyzer (Applied Biosystems). The sequence fragments were assembled with Bioedit v. 7.0 [70] and manually adjusted. Sequence alignment and selection of the substitution model We used a dataset of 2 loci of 82 specimens representing 25 species of the Xanthoparmelia pulla group and 6 specimens as outgroup. The sequences were mainly generated in this study (140 sequences) and the others taken from our previous studies [33,41]. The outgroup selection was based on previous phylogenetic studies [41]. GenBank accession numbers and sources of the specimens are listed in Table 2. The two loci were aligned separately with Muscle 3.6 [71] and the ambiguous positions were removed manually. The general time reversible model including estimation of invariant sites (GTR+I+G) was selected by jmodeltest v [72] as the most appropriate nucleotide substitution model for both loci. Phylogenetic Analyses Potential conflict between the two loci was assessed by comparison of the ML analyses obtained with Garli 0.96 [73] for each locus, using 100 pseudoreplicates for the bootstrap analyses. The phylogenetic analyses of the combined matrix were done using maximum likelihood (ML) and a Bayesian approach. ML analysis was performed using Garli 0.96 [73] with default settings and 100 replicates for the bootstrap analyses. The Bayesian analysis was performed using MrBayes [74] using the GTR+I+G model, and the data set partitioned into nu ITS and nu LSU. Each partition was allowed to have its own parameter values [75]. Heating of chains was set to 0.2, with 5 PLoS ONE 9 June 2012 Volume 7 Issue 6 e39683

10 million generations sampled every 500th tree. The first 1000 trees were discarded as burn in. We used AWTY [76] to compare splits frequencies in the different runs and to plot cumulative split frequencies to insure that stationarity was reached. Of the remaining trees (9000 from each of the parallel runs) a majority rule consensus tree with average branch lengths and posterior probabilities was calculated using the sumt option of MrBayes. Clades with bootstrap support equal or above 70 % under ML and/or posterior probabilities $0.95 in the Bayesian analysis were considered as strongly supported. Phylogenetic trees were visualized using the program FigTree [77]. Calibration of nodes and dating analysis The ages of the X. pulla group and its major clades were estimated by a divergence time analysis based on a calibrated phylogeny of the parmelioid lichens [78]. We used a matrix of two loci (nu ITS and LSU) with a proportional number of samples of each parmelioid clade to have a representative tree and trend in speciation through time [79]. The matrix included 299 specimens of parmelioid lichens and 3 specimens of the genus Usnea (as outgroup); 62 new sequences of Xanthoparmelia species outside the X. pulla group were included. GenBank accession numbers with the specimens of the dating analysis are listed in Table S1. The sequence alignment, selection of the nucleotide substitution model, and phylogenetic analyses were done using the same procedures used for the Xanthoparmelia pulla dataset (see above). The divergence time analyses were performed using BEAST v [80]. We used a starting tree obtained from a ML analysis using Garli [73] of the concatenated dataset, calculated an ultrametric tree using nonparametric rate smoothing (NPRS) implemented in TreeEdit v.10a10 [81]. The age of the crown node of the parmelioid lichens was calibrated at 60 Ma, following Amo de Paz et al. [78]. The starting tree was topologically congruent with the parmelioid phylogeny presented in Crespo et al. [32]. For the divergence time analyses we used two points of calibration: the age of the crown node of the parmelioid lichens set at Ma ( Ma) [78], and the age of the crown node of the genus Parmelia (dated from the fossil Parmelia ambra from the Dominican amber, Ma, [82] as discussed previously) [78]. The BEAST analysis was performed using the GTR+I+G substitution model, a Birth-Death process tree prior, and a relaxed clock model (uncorrelated lognormal) for the concatenated dataset. Calibration points were defined as prior distribution, minimal ages and calibrated with a lognormal distributions: 1) the parmelioid crown node at uniform distribution between Ma; 2) the Parmelia crown node at log-normal mean = 2.77, offset = 14, lognormal standard deviation = 0.5. The analysis was run for 10 million generations, with parameter values sampled every 1000 generation. We checked the stationary plateau with Tracer v. 1.5 [83]. We discarded 10% of the initial trees as burn in and the References 1. Wiens JJ, Penkrot TA (2002) Delimiting species using DNA and morphological variation and discordant species limits in spiny lizards (Sceloporus). Systematic Biology 51: Cracraft J (1983) Species concepts and speciation analysis. Current Ornithology 1: Mayr E (1963) Animal species and evolution. Cambridge, MA: Harvard University Press. 4. Arguello A, Del Prado R, Cubas P, Crespo A (2007) Parmelina quercina (Parmeliaceae, Lecanorales) includes four phylogenetically supported morphospecies. Biological Journal of the Linnean Society 91: Leavitt SD, Fankhauser JD, Leavitt DH, Porter LD, Johnson LA, et al. (2011) Complex patterns of speciation in cosmopolitan rock posy lichens consensus tree was calculated using Tree Annotator v [80]. The results were visualized with FigTree v [84]. Ages of the X. pulla clades were estimated for nodes with more than 0.95 of posterior probability in the BEAST runs and in the previous Bayesian analysis. Ancestral range reconstructions The biogeographical analysis to reconstruct the ancestral area was performed using an indirect Bayesian approach to character state reconstruction [85] implemented in SIMMAP v1.5 [86] following [87]. This analysis integrates the combination of the uncertainty in the tree, branch lengths and the substitution models using Markov chain Monte Carlo. We treated the biogeographic regions as discrete characters. The major areas in which X. pulla species are distributed were categorised broadly into five areas: California, Mediterranean basin (including Macaronesia), South America,, and Australia. Presence/absence was coded as binary states and each area was given equal probability. We performed the ancestral state reconstruction analysis on a subsample of 1000 trees derived from the MrBayes tree sampling of the Xanthoparmelia pulla group. We also performed ancestral range reconstruction analysis using dispersal-extinction-cladogenesis (DEC) implemented in Lagrange program [88]. The results were inconclusive due to the lack of confidence in parts of the X. pulla phylogeny and hence the results are not included in this paper. Supporting Information Figure S1 Chronogram of parmelioid lichens focusing in Xanthoparmelia pulla group. Calibration points: A, inferred age of radiation of parmelioid lichens and B, age of the Parmelia fossil. The Xanthoparmelia pulla group is highlighted by a box and the dated clades are indicated by a branch in bold. (TIF) Table S1 GenBank accession numbers of parmelioid lichens (except X. pulla group, see Table 2) used for divergence time analysis. New sequences are in bold. (DOC) Acknowledgments We thank L. Mucina (Perth) for his help in the organization of the field work in. GAP thanks W. Quilot (Valparaiso) and C. Rubio (Valparaiso) and C. Cargill (Canberra) who kindly facilitated the visits to the Universidad de Valparaiso (Chile) and ANU and CANBR (Australia). Author Contributions Conceived and designed the experiments: HTL AC PC. Performed the experiments: GAP. Analyzed the data: HTL PC GAP. Contributed reagents/materials/analysis tools: JAE GAP. Wrote the paper: HTL AC PC GAP. Discovering and delimiting cryptic fungal species in the lichen-forming Rhizoplaca melanophthalma species-complex (Lecanoraceae, Ascomycota). Molecular Phylogenetics and Evolution 59: Leavitt SD, Johnson L, St Clair LL (2011) Species delimitation and evolution in morphologially and chmically diverse communities of the lichen-forming genus Xanthoparmelia (Parmeliaceae, Ascomycota) in western North America. American Journal of Botany 98: Leavitt SD, Johnson LA, Goward T, St Clair LL (2011) Species delimitation in taxonomically difficult lichen-forming fungi: An example from morphologically and chemically diverse Xanthoparmelia (Parmeliaceae) in North America. Molecular Phylogenetics and Evolution 60: PLoS ONE 10 June 2012 Volume 7 Issue 6 e39683

11 8. Lumbsch HT, Leavitt SD (2011) Goodbye morphology? A paradigm shift in the delimitation of species in lichenized fungi Fungal Diversity 50: Elix JA, Corush J, Lumbsch HT (2009) Triterpene chemosyndromes and subtle morphological characters characterise lineages in the Physcia aipolia group in Australia (Ascomycota). Systematics and Biodiversity 7: Wirtz N, Printzen C, Lumbsch HT (2008) The delimitation of Antarctic and bipolar species of neuropogonoid Usnea (Ascomycota, Lecanorales): a cohesion approach of species recognition for the Usnea perpusilla complex. Mycological Research 112: Cubero OF, Crespo A, Esslinger TL, Lumbsch HT (2004) Molecular phylogeny of the genus Physconia (Ascomycota, Lecanorales) inferred from a Bayesian analysis of nuclear ITS rdna sequences. Mycological Research 108: Molina MD, Crespo A, Blanco O, Lumbsch HT, Hawksworth DL (2004) Phylogenetic relationships and species concepts in Parmelia s.str. (Parmeliaceae) inferred from nuclear ITS rdna and beta-tubulin sequences. Lichenologist 36: Crespo A, Molina MC, Blanco O, Schroeter B, Sancho LG, et al. (2002) rdna ITS and beta-tubulin gene sequence analyses reveal two monophyletic groups within the cosmopolitan lichen Parmelia saxatilis. Mycological Research 106: Divakar PK, Figueras G, Hladun NL, Crespo A (2010) Molecular phylogenetic studies reveal an undescribed species within the North American concept of Melanelixia glabra (Parmeliaceae). Fungal Diversity 42: Hodkinson BP, Lendemer JC (2011) Molecular analyses reveal semi-cryptic species in Xanthoparmelia tasmanica. Bibliotheca Lichenologica 106: Thell A, Elix JA, Søchting U (2009) Xanthoparmelia lineola s. l. in Australia and North America. Bibliotheca Lichenologica 99: Velmala S, Myllys L, Halonen P, Goward T, Ahti T (2009) Molecular data show that Bryoria fremontii and B. tortuosa (Parmeliaceae) are conspecific. Lichenologist 41: Otálora MAG, Martínez I, Aragón G, Molina MC (2010) Phylogeography and divergence date estimates of a lichen species complex with a disjunct distribution pattern. American Journal of Botany 97: Otalora MAG, Martinez I, Molina MC, Aragon G, Lutzoni F (2008) Phylogenetic relationships and taxonomy of the Leptogium lichenoides group (Collemataceae, Ascomycota) in Europe. Taxon 57: Sérusiaux E, Villarreal A JC, Wheeler T, Goffinet B (2011) Recent origin, active speciation and dispersal for the lichen genus Nephroma (Peltigerales) in Macaronesia. Journal of Biogeography: no no. 21. Wedin M, Westberg M, Crewe AT, Tehler A, Purvis OW (2009) Species delimitation and evolution of metal bioaccumulation in the lichenized Acarospora smaragdula (Ascomycota, Fungi) complex. Cladistics 25: Kroken S, Taylor JW (2001) A gene genealogical approach to recognize phylogenetic species boundaries in the lichenized fungus Letharia. Mycologia 93: Crespo A, Lumbsch HT (2010) Cryptic species in lichen-forming fungi. IMA Fungus 1: Divakar PK, Blanco O, Hawksworth DL, Crespo A (2005) Molecular phylogenetic studies on the Parmotrema reticulatum (syn. Rimelia reticulata) complex, including the confirmation of P. pseudoreticulatum as a distinct species. Lichenologist 37: Baloch E, Grube M (2009) Pronounced genetic diversity in tropical epiphyllous lichen fungi. Molecular Ecology 18: Högnabba F, Wedin M (2003) Molecular phylogeny of the Sphaerophorus globosus species complex. Cladistics 19: Fenchel T, Finlay BJ (2004) The ubiquity of small species: Patterns of local and global diversity. Bioscience 54: Finlay BJ (2002) Global dispersal of free-living microbial eukaryote species. Science 296: Lumbsch HT, Buchanan PK, May TW, Mueller GM (2008) Phylogeography and biogeography of fungi. Mycological Research 112: Crespo A, Lumbsch HT, Mattsson JE, Blanco O, Divakar PK, et al. (2007) Testing morphology-based hypotheses of phylogenetic relationships in Parmeliaceae (Ascomycota) using three ribosomal markers and the nuclear RPB1 gene. Molecular Phylogenetics and Evolution 44: Hale BW, DePriest PT (1999) Mason E. Hale s list of epithets in the parmelioid genera. Bryologist 102: Crespo A, Kauff F, Divakar PK, Amo G, Arguello A, et al. (2010) Phylogenetic generic classification of parmelioid lichens (Parmeliaceae, Ascomycota) based on molecular, morphological and chemical evidence. Taxon 59: Blanco O, Crespo A, Elix JA, Hawksworth DL, Lumbsch HT (2004) A molecular phylogeny and a new classification of parmelioid lichens containing Xanthoparmelia-type lichenan (Ascomycota: Lecanorales). Taxon 53: Del Prado R, Ferencova Z, Armas-Crespo V, De Paz GA, Cubas P, et al. (2007) The arachiform vacuolar body: an overlooked shared character in the ascospores of a large monophyletic group within Parmeliaceae (Xanthoparmelia clade, Lecanorales). Mycological Research 111: Blanco O, Crespo A, Ree RH, Lumbsch HT (2006) Major clades of parmeliold lichens (Parmeliaceae, Ascomycota) and the evolution of their morphological and chemical diversity. Molecular Phylogenetics and Evolution 39: Lumbsch HT, Hipp AL, Divakar PK, Blanco O, Crespo A (2008) Accelerated evolutionary rates in tropical and oceanic parmelioid lichens (Ascomycota). BMC Evolutionary Biology Esslinger TL (1978) A new status for the brown Parmeliae. Mycotaxon 7: Esslinger TL (2000) Notes on the brown-colored species of Parmeliaceae (lichenized Ascomycota) in southern Africa. Bryologist 103: Elix JA (2003) The lichen genus Paraparmelia, a synonym of Xanthoparmelia (Ascomycota, Parmeliaceae). Mycotaxon 87: Amo de Paz G, Lumbsch HT, Cubas P, Elix JA, Crespo A (2010) The morphologically deviating genera Omphalodiella and Placoparmelia belong to Xanthoparmelia (Parmeliaceae). The Bryologist 113: Amo de Paz G, Lumbsch HT, Cubas P, Elix JA, Crespo A (2010) The genus Karoowia (Parmeliaceae, Ascomycota) includes unrelated clades nested within Xanthoparmelia. Australian Systematic Botany 23: Thell A, Feuerer T, Elix JA, Kärnefelt I (2006) A contribution to the phylogeny and taxonomy of Xanthoparmelia (Ascomycota, Parmeliaceae). Journal of the Hattori Botanical Laboratory 100: Hawksworth DL, Crespo A (2002) Proposal to conserve the name Xanthoparmelia against Chondropsis nom. cons. (Parmeliaceae). Taxon 51: Esslinger TL (1977) A chemosystematic revision of the brown Parmeliae. Jour Hattori Bot Lab 42: Crespo A, Perez-Ortega S (2009) Cryptic species and species pairs in lichens: A discussion on the relationship between molecular phylogenies and morphological characters. Anales del Jardín Botánico de Madrid 66: Culberson CF, Culberson WL, Esslinger TL (1977) Chemosyndromic variation in the Parmelia pulla group. Bryologist 80: Elix JA (2002) Chemical variations of the lichen Neofuscelia pulla (Ascomycotina: Parmeliaceae) sensu Esslinger. Australasian Lichenology 51: Elix JA, Corush J, Lumbsch HT (2009) Triterpene chemosyndromes and subtle morphological characters characterise lineages in the Physcia aipolia group in Australia (Ascomycota). Systematics and Biodiversity 7: Richardson JE, Weitz FM, Fay MF, Cronk QCB, Linder HP, et al. (2001) Rapid and recent origin of species richness in the Cape flora of. Nature 412: Lumbsch HT, Schmitt I, Barker D, Pagel M (2006) Evolution of micromorphological and chemical characters in the lichen-forming fungal family Pertusariaceae. Biological Journal of the Linnean Society 89: Schmitt I, Lumbsch HT (2004) Molecular phylogeny of the Pertusariaceae supports secondary chemistry as an important systematic character set in lichenforming ascomycetes. Molecular Phylogenetics and Evolution 33: Lumbsch HT, Schmitt I (2002) Molecular data shake the Pertusariaceae tree into order. Lichenology 1: Sanmartin I, Ronquist F (2004) Southern Hemisphere biogeography inferred by event-based models: Plant versus animal patterns. Systematic Biology 53: Sanmartin I, Wanntorp L, Winkworth RC (2007) West Wind Drift revisited: testing for directional dispersal in the Southern Hemisphere using event-based tree fitting. Journal of Biogeography 34: Leavitt SD, Johnson L, St Clair LL (2011) Species delimitation and evolution in morphologically and chemically diverse communities of the lichen-forming genus Xanthoparmelia (Parmeliaceae, Ascomycota) in western North America. American Journal of Botany 98: Leavitt SD, Johnson LA, Goward T, St. Clair LL (2011) Species delimitation in taxonomically difficult lichen-forming fungi: an example from morphologically and chemically diverse Xanthoparmelia (Parmeliaceae) in North America. Molecular Phylogenetics and Evolution in press: xx xx. 57. Wirtz N, Printzen C, Lumbsch HT (2012) Using haplotype networks, estimation of gene flow and phenotypic characters to understand species delimitation in fungi of a predominantly Antarctic Usnea group (Ascomycota, Parmeliaceae). Organisms Diversity & Evolution 12: Rivas Plata E, Lumbsch HT (2011) Parallel evolution and phenotypic divergence in lichenized fungi: a case study in the lichen-forming fungal family Graphidaceae (Ascomycota: Lecanoromycetes: Ostropales). Molecular Phylogenetics and Evolution 61: Elix JA (1999) New species of Neofuscelia (Lichenized Ascomycotina, Parmeliaceae) from the Southern Hemisphere. Mycotaxon 71: Elix JA (1994) Neofuscelia. Flora of Australia 55: Culberson CF (1972) Improved conditions and new data for the identification of lichen products by a standardized thin-layer chromatographic method. Journal of Chromatography 72: Culberson CF, Culberson WL, Johnson A (1981) A standardized TLC analysis of ß-orcinol depsidones. Bryologist 84: Culberson CF, Johnson A (1982) Substitution of methyl tert.-butyl ether for diethyl ether in standardized thin-layer chromatographic method for lichen products. Journal of Chromatography 238: Orange A, James PW, White FJ (2001) Microchemical Methods for the Identification of Lichens: British Lichen Society. 101 p. 65. Feige GB, Lumbsch HT, Huneck S, Elix JA (1993) Identification of lichen substance by a standardized high-performance liquid-chromatographic method. Journal of Chromatography 646: Crespo A, Blanco O, Hawksworth DL (2001) The potential of mitochondrial DNA for establishing phylogeny and stabilising generic concepts in the parmelioid lichens. Taxon 50: Myllys L, Lohtander K, Källersjö M, Tehler A (1999) Sequence insertions and ITS data provide congruent information on Roccella canariensis and R. tuberculata (Arthoniales, Euascomycetes) phylogeny. Molecular Phylogenetics and Evolution 12: PLoS ONE 11 June 2012 Volume 7 Issue 6 e39683