Testing morphology-based hypotheses of phylogenetic relationships in Parmeliaceae (Ascomycota) using three ribosomal markers and the nuclear RPB1 gene

Global phylogeny of the family Gomphillaceae (Ascomycota, Graphidales) sheds light on the origin, diversification and endemism in foliicolous lineages

Foliicolous lichens grow on living leaves of vascular plants. They are mostly found in tropical to subtropical or temperate rainforests. Many phenotype-based species are considered as pantropical or even sub-cosmopolitan, either attributed to old ages, having existed prior to continental breakups or long-distance dispersal. We built a much expanded, global phylogeny of Gomphillaceae, the most diverse group of leaf-dwelling lichenised fungi. Our sampling encompassed six major biodiversity hotspots: MIOI (Madagascar and the Indian Ocean Islands), the Caribbean, New Caledonia, the Colombian Chocó, Mesoamerica and the Atlantic coast of Brazil. It was based on multilocus sequence data (mtSSU rDNA, nuLSU rDNA and RPB1), including 2207 sequences of 1256 specimens. Species delimitation methods combined with a phenotype matrix identified 473 putative species. Amongst these, 104 are confirmed as described, 213 are classified as cryptic or near cryptic (hidden diversity), 100 represent new species to science (identified on the basis of phenotype) and 56 remain unidentified. Amongst the 104 species with a valid name, 40.5% are distributed across 2-5 continents (lichenogeographical regions) by applying the phenotype-based species concept. However, using the integrative approach to delineate species, this estimate is reduced to 9%. We estimate the global species richness of Gomphillaceae at 1,861-2,356 species. The timing of species-level divergences suggests that the current distribution of foliicolous lichens is shaped more by long-distance dispersal and rapid diversification than by vicariance. The origin of the family and major clades appears to be in the Neotropics, with subsequent numerous dispersal events. Our results support the separation of three major lineages, corresponding to the former families Asterothyriaceae, Gomphillaceae s.str. and Solorinellaceae, which should be recognised at the subfamily level.

The First Miniature, Small Foliose, Brown Xanthoparmelia in the Northern Hemisphere

The genus Xanthoparmelia includes several subcrustose, squamulose, small foliose, and small subfruticose species, primarily in the Southern Hemisphere. Here, we report on the first small foliose species lacking usnic acid in the genus occurring in the Holarctic. The species has been previously known as Lecanora olivascens Nyl., but subsequent studies of the morphology, secondary chemistry, and molecular data of the nuITS rDNA indicate that this species instead belongs to Xanthoparmelia. Consequently, the new combination Xanthoparmelia olivascens (Nyl.) V.J. Rico and G. Amo is proposed, and an epitype is designated here. We discuss the unique presence of a subcrustose Xanthoparmelia species lacking cortical usnic acid in the Northern Hemisphere. This species fits phylogenetically into a clade that was previously only known from the Southern Hemisphere, and hence represents another example of N-S disjunction in lichenized fungi.

Mycobiont-specific primers facilitate the amplification of mitochondrial small subunit ribosomal DNA: a focus on the lichenized fungal genus Melanelia (Ascomycota, Parmeliaceae) in Iceland

The fungal mitochondrial small subunit (mtSSU) ribosomal DNA is one of the most commonly used loci for phylogenetic analysis of lichen-forming fungi, but their primer specificity to mycobionts has not been evaluated. The current study aimed to design mycobiont-specific mtSSU primers and highlights their utility with an example from the saxicolous lichen-forming fungal genus Melanelia Essl. in Iceland. The study found a 12.5% success rate (3 out of 24 specimens with good-quality mycobiont mtSSU sequences) using universal primers (i.e. mrSSU1 and mrSSU3R), not including off-target amplification of environmental fungi, e.g. Cladophialophoracarrionii and Lichenotheliaconvexa. New mycobiont-specific primers (mt-SSU-581-5' and mt-SSU-1345-3') were designed by targeting mycobiont-specific nucleotide sites in comparison with environmental fungal sequences, and assessed for mycobiont primer specificity using in silico PCR. The new mycobiont-specific mtSSU primers had a success rate of 91.7% (22 out of 24 specimens with good-quality mycobiont mtSSU sequences) on the studied Melanelia specimens. Additional testing confirmed the specificity and yielded amplicons from 79 specimens of other Parmeliaceae mycobiont lineages. This study highlights the effectiveness of designing mycobiont-specific primers for studies on lichen identification, barcoding and phylogenetics.

Diversity and Host Relationships of the Mycoparasite Sepedonium (Hypocreales, Ascomycota) in Temperate Central Chile

We present the first major survey of regional diversity, distribution and host-association of Sepedonium. Whereas the rather scarce worldwide records of this mycoparasitic fungus suggested no specific distribution pattern of most species before, we provide new evidence of endemic and specific host-parasite guilds of Sepedonium in Southern South America, including the description of a new species. The corresponding inventory was performed in temperate central Chile. The regional landscape, a mosaic of exotic timber plantations and remnants of native Nothofagus forests, facilitates a unique combination of endemic and adventitious Boletales hosts. During a two-year survey, 35 Sepedonium strains were isolated and cultured from infected basidiomata of allochthonous Chalciporus piperatus, Paxillus involutus, Rhizopogon spp. and Suillus spp., as well as from the native Boletus loyita, B. loyo, B. putidus and Gastroboletus valdivianus. Taxonomic diagnosis included morphology of conidia and conidiophores, sequences of ITS, RPB2 and EF1 molecular markers and characteristics of in vitro cultures. Phylogenetic reconstructions were performed using Bayesian methods. Four Sepedonium species could be identified and characterized, viz.: S. ampullosporum, S. chrysospermum, S. laevigatum and the newly described species S. loyorum. The most frequent species on introduced Boletales was S. ampullosporum, followed by S. chrysospermum and S. laevigatum. S. loyorum sp. nov. was found exclusively on native boletacean hosts, separated from its closest relative S. chalcipori by micromorphological and molecular attributes. Species descriptions and identification keys are provided. Ecological and biogeographical aspects of endemic and allochthonous symbiotic units consisting of mycoparasite, ectomycorrhizal fungal host and respective mycorrhizal tree are discussed.

Using target enrichment sequencing to study the higher-level phylogeny of the largest lichen-forming fungi family: Parmeliaceae (Ascomycota)

Parmeliaceae is the largest family of lichen-forming fungi with a worldwide distribution. We used a target enrichment data set and a qualitative selection method for 250 out of 350 genes to infer the phylogeny of the major clades in this family including 81 taxa, with both subfamilies and all seven major clades previously recognized in the subfamily Parmelioideae. The reduced genome-scale data set was analyzed using concatenated-based Bayesian inference and two different Maximum Likelihood analyses, and a coalescent-based species tree method. The resulting topology was strongly supported with the majority of nodes being fully supported in all three concatenated-based analyses. The two subfamilies and each of the seven major clades in Parmelioideae were strongly supported as monophyletic. In addition, most backbone relationships in the topology were recovered with high nodal support. The genus Parmotrema was found to be polyphyletic and consequently, it is suggested to accept the genus Crespoa to accommodate the species previously placed in Parmotrema subgen. Crespoa. This study demonstrates the power of reduced genome-scale data sets to resolve phylogenetic relationships with high support. Due to lower costs, target enrichment methods provide a promising avenue for phylogenetic studies including larger taxonomic/specimen sampling than whole genome data would allow.

A new species concept for the clinically relevant Mucor circinelloides complex

Mucor species are common soil fungi but also known as agents of human infections (mucormycosis) and used in food production and biotechnology. Mucor circinelloides is the Mucor species that is most frequently isolated from clinical sources. The taxonomy of Mucor circinelloides and its close relatives (Mucor circinelloides complex – MCC) is still based on morphology and mating behaviour. The aim of the present study was a revised taxonomy of the MCC using a polyphasic approach. Using a set of 100 strains molecular phylogenetic analysis of five markers (ITS, rpb1, tsr1, mcm7, and cfs, introduced here) were performed, combined with phenotypic studies, mating tests and the determination of the maximum growth temperatures. The multi-locus analyses revealed 16 phylogenetic species of which 14 showed distinct phenotypical traits and were recognised as discrete species. Five of these species are introduced as novel taxa: M. amethystinus sp. nov., M. atramentarius sp. nov., M. variicolumellatus sp. nov., M. pseudocircinelloides sp. nov., and M. pseudolusitanicus sp. nov. The former formae of M. circinelloides represent one or two separate species. In the MCC, the simple presence of well-shaped zygospores only indicates a close relation of both strains, but not necessarily conspecificity. Seven species of the MCC have been implemented in human infection: M. circinelloides, M. griseocyanus, M. janssenii, M. lusitanicus, M. ramosissimus, M. variicolumellatus, and M. velutinosus.

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha-proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein-coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome-wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad-scale genome streamlining via loss of protein-coding genes as well as noncoding, parasitic DNA elements.

Lichen chemistry is concordant with multilocus gene genealogy in the genus Cetrelia (Parmeliaceae, Ascomycota)

The lichen genus Cetrelia represents a taxonomically interesting case where morphologically almost uniform populations differ considerably from each other chemically. Similar variation is not uncommon among lichenized fungi, but it is disputable whether such populations should be considered entities at the species level. Species boundaries in Cetrelia are traditionally delimited either as solely based on morphology or as combinations of morpho- and chemotypes. A dataset of four nuclear markers (ITS, IGS, Mcm7, RPB1) from 62 specimens, representing ten Cetrelia species, was analysed within Bayesian and maximum likelihood frameworks. Analyses recovered a well-resolved phylogeny where the traditional species generally were monophyletic, with the exception of Cetrelia chicitae and Cetrelia pseudolivetorum. Species delimitation analyses supported the distinction of 15 groups within the studied Cetrelia taxa, dividing three traditionally identified species into some species candidates. Chemotypes, distinguished according to the major medullary substance, clearly correlated with clades recovered within Cetrelia, while samples with the same reproductive mode were dispersed throughout the phylogenetic tree. Consequently, delimiting Cetrelia species based only on reproductive morphology is not supported phylogenetically. Character analyses suggest that chemical characters have been more consistent compared to reproductive mode and indicate that metabolite evolution in Cetrelia towards more complex substances is probable.

Thorsten Lumbsch. Population genomic analyses of RAD sequences resolves the phylogenetic relationship of the lichen-forming fungal species Usneaantarctica and Usneaaurantiacoatra

Neuropogonoid species in the lichen-forming fungal genus Usnea exhibit great morphological variation that can be misleading for delimitation of species. We specifically focused on the species delimitation of two closely-related, predominantly Antarctic species differing in the reproductive mode and representing a so-called species pair: the asexual U.antarctica and the sexual U.aurantiacoatra. Previous studies have revealed contradicting results. While multi-locus studies based on DNA sequence data provided evidence that these two taxa might be conspecific, microsatellite data suggested they represent distinct lineages. By using RADseq, we generated thousands of homologous markers to build a robust phylogeny of the two species. Furthermore, we successfully implemented these data in fine-scale population genomic analyses such as DAPC and fineRADstructure. Both Usnea species are readily delimited in phylogenetic inferences and, therefore, the hypothesis that both species are conspecific was rejected. Population genomic analyses also strongly confirmed separated genomes and, additionally, showed different levels of co-ancestry and substructure within each species. Lower co-ancestry in the asexual U.antarctica than in the sexual U.aurantiacoatra may be derived from a wider distributional range of the former species. Our results demonstrate the utility of this RADseq method in tracing population dynamics of lichens in future analyses.

Database establishment for the secondary fungal DNA barcode translational elongation factor 1α (TEF1α) 1

With new or emerging fungal infections, human and animal fungal pathogens are a growing threat worldwide. Current diagnostic tools are slow, non-specific at the species and subspecies levels, and require specific morphological expertise to accurately identify pathogens from pure cultures. DNA barcodes are easily amplified, universal, short species-specific DNA sequences, which enable rapid identification by comparison with a well-curated reference sequence collection. The primary fungal DNA barcode, ITS region, was introduced in 2012 and is now routinely used in diagnostic laboratories. However, the ITS region only accurately identifies around 75% of all medically relevant fungal species, which has prompted the development of a secondary barcode to increase the resolution power and suitability of DNA barcoding for fungal disease diagnostics. The translational elongation factor 1α (TEF1α) was selected in 2015 as a secondary fungal DNA barcode, but it has not been implemented into practice, due to the absence of a reference database. Here, we have established a quality-controlled reference database for the secondary barcode that together with the ISHAM-ITS database, forms the ISHAM barcode database, available online at http://its.mycologylab.org/ . We encourage the mycology community for active contributions.

The complete mitochondrial genomes of five lichenized fungi in the genus Usnea (Ascomycota: Parmeliaceae)

Known colloquially as 'Old Man's Beard', Usnea is a genus of lichenized Ascomycete fungi characterized by having a fruticose growth form and cartilaginous central axis. The complete mitochondrial genomes of Usnea halei, U. mutabilis, U. subfusca, U. subgracilis, and U. subscabrosa were sequenced using Illumina data and then assembled de novo. These mitogenomes ranged in size from 52,486 bp (U. subfusca) to 94,464 bp (U. subgracilis). All were characterized by having high levels of intronic and intergenic variation, such as ORFs that encode proteins with homology to two homing endonuclease types, LAGLIDADG and GIY-YIG. Genes annotated within these mitogenomes include 14 protein-coding genes, the large and small ribosomal subunits (LSU and SSU), and 23-26 tRNAs. Notably, the atp9 gene was absent from each genome. Genomic synteny was highly conserved across the five species. Five conserved mitochondrial genes (nad2, nad4, cox1, cox2, and cox3) were used to infer a best estimate maximum likelihood phylogeny among these five Usnea and other relatives, which yielded relationships consistent with prior published phylogenies.

Lichenized Fungi and the Evolution of Symbiotic Organization

Lichen symbioses comprise a fascinating relationship between algae and fungi. The lichen symbiotic lifestyle evolved early in the evolution of ascomycetes and is also known from a few basidiomycetes. The ascomycete lineages have diversified in the lichenized stage to give rise to a tremendous variety of morphologies. Their thalli are often internally complex and stratified for optimized integration of algal and fungal metabolisms. Thalli are frequently colonized by specific nonlichenized fungi and occasionally also by other lichens. Microscopy has revealed various ways these fungi interact with their hosts. Besides the morphologically recognizable diversity of the lichen mycobionts and lichenicolous (lichen-inhabiting) fungi, many other microorganisms including other fungi and bacterial communities are now detected in lichens by culture-dependent and culture-independent approaches. The application of multi-omics approaches, refined microscopic techniques, and physiological studies has added to our knowledge of lichens, not only about the taxa involved in the lichen interactions, but also about their functions.

Evaluation of the Antioxidant Capacities and Cytotoxic Effects of Ten Parmeliaceae Lichen Species

Parmeliaceae represents the largest and widespread family of lichens and includes species that attract much interest regarding pharmacological activities, due to their production of unique secondary metabolites. The current work aimed to investigate the in vitro antioxidant and cytotoxic activities of the methanol extracts of ten Parmeliaceae species, collected in different continents. Methanol extraction afforded high phenolic content in the extracts. The antioxidant activity displayed by lichens was evaluated through chemical assays, such as the ORAC (Oxygen Radical Absorbance Capacity) and 1,1-diphenyl-2-picrylhydrazyl (DPPH) radical scavenging activities and the ferric reducing antioxidant power (FRAP). A moderately positive correlation was found between the phenolic content and the antioxidant properties for all the species: R: 0.7430 versus ORAC values, R: 0.7457 versus DPPH scavenging capacity, and R: 0.7056 versus FRAP reducing power. The methanol extract of Flavoparmelia euplecta exhibited the highest ORAC value, the extract of Myelochroa irrugans showed the maximum DPPH scavenging capacity, and Hypotrachyna cirrhata methanol extract demonstrated the highest reducing power. Further, the cytotoxic activity of the ten species was investigated on the human cancer cell lines HepG2 and MCF-7; Myelochroa irrugans exhibited the highest anticancer potential. The pharmacological activities shown here could be attributed to their phytochemical constituents.

Limitations of Species Delimitation Based on Phylogenetic Analyses: A Case Study in the Hypogymnia hypotrypa Group (Parmeliaceae, Ascomycota)

Delimiting species boundaries among closely related lineages often requires a range of independent data sets and analytical approaches. Similar to other organismal groups, robust species circumscriptions in fungi are increasingly investigated within an empirical framework. Here we attempt to delimit species boundaries in a closely related clade of lichen-forming fungi endemic to Asia, the Hypogymnia hypotrypa group (Parmeliaceae). In the current classification, the Hypogymnia hypotrypa group includes two species: H. hypotrypa and H. flavida, which are separated based on distinctive reproductive modes, the former producing soredia but absent in the latter. We reexamined the relationship between these two species using phenotypic characters and molecular sequence data (ITS, GPD, and MCM7 sequences) to address species boundaries in this group. In addition to morphological investigations, we used Bayesian clustering to identify potential genetic groups in the H. hypotrypa/H. flavida clade. We also used a variety of empirical, sequence-based species delimitation approaches, including: the “Automatic Barcode Gap Discovery” (ABGD), the Poisson tree process model (PTP), the General Mixed Yule Coalescent (GMYC), and the multispecies coalescent approach BPP. Different species delimitation scenarios were compared using Bayes factors delimitation analysis, in addition to comparisons of pairwise genetic distances, pairwise fixation indices (F_ST). The majority of the species delimitation analyses implemented in this study failed to support H. hypotrypa and H. flavida as distinct lineages, as did the Bayesian clustering analysis. However, strong support for the evolutionary independence of H. hypotrypa and H. flavida was inferred using BPP and further supported by Bayes factor delimitation. In spite of rigorous morphological comparisons and a wide range of sequence-based approaches to delimit species, species boundaries in the H. hypotrypa group remain uncertain. This study reveals the potential limitations of relying on distinct reproductive strategies as diagnostic taxonomic characters for Hypogymnia and also the challenges of using popular sequence-based species delimitation methods in groups with recent diversification histories.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

A Tale of Two Hyper-diversities: Diversification dynamics of the two largest families of lichenized fungi

Renewed interests in macroevolutionary dynamics have led to the proliferation of studies on diversification processes in large taxonomic groups, such as angiosperms, mammals, and birds. However, such a study has yet to be conducted in lichenized fungi--an extremely successful and diverse group of fungi. Analysing the most comprehensive time-calibrated phylogenies with a new analytical method, we illustrated drastically different diversification dynamics between two hyper-diverse families of lichenized fungi, Graphidaceae and Parmeliaceae, which represent more than a fourth of the total species diversity of lichenized fungi. Despite adopting a similar nutrition mode and having a similar number of species, Graphidaceae exhibited a lower speciation rate, while Parmeliaceae showed a sharp increase in speciation rate that corresponded with the aridification during the Oligocene-Miocene transition, suggesting their adaptive radiation into a novel arid habitat.

Coalescent-based species delimitation approach uncovers high cryptic diversity in the cosmopolitan lichen-forming fungal genus Protoparmelia (Lecanorales, Ascomycota)

Species recognition in lichen-forming fungi has been a challenge because of unsettled species concepts, few taxonomically relevant traits, and limitations of traditionally used morphological and chemical characters for identifying closely related species. Here we analyze species diversity in the cosmopolitan genus Protoparmelia s.l. The ~25 described species in this group occur across diverse habitats from the boreal -arctic/alpine to the tropics, but their relationship to each other remains unexplored. In this study, we inferred the phylogeny of 18 species currently assigned to this genus based on 160 specimens and six markers: mtSSU, nuLSU, ITS, RPB1, MCM7, and TSR1. We assessed the circumscription of species-level lineages in Protoparmelia s. str. using two coalescent-based species delimitation methods – BP&P and spedeSTEM. Our results suggest the presence of a tropical and an extra-tropical lineage, and eleven previously unrecognized distinct species-level lineages in Protoparmelia s. str. Several cryptic lineages were discovered as compared to phenotype-based species delimitation. Many of the putative species are supported by geographic evidence.

DNA barcoding of fungi causing infections in humans and animals

Correct species identification is becoming increasingly important in clinical diagnostics. Till now, many mycological laboratories rely on conventional phenotypic identification. But this is slow and strongly operator-dependent. Therefore, to improve the quality of pathogen identification, rapid, reliable, and objective identification methods are essential. One of the most encouraging approaches is molecular barcoding using the internal transcribed spacer (ITS) of the rDNA, which is rapid, easily achievable, accurate, and applicable directly from clinical specimens. It relies on the comparison of a single ITS sequence with a curated reference database. The International Society for Human and Animal Mycology (ISHAM) working group for DNA barcoding has recently established such a database, focusing on the majority of human and animal pathogenic fungi (ISHAM-ITS, freely accessible at http://www.isham.org/ or directly from http://its.mycologylab.org). For some fungi the use of secondary barcodes may be necessary.

International Society of Human and Animal Mycology (ISHAM)-ITS reference DNA barcoding database--the quality controlled standard tool for routine identification of human and animal pathogenic fungi

Human and animal fungal pathogens are a growing threat worldwide leading to emerging infections and creating new risks for established ones. There is a growing need for a rapid and accurate identification of pathogens to enable early diagnosis and targeted antifungal therapy. Morphological and biochemical identification methods are time-consuming and require trained experts. Alternatively, molecular methods, such as DNA barcoding, a powerful and easy tool for rapid monophasic identification, offer a practical approach for species identification and less demanding in terms of taxonomical expertise. However, its wide-spread use is still limited by a lack of quality-controlled reference databases and the evolving recognition and definition of new fungal species/complexes. An international consortium of medical mycology laboratories was formed aiming to establish a quality controlled ITS database under the umbrella of the ISHAM working group on "DNA barcoding of human and animal pathogenic fungi." A new database, containing 2800 ITS sequences representing 421 fungal species, providing the medical community with a freely accessible tool at http://www.isham.org/ and http://its.mycologylab.org/ to rapidly and reliably identify most agents of mycoses, was established. The generated sequences included in the new database were used to evaluate the variation and overall utility of the ITS region for the identification of pathogenic fungi at intra-and interspecies level. The average intraspecies variation ranged from 0 to 2.25%. This highlighted selected pathogenic fungal species, such as the dermatophytes and emerging yeast, for which additional molecular methods/genetic markers are required for their reliable identification from clinical and veterinary specimens.

A multigene phylogenetic synthesis for the class Lecanoromycetes (Ascomycota): 1307 fungi representing 1139 infrageneric taxa, 317 genera and 66 families

The Lecanoromycetes is the largest class of lichenized Fungi, and one of the most species-rich classes in the kingdom. Here we provide a multigene phylogenetic synthesis (using three ribosomal RNA-coding and two protein-coding genes) of the Lecanoromycetes based on 642 newly generated and 3329 publicly available sequences representing 1139 taxa, 317 genera, 66 families, 17 orders and five subclasses (four currently recognized: Acarosporomycetidae, Lecanoromycetidae, Ostropomycetidae, Umbilicariomycetidae; and one provisionarily recognized, 'Candelariomycetidae'). Maximum likelihood phylogenetic analyses on four multigene datasets assembled using a cumulative supermatrix approach with a progressively higher number of species and missing data (5-gene, 5+4-gene, 5+4+3-gene and 5+4+3+2-gene datasets) show that the current classification includes non-monophyletic taxa at various ranks, which need to be recircumscribed and require revisionary treatments based on denser taxon sampling and more loci. Two newly circumscribed orders (Arctomiales and Hymeneliales in the Ostropomycetidae) and three families (Ramboldiaceae and Psilolechiaceae in the Lecanorales, and Strangosporaceae in the Lecanoromycetes inc. sed.) are introduced. The potential resurrection of the families Eigleraceae and Lopadiaceae is considered here to alleviate phylogenetic and classification disparities. An overview of the photobionts associated with the main fungal lineages in the Lecanoromycetes based on available published records is provided. A revised schematic classification at the family level in the phylogenetic context of widely accepted and newly revealed relationships across Lecanoromycetes is included. The cumulative addition of taxa with an increasing amount of missing data (i.e., a cumulative supermatrix approach, starting with taxa for which sequences were available for all five targeted genes and ending with the addition of taxa for which only two genes have been sequenced) revealed relatively stable relationships for many families and orders. However, the increasing number of taxa without the addition of more loci also resulted in an expected substantial loss of phylogenetic resolving power and support (especially for deep phylogenetic relationships), potentially including the misplacements of several taxa. Future phylogenetic analyses should include additional single copy protein-coding markers in order to improve the tree of the Lecanoromycetes. As part of this study, a new module ("Hypha") of the freely available Mesquite software was developed to compare and display the internodal support values derived from this cumulative supermatrix approach.

Host switching promotes diversity in host-specialized mycoparasitic fungi: uncoupled evolution in the Biatoropsis-usnea system

Fungal mycoparasitism-fungi parasitizing other fungi-is a common lifestyle in some basal lineages of the basidiomycetes, particularly within the Tremellales. Relatively nonaggressive mycoparasitic fungi of this group are in general highly host specific, suggesting cospeciation as a plausible speciation mode in these associations. Species delimitation in the Tremellales is often challenging because morphological characters are scant. Host specificity is therefore a great aid to discriminate between species but appropriate species delimitation methods that account for actual diversity are needed to identify both specialist and generalist taxa and avoid inflating or underestimating diversity. We use the Biatoropsis-Usnea system to study factors inducing parasite diversification. We employ morphological, ecological, and molecular data-methods including genealogical concordance phylogenetic species recognition (GCPSR) and the general mixed Yule-coalescent (GMYC) model-to assess the diversity of fungi currently assigned to Biatoropsis usnearum. The degree of cospeciation in this association is assessed with two cophylogeny analysis tools (ParaFit and Jane 4.0). Biatoropsis constitutes a species complex formed by at least seven different independent lineages and host switching is a prominent force driving speciation, particularly in host specialists. Combining ITS and nLSU is recommended as barcode system in tremellalean fungi.

Pleistocene speciation in North American lichenized fungi and the impact of alternative species circumscriptions and rates of molecular evolution on divergence estimates

Pleistocene climatic fluctuations influenced patterns of genetic variation and promoted speciation across a wide range of species groups. Lichens are commonly found in habitats that were directly impacted by glacial cycles; however, the role of Pleistocene climate in driving speciation in most lichen symbionts remains unclear. This uncertainty is due in part to limitations in our ability to accurately recognize independently evolving lichen-forming fungal lineages and a lack of relevant fossil calibrations. Using a coalescent-based species tree approach, we estimated divergence times for two sister clades in the genus Xanthoparmelia (Parmeliaceae) restricted to western North America. We assessed the influence of two different species circumscription scenarios and various locus-specific rates of molecular evolution on divergence estimates. Species circumscriptions were validated using the program BP&P. although speciation was generally supported in both scenarios, divergence times differed between traditional species circumscriptions and those based on genetic data, with more recent estimates resulting from the former. Similarly, rates of evolution for different loci resulted in variable divergence time estimates. However, our results unambiguously indicate that diversification in the sampled Xanthoparmelia clades occurred during the Pleistocene. Our study highlights the potential impact of ambiguous species circumscriptions and uncertain rates of molecular evolution on estimating divergence times within a multilocus species tree framework.

Understanding phenotypical character evolution in parmelioid lichenized fungi (Parmeliaceae, Ascomycota)

Parmelioid lichens form a species-rich group of predominantly foliose and fruticose lichenized fungi encompassing a broad range of morphological and chemical diversity. Using a multilocus approach, we reconstructed a phylogeny including 323 OTUs of parmelioid lichens and employed ancestral character reconstruction methods to understand the phenotypical evolution within this speciose group of lichen-forming fungi. Specifically, we were interested in the evolution of growth form, epicortex structure, and cortical chemistry. Since previous studies have shown that results may differ depending on the reconstruction method used, here we employed both maximum-parsimony and maximum-likelihood approaches to reconstruct ancestral character states. We have also implemented binary and multistate coding of characters and performed parallel analyses with both coding types to assess for potential coding-based biases. We reconstructed the ancestral states for nine well-supported major clades in the parmelioid group, two higher-level sister groups and the ancestral character state for all parmelioid lichens. We found that different methods for coding phenotypical characters and different ancestral character state reconstruction methods mostly resulted in identical reconstructions but yield conflicting inferences of ancestral states, in some cases. However, we found support for the ancestor of parmelioid lichens having been a foliose lichen with a non-pored epicortex and pseudocyphellae. Our data suggest that some traits exhibit patterns of evolution consistent with adaptive radiation.

Neogene-dominated diversification in neotropical montane lichens: dating divergence events in the lichen-forming fungal genus Oropogon (Parmeliaceae)

PREMISE OF THE STUDY

Diversification in neotropical regions has been attributed to both Tertiary geological events and Pleistocene climatic fluctuations. However, the timing and processes driving speciation in these regions remain unexplored in many important groups. Here, we address the timing of diversification in the neotropical lichenized fungal genus Oropogon (Ascomycota) and assess traditional species boundaries.

METHODS

We analyzed sequence data from three loci to assess phenotypically circumscribed Oropogon species from the Oaxacan Highlands, Mexico. We provide a comparison of dated divergence estimates between concatenated gene trees and a calibrated multilocus species-tree using substitution rates for two DNA regions. We also compare estimates from a data set excluding ambiguously aligned regions and a data set including the hyper-variable regions in two ribosomal markers.

KEY RESULTS

Phylogenetic reconstructions were characterized by well-supported monophyletic clades corresponding to traditionally circumscribed species, with the exception of a single taxon. Divergence estimates indicate that most diversification of the sampled Oropogon species occurred throughout the Oligocene and Miocene, although diversification of a single closely related clade appears to have occurred during the late Pliocene and into the Pleistocene. Divergence estimates calculated from a data set with ambiguously aligned regions removed were much more recent than those from the full data set.

CONCLUSIONS

Overall, our analyses place the majority of divergence events of Oropogon species from the Oaxacan Highlands within the Neogene and provide strong evidence that climatic changes during the Pleistocene were not a major factor driving speciation in the lichenized genus Oropogon in neotropical highlands.

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

Background

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.

Results

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.

Conclusions

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.

Transoceanic dispersal and subsequent diversification on separate continents shaped diversity of the Xanthoparmelia pulla group (Ascomycota)

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. South Africa 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.

Origin and diversification of major clades in parmelioid lichens (Parmeliaceae, Ascomycota) during the Paleogene inferred by Bayesian analysis

There is a long-standing debate on the extent of vicariance and long-distance dispersal events to explain the current distribution of organisms, especially in those with small diaspores potentially prone to long-distance dispersal. Age estimates of clades play a crucial role in evaluating the impact of these processes. The aim of this study is to understand the evolutionary history of the largest clade of macrolichens, the parmelioid lichens (Parmeliaceae, Lecanoromycetes, Ascomycota) by dating the origin of the group and its major lineages. They have a worldwide distribution with centers of distribution in the Neo- and Paleotropics, and semi-arid subtropical regions of the Southern Hemisphere. Phylogenetic analyses were performed using DNA sequences of nuLSU and mtSSU rDNA, and the protein-coding RPB1 gene. The three DNA regions had different evolutionary rates: RPB1 gave a rate two to four times higher than nuLSU and mtSSU. Divergence times of the major clades were estimated with partitioned BEAST analyses allowing different rates for each DNA region and using a relaxed clock model. Three calibrations points were used to date the tree: an inferred age at the stem of Lecanoromycetes, and two dated fossils: Parmelia in the parmelioid group, and Alectoria. Palaeoclimatic conditions and the palaeogeological area cladogram were compared to the dated phylogeny of parmelioid. The parmelioid group diversified around the K/T boundary, and the major clades diverged during the Eocene and Oligocene. The radiation of the genera occurred through globally changing climatic condition of the early Oligocene, Miocene and early Pliocene. The estimated divergence times are consistent with long-distance dispersal events being the major factor to explain the biogeographical distribution patterns of Southern Hemisphere parmelioids, especially for Africa-Australia disjunctions, because the sequential break-up of Gondwana started much earlier than the origin of these clades. However, our data cannot reject vicariance to explain South America-Australia disjunctions.

Phylogeny and character evolution in the jelly fungi (Tremellomycetes, Basidiomycota, Fungi)

The Tremellomycetes (Agaricomycotina, Basidiomycota, Fungi) are a nutritionally heterogeneous group comprising saprotrophs, animal parasites, and fungicolous species (fungal-inhabiting, including lichen-inhabiting). The relationships of many species, particularly those with a lichenicolous habit, have never been investigated by molecular methods. We present a phylogeny of the Tremellomycetes based on three nuclear DNA ribosomal markers (nSSU, 5.8S and nLSU), representing all main taxonomic groups and life forms, including lichenicolous taxa. The Cystofilobasidiales, Filobasidiales, Holtermanniales, and Tremellales (including the Trichosporonales) are recovered as monophyletic, but this is not the case for the Tremellomycetes. We suggest, however, that the Cystofilobasidiales tentatively continue to be included in the Tremellomycetes. As currently circumscribed, the Filobasidiaceae, Sirobasidiaceae, Syzygosporaceae and Tremellaceae are non-monophyletic. Cuniculitremaceae, Sirobasidiaceae and Tetragoniomycetaceae are nested within Tremellaceae. The lichenicolous species currently included within the Tremellomycetes belong in this group, distributed across the Filobasidiales and Tremellales. Lichen-inhabiting taxa do not form a monophyletic group; they are distributed in several clades and sometimes intermixed with taxa of other nutritional habits. Character state reconstruction indicates that two morphological traits claimed to characterize groups in the Tremellomycetes (the basidium habit and basidium septation) are highly homoplastic. Comparative phylogenetic methods suggest that the transitions between single and catenulate basidia in the Tremellales are consistent with a punctuational model of evolution whereas basidium septation is likely to have evolved under a graduational model in the clade comprising the Holtermanniales, Filobasidiales, and Tremellales.

Species delimitation and evolution in morphologically and chemically diverse communities of the lichen-forming genus Xanthoparmelia (Parmeliaceae, Ascomycota) in western North America

PREMISE OF THE STUDY

Accurate species delimitation is important for understanding the diversification of biota and has critical implications for ecological and conservation studies. However, a growing body of evidence indicates that morphology-based species circumspection in lichenized fungi misrepresents fungal diversity. The foliose lichen genus Xanthoparmelia includes over 800 species displaying a complex array of morphological and secondary metabolite diversity.

METHODS

We used a multifaceted approach, applying phylogenetic, population genetic, and genealogical analyses to delimit species in a single well-supported monophyletic clade containing 10 morphologically and chemically diverse Xanthoparmelia species in western North America. Sequence data from four ribosomal and two low-copy, protein-coding markers, along with chemical and morphological data were used to assess species diversity.

KEY RESULTS

We found that traditionally circumscribed species are not supported by molecular data. Rather, all sampled taxa were better represented by three polymorphic population clusters. Our results suggest that secondary metabolite variation may have limited utility in diagnosing lineages within this group, while identified populations clusters did not reflect major phylogeographic or ecological patterns.

CONCLUSIONS

In contrast to studies revealing previously undiscovered fungal lineages masked within lichen species circumscribed by traditional morphological and chemical concepts, the present study suggests that species diversity has been overestimated in the species-rich genus Xanthoparmelia.

Thallus anatomy of Canoparmelia texana (Parmeliaceae, lichenized Ascomycota)

Conventional techniques for structural studies under light microscope and scanning electron microscope were employed to describe the histology of thallus in Canoparmelia texana, a lichen with wide distribution in open environments and fairly common in the cerrados and urban areas of Brazil. This study describes a new type of cortical organization for the family Parmeliaceae, in C. texana the upper cortex is lacunar, showing a large quantity of small intercellular spaces or lacunae. The anatomical features including medulla thickness, hyphal orientation pattern, rhizines thickness and crystals inclusions have an important adaptive role for the success of C. texana in Brazilian cerrado, a savanna-like ecosystem.